Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,705,971 C→A 20.7% intergenic (‑181/+127) 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,9710CA20.7% 261.4 / 62.9 107intergenic (‑181/+127)dppB/dppAdipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor
Reads supporting (aligned to +/- strand):  ref base C (35/50);  new base A (18/4);  total (53/54)
Fisher's exact test for biased strand distribution p-value = 7.35e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.17e-02

Only 100 of 115 total aligned reads displayed.

GGAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCT  >  NC_000913/3705874‑3706070
                                                                                                 |                                                                                                       
ggAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcgg                                                                                                       >  2:433167/1‑100 (MQ=255)
 gAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcggc                                                                                                      >  1:527612/1‑100 (MQ=255)
  gcataaacaaagcgcactttgtcaataatctgaaaccccggcacatggccggggcttCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATgcggcg                                                                                                     <  1:652168‑M2/43‑1 (MQ=255)
  aaTGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcggcg                                                                                                     <  1:2154949/100‑1 (MQ=255)
     gAAGTATTGCTGCCGGATGGCGCTGCACGCAGGGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTa                                                                                                  <  2:2049417/100‑1 (MQ=255)
      aaGTAGTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTaa                                                                                                 <  2:1462738/100‑1 (MQ=255)
       aGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAAc                                                                                                <  2:2163682/100‑1 (MQ=255)
        gTATTGCTGCCGGATGGCGCTGCACGCAGAGGCCCAGCGCCACCGGGTATACCGATTGTTGACCAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACg                                                                                               >  2:351631/1‑100 (MQ=255)
        gTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACg                                                                                               <  2:2015026/100‑1 (MQ=255)
           ttGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCt                                                                                            <  1:2220446/100‑1 (MQ=255)
                  cGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCgg                                                                                     <  2:2386273/100‑1 (MQ=255)
                  cGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCgg                                                                                     <  2:1965785/100‑1 (MQ=255)
                  cGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCg                                                                                      >  1:2083787/1‑99 (MQ=255)
                   ggtgggCTCTGCACGCAGTGGCACAGCGCGACGGCGTATACAGATTGTTGACAAAGTGCGCTTGGTTTATGCTGGATAAGGCGTTAACGCCCGATCCGGc                                                                                    <  2:1485043/96‑1 (MQ=18)
                   ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGc                                                                                    <  1:1139263/100‑1 (MQ=255)
                   ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGc                                                                                    >  1:2031096/1‑100 (MQ=255)
                   ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGc                                                                                    <  1:1659655/100‑1 (MQ=255)
                   ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGc                                                                                    >  1:533609/1‑100 (MQ=255)
                   ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGccgcatacggc                                                                                    <  2:868053/100‑10 (MQ=37)
                   ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGc                                                                                    <  1:2265040/100‑1 (MQ=37)
                    gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                   <  2:257273/100‑1 (MQ=255)
                    gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                   <  1:1026546/100‑1 (MQ=255)
                    gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                   <  2:2246673/100‑1 (MQ=255)
                    gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGca                                                                                   <  1:1321008/100‑2 (MQ=37)
                     aTGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCt                                                                                  <  2:1049687/100‑1 (MQ=255)
                      tGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCGTATCCGGCCTa                                                                                 <  2:583207/100‑1 (MQ=255)
                        gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑a                                                                              >  2:1289056/1‑100 (MQ=255)
                        gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑a                                                                              >  1:2233800/1‑100 (MQ=255)
                            tGCACGCAGAGGCACAGCGCCACCGGGTATACAGGTTGTTGACAAAGGGTGCTTTGTTTCTGCTGGATGCGGCGTTAACGCCTTAGCCGCCCTAC‑a                                                                              >  2:1420308/1‑96 (MQ=37)
                               aCGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGt                                                                       >  1:1906890/1‑100 (MQ=255)
                                 gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGc                                                                     <  1:1757996/100‑1 (MQ=255)
                                 gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGc                                                                     <  2:1625146/100‑1 (MQ=255)
                                  cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCa                                                                    <  2:95124/100‑1 (MQ=255)
                                  cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCa                                                                    >  2:1598712/1‑100 (MQ=255)
                                   agagGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGc                                                                     <  1:970509/98‑1 (MQ=255)
                                    gagGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCaaa                                                                  >  2:367193/1‑100 (MQ=255)
                                     agGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAAt                                                                 >  2:910870/1‑100 (MQ=255)
                                       gCACAGCGCCACCGGGTATACAGATTGTTGACCAAGTGCGCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTc                                                               <  2:1672339/100‑1 (MQ=255)
                                       gCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTc                                                               <  1:649986/100‑1 (MQ=255)
                                       gCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTc                                                               <  1:531135/100‑1 (MQ=255)
                                        cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCa                                                              >  2:1583913/1‑100 (MQ=255)
                                        cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCa                                                              <  2:2132931/100‑1 (MQ=255)
                                         acaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCaa                                                             <  1:595832/100‑1 (MQ=255)
                                          caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAAt                                                            >  1:1205313/1‑100 (MQ=255)
                                            gcgcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATAc                                                          >  1:791972/1‑100 (MQ=255)
                                               ccACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACAtt                                                       >  2:2019038/1‑100 (MQ=255)
                                                   cGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGc                                                     <  2:2011285/98‑1 (MQ=255)
                                                   cGGGTATACAGATTGTTCACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCa                                                    <  2:2168621/99‑1 (MQ=255)
                                                     ggTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑Gaa                                                >  2:1997370/1‑100 (MQ=255)
                                                        ataCGGATTGTTGACAAAGTGCGCGTTGTTTAAGCTGGATGCGGCGTTAACGCCTTATCGGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑g                                                  <  2:2214437/95‑1 (MQ=14)
                                                         taCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAtttt                                            >  2:1234753/1‑100 (MQ=255)
                                                          aCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAATTTTc                                           >  2:336257/1‑100 (MQ=255)
                                                           cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAATTTTC‑‑g                                        <  1:1540601/100‑1 (MQ=255)
                                                             gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAATTTTC‑‑GTa                                      <  2:359398/100‑1 (MQ=255)
                                                              aTTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAATTTTC‑‑GTAg                                     <  1:913858/100‑1 (MQ=255)
                                                               gcataaaCAAAGCGCACTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAATTTTC‑‑GTAgg                                    >  1:1159376/7‑100 (MQ=17)
                                                               gcataaaCAAAGCGCACTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAATTTTC‑‑GTAgg                                    >  1:224687/7‑100 (MQ=18)
                                                               gcataaaCAAAGCGCACTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAATTTTC‑‑GTAgg                                    >  1:1556835/7‑100 (MQ=18)
                                                                tgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAATTTTC‑‑GTAGGc                                   >  1:1266934/1‑100 (MQ=255)
                                                                tgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAATTTTC‑‑GTAGGc                                   <  2:1541466/100‑1 (MQ=255)
                                                                  ttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAATTTTC‑‑GTAGGcc                                  <  1:1978294/99‑1 (MQ=255)
                                                                   tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAATTTTC‑‑GTAGGCCTg                                >  1:389865/1‑100 (MQ=255)
                                                                   tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAATTTTC‑‑GTAGGCCTg                                <  1:1234760/100‑1 (MQ=255)
                                                                      cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAATTTTC‑‑GTAGGCCTGATa                             <  1:571042/100‑1 (MQ=255)
                                                                       aaatggcgcTTTATATAGGCCGGATGCGGCGTTAACGCCCTATCAGGCCTGC‑AAAAGCGTGCATTTTCACTACATTGCA‑TAATTTTC‑‑GTAGTCCTGATa                             <  2:1548107/94‑1 (MQ=16)
                                                                          gTGCGCTTTGTTTATGCCGGATGCGGCTTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAATTTTC‑‑GTAGGCCTGATAAGCg                         <  1:992402/100‑1 (MQ=255)
                                                                           tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTCC‑AAACTCGTGCAAATTCAATACATTGCA‑GAATTTTC‑‑GTAGGCCTGATAAGCgg                        <  2:1544585/100‑1 (MQ=255)
                                                                           tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAATTTTC‑‑GTAGGCCTGATAAGCgg                        <  2:1472279/100‑1 (MQ=255)
                                                                            gcgcTTTGTTTATGTCGGATGCGGCGTTAACGCCTTATCCGGGCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAATTTTC‑‑GTAGGCCTGATAAGCGGa                       <  2:851541/100‑1 (MQ=255)
                                                                               gattttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTGCA‑CCATCCTC‑‑GTAGGCCTGATAAGCGTAgcg                    >  1:1437372/11‑100 (MQ=18)
                                                                                 ttttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTTTGCAAATTCAATATATTGC‑‑‑AATCTCCGTGTAGGCCTGATAAGCGTAGCGCa                  >  1:1368646/9‑100 (MQ=18)
                                                                                 ttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCa                  <  1:1472132/100‑1 (MQ=255)
                                                                                   gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATc                <  2:1065081/100‑1 (MQ=255)
                                                                                   gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATc                >  2:797622/1‑100 (MQ=255)
                                                                                    tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCa               <  1:1806524/100‑1 (MQ=255)
                                                                                     gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCGTAC‑GAAATCGTGCTAATTCAATATATTGTG‑GGCTCCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAg              >  1:61560/5‑100 (MQ=21)
                                                                                     gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑‑AAATCATGCAAATCCAAAATATTGCA‑GGAGTTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg             >  1:1017640/5‑100 (MQ=21)
                                                                                      taggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGCA‑GAAACCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg             >  1:1923143/4‑100 (MQ=17)
                                                                                      taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAGCATGAAAATTCAATATATTGCA‑GGAGCTGC‑‑GTAGGCCTGATAGGCGTAGCGCATCAgg             >  1:2351496/4‑100 (MQ=18)
                                                                                      taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑‑AAATCATGCAAATCCAAAATATTGCA‑GGAGTTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc            >  1:2073405/4‑100 (MQ=21)
                                                                                       aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAGTTCGTGCAAATTCAATAAATTGCA‑‑ATATGAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGcg           >  1:1094191/3‑99 (MQ=21)
                                                                                       aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAGCATGAAAATTCAATATATTGCA‑GGAGCTGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc            >  2:1311126/3‑100 (MQ=18)
                                                                                       aggCCGGATAAGGCGTTAACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGc            >  1:2300622/3‑100 (MQ=18)
                                                                                       aTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGc            >  1:1154672/1‑100 (MQ=255)
                                                                                        ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCGTGCAAATTCAATACATTGCATGGGCCAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa           >  1:181213/2‑100 (MQ=17)
                                                                                        ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑‑AAATCATGCAAATCCAAAATATTGCA‑GGAGTTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAt          >  1:1788110/2‑100 (MQ=21)
                                                                                        ggCCGGATAAGGCGTTAACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCa           >  2:2323033/2‑100 (MQ=18)
                                                                                         gCCGGATGCGGCGTGAACGCCTTATCCGGCCTAT‑AAAGTCGTATAAATTCAATGTATTGCA‑GTTTATTG‑‑GTAGGCCTGATGCGCTGCGCTTATCAGGcct          <  1:464619/100‑3 (MQ=21)
                                                                                         gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCA‑GGACCAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag          >  1:1884715/1‑99 (MQ=18)
                                                                                         gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATTGTGTAAATTCAATAAATTGAT‑CGATCCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag          >  2:2153720/1‑99 (MQ=18)
                                                                                         gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAGCATGAAAATTCAATATATTGCA‑GGAGCTGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag          >  1:1336598/1‑99 (MQ=18)
                                                                                          ccGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtt         >  1:150657/1‑100 (MQ=255)
                                                                                          ccGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtt         <  1:802835/100‑1 (MQ=255)
                                                                                           cGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAttt        >  1:272095/1‑100 (MQ=255)
                                                                                             gATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAttttt      <  1:1981795/100‑1 (MQ=255)
                                                                                               tGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgc    <  1:1555112/100‑1 (MQ=255)
                                                                                               tGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgc    >  1:1157620/1‑100 (MQ=255)
                                                                                                gcggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcc   >  1:879450/1‑100 (MQ=255)
                                                                                                gcggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcc   >  1:64634/1‑100 (MQ=255)
                                                                                                 cggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct  >  1:451735/1‑100 (MQ=255)
                                                                                                 |                                                                                                       
GGAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCT  >  NC_000913/3705874‑3706070

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