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.

GGAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑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)
                        gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACa                                                                              >  2:1289056/1‑100 (MQ=255)
                        gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACa                                                                              >  1:2233800/1‑100 (MQ=255)
                            tGCACGCAGAGGCACAGCGCCACCGGGTATACAGGTTGTTGACAAAGGGTGCTTTGTTTCTGCTGGATGCGGCGTTAACGCCTTAGCCGCCCTACa                                                                              >  2:1420308/1‑96 (MQ=37)
                               aCGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGt                                                                       >  1:1906890/1‑100 (MQ=255)
                                 gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGc                                                                     <  1:1757996/100‑1 (MQ=255)
                                 gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGc                                                                     <  2:1625146/100‑1 (MQ=255)
                                  cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCa                                                                    <  2:95124/100‑1 (MQ=255)
                                  cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCa                                                                    >  2:1598712/1‑100 (MQ=255)
                                   agagGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGc                                                                     <  1:970509/98‑1 (MQ=255)
                                    gagGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCaaa                                                                  >  2:367193/1‑100 (MQ=255)
                                     agGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAAt                                                                 >  2:910870/1‑100 (MQ=255)
                                       gCACAGCGCCACCGGGTATACAGATTGTTGACCAAGTGCGCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTc                                                               <  2:1672339/100‑1 (MQ=255)
                                       gCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTc                                                               <  1:649986/100‑1 (MQ=255)
                                       gCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTc                                                               <  1:531135/100‑1 (MQ=255)
                                        cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCa                                                              >  2:1583913/1‑100 (MQ=255)
                                        cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCa                                                              <  2:2132931/100‑1 (MQ=255)
                                         acaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCaa                                                             <  1:595832/100‑1 (MQ=255)
                                          caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAAt                                                            >  1:1205313/1‑100 (MQ=255)
                                            gcgcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATAc                                                          >  1:791972/1‑100 (MQ=255)
                                               ccACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACAtt                                                       >  2:2019038/1‑100 (MQ=255)
                                                   cGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGc                                                     <  2:2011285/98‑1 (MQ=255)
                                                   cGGGTATACAGATTGTTCACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCa                                                    <  2:2168621/99‑1 (MQ=255)
                                                     ggTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGaa                                                >  2:1997370/1‑100 (MQ=255)
                                                        ataCGGATTGTTGACAAAGTGCGCGTTGTTTAAGCTGGATGCGGCGTTAACGCCTTATCGGGCCTACAAAATCGTGCAAATTCAATACATTGCAg                                                  <  2:2214437/95‑1 (MQ=14)
                                                         taCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAtttt                                            >  2:1234753/1‑100 (MQ=255)
                                                          aCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTc                                           >  2:336257/1‑100 (MQ=255)
                                                           cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑g                                        <  1:1540601/100‑1 (MQ=255)
                                                             gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTa                                      <  2:359398/100‑1 (MQ=255)
                                                              aTTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAg                                     <  1:913858/100‑1 (MQ=255)
                                                               gcataaaCAAAGCGCACTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAgg                                    >  1:1159376/7‑100 (MQ=17)
                                                               gcataaaCAAAGCGCACTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAgg                                    >  1:224687/7‑100 (MQ=18)
                                                               gcataaaCAAAGCGCACTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAgg                                    >  1:1556835/7‑100 (MQ=18)
                                                                tgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGc                                   >  1:1266934/1‑100 (MQ=255)
                                                                tgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGc                                   <  2:1541466/100‑1 (MQ=255)
                                                                  ttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGcc                                  <  1:1978294/99‑1 (MQ=255)
                                                                   tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTg                                >  1:389865/1‑100 (MQ=255)
                                                                   tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTg                                <  1:1234760/100‑1 (MQ=255)
                                                                      cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATa                             <  1:571042/100‑1 (MQ=255)
                                                                       aaatggcgcTTTATATAGGCCGGATGCGGCGTTAACGCCCTATCAGGCCTGCAAAAGCGTGCATTTTCACTACATTGCATAATTTTC‑‑GTAGTCCTGATa                             <  2:1548107/94‑1 (MQ=16)
                                                                          gTGCGCTTTGTTTATGCCGGATGCGGCTTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCg                         <  1:992402/100‑1 (MQ=255)
                                                                           tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTCCAAACTCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCgg                        <  2:1544585/100‑1 (MQ=255)
                                                                           tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCgg                        <  2:1472279/100‑1 (MQ=255)
                                                                            gcgcTTTGTTTATGTCGGATGCGGCGTTAACGCCTTATCCGGGCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGa                       <  2:851541/100‑1 (MQ=255)
                                                                               gattttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAATCGCTTAAATTCAATATATTGCACCATCCTC‑‑GTAGGCCTGATAAGCGTAgcg                    >  1:1437372/11‑100 (MQ=18)
                                                                                 ttttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACTTTGCAAATTCAATATATTGCAATCTCCGTGTAGGCCTGATAAGCGTAGCGCa                  >  1:1368646/9‑100 (MQ=18)
                                                                                 ttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCa                  <  1:1472132/100‑1 (MQ=255)
                                                                                   gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATc                <  2:1065081/100‑1 (MQ=255)
                                                                                   gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATc                >  2:797622/1‑100 (MQ=255)
                                                                                    tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCa               <  1:1806524/100‑1 (MQ=255)
                                                                                     gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCGTACGAAATCGTGCTAATTCAATATATTGTGGGCTCCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAg              >  1:61560/5‑100 (MQ=21)
                                                                                     gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAATCATGCAAATCCAAAATATTGCAGGAGTTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg             >  1:1017640/5‑100 (MQ=21)
                                                                                      taggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg             >  1:1923143/4‑100 (MQ=17)
                                                                                      taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGC‑‑GTAGGCCTGATAGGCGTAGCGCATCAgg             >  1:2351496/4‑100 (MQ=18)
                                                                                      taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAATCATGCAAATCCAAAATATTGCAGGAGTTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc            >  1:2073405/4‑100 (MQ=21)
                                                                                       aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGcg           >  1:1094191/3‑99 (MQ=21)
                                                                                       aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc            >  2:1311126/3‑100 (MQ=18)
                                                                                       aggCCGGATAAGGCGTTAACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGc            >  1:2300622/3‑100 (MQ=18)
                                                                                       aTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGc            >  1:1154672/1‑100 (MQ=255)
                                                                                        ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCATGGGCCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa           >  1:181213/2‑100 (MQ=17)
                                                                                        ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAATCATGCAAATCCAAAATATTGCAGGAGTTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAt          >  1:1788110/2‑100 (MQ=21)
                                                                                        ggCCGGATAAGGCGTTAACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCa           >  2:2323033/2‑100 (MQ=18)
                                                                                         gCCGGATGCGGCGTGAACGCCTTATCCGGCCTATAAAGTCGTATAAATTCAATGTATTGCAGTTTATTG‑‑GTAGGCCTGATGCGCTGCGCTTATCAGGcct          <  1:464619/100‑3 (MQ=21)
                                                                                         gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCAGGACCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag          >  1:1884715/1‑99 (MQ=18)
                                                                                         gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAATTGTGTAAATTCAATAAATTGATCGATCCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag          >  2:2153720/1‑99 (MQ=18)
                                                                                         gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag          >  1:1336598/1‑99 (MQ=18)
                                                                                          ccGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtt         >  1:150657/1‑100 (MQ=255)
                                                                                          ccGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtt         <  1:802835/100‑1 (MQ=255)
                                                                                           cGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAttt        >  1:272095/1‑100 (MQ=255)
                                                                                             gATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAttttt      <  1:1981795/100‑1 (MQ=255)
                                                                                               tGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgc    <  1:1555112/100‑1 (MQ=255)
                                                                                               tGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgc    >  1:1157620/1‑100 (MQ=255)
                                                                                                gcggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcc   >  1:879450/1‑100 (MQ=255)
                                                                                                gcggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcc   >  1:64634/1‑100 (MQ=255)
                                                                                                 cggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct  >  1:451735/1‑100 (MQ=255)
                                                                                                 |                                                                                                       
GGAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑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: