Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,886,756 G→A 6.7% intergenic (+34/‑72) yidC → / → mnmE membrane protein insertase/tRNA U34 5‑methylaminomethyl‑2‑thiouridine modification GTPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,886,7560GA 6.7% 379.7 / 12.5 119intergenic (+34/‑72)yidC/mnmEmembrane protein insertase/tRNA U34 5‑methylaminomethyl‑2‑thiouridine modification GTPase
Reads supporting (aligned to +/- strand):  ref base G (54/57);  new base A (4/4);  total (58/61)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.61e-01

Only 100 of 121 total aligned reads displayed.

GGTCTGGTGCTGTACTATATCGTCAGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACGTGGCGGCGTTGGCATCCTGCGCATCTCCGG  >  NC_000913/3886609‑3886904
                                                                                                                                                   |                                                                                                                                                     
ggTCTGGTGCTGTACTATATCGTCAGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGt                                                                                                                                                      <  3:388573/149‑1 (MQ=255)
ggTCTGGTGCTGTACTATATCGTCAGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTAAGt                                                                                                                                                      <  3:454754/149‑1 (MQ=255)
 gTCTGGTGCTGTACTATATCGTCAGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGtt                                                                                                                                                     >  4:816407/1‑149 (MQ=255)
     ggTGCTGTACTATATCGTCAGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGAc                                                                                                                                                  <  2:516066/148‑1 (MQ=255)
       tGCTGTACTATATCGTCAGCAACCTGGTAACCATTATTCACCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGACGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAATGGCGGTCAACTGACCg                                                                                                                                                <  2:433083/148‑1 (MQ=255)
        gCTGTACTATATCGTCAGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑cc                                                                                                                                             >  4:522194/1‑149 (MQ=255)
        gCTGTACTATATCGTCAGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑c                                                                                                                                              >  1:4399/1‑148 (MQ=255)
         cTGTACTATATCGTCAGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAACTGACCG‑ccc                                                                                                                                            <  4:5846/149‑2 (MQ=255)
           gTACTATATCGTCAGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCtt                                                                                                                                           >  2:339948/1‑148 (MQ=255)
           gTACTATATCGTCAGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCtt                                                                                                                                           >  2:567746/1‑148 (MQ=255)
            tACTATATCGTCAGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCtttt                                                                                                                                         <  4:760722/149‑1 (MQ=255)
            tACTATATCGTCAGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCtttt                                                                                                                                         <  4:186188/149‑1 (MQ=255)
             aCTATATCGTCAGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCtttt                                                                                                                                         >  2:416274/1‑148 (MQ=255)
               tataTCGTCAGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCtttttt                                                                                                                                       <  2:818154/148‑1 (MQ=255)
                atatCGTCAGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGATCG‑CCTTttttc                                                                                                                                      >  2:481568/1‑148 (MQ=255)
                 tatCGTCAGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGTGCTGTCAGTTGACCG‑CCTTttttct                                                                                                                                     <  2:639736/148‑1 (MQ=255)
                 tatCGTCAGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTttttct                                                                                                                                     <  1:194379/148‑1 (MQ=255)
                    cGTCAGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTttttcttttc                                                                                                                                 <  3:453429/149‑1 (MQ=255)
                      tCAGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTACATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCg                                                                                                                                >  1:118960/1‑148 (MQ=255)
                       cAGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGt                                                                                                                               <  2:24824/148‑1 (MQ=255)
                        aGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTa                                                                                                                              <  1:128412/148‑1 (MQ=255)
                        aGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAg                                                                                                                             >  3:136977/1‑149 (MQ=255)
                           aaCCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTttttct                                                                                                                                     <  2:183503/138‑1 (MQ=255)
                           aaCCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAggg                                                                                                                           <  2:4399/148‑1 (MQ=255)
                             ccTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCg                                                                                                                         <  2:948301/148‑1 (MQ=255)
                              cTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCgg                                                                                                                        <  1:128023/148‑1 (MQ=255)
                                  tAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTa                                                                                                                              <  1:766962/138‑1 (MQ=255)
                                  tAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATaa                                                                                                                    <  1:552306/148‑1 (MQ=255)
                                  tAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAg                                                                                                                   <  3:522194/149‑1 (MQ=255)
                                       attattCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATaa                                                                                                                    <  1:374751/143‑1 (MQ=255)
                                       attattCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATaa                                                                                                                    >  2:374751/1‑143 (MQ=255)
                                        ttattCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACcg                                                                                                              >  2:13364/1‑148 (MQ=255)
                                            tCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCa                                                                                                          >  1:678335/1‑148 (MQ=255)
                                             cagcagcagCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCAt                                                                                                         >  1:688915/1‑148 (MQ=255)
                                               gcagcagcTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCg                                                                                                      >  4:118826/1‑149 (MQ=255)
                                                  gcagcTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCca                                                                                                   <  3:129337/149‑1 (MQ=255)
                                                    agcTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCcac                                                                                                  <  1:13364/148‑1 (MQ=255)
                                                    agcTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCAAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCcac                                                                                                  <  2:414240/148‑1 (MQ=255)
                                                      cTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCcacac                                                                                                >  1:344257/1‑148 (MQ=255)
                                                       tGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATaa                                                                                                                    >  4:111514/1‑127 (MQ=255)
                                                            tACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCa                                                                                          <  1:618545/148‑1 (MQ=255)
                                                             aCCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAAc                                                                                        <  3:21284/149‑1 (MQ=255)
                                                                  ggTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGaa                                                                                   <  3:534093/149‑1 (MQ=255)
                                                                  ggTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGaa                                                                                   <  4:149292/149‑1 (MQ=255)
                                                                    tCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACa                                                                                       <  1:851325/143‑1 (MQ=255)
                                                                    tCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACa                                                                                       >  2:851325/1‑143 (MQ=255)
                                                                      tGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACa                                                                                       <  4:253356/141‑1 (MQ=255)
                                                                      tGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACa                                                                                       >  3:253356/1‑141 (MQ=255)
                                                                      tGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAAcatc                                                                               >  3:288506/1‑149 (MQ=255)
                                                                      tGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAAcat                                                                                <  1:527255/148‑1 (MQ=255)
                                                                        gAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAAcatcat                                                                             <  3:121171/149‑1 (MQ=255)
                                                                         aaaaaCGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGt                                                                                                                                                      <  2:841385/76‑1 (MQ=255)
                                                                         aaaaaCGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGt                                                                                                                                                      >  1:841385/1‑76 (MQ=255)
                                                                          aaaaCGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATg                                                                            >  2:70017/1‑148 (MQ=255)
                                                                          aaaaCGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGa                                                                           <  4:166402/149‑1 (MQ=255)
                                                                            aaCGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAg                                                                          <  2:351111/148‑1 (MQ=255)
                                                                               gTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGAt                                                                      >  3:592776/1‑149 (MQ=255)
                                                                                 ggCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTCTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATaa                                                                    <  3:249948/149‑1 (MQ=255)
                                                                                  gCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAAt                                                                   <  4:326240/149‑1 (MQ=255)
                                                                                    cTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGa                                                                 <  4:266164/149‑1 (MQ=255)
                                                                                      gCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGac                                                                <  1:202189/148‑1 (MQ=255)
                                                                                      gCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGac                                                                >  1:82215/1‑148 (MQ=255)
                                                                                              gcgAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTa                                                       <  4:329003/149‑1 (MQ=255)
                                                                                               cgAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGAACAGCGCATCCGCCACACAAAGCGACAGGAACATCATGAGCGATAATGACACTATCGTa                                                       >  2:273872/1‑148 (MQ=255)
                                                                                                      aaaaaaTCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAgg                                                >  1:748668/1‑148 (MQ=255)
                                                                                                      aaaaaaTCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAgg                                                >  1:925712/1‑148 (MQ=255)
                                                                                                      aaaaaaGCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAgg                                                >  2:531036/1‑148 (MQ=255)
                                                                                                         aaaTCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCa                                             >  2:192690/1‑148 (MQ=255)
                                                                                                               tGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTcc                                      >  4:287060/1‑149 (MQ=255)
                                                                                                                gATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTcc                                      <  2:688915/148‑1 (MQ=255)
                                                                                                                gATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTcc                                      >  2:835370/1‑148 (MQ=255)
                                                                                                                      gTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACg                                <  1:914071/148‑1 (MQ=255)
                                                                                                                       tGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACGTg                              >  4:752110/1‑149 (MQ=255)
                                                                                                                         aGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACGTggc                            <  4:59395/149‑1 (MQ=255)
                                                                                                                               cgttAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACGTGGCGGCGtt                      <  3:634976/146‑1 (MQ=255)
                                                                                                                               cGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACGTGGCGGCGt                       <  1:413776/148‑1 (MQ=255)
                                                                                                                               cGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACGTGGCGGCGt                       <  1:436481/148‑1 (MQ=255)
                                                                                                                               cGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACGTGGCGGCGt                       >  1:850744/1‑148 (MQ=255)
                                                                                                                                gCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACGTGGCGGCGTTg                     <  3:259648/149‑1 (MQ=255)
                                                                                                                                 cTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACGTGGCGGCGTTg                     <  1:70017/148‑1 (MQ=255)
                                                                                                                                  tAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACGTGGCGGCGTTGGc                   >  3:306116/1‑149 (MQ=255)
                                                                                                                                   aaaaTAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACGTGGCGGCGTTGGc                   >  2:241480/1‑148 (MQ=255)
                                                                                                                                   aaaaTAAGCGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATGCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACATGGCGGCGTTGGc                   >  2:736976/1‑148 (MQ=255)
                                                                                                                                      aTAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACGTGGCGGCGTTGGCATcc               >  3:724458/1‑149 (MQ=255)
                                                                                                                                       tAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACGTGGCGGCGTTGGCATcc               >  2:224758/1‑148 (MQ=255)
                                                                                                                                        aaaggCGGTCAACTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCg                                     <  1:156895/122‑1 (MQ=255)
                                                                                                                                        aaGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACGTGGCGGCGTTGGCATCCt              >  1:810668/1‑148 (MQ=255)
                                                                                                                                         aGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACGTGGCGGCGTTGGCATCCTgc            <  3:939154/149‑1 (MQ=255)
                                                                                                                                          gggCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAg                                                                          >  4:248188/1‑86 (MQ=255)
                                                                                                                                          gggCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAg                                                                          <  3:248188/86‑1 (MQ=255)
                                                                                                                                          gggCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACGTGGCGGCGTTGGCATCCTgcg           >  3:143734/1‑149 (MQ=255)
                                                                                                                                          gggCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACGTGGCGGCGTTGGCATCCTgcg           >  4:813539/1‑149 (MQ=255)
                                                                                                                                          aggCGGTCAACTGACCGCCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACGTGGCGGCGTTGGCATCCTg             >  2:659690/2‑148 (MQ=255)
                                                                                                                                           ggCGGTCAACTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTAt                                                           >  2:543629/1‑100 (MQ=255)
                                                                                                                                           ggCGGTCAACTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTAt                                                           <  1:543629/100‑1 (MQ=255)
                                                                                                                                              ggTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACGTGGCGGCGTTGGCATCCTGCGCAtc       <  4:307471/149‑1 (MQ=255)
                                                                                                                                                 cAGTTTACCG‑CCTTTTTTCGTTTCGTAGGGCGGATAAGCACCGCGCCTCCGCCACACAAAGCAACATGGACATCCTGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACGTGGCGGCGTTGGCATCCTGCGCAtctc     >  2:376492/1‑148 (MQ=255)
                                                                                                                                                  aGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGTGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACGTGGCGGCGTTGGCATCCTGCGCATCTCCg   >  4:911905/1‑149 (MQ=255)
                                                                                                                                                   gTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCCTCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACGTGGCGGCGTTGGCATCCTGCGCATCTCCgg  <  3:816407/149‑1 (MQ=255)
                                                                                                                                                   gTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACGTGGCGGCGTTGGCATCCTGCGCATCTCCgg  >  4:929631/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GGTCTGGTGCTGTACTATATCGTCAGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCG‑CCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACGTGGCGGCGTTGGCATCCTGCGCATCTCCGG  >  NC_000913/3886609‑3886904

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑