Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,705,985 T→C 30.7% intergenic (‑195/+113) 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,705,9850TC30.7% 441.0 / 239.7 287intergenic (‑195/+113)dppB/dppAdipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor
Reads supporting (aligned to +/- strand):  ref base T (90/107);  new base C (77/10);  total (169/118)
Fisher's exact test for biased strand distribution p-value = 2.62e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.66e-03

Only 100 of 291 total aligned reads displayed.

GCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTT  >  NC_000913/3705887‑3706084
                                                                                                  |                                                                                                         
gctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTa                                                                                                          >  1:2410711/1‑100 (MQ=255)
  tgcCGGATGGCGCTGCACGCAGAGGCATAGCGTCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATc                                                                                                        <  1:522855/100‑1 (MQ=255)
  tgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATc                                                                                                        <  1:2067364/100‑1 (MQ=255)
     cGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCgg                                                                                                     <  1:1515046/100‑1 (MQ=255)
     cGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCACATTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCgg                                                                                                     >  1:2094199/1‑100 (MQ=255)
      ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTGATCCGGc                                                                                                    <  2:1118786/100‑1 (MQ=255)
      ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGc                                                                                                    <  2:212884/100‑1 (MQ=37)
        aTGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                                   <  1:3780093/99‑1 (MQ=255)
        aTGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGcat                                                                                                  <  1:1762180/100‑3 (MQ=37)
         tGGCGCTTCACGTTTAGGTACAGCGCTATCGTGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGCAT‑‑a                                                                                               <  1:2349052/100‑1 (MQ=12)
         tGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑a                                                                                               <  2:46708/100‑1 (MQ=255)
           gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGCAT‑‑aaa                                                                                             <  2:2768549/100‑3 (MQ=25)
           gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGCAT‑‑aaa                                                                                             <  1:2142847/100‑3 (MQ=25)
            cgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAATGAGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑aa                                                                                           >  1:992641/1‑100 (MQ=255)
             gcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑aaa                                                                                          >  2:1200600/1‑100 (MQ=255)
                gCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCATTATCCGGCCT‑‑AC‑AAAATc                                                                                       >  2:3172248/1‑100 (MQ=255)
                 cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCg                                                                                      <  1:415612/100‑1 (MQ=255)
                  aCGCAGAGGCACAGGGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATACGGCCT‑‑CC‑AAAATCGt                                                                                     >  2:2862386/1‑100 (MQ=255)
                  aCGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGt                                                                                     >  2:2514592/1‑100 (MQ=255)
                     cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCa                                                                                  >  2:166820/1‑100 (MQ=255)
                      agagGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCaa                                                                                 >  2:3518179/1‑100 (MQ=255)
                         ggCACAGCGCCACCGGGTATACAGATTGTTTACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAAtt                                                                              <  1:1107384/100‑1 (MQ=255)
                          gCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTc                                                                             >  2:39680/1‑100 (MQ=255)
                          gCACAGCGCAACCGGGTATACAGATTGTTTACAAAGTGTGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAAtcc                                                                             <  1:2655158/100‑3 (MQ=18)
                           cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCa                                                                            >  2:1227680/1‑100 (MQ=255)
                             caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAAt                                                                          <  2:694770/100‑1 (MQ=255)
                              aGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATa                                                                         <  2:1793287/100‑1 (MQ=255)
                                 gcCACCGTGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACAt                                                                      <  1:1936597/100‑1 (MQ=255)
                                 gcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACAt                                                                      >  2:3005950/1‑100 (MQ=255)
                                 gcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACAt                                                                      >  1:781340/1‑100 (MQ=255)
                                  ccACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACAtt                                                                     <  1:3307646/100‑1 (MQ=255)
                                   cACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTg                                                                    >  1:2690150/1‑100 (MQ=255)
                                     ccGGTTATACAGATTGTTGACAAATTTCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGc                                                                   <  1:3398284/99‑1 (MQ=255)
                                           ataCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑Ga‑a                                                             >  1:2722163/1‑97 (MQ=255)
                                           ataCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑At‑tt                                                         >  2:1985342/1‑100 (MQ=255)
                                              cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCg                                                      <  2:1826718/100‑1 (MQ=255)
                                                  gcataaaCAAAGCGCACTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAgg                                                  >  2:3770045/7‑100 (MQ=18)
                                                  gcataaaCAAAGCGCACTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAgg                                                  >  2:3812832/7‑100 (MQ=18)
                                                    gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGcc                                                <  1:291382/100‑1 (MQ=255)
                                                    gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AA‑TTTCGTAGGcc                                                <  2:2044688/100‑1 (MQ=255)
                                                     ttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCt                                               >  2:1525194/1‑100 (MQ=255)
                                                        aCAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGAt                                            <  1:3460860/100‑1 (MQ=255)
                                                         cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCtagat                                           >  1:3271586/1‑96 (MQ=255)
                                                           aaGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAg                                         <  1:3469475/100‑1 (MQ=255)
                                                             gTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCg                                       <  1:246085/100‑1 (MQ=255)
                                                              tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCgg                                      <  1:1694237/100‑1 (MQ=255)
                                                                 gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAgc                                   >  2:2706405/1‑100 (MQ=255)
                                                                 gcTTTGTTTATGCCGGATGCGGAGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAgc                                   >  2:854370/1‑100 (MQ=255)
                                                                 gcTTTGTATATACCTGATGCTGCGTTAACACCTTATTCGGCCA‑‑AC‑AAAATAGTGCAAATTCATTATATTGCA‑GA‑AT‑TTTCGTAGTCCTGATCAGCgtata                                   <  2:3336612/100‑5 (MQ=255)
                                                                  cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAgcg                                  <  2:581198/100‑1 (MQ=255)
                                                                  cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAgcg                                  >  1:1099675/1‑100 (MQ=255)
                                                                  cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAgcg                                  >  1:1172899/1‑100 (MQ=255)
                                                                    ttttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACAAAAATCTTGCCAATTCAATATATTGCA‑GG‑AC‑CAT‑GTAGGCCTGATAAGCGTAGCGCa                                >  1:3615220/9‑100 (MQ=18)
                                                                     tGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCAt                               >  2:1263790/1‑100 (MQ=255)
                                                                      gTTTATGCCGGATGCGGCGTTAACGCCTTCTCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATc                              >  2:3658973/1‑100 (MQ=255)
                                                                      gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑ATTCGTAGGCCTGATAAGCGGAGCGCATc                              >  1:2634574/1‑100 (MQ=255)
                                                                       tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCa                             <  2:3100519/100‑1 (MQ=255)
                                                                       tgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑AC‑AAAAGCATGCAAATTCAATATATTGCA‑GA‑GT‑TCACGTAGGCCTGATAAGCGTAGCGCATCa                             >  1:691795/6‑100 (MQ=14)
                                                                        ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATAGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAg                            <  2:2414482/100‑1 (MQ=255)
                                                                        gtaggCCGGATAAGGCGTTAACGCCGCATCCGGCAT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAg                            >  1:3028124/5‑100 (MQ=18)
                                                                         taggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑AC‑AAGATCGTGCAAATTCAATATATTGCA‑TG‑AA‑CCACGTAGGCCTGATAAACGTAGAGCATCaag                           >  1:1678594/4‑98 (MQ=14)
                                                                          aggCCGGATAAGGCGTTTACGCCGCATCTGGCCT‑‑AC‑AAATTCTTCCAAATTCAATATATTGCA‑GA‑AA‑TCATGTAGGCCTGATAAGCGTAGCGCATCAGGc                          >  2:2490923/3‑100 (MQ=18)
                                                                          aggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑AC‑AAAACTGTGCAAATTCAATATATTGCA‑GG‑AA‑ACACGTAGGCCTGATAAGCGAAGCGCATCAGGc                          >  1:225326/3‑100 (MQ=18)
                                                                          aggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑AC‑AAAACCTTGCTAATTCAATATATTGCA‑GG‑GA‑CTATGTAGGCCTGATAAGCATAGCGCATCAGGc                          >  1:3165231/3‑100 (MQ=18)
                                                                          aggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑AC‑‑AAACCGCGCAAATTCAATATATTGCG‑GAGAA‑AAT‑GTAGGCCTGATAAGCGTAGCGCATCAGGct                         >  2:2174068/3‑99 (MQ=18)
                                                                           tGCCGGATGCGGCGTTAACGCCTTATACGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAAGCCTGATAAGCGGAGCGCATCAGGCa                         >  2:290824/1‑100 (MQ=255)
                                                                           ggCCGGATAAGGCGTTCGCGCCGCATCCGGCCT‑‑AC‑AAAAGTTTGCAAATTCAATAAATTGCA‑GA‑AT‑TCATGTAGGCCTGATAAGCGAAGCGCATCAGGCa                         >  2:3572191/2‑100 (MQ=18)
                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACAAAAATCTTGCCAATTCAATATATTGCA‑GA‑A‑‑CCAAGTAGGCCTGATAAGCGTAGCGCATCAGGCa                         >  1:858175/2‑100 (MQ=17)
                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑AC‑AAAATTGTGTAAATTCAATAAATTGAT‑CG‑AT‑CCACGTAGGCCTGATAAGCGTAGCGCATCAGGCa                         >  1:2602120/2‑100 (MQ=21)
                                                                            gCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCAt                        >  1:1455443/1‑100 (MQ=255)
                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACAAAAATCTTGCCAATTCAATATATTGCA‑GG‑AC‑CAT‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag                        >  2:2806682/1‑99 (MQ=18)
                                                                             ccggataaggcgTTCACGCCGCATCCGGCCT‑‑AC‑AAAATCGCTTAAATTCAATATATTGCA‑CA‑AT‑CCTCGTAGGCCTGATAAGCGTAGCGCATCAGGCaat                       >  1:1647455/9‑98 (MQ=18)
                                                                             ccGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCAtt                       >  2:3260231/1‑100 (MQ=255)
                                                                             ccGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCAtt                       >  1:3102731/1‑100 (MQ=255)
                                                                             ccGGATAAGGCGTTTACGCCGCATCCGGCATAAAC‑AAAATCGTGCTAATTCAATATATTGCA‑GA‑AA‑CCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCa                         >  1:1961409/1‑100 (MQ=17)
                                                                              cggataaggcgTTTACGCCGCATCCGGCCT‑‑AC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑GA‑ATCCCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAtttg                     >  1:2950735/8‑99 (MQ=18)
                                                                              cGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCAttt                      >  1:60220/1‑100 (MQ=255)
                                                                               ggataaggcgTTTACGCCGCATCCGGCCT‑‑AC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑GA‑ATCCCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAtttgc                    >  2:2218687/7‑98 (MQ=21)
                                                                               ggataaggcgTTTACGCCGCATCCGGCCT‑‑AC‑AAAAGCGTGCAAATTCAATACATTGCATGG‑GC‑CAT‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAattt                     >  2:2089127/7‑100 (MQ=18)
                                                                               ggataaggcgTTCACGCCGCATCCGGCCT‑‑ACAAAAATCTTGCCAATTCAATATATTGCA‑GG‑AC‑CAT‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                     >  1:1605903/7‑100 (MQ=18)
                                                                               ggataaggcgTTCACGCCGCATCCGGCCT‑‑AC‑GAAATCGTGCTAATTCAATATATTGTG‑GG‑CT‑CCACGTAGGCCTGATAAGCGTAGCGCATCAGGCagt                       >  1:2016312/7‑96 (MQ=18)
                                                                               ggataaggcgTTCACGCCGCATCCGGCCT‑‑AC‑AAAAACATGCAAATTCAATAGATTGCA‑GA‑GA‑TTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAattt                     >  1:2278495/7‑100 (MQ=17)
                                                                               ggataaggcgTTAACGCCGCATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCAtttt                     >  2:534952/7‑100 (MQ=18)
                                                                                gataaggcgTTCACGCCGCATCCGGCCT‑‑AC‑‑AAATCATGCAAATCCAAAATATTGCA‑GG‑AG‑TTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAtttatg                   >  1:1355078/6‑97 (MQ=21)
                                                                                 aTGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                   <  1:1875461/100‑1 (MQ=255)
                                                                                 aTGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                   <  1:1197727/100‑1 (MQ=255)
                                                                                 aTGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                   <  1:2571630/100‑1 (MQ=255)
                                                                                    cggcgTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct                >  2:2382609/1‑100 (MQ=255)
                                                                                    cggcgTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct                <  1:2861154/100‑1 (MQ=255)
                                                                                    cggcgTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcc                 <  2:3319050/99‑1 (MQ=255)
                                                                                      gcgTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgccctt              <  2:1021377/100‑1 (MQ=255)
                                                                                        gTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAGTTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgc            <  1:1265379/100‑1 (MQ=255)
                                                                                           aaCGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCAt         >  2:1431861/1‑100 (MQ=255)
                                                                                            aCGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc        <  1:2239539/100‑1 (MQ=255)
                                                                                             cGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGTA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCa       <  1:46533/100‑1 (MQ=255)
                                                                                             cGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATAAGGCATTTTTGCCTTTGCCATCa       >  1:1015696/1‑100 (MQ=255)
                                                                                              gCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAg      >  1:2389510/1‑100 (MQ=255)
                                                                                                cTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAATc    <  2:261438/100‑1 (MQ=255)
                                                                                                 ttATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCt   <  1:1551482/100‑1 (MQ=255)
                                                                                                  tATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TAGCGTAGGCCTGATAAGCGGAGCGCATCAGGCAGTTTTGCCTTTGCCATCAGTGtt  >  2:946154/1‑100 (MQ=17)
                                                                                                  |                                                                                                         
GCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTT  >  NC_000913/3705887‑3706084

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: ‑