Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,705,985 T→C 27.6% 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,9850TC27.6% 391.8 / 167.2 218intergenic (‑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 (74/84);  new base C (56/4);  total (130/88)
Fisher's exact test for biased strand distribution p-value = 4.13e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.74e-02

Only 100 of 220 total aligned reads displayed.

GCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTT  >  NC_000913/3705887‑3706084
                                                                                                  |                                                                                                       
gctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTa                                                                                                        >  1:1067817/1‑100 (MQ=255)
gctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTa                                                                                                        >  1:1892991/1‑100 (MQ=255)
gctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTa                                                                                                        >  2:1069146/1‑100 (MQ=255)
  tgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATc                                                                                                      >  2:2914380/1‑100 (MQ=255)
  tgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATc                                                                                                      >  1:3725059/1‑100 (MQ=255)
   gcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAGCGCCTTATcc                                                                                                     >  1:889605/1‑100 (MQ=255)
   gcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATcc                                                                                                     <  2:1130074/100‑1 (MQ=255)
   gcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATcc                                                                                                     >  2:3550257/1‑100 (MQ=255)
    ccGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGTTATACAGATTGTTGACAAAGTGCGCTTTGTTTATTCCGGATGCGGCGTTAACGCCTTATCCg                                                                                                    <  1:577404/100‑1 (MQ=255)
     cGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTTCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCgg                                                                                                   <  1:2074348/100‑1 (MQ=255)
     cGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATccag                                                                                                   <  1:1556217/100‑3 (MQ=255)
        aTGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTACCGCCTTATCCGGcc                                                                                                 >  2:815701/1‑99 (MQ=255)
         tGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTa                                                                                               >  1:2546238/1‑100 (MQ=255)
         tGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTa                                                                                               >  1:1480729/1‑100 (MQ=255)
         tGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTa                                                                                               <  1:236916/100‑1 (MQ=255)
         tGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTa                                                                                               <  1:1363576/100‑1 (MQ=255)
            cgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACCAAGCGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAcca                                                                                            >  1:2996871/1‑98 (MQ=255)
            cgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGcataaac                                                                                            <  1:319447/100‑7 (MQ=25)
             gcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGcataaaca                                                                                           <  1:764793/100‑8 (MQ=25)
              cTGCACGCAGAGGCACAGCGCCACCGGGTATACATATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑aaaa                                                                                         <  2:772281/100‑1 (MQ=255)
              cTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑aaaa                                                                                         <  2:1986667/100‑1 (MQ=255)
                gCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGGCTTATCCGGCCTAC‑AAAATc                                                                                       >  1:528286/1‑100 (MQ=255)
                gCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATc                                                                                       <  2:615006/100‑1 (MQ=255)
                 cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCg                                                                                      <  1:5966/100‑1 (MQ=255)
                 cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCg                                                                                      <  2:1395981/100‑1 (MQ=255)
                  aCGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGt                                                                                     <  2:3055300/100‑1 (MQ=255)
                     cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGGCTTATCCGGCCTAC‑AAAATCGTGCa                                                                                  >  1:452354/1‑100 (MQ=255)
                       gagGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCaaa                                                                                >  2:2328204/1‑100 (MQ=255)
                         ggCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAAtt                                                                              >  2:3724192/1‑100 (MQ=255)
                           cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCa                                                                            >  2:2164254/1‑100 (MQ=255)
                           cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGATAACGACTTATCCGCCCTAC‑AAAATCGTGCAAATTCa                                                                            >  1:2563651/1‑100 (MQ=255)
                             caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAAt                                                                          >  2:1048401/1‑100 (MQ=255)
                             caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAAt                                                                          <  2:2123534/100‑1 (MQ=255)
                             caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAAt                                                                          >  1:924750/1‑100 (MQ=255)
                             caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAAt                                                                          <  1:2869746/100‑1 (MQ=255)
                              aGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAAtn                                                                         <  1:1120755/100‑2 (MQ=255)
                                   cACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATTCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACAtt                                                                     <  1:2665355/99‑1 (MQ=255)
                                   cACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTg                                                                    >  2:1791100/1‑100 (MQ=255)
                                     ccGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCa                                                                  <  1:1725242/100‑1 (MQ=255)
                                      cGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACAtt                                                                     <  1:2555084/96‑1 (MQ=255)
                                      cGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑g                                                                <  1:329690/100‑1 (MQ=255)
                                       gggTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑Ga                                                               >  2:3093577/1‑100 (MQ=255)
                                              cAGATTGTTGACAAAGTGCACTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑TTCg                                                      <  1:2145348/100‑1 (MQ=255)
                                                  ttgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑TTCGTAgg                                                  <  1:954639/100‑1 (MQ=255)
                                                  gcataaaCAAAGCGCACTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑TTCGTAgg                                                  >  1:3457812/7‑100 (MQ=17)
                                                   tgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑TTCGTAGGc                                                 <  2:1654429/100‑1 (MQ=255)
                                                    gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑TTCGTAGGcc                                                >  1:1369906/1‑100 (MQ=255)
                                                     ttgACAAATGGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑TTCGTAGGCCt                                               <  1:1179216/100‑1 (MQ=255)
                                                      tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑TTCGTAGGCCTg                                              <  2:1226554/100‑1 (MQ=255)
                                                      tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑TTCGTAGGCCTg                                              >  2:506837/1‑100 (MQ=255)
                                                       gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑TTCGTAGGCCTGa                                             >  1:1846862/1‑100 (MQ=255)
                                                         cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTTCA‑GA‑ATT‑TTCGTAGGCCTGATa                                           <  1:2203444/100‑1 (MQ=255)
                                                         cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑TTCGTAGGCCTGATa                                           <  2:262355/100‑1 (MQ=255)
                                                         cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑TTCGTAGGCCTGATa                                           <  1:180243/100‑1 (MQ=255)
                                                           aaGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑TTCGTAGGCCTGATAAg                                         >  1:3255959/1‑100 (MQ=255)
                                                            aGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑TTCGTAGGCCTGATAAGc                                        >  2:989331/1‑100 (MQ=255)
                                                               gcgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑TTCGTAGGCCTGATAAGCGGa                                     >  1:1695818/1‑100 (MQ=255)
                                                                 gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑TTCGTAGGCCTGATAAGCGGAgc                                   <  2:1180077/100‑1 (MQ=255)
                                                                 gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑TTCGTAGGCCTGATAAGCGGAgc                                   >  1:1953889/1‑100 (MQ=255)
                                                                 gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑TTCGTAGGCCTGATAAGCGGAgc                                   <  1:1224607/100‑1 (MQ=255)
                                                                  cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑TTCGTAGGCCTGATAAGCGGAgcg                                  >  2:3353863/1‑100 (MQ=255)
                                                                  cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑TTCGTAGGCCTGATAAGCGGAgcg                                  <  1:3030356/100‑1 (MQ=255)
                                                                  cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑TTCGTAGGCCTGATAAGCGGAgcg                                  >  1:1297646/1‑100 (MQ=255)
                                                                  cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑TTCGTAGGCCTGATAAGCGGAgcg                                  <  2:1907433/100‑1 (MQ=255)
                                                                    ttttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAACCTGCAAATTCAATATATTGCA‑GG‑ATA‑CGTGTAGGCCTGATAAGCGTAGCGCa                                >  2:882405/9‑100 (MQ=17)
                                                                     tGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑TTCGTAGGCCTGATAAGCGGAGCGCAt                               <  2:2250984/100‑1 (MQ=255)
                                                                      ttgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGCA‑GA‑AAC‑CTTGTAGGCCTGATAAGCGTAGCGCATc                              >  1:1135662/7‑100 (MQ=18)
                                                                       tttATGCCGGACGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑TTCGTAGGCCTGATAAGCGGAGCGCATCa                             <  2:3401971/100‑1 (MQ=255)
                                                                        gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAGATCGTGCAAATTCAATATATTGCA‑TG‑AGC‑CACGTAGGCCTGATAAGCGTAGCGCATCAg                            >  2:1549921/5‑100 (MQ=17)
                                                                        gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTGTGCAAATTCAATATATTGCA‑GG‑AAA‑CACGTAGGCCTGATAAGCGAAGCGCATCAg                            >  1:2117020/5‑100 (MQ=18)
                                                                        gtaggCCGGATAAGGCGTTAACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑TTCGTAGGCCTGATAAGCGGAGCGCATCAg                            >  1:3153454/5‑100 (MQ=17)
                                                                          aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCGTGCAAATTCAATACATTGCA‑TG‑GGC‑CATGTAGGCCTGATAAGCGTAGCGCATCAGGc                          >  1:3201060/3‑100 (MQ=17)
                                                                          aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTTTGCAAATTCAATATATTGCA‑AT‑CTC‑CGTGTAGGCCTGATAAGCGTAGCGCATCAGGc                          >  1:722526/3‑100 (MQ=18)
                                                                          aggCCGGATAAGGCGTTAACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGc                          >  2:618701/3‑100 (MQ=18)
                                                                           tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑TTCGTAGGCATGATAAGCGGAGCGCATCAGGc                          >  2:2834524/1‑99 (MQ=255)
                                                                           ggCCGGATAAGGCGTTTACGCCGCATCTGGCCTAC‑AAATTCTTCCAAATTCAATATATTGCA‑GA‑AAT‑CATGTAGGCCTGATAAGCGTAGCGCATCAGGCa                         >  2:1260548/2‑100 (MQ=17)
                                                                           ggCCGGATAAGGCGTTAACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCa                         >  2:3661315/2‑100 (MQ=18)
                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAGCACGCAAACTCAATATATTGCA‑GAGATC‑AT‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa                        >  2:3302709/1‑99 (MQ=17)
                                                                             ccggataaggcgTTCACGCCGCATCCGGCATAC‑AAAACTGTGCAAATTCAATATATTGCA‑GG‑AAA‑CACGTAGGCCTGATAAGCGAAGCGCATCAGGCagt                       >  2:1862509/9‑98 (MQ=18)
                                                                              cggctaaggcgTTTACGCCGCAACCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑GA‑ATCCCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAtttg                     >  1:1918263/8‑99 (MQ=18)
                                                                               ggataaggcgTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCA‑GG‑ACC‑AT‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                     >  2:246100/7‑100 (MQ=18)
                                                                               ggataaggcgTTCACGCCGCATCCGACCCAC‑AAAGTCTTGCAAAATCAATAAATTGCAGGA‑ACG‑AT‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaacta                     >  2:3147289/7‑96 (MQ=25)
                                                                               ggataaggcgTTCACGCCGCATCCAGCCTAC‑AAAATTGTGCAAAATCAATGGATTGCA‑CA‑GCT‑AGCGTAGGCCTGATAAGCGTAGCGCATCAGGCagatt                     >  1:2434254/7‑96 (MQ=25)
                                                                                 aTGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                   <  2:70289/100‑1 (MQ=255)
                                                                                 aTGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                   <  1:1032340/100‑1 (MQ=255)
                                                                                  taaggcgTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTGCA‑GA‑GTT‑CACGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTtttgc                  >  2:2685643/4‑100 (MQ=17)
                                                                                   gcggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcc                 <  2:3347420/100‑1 (MQ=255)
                                                                                    cggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct                >  2:151445/1‑100 (MQ=255)
                                                                                     ggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTACAAATTCAATACATTGCA‑GA‑ATT‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctt               >  1:2504716/1‑100 (MQ=255)
                                                                                      gcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgccttt              <  2:1597117/100‑1 (MQ=255)
                                                                                       cgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg             >  2:3285306/1‑100 (MQ=255)
                                                                                        gTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑TTTGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgc            <  1:1097959/100‑1 (MQ=255)
                                                                                        gTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgc            <  1:541019/100‑1 (MQ=255)
                                                                                         ttAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgcc           <  2:2016149/100‑1 (MQ=255)
                                                                                            aCGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc        <  2:3495213/100‑1 (MQ=255)
                                                                                             cGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc        >  1:3679804/1‑99 (MQ=255)
                                                                                              gCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑GTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAg      >  1:3790582/1‑100 (MQ=255)
                                                                                               ccTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGt     <  1:2790748/100‑1 (MQ=255)
                                                                                               ccTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGt     >  1:3615142/1‑100 (MQ=255)
                                                                                               ccTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATC‑TTCGTGGCCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGt     >  1:164981/1‑100 (MQ=255)
                                                                                                  |                                                                                                       
GCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTT  >  NC_000913/3705887‑3706084

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