Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,705,985 T→C 30.7% intergenic (‑195/+113) dppB ← / ← dppA dipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,705,9850TC30.7% 394.1 / 211.9 239intergenic (‑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 (78/87);  new base C (68/5);  total (146/93)
Fisher's exact test for biased strand distribution p-value = 8.77e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.71e-02

Only 100 of 243 total aligned reads displayed.

GCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCT  >  NC_000913/3705887‑3706083
                                                                                                  |                                                                                                        
gctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTa                                                                                                         <  1:1195709/100‑1 (MQ=255)
 ctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTAt                                                                                                        <  2:1429576/100‑1 (MQ=255)
     cGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCgg                                                                                                    >  1:433704/1‑100 (MQ=255)
       gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                                  <  1:3785481/100‑1 (MQ=255)
       gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                                  <  1:283033/100‑1 (MQ=255)
          ggCGCTGTACGCAGAGGCACAGCGCAAGCGATTATACAGATTGTGGAAAAAGTGCGCGTTGTTTATGCCTTATAAGGCGTTAACGCCGCATCCGGca                                                                                                  <  1:903820/97‑2 (MQ=16)
           gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACa                                                                                              >  1:2980986/1‑100 (MQ=255)
                 cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAGAAAATCg                                                                                        <  2:2023992/100‑1 (MQ=255)
                 cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCg                                                                                        >  1:1866644/1‑100 (MQ=255)
                  aCGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGt                                                                                       >  2:588420/1‑100 (MQ=255)
                   cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTg                                                                                      <  2:1831757/100‑1 (MQ=255)
                   cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTg                                                                                      >  2:452989/1‑100 (MQ=255)
                   cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCAGGCCTACAAAATCGTg                                                                                      >  2:292574/1‑100 (MQ=255)
                    gatctggCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGc                                                                                     <  1:100918/95‑1 (MQ=255)
                      agagGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCaa                                                                                   >  1:1928904/1‑100 (MQ=255)
                       gagGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCaaa                                                                                  >  1:3690706/1‑100 (MQ=255)
                       gagGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCaaa                                                                                  >  2:73991/1‑100 (MQ=255)
                        agGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAAt                                                                                 <  2:428649/100‑1 (MQ=255)
                         ggCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCagatc                                                                                >  2:100920/1‑96 (MQ=255)
                          gCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTc                                                                               <  1:2172662/100‑1 (MQ=255)
                          gCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGGTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTc                                                                               <  1:696068/100‑1 (MQ=255)
                             caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCTGGCCTACAAAATCGTGCAAATTCAAt                                                                            <  2:2348666/100‑1 (MQ=255)
                             caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAAt                                                                            <  2:3623243/100‑1 (MQ=255)
                                    aCCGGGGATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTAGCCGGCCTACAAAATCGTGCAAATTCAATACATTGc                                                                     <  2:1704188/100‑1 (MQ=255)
                                     ccGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCa                                                                    <  2:819002/100‑1 (MQ=255)
                                      cGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAg                                                                  <  2:689259/100‑1 (MQ=255)
                                      cGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAg                                                                  >  2:1595785/1‑100 (MQ=255)
                                           ataCAGATTGTTGCACAAGTGCTCTGTGTTTATTCCGTATGCGGCGTTTACGCCGTATCCTGCCTACAAAATCGTGAAAATTCAATGCATTGCAG‑‑‑a                                                              <  1:3735098/96‑1 (MQ=12)
                                             aCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTc                                                      >  1:1134631/1‑100 (MQ=255)
                                               aGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTCGt                                                    >  1:238974/1‑100 (MQ=255)
                                                 aTTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTCGTAg                                                  <  1:2145517/100‑1 (MQ=255)
                                                 aTTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTCGTAg                                                  <  1:604568/100‑1 (MQ=255)
                                                   tgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTCGTAgcc                                                >  1:2342165/1‑98 (MQ=255)
                                                   cataaaCAAAGCGCACTTTGGTTATGCAGCATGCGGCGCTAACACCGTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTCGTa                                                   >  1:2492197/6‑97 (MQ=17)
                                                      tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTCGTAGGCCTg                                             >  2:1555601/1‑100 (MQ=255)
                                                       gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTCGTAGGCCTGa                                            >  2:765150/1‑100 (MQ=255)
                                                        aCAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTCGTAGGCCTGAt                                           >  2:1519232/1‑100 (MQ=255)
                                                        aCAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTACAG‑‑‑AATTTTCGTAGGCCTGAt                                           <  2:1123638/100‑1 (MQ=255)
                                                          aaaGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTCGTAGGCCTGATaa                                         <  2:1208219/100‑1 (MQ=255)
                                                             gTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTCGTAGGCCTGATAAGCg                                      <  1:625695/100‑1 (MQ=255)
                                                                 gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTCGTAGGCCTGATAAGCGGAgc                                  <  1:1257887/100‑1 (MQ=255)
                                                                  cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTCGTAGGCCTGATAAGCGGAgcg                                 >  1:3477021/1‑100 (MQ=255)
                                                                  cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTCGTAGGCCTGATAAGCGGAgcg                                 >  2:3190479/1‑100 (MQ=255)
                                                                     tttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAG‑‑‑GGACTATGTAGGCCTGATAAGCATAGCGCAt                              >  2:2718370/8‑100 (MQ=18)
                                                                     tttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAG‑‑‑AGATTATGTAGGCCTGATAAGCGTAGCGCAt                              >  2:297290/8‑100 (MQ=17)
                                                                      ttgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGCAG‑‑‑AGTTCACGTAGGCCTGATAAGCGTAGCGCATc                             >  1:1985931/7‑100 (MQ=17)
                                                                       tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTCGTAGGCCTGATAAGCGGAGCGCATCa                            <  2:2686460/100‑1 (MQ=255)
                                                                       tctaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAATCATGCAAATCCAAAATATTGCAG‑‑‑GAGTTATGTAGGCCTGATAAGCGTAGCGCATCAg                           <  1:2876362/95‑1 (MQ=21)
                                                                        ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTCGTAGGCCTGATAAGCGGAGCGCATCAg                           >  2:3411144/1‑100 (MQ=255)
                                                                        ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTCGTAGGCCTGATAAGCGGAGCGCATCAg                           <  1:3191798/100‑1 (MQ=255)
                                                                        gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCATACAAAAACCTGCAAATTCAATATATTGCAG‑‑‑GATACGTGTAGGCCTGATAAGCGTAGCGCATCAg                           >  2:3044992/5‑100 (MQ=18)
                                                                         taggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAG‑‑‑AAACCTTGTAGGCCTGATAAGCGTAGCGCATCAgg                          >  2:3749263/4‑100 (MQ=18)
                                                                          aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAATTCTGCAAATTCAACAAATTGCG‑‑‑AATCCCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCa                        >  1:3093945/3‑100 (MQ=18)
                                                                          aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTACAG‑‑‑AAATCTTGTAGGCCTGATAAGCGTAGCGCATCACGc                         >  1:2711244/3‑100 (MQ=14)
                                                                          aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCA‑‑‑ATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGcg                        >  2:3844775/3‑99 (MQ=18)
                                                                          aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAATCATGCAAATCCAAAATATTGCAG‑‑‑GAGTTATGTAGGCCTGATAAGCGTAGCGCATCAgaga                        >  2:2876384/3‑97 (MQ=21)
                                                                          aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAATCATGCAAATCCAAAATATTGCAG‑‑‑GAGTTATGTAGGCCTGATAAGCGTAGCGCATCAGGCa                        >  2:2378803/3‑100 (MQ=21)
                                                                          aggCCGGATAAGGCGTTCACGCCGCATCCGGCATACAAAATTGTGTAAATTCAATAAATTGATC‑‑‑GATCCACGTAGGCCTGATAAGCGTAGCGCATCAGGc                         >  2:3519958/3‑100 (MQ=21)
                                                                           ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAATTCTGCAAATTCAACAAATTGCG‑‑‑AATCCCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAt                       >  2:3805789/2‑100 (MQ=18)
                                                                           ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAG‑‑‑AAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCa                        >  2:3286393/2‑100 (MQ=17)
                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAG‑‑‑GAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCa                        >  1:951293/2‑100 (MQ=17)
                                                                            gCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAATTCTGCAAATTCAACAAATTGCG‑‑‑AATCCCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAtt                      >  2:1277516/1‑100 (MQ=18)
                                                                             ccggataaggcgTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAG‑‑‑TTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCaat                      >  2:2003096/9‑98 (MQ=25)
                                                                             ccggataaggcgTTCACGCCGCATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAA‑‑‑GATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCacg                      >  1:3401505/9‑98 (MQ=17)
                                                                             ccggataaggcgTTCACGCCGCATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAG‑‑‑GAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCagt                      >  2:3218106/9‑98 (MQ=17)
                                                                             ccggataaggcgTTCACGCCGCATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAG‑‑‑GAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCagt                      >  2:2641618/9‑98 (MQ=18)
                                                                             ccggataaggcgTTCACGCCGCATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAG‑‑‑GGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCagc                      >  1:2696896/9‑98 (MQ=18)
                                                                             ccggataaggcgTTCACGCCGCATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAG‑‑‑AGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCaat                      >  2:3697232/9‑98 (MQ=17)
                                                                             ccGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCAtt                      >  2:1906813/1‑100 (MQ=255)
                                                                             ccGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCAtt                      >  1:3546130/1‑100 (MQ=255)
                                                                               ngataaggcgTTCACGCCGCATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGG‑‑‑GACTATGTAGGCCTGATAAGCATAGCGCATCAGGCagctt                    >  2:3145180/7‑96 (MQ=18)
                                                                               ggataaggcgTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGCAG‑‑‑AGTTCACGTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                    >  2:3307601/7‑100 (MQ=17)
                                                                               ggataaggcgTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGCAG‑‑‑AGTTCACGTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                    >  2:934387/7‑100 (MQ=17)
                                                                               ggataaggcgTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGCAG‑‑‑AGTTCACGTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                    >  2:1673375/7‑100 (MQ=17)
                                                                               ggataaggcgTTTACGCCGCATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTT‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaatt                     >  2:929691/7‑97 (MQ=21)
                                                                               ggataaggcgTTCACGCCGCATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAA‑‑‑GATTTTCGTAGGCCTGATAAGCGTAGCGAATCAGGCacgc                     >  2:2779062/7‑96 (MQ=17)
                                                                               ggataaggcgTTCACGCCGCATCCGGCCTACAAATCATGCAAATCCAAAATATTGCAG‑‑‑GAGTTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAtttat                   >  2:2166157/7‑98 (MQ=21)
                                                                               ggataaggcgTTCACGCCGCATCCGACCCACAAAGTCTTGCAAAATCAATAAATTGCAG‑‑‑GAACGATGTAGGCCTGATAAGCGTAGCGCATCAGGCaacta                    >  1:1013021/7‑96 (MQ=21)
                                                                               ggATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCAtttt                    <  2:1220856/100‑1 (MQ=255)
                                                                               ggATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCAtttt                    >  1:1435319/1‑100 (MQ=255)
                                                                                gataaggcgTTCACGCCGCATCCGGCCTACAAATCCTTGCGAATTCAATATATTGTATG‑‑‑GAAATTC‑‑AGGCCTGATAAGCGTAGCGCATCAGGCttttcac                  >  1:2805802/6‑97 (MQ=37)
                                                                                gataaggcgTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAG‑‑‑TTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCaatctt                   >  2:637052/6‑95 (MQ=25)
                                                                                gataaggcgTTCACGCCGCATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAG‑‑‑GAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAgtttt                   >  2:436250/6‑100 (MQ=18)
                                                                                 aTGCGGCGTTAACGCTTTATCCGGCCTACAAAATCGTGCAAATCCAATACATTGCAG‑‑‑AATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                  <  2:491520/100‑1 (MQ=255)
                                                                                   gcggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTCGTAGGCCTGATACGCGGAGCGCATCAGCCAttttttc                 >  2:663766/1‑97 (MQ=255)
                                                                                    cggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCAttttttttt               >  1:6669/1‑96 (MQ=255)
                                                                                    cggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct               <  2:2436025/100‑1 (MQ=255)
                                                                                    cggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct               >  2:178452/1‑100 (MQ=255)
                                                                                       cgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg            <  1:1461450/100‑1 (MQ=255)
                                                                                       cgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg            <  1:1056723/100‑1 (MQ=255)
                                                                                        gTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgc           <  1:314861/100‑1 (MQ=255)
                                                                                         cttccgatctATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgcc          <  1:2592890/91‑1 (MQ=255)
                                                                                          tAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCa         >  2:2065443/1‑100 (MQ=255)
                                                                                          tAACGCCTTATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCa         <  1:1558964/100‑1 (MQ=255)
                                                                                            aCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc       <  2:3528349/100‑1 (MQ=255)
                                                                                             cGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGTTTTTTCTATCa      >  1:1495488/1‑100 (MQ=17)
                                                                                             cGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc       <  1:515611/99‑1 (MQ=255)
                                                                                              gCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAg     >  1:202999/1‑100 (MQ=255)
                                                                                              gCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAg     <  1:1743347/100‑1 (MQ=255)
                                                                                               ncTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGt    >  1:823598/2‑100 (MQ=255)
                                                                                                  |                                                                                                        
GCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCT  >  NC_000913/3705887‑3706083

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