Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,705,984 T→G 28.0% intergenic (‑194/+114) 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,9840TG28.0% 439.0 / 201.8 218intergenic (‑194/+114)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/83);  new base G (56/5);  total (130/88)
Fisher's exact test for biased strand distribution p-value = 2.04e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.08e-01

Only 100 of 220 total aligned reads displayed.

TTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGT  >  NC_000913/3705885‑3706081
                                                                                                   |                                                                                                     
  gctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTa                                                                                                     >  1:1067817/1‑100 (MQ=255)
  gctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTa                                                                                                     >  2:1069146/1‑100 (MQ=255)
  gctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTa                                                                                                     >  1:1892991/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:3550257/1‑100 (MQ=255)
     gcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATcc                                                                                                  >  1:562767/1‑100 (MQ=255)
     gcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATcc                                                                                                  <  2:1130074/100‑1 (MQ=255)
      ccGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGTTATACAGATTGTTGACAAAGTGCGCTTTGTTTATTCCGGATGCGGCGTTAACGCCTTATCCg                                                                                                 <  1:577404/100‑1 (MQ=255)
       cGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTTCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCgg                                                                                                <  1:2074348/100‑1 (MQ=255)
       cGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATccag                                                                                                <  1:1556217/100‑3 (MQ=255)
        ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGc                                                                                               >  2:2198290/1‑100 (MQ=255)
          aTGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTACCGCCTTATCCGGcc                                                                                              >  2:815701/1‑99 (MQ=255)
           tGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTa                                                                                            <  1:1363576/100‑1 (MQ=255)
           tGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTa                                                                                            >  1:1480729/1‑100 (MQ=255)
           tGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTa                                                                                            >  1:2546238/1‑100 (MQ=255)
           tGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTa                                                                                            <  1:236916/100‑1 (MQ=255)
              cgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACCAAGCGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAcca                                                                                         >  1:2996871/1‑98 (MQ=255)
              cgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGcataaac                                                                                         <  1:319447/100‑7 (MQ=25)
              cgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGAGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑aa                                                                                        >  2:198668/1‑100 (MQ=255)
               gcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGcataaaca                                                                                        <  1:764793/100‑8 (MQ=25)
                cTGCACGCAGAGGCACAGCGCCACCGGGTATACATATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑aaaa                                                                                      <  2:772281/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)
                             cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGATAACGACTTATCCGCCCTAC‑AAAATCGTGCAAATTCa                                                                         >  1:2563651/1‑100 (MQ=255)
                               caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAAt                                                                       >  1:924750/1‑100 (MQ=255)
                               caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAAt                                                                       <  1:2869746/100‑1 (MQ=255)
                               caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAAt                                                                       >  2:2091380/1‑100 (MQ=255)
                               caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAAt                                                                       >  2:1048401/1‑100 (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                               >  1:1297646/1‑100 (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                               <  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:2151747/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                <  1:1032340/100‑1 (MQ=255)
                                                                                   aTGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                <  2:70289/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)
                                                                                              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)
                                                                                                   |                                                                                                     
TTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGT  >  NC_000913/3705885‑3706081

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