Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,706,048 G→T 25.9% 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,0480GT25.9% 325.6 / 100.9 158intergenic (‑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 (63/50);  new base T (35/3);  total (103/55)
Fisher's exact test for biased strand distribution p-value = 2.40e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.66e-01

Only 100 of 163 total aligned reads displayed.

TGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCAACC  >  NC_000913/3705949‑3706146
                                                                                                         |                                                                                                  
 gcgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑TTC‑‑GTAGGCCTGATCAGCGGa                                                                                                   >  2:1217667/1‑100 (MQ=255)
  cgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑TTC‑‑GTAGGCCTGATAAGCGGAg                                                                                                  >  1:587145/1‑100 (MQ=255)
    cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑TTC‑‑GTAGGCCTGATAAGCGGAgcg                                                                                                <  1:3040692/100‑1 (MQ=255)
    cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑TTC‑‑GTAGGCCTGATAAGCGGAgc                                                                                                 <  1:405163/99‑1 (MQ=255)
     tttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑TTC‑‑GTAGGCCTGATAAGCGGAgcgc                                                                                               >  1:746255/1‑100 (MQ=255)
       tttgtaggCCGGATAAGTCGTTCACTCCGCATCCGGCCTACAAAAGCACGCAAATTCAATATATTACA‑GAGA‑T‑‑TGC‑‑GTAGGCCTGATAAGCGTAg                                                                                                  >  2:1462092/8‑95 (MQ=11)
       tttgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCA‑ATTC‑T‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCAt                                                                                             >  1:2836239/8‑100 (MQ=18)
        ttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACTTTGCAAATTCAATATATTGC‑‑‑AAT‑C‑‑TCCGTGTAGGCCTGATAAGCGTAGCGCATc                                                                                            >  2:2399942/7‑100 (MQ=18)
        gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCTTGCAAATTCAATACATTGCA‑GAAA‑T‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATc                                                                                            >  1:1652348/1‑100 (MQ=255)
        gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GGAT‑T‑‑TTC‑‑GTAGGCCTGATCAGCGGAGCGCATc                                                                                            <  1:2044564/100‑1 (MQ=255)
         tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGAATCa                                                                                           >  1:577411/1‑100 (MQ=255)
         tttATGCCGGATGCGGCGTTAACGACTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑ATC‑‑GTAGGCCTGATAAACGGAGCGCATCa                                                                                           >  2:2268461/1‑100 (MQ=255)
         tgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCA‑TGGG‑C‑‑CAT‑‑GTAGGCCTGATAAGCGTAGCGCATCa                                                                                           >  1:1815373/6‑100 (MQ=17)
          gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCA‑GGAG‑C‑‑TGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAg                                                                                          >  2:3020255/5‑100 (MQ=18)
          gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACTTTGCAAATTCAATATATTGC‑‑‑AAT‑C‑‑TCCGTGTAGGCCTGATAAGCGTAGCGCATCAg                                                                                          <  1:2399970/96‑1 (MQ=18)
           tATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATCCAATACATTGCA‑GAAT‑T‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATAAgg                                                                                         >  2:1289242/1‑100 (MQ=255)
             ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCA‑GAAA‑C‑‑CCT‑‑GTAGGCCGGATAAGCGTAGCGCATCAGGCa                                                                                       >  2:2050239/2‑100 (MQ=17)
             ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGCA‑GAGT‑T‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                                       >  1:85464/2‑100 (MQ=14)
             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCA‑GGGA‑C‑‑TAT‑‑GTAGACCTGATAAGCATAGCGCATCAGGCa                                                                                       >  1:1212855/2‑100 (MQ=18)
             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCA‑GAGA‑C‑‑TAT‑‑GTAGGCCTGATAAGCATAGCGCATCAGGCa                                                                                       >  1:401355/2‑100 (MQ=17)
             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCA‑GGAT‑A‑‑CGT‑‑ATAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                                       >  1:321438/2‑100 (MQ=17)
              gCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt                                                                                      >  1:370551/1‑100 (MQ=255)
              gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACACAAACATGCAAATTCAATAGATTGCA‑GAGA‑T‑‑TAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa                                                                                      >  1:25540/1‑99 (MQ=17)
              gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCA‑‑ATA‑T‑‑GAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGcgat                                                                                     >  1:494082/1‑97 (MQ=21)
               ccggataaggcgTTTACGCCGCATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCA‑‑ATT‑CAACTT‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaa                                                                                      >  1:449408/9‑99 (MQ=21)
               ccggataaggcgTTCACGCCGCATCCGGCCTACAAGATCGTGCAAATTCAATATATTGCA‑TGAG‑C‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagt                                                                                     >  1:146806/9‑98 (MQ=17)
               ccggataaggcgTTCACGCCGCATCCGGCCTACAAAACTTTGCAAATTCAATATATTGC‑‑‑AAT‑C‑‑TCCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCaat                                                                                     >  2:257069/9‑98 (MQ=18)
               ccGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtt                                                                                     >  1:2028258/1‑100 (MQ=255)
               ccGGATCCGGCGTTAACGCCTTATCCAGCCTACAAAATTGTGCAAAATCAATTGATTGCA‑CAGC‑T‑‑AGC‑‑GTAGGCCTGATGCGCTACGCTCATCAGGc                                                                                        <  2:1483766/97‑1 (MQ=18)
                 ggataaggcgTTTACGCCGCATCCGGCCTACAAATTC‑TGCAAATTCAACAAATTGC‑‑GAATCC‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtttgc                                                                                  >  1:2576818/7‑98 (MQ=18)
                 ggataaggcgTTCACGCCGCATCCGGCCTACAAGATCGTGCAAATTCAATATATTGCA‑TGAG‑C‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                                                   >  1:131475/7‑100 (MQ=17)
                 ggataaggcgTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCA‑GTTT‑A‑‑TTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaatct                                                                                   >  2:756655/7‑96 (MQ=25)
                 ggataaggcgTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCA‑GTTT‑A‑‑TTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaatct                                                                                   >  2:1708710/7‑96 (MQ=25)
                 ggataaggcgTTCACGCCGCATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGAC‑T‑‑AT‑‑‑GTAGGCCTGATAAGCATAGCGCATCAGGCagctt                                                                                   >  2:1239319/7‑96 (MQ=21)
                 ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAATCATGCAAATCCAAAATATTGCA‑GGAG‑T‑‑TAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtttat                                                                                  >  1:1066030/7‑98 (MQ=21)
                 ggataaggcgTTCACGCCGCATCCGGCATAC‑AAATCATGCAAATCCAAAATATTGCA‑GGAG‑A‑‑TAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGACAttta                                                                                   >  1:2756742/7‑98 (MQ=18)
                   aTGGGGCGTTATCGCTTTATAGGGCGAAAAAAATCGTGGAAATTCACTACATTGCA‑GAAT‑T‑‑TTC‑‑GGAGGGCTAATAAGTGGAGCGCATCAGGCATTtttg                                                                                 <  2:1974961/100‑1 (MQ=16)
                   aTGCGGCGTTAACGCCTTATCCGGCCTACAAATTCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                                                                                 <  2:89473/100‑1 (MQ=255)
                     gcggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcc                                                                               >  1:895262/1‑100 (MQ=255)
                      cggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct                                                                              <  1:2913962/100‑1 (MQ=255)
                         cgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg                                                                           >  2:636789/1‑100 (MQ=255)
                           ttAACGCCTTATCCGGCCTACAAAATCGTGTAAATTCAATACATTGCA‑GAAT‑T‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg                                                                           <  2:2555109/98‑1 (MQ=255)
                           ttAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgcc                                                                         <  2:475307/100‑1 (MQ=255)
                            tAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgcc                                                                         >  2:1710077/1‑99 (MQ=255)
                            tAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCa                                                                        <  2:1381428/100‑1 (MQ=255)
                              aCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc                                                                      <  1:2420009/100‑1 (MQ=255)
                               cGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCTATCa                                                                     >  1:2333888/1‑100 (MQ=255)
                                gCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAg                                                                    <  1:1116001/100‑1 (MQ=255)
                                 ccTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTTCCTTTGCCATCatt                                                                   >  2:3045457/1‑98 (MQ=255)
                                       ccGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTa                                                             >  1:1980526/1‑100 (MQ=255)
                                       ccGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTa                                                             <  1:2202677/100‑1 (MQ=255)
                                          gCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑TTC‑‑TTAGGCCTGATAAGCGTAGCGCATCATGCATTTTTGCCTTTGCCATCAGTCtt                                                                <  2:2656694/94‑1 (MQ=255)
                                           ccTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGc                                                         <  2:1831134/100‑1 (MQ=255)
                                           ccTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGc                                                         <  2:1626755/100‑1 (MQ=255)
                                               cAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCtaga                                                     >  1:2102142/1‑97 (MQ=255)
                                                   aTCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAatt                                                 >  2:919909/1‑100 (MQ=255)
                                                    tCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAatta                                                <  2:388500/100‑1 (MQ=255)
                                                     cGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAattat                                               >  1:1187900/1‑100 (MQ=255)
                                                      gTGCAAATTCAATACATTGCA‑GAAT‑T‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAattatt                                              >  2:2082608/1‑100 (MQ=255)
                                                      gTGCAAATTCAATACATTGCA‑GAAT‑T‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAattatt                                              >  2:2680400/1‑100 (MQ=255)
                                                       tGCAAATTCAATACATTGCA‑GAAT‑T‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCTCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTc                                             <  1:997315/100‑1 (MQ=255)
                                                       tGCAAATTCAATACATTGCA‑GAAT‑T‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGTCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTc                                             >  1:2595544/1‑100 (MQ=255)
                                                       tGCAAATTCAATACATTGCA‑GAAT‑T‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTc                                             <  2:2148835/100‑1 (MQ=255)
                                                       tGCAAATTCAATACATTGCA‑GAAT‑T‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTc                                             <  1:1549065/100‑1 (MQ=255)
                                                        gCAAATTCAATACATTGCA‑GAAT‑T‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTc                                             >  2:1093481/1‑99 (MQ=255)
                                                        gCAAATTCAATACATTGCA‑GAAT‑T‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCg                                            >  2:1334593/1‑100 (MQ=255)
                                                        gCAAATTCAATACATTGCA‑GAAT‑T‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCg                                            >  2:2894226/1‑100 (MQ=255)
                                                        gCAAATTCAATACATTGCA‑GAAT‑T‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCg                                            <  1:1476580/100‑1 (MQ=255)
                                                          aaaTTCAATACATTGCA‑GAAG‑T‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCAGCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGAt                                          <  2:2741684/100‑1 (MQ=255)
                                                             ttCAATACATTGCA‑GAAT‑T‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATaga                                       <  2:1690613/100‑1 (MQ=255)
                                                              tCAATACATTGCA‑GAAT‑T‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTATATGGCTTTTAATTATTCGATagag                                      >  1:632627/1‑100 (MQ=255)
                                                                 aTACATTGCA‑GAAT‑T‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACg                                   >  2:1666298/1‑100 (MQ=255)
                                                                    cATTGCA‑GAAT‑T‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGttt                                >  1:68886/1‑100 (MQ=255)
                                                                     aTTGCA‑GAAT‑T‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGtttt                               >  2:1113891/1‑100 (MQ=255)
                                                                        gCA‑GAAT‑T‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGa                            <  1:259708/100‑1 (MQ=255)
                                                                        gCA‑GAAT‑T‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGa                            <  1:260192/100‑1 (MQ=255)
                                                                         cA‑GAAT‑T‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGaa                           >  1:276113/1‑100 (MQ=255)
                                                                          a‑GAAT‑T‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAg                          <  2:2333861/100‑1 (MQ=255)
                                                                          a‑GAAT‑T‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAg                          <  2:2391597/100‑1 (MQ=255)
                                                                            gAAT‑T‑‑TTC‑‑GTAGGCCTGATAAGCTGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGt                         <  2:68886/100‑1 (MQ=255)
                                                                            gAAT‑T‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGt                         <  2:2233571/100‑1 (MQ=255)
                                                                              aT‑T‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGa                       >  2:2302553/1‑100 (MQ=255)
                                                                              aT‑T‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGa                       >  1:2223029/1‑100 (MQ=255)
                                                                              aT‑T‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGa                       >  2:1879276/1‑100 (MQ=255)
                                                                                     tC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGtt                   <  2:2521330/100‑1 (MQ=255)
                                                                                           agtccTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCTATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGcc               <  1:429393/97‑1 (MQ=255)
                                                                                           aGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGcc               <  2:950906/100‑1 (MQ=255)
                                                                                           aGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGcc               <  2:701614/100‑1 (MQ=255)
                                                                                           aGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGcc               <  2:468000/100‑1 (MQ=255)
                                                                                            ggCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCt              >  2:1048013/1‑100 (MQ=255)
                                                                                               cTGCTAAGCGCAGCCTATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAAt           <  2:1944670/100‑1 (MQ=39)
                                                                                                tGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAAt           >  2:1132205/1‑99 (MQ=255)
                                                                                                tGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATg          >  1:2870006/1‑100 (MQ=255)
                                                                                                    aaGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATc      <  2:284376/100‑1 (MQ=255)
                                                                                                     aGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCa     >  1:2970680/1‑100 (MQ=255)
                                                                                                     aGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCa     >  2:2172312/1‑100 (MQ=255)
                                                                                                       cGGAGCGCATCAGGCATTTTTGCCTTTTCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCAAc   >  2:3008719/1‑100 (MQ=255)
                                                                                                        ggAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCAAcc  >  2:374413/1‑100 (MQ=255)
                                                                                                        ggAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCAAcc  <  1:1816693/100‑1 (MQ=255)
                                                                                                         gAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCAAcc  <  1:1786123/99‑1 (MQ=255)
                                                                                                         |                                                                                                  
TGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCAACC  >  NC_000913/3705949‑3706146

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