Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,706,048 G→T 37.0% intergenic (‑258/+50) dppB ← / ← dppA dipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,706,0480GT37.0% 489.8 / inf 328intergenic (‑258/+50)dppB/dppAdipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor
Reads supporting (aligned to +/- strand):  ref base G (90/111);  new base T (113/3);  total (214/114)
Fisher's exact test for biased strand distribution p-value = 4.31e-25
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.10e-01

Only 100 of 338 total aligned reads displayed.

TGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCAACCA  >  NC_000913/3705949‑3706147
                                                                                                            |                                                                                                   
tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCgg                                                                                                     <  1:2174010/100‑1 (MQ=255)
   gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAgc                                                                                                  >  1:3912159/1‑100 (MQ=255)
   gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAAGCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAg                                                                                                   <  1:1042068/99‑1 (MQ=255)
       tttgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTGCA‑GA‑GT‑T‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCAt                                                                                              >  1:2618801/8‑100 (MQ=17)
       tttgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTGCA‑GA‑GT‑T‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCAt                                                                                              >  1:2389468/8‑100 (MQ=17)
       tGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCAt                                                                                              >  2:1606354/1‑100 (MQ=255)
       tGTCTATGCCAGAGGGGGCCTTAACGGCGTATCCGGCCTAC‑AAGAGCGTGCAAATTCAATTCGTTGCA‑GA‑AG‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGTgc                                                                                                  <  2:202789/96‑1 (MQ=12)
        ttgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑GA‑AT‑CC‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCa                                                                                            >  1:3675378/7‑100 (MQ=18)
        gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATc                                                                                             >  1:2118781/1‑100 (MQ=255)
         tgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑GA‑ATCC‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAg                                                                                           >  1:1495145/6‑100 (MQ=18)
         tgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTGCA‑GA‑GT‑T‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCa                                                                                            >  2:2963406/6‑100 (MQ=17)
          gaaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAACAGTGCAAATTCAATATATTGCA‑AT‑CT‑T‑‑‑CTG‑‑GTAGGCCTGATAAGCATCGCGCATCAg                                                                                           >  1:2830811/5‑100 (MQ=18)
           tATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAg                                                                                           <  2:1484001/99‑1 (MQ=255)
            aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAACAGTGCAAATTCAATATATTGCA‑AT‑CT‑T‑‑‑CTG‑‑GTAGGCCTGATAAGCATCGCGCATCAGGc                                                                                         >  2:2507734/3‑100 (MQ=18)
            aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTGCA‑CC‑AT‑C‑‑‑CTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                                                                                         >  1:1359021/3‑100 (MQ=21)
            aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAACCTGCAAATTCAATATATTGCA‑GG‑AT‑A‑‑‑CGT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                                                                                         >  1:2197207/3‑100 (MQ=17)
            aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑‑AAATCATGCAAATCCAAAATATTGCA‑GG‑AG‑T‑‑‑TAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                                        >  2:366725/3‑100 (MQ=21)
             tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCa                                                                                        <  2:3006328/100‑1 (MQ=255)
             ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGCA‑GA‑AA‑C‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                                        >  1:3237597/2‑100 (MQ=17)
             ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAACATATTGCA‑AT‑TC‑T‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGct                                                                                        >  1:3112745/2‑99 (MQ=18)
             ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCGTGCAAATTCAATGCATTGCA‑CG‑GG‑C‑‑‑CAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                                        >  2:1137533/2‑100 (MQ=17)
             ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTGCA‑GA‑GT‑T‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                                        >  1:75441/2‑100 (MQ=17)
             ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAATGTGCAAATTCAATAAATTGC‑‑‑A‑AT‑TCAACTT‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGc                                                                                         >  1:3585531/2‑100 (MQ=21)
             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTTTGCAAATTCAATATATTGCA‑AT‑CT‑C‑‑‑CGT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                                        >  2:1888515/2‑100 (MQ=18)
             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTTTGCAAATTCAATATATTGC‑‑‑A‑AT‑C‑‑‑TCCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                                        >  2:3855021/2‑100 (MQ=18)
             ggCCGGATAAGGCGTTCACGCCGCATCCGGCATAC‑AAGTTCGTGCAAATTCAATAAATTGCA‑‑A‑TA‑T‑‑‑GAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGcga                                                                                       >  1:3519790/2‑98 (MQ=21)
             gcccggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAAGCACGCAAACTCAATATATTGCA‑GAGAT‑C‑‑‑AT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                                        >  1:3893438/11‑100 (MQ=17)
              gCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAggaa                                                                                        >  2:2838325/1‑97 (MQ=255)
              gCCGGATAAGGCGTTTACGCCGCATCTGGCCTAC‑AAATTCTTCCAAATTCAATATATTGCA‑GA‑AA‑T‑‑‑CAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa                                                                                       >  1:1816695/1‑99 (MQ=17)
              gCCGGATAAGGCGTTTACGCCGCATCTGGCCTAC‑AAATTCTTCCAAATTCAATATATTGCA‑GA‑AA‑T‑‑‑CAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa                                                                                       >  2:4136408/1‑99 (MQ=18)
              gCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑GA‑AT‑CC‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtt                                                                                      >  1:3365744/1‑100 (MQ=18)
              gCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAACATATTGCA‑‑A‑TTCT‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGctg                                                                                       >  2:1423190/1‑98 (MQ=18)
              gCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAATGTGCAAATTCAATAAATTGC‑‑‑A‑AT‑TCAACTT‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCa                                                                                        >  1:2267381/1‑100 (MQ=21)
               cctgataaggcgTTTACGCCTTATCAGGCCTAC‑AAGATCGTGCAAATTCAACATATTGC‑‑CA‑CT‑C‑‑‑ACCCAGTAGGCCTGATAAGCGCAGCGCATCAGGCaa                                                                                       >  2:2651227/9‑99 (MQ=21)
               ccggataaggcgTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGCA‑GA‑AA‑C‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagt                                                                                      >  1:2446499/9‑98 (MQ=18)
               ccggataaggcgTTTACGCCGCATCCGGCCTAC‑AAAAGCGTGCAAATTCAATACATTGCATGG‑GC‑C‑‑‑AT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaat                                                                                      >  1:1072803/9‑98 (MQ=17)
               ccggataaggcgTTTACGCCGCATCCGGCCTAC‑AAAAATGTGCAAATTCAATAAATTGC‑‑‑A‑AT‑TCAACTT‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaa                                                                                       >  1:92542/9‑99 (MQ=21)
               ccggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTGCA‑CC‑AT‑C‑‑‑CTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaat                                                                                      >  2:914270/9‑98 (MQ=21)
               ccggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAATCATGAAAATTCAATACATTGCA‑AG‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCacg                                                                                      >  1:4028903/9‑98 (MQ=17)
               ccggataaggcgTTCACGCCGCATCCGGCCTAC‑‑AAATCATGCAAATCCAAAATATTGCA‑GG‑AG‑T‑‑‑TAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAttt                                                                                     >  1:2554907/9‑100 (MQ=21)
               ccGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtt                                                                                      >  2:1190945/1‑100 (MQ=255)
                cggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAAGCACGCAAACTCAATATATTGCA‑GAGAT‑C‑‑‑AT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaatt                                                                                     >  2:2109531/8‑97 (MQ=17)
                cGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAttt                                                                                     <  2:3259963/100‑1 (MQ=255)
                 ggataaggcgTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTGCA‑GA‑GT‑T‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                                                    >  1:3278712/7‑100 (MQ=17)
                 ggataaggcgTTTACGCCGCATCCGGCCTAC‑AAAAATGTGCAAATTCAATAAATTGC‑‑‑A‑AT‑TCAACTT‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaatt                                                                                     >  1:3842269/7‑97 (MQ=21)
                 ggataaggcgTTCACGCCTTATCCGGCCTAC‑AAAATTGTGTAAATTCAATAAATTGAT‑CG‑AT‑C‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                                                    >  1:621797/7‑100 (MQ=21)
                 ggataaggcgTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCA‑GG‑AC‑C‑‑‑AT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                                                    >  1:3784957/7‑100 (MQ=18)
                 ggataaggcgTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCA‑GG‑AC‑C‑‑‑AT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                                                    >  1:860604/7‑100 (MQ=18)
                   aTGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                                                                                  <  1:116989/100‑1 (MQ=255)
                   aTGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                                                                                  <  2:396156/100‑1 (MQ=255)
                        gcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgccttt                                                                             >  1:1514490/1‑100 (MQ=255)
                         cgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg                                                                            <  1:2731698/100‑1 (MQ=255)
                         cgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg                                                                            <  2:3091136/100‑1 (MQ=255)
                           ttAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATCCATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgcc                                                                          <  2:2698226/100‑1 (MQ=255)
                             aaCGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCAt                                                                        >  2:3610401/1‑100 (MQ=255)
                              aCGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc                                                                       <  1:448589/100‑1 (MQ=255)
                              aCGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc                                                                       <  2:1332259/100‑1 (MQ=255)
                                  cTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTc                                                                   >  1:131155/1‑100 (MQ=255)
                                       ccGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCaagatcgga                                                              >  1:2638605/1‑92 (MQ=255)
                                       ccGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTa                                                              <  1:672771/100‑1 (MQ=255)
                                       ccGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTa                                                              <  1:1698568/100‑1 (MQ=255)
                                         ggCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATg                                                            <  2:257182/100‑1 (MQ=255)
                                            cTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCt                                                         >  2:3146875/1‑100 (MQ=255)
                                              aC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCttt                                                       <  2:3782099/100‑1 (MQ=255)
                                                 aaaaTCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTa                                                     <  2:4046120/100‑1 (MQ=255)
                                                  aaaTCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTaa                                                    <  1:4120659/100‑1 (MQ=255)
                                                   aaTCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAat                                                   <  1:2975147/100‑1 (MQ=255)
                                                      cGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAattat                                                <  1:3955061/100‑1 (MQ=255)
                                                       gTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAattatt                                               <  2:549819/100‑1 (MQ=255)
                                                       gTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAattatt                                               >  2:1069108/1‑100 (MQ=255)
                                                        tGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTc                                              >  2:3494064/1‑100 (MQ=255)
                                                        tGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTc                                              <  1:1606358/100‑1 (MQ=255)
                                                        tGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTc                                              <  1:732078/100‑1 (MQ=255)
                                                         gCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCg                                             <  2:3975158/100‑1 (MQ=255)
                                                            aaTTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATa                                          <  1:1903744/100‑1 (MQ=255)
                                                              ttCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATaga                                        >  2:829143/1‑100 (MQ=255)
                                                                       ttGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTGTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTc                               >  1:2940773/1‑100 (MQ=255)
                                                                       ttGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTc                               >  2:223491/1‑100 (MQ=255)
                                                                        tGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGtttt                                <  2:3441592/98‑1 (MQ=255)
                                                                         gCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGTATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGa                             <  2:1772733/100‑1 (MQ=255)
                                                                         gCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGa                             <  2:1719885/100‑1 (MQ=255)
                                                                         gCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGa                             >  1:3726310/1‑100 (MQ=255)
                                                                         gCA‑GA‑AT‑G‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTGGCCATCAGTCTTGTATGGCGTTTAATTATTCGATAGAGACGTTTTCGa                             <  2:2771372/100‑1 (MQ=255)
                                                                           a‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAg                           <  1:506221/100‑1 (MQ=255)
                                                                             gA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGt                          <  2:4184273/100‑1 (MQ=255)
                                                                                aT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCCGAGCACATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTg                         <  2:505188/99‑1 (MQ=255)
                                                                                             tAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGc                 <  1:3610598/100‑1 (MQ=255)
                                                                                              aGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGcc                <  2:1681876/100‑1 (MQ=255)
                                                                                               ggCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCt               <  1:4058084/100‑1 (MQ=255)
                                                                                                gCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTa              <  2:2452789/100‑1 (MQ=255)
                                                                                                   tGATAAGCGGAGCGCATCAGGCATTTTTTTTTTTTTTTTCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATg           >  2:1278312/1‑100 (MQ=255)
                                                                                                   tGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATg           <  1:490018/100‑1 (MQ=255)
                                                                                                    gATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCCGTCTTGTATGGCTTTTAATTATTCCATAGAGACGTTTTCGAAGTGATGTTTGCGTAATgg          >  2:1359037/1‑100 (MQ=255)
                                                                                                     aTAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGa         <  1:3417110/100‑1 (MQ=255)
                                                                                                       aaGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATc       <  2:664561/100‑1 (MQ=255)
                                                                                                        aGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGAt        >  2:3758888/1‑98 (MQ=255)
                                                                                                        aGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCa      <  2:2165512/100‑1 (MQ=255)
                                                                                                         gCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCaa     >  1:1206760/1‑100 (MQ=255)
                                                                                                         gCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCaa     <  2:3369003/100‑1 (MQ=255)
                                                                                                            gAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCAACca  <  1:1937218/100‑1 (MQ=255)
                                                                                                            |                                                                                                   
TGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCAACCA  >  NC_000913/3705949‑3706147

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑