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.

GCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑G‑‑‑A‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCT  >  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)
                                      cGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑g                                                                  <  2:689259/100‑1 (MQ=255)
                                      cGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑g                                                                  >  2:1595785/1‑100 (MQ=255)
                                           ataCAGATTGTTGCACAAGTGCTCTGTGTTTATTCCGTATGCGGCGTTTACGCCGTATCCTGCCTACAAAATCGTGAAAATTCAATGCATTGCA‑G‑‑‑a                                                              <  1:3735098/96‑1 (MQ=12)
                                             aCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑G‑‑‑A‑AT‑TTTc                                                      >  1:1134631/1‑100 (MQ=255)
                                               aGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑G‑‑‑A‑AT‑TTTCGt                                                    >  1:238974/1‑100 (MQ=255)
                                                 aTTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑G‑‑‑A‑AT‑TTTCGTAg                                                  <  1:2145517/100‑1 (MQ=255)
                                                 aTTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑G‑‑‑A‑AT‑TTTCGTAg                                                  <  1:604568/100‑1 (MQ=255)
                                                   tgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑G‑‑‑A‑AT‑TTTCGTAgcc                                                >  1:2342165/1‑98 (MQ=255)
                                                   cataaaCAAAGCGCACTTTGGTTATGCAGCATGCGGCGCTAACACCGTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑G‑‑‑A‑AT‑TTTCGTa                                                   >  1:2492197/6‑97 (MQ=17)
                                                      tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑G‑‑‑A‑AT‑TTTCGTAGGCCTg                                             >  2:1555601/1‑100 (MQ=255)
                                                       gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑G‑‑‑A‑AT‑TTTCGTAGGCCTGa                                            >  2:765150/1‑100 (MQ=255)
                                                        aCAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑G‑‑‑A‑AT‑TTTCGTAGGCCTGAt                                           >  2:1519232/1‑100 (MQ=255)
                                                        aCAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTACA‑G‑‑‑A‑AT‑TTTCGTAGGCCTGAt                                           <  2:1123638/100‑1 (MQ=255)
                                                          aaaGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑G‑‑‑A‑AT‑TTTCGTAGGCCTGATaa                                         <  2:1208219/100‑1 (MQ=255)
                                                             gTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑G‑‑‑A‑AT‑TTTCGTAGGCCTGATAAGCg                                      <  1:625695/100‑1 (MQ=255)
                                                                 gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑G‑‑‑A‑AT‑TTTCGTAGGCCTGATAAGCGGAgc                                  <  1:1257887/100‑1 (MQ=255)
                                                                  cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑G‑‑‑A‑AT‑TTTCGTAGGCCTGATAAGCGGAgcg                                 >  1:3477021/1‑100 (MQ=255)
                                                                  cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑G‑‑‑A‑AT‑TTTCGTAGGCCTGATAAGCGGAgcg                                 >  2:3190479/1‑100 (MQ=255)
                                                                     tttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCA‑G‑‑‑G‑GA‑CTATGTAGGCCTGATAAGCATAGCGCAt                              >  2:2718370/8‑100 (MQ=18)
                                                                     tttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCA‑G‑‑‑AGAT‑TAT‑GTAGGCCTGATAAGCGTAGCGCAt                              >  2:297290/8‑100 (MQ=17)
                                                                      ttgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGCA‑G‑‑‑A‑GT‑TCACGTAGGCCTGATAAGCGTAGCGCATc                             >  1:1985931/7‑100 (MQ=17)
                                                                       tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑G‑‑‑A‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCa                            <  2:2686460/100‑1 (MQ=255)
                                                                       tctaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAATCATGCAAATCCAAAATATTGCA‑G‑‑‑G‑AG‑TTATGTAGGCCTGATAAGCGTAGCGCATCAg                           <  1:2876362/95‑1 (MQ=21)
                                                                        ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑G‑‑‑A‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAg                           >  2:3411144/1‑100 (MQ=255)
                                                                        ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑G‑‑‑A‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAg                           <  1:3191798/100‑1 (MQ=255)
                                                                        gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCATACAAAAACCTGCAAATTCAATATATTGCA‑G‑‑‑G‑AT‑ACGTGTAGGCCTGATAAGCGTAGCGCATCAg                           >  2:3044992/5‑100 (MQ=18)
                                                                         taggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCA‑G‑‑‑A‑AA‑CCTTGTAGGCCTGATAAGCGTAGCGCATCAgg                          >  2:3749263/4‑100 (MQ=18)
                                                                          aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAATTC‑TGCAAATTCAACAAATTGC‑‑G‑‑‑A‑ATCCCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCa                        >  1:3093945/3‑100 (MQ=18)
                                                                          aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTACA‑G‑‑‑A‑AA‑TCTTGTAGGCCTGATAAGCGTAGCGCATCACGc                         >  1:2711244/3‑100 (MQ=14)
                                                                          aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCA‑‑‑‑‑A‑TA‑TGACGTAGGCCTGATAAGCGTAGCGCATCAGGcg                        >  2:3844775/3‑99 (MQ=18)
                                                                          aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAATCATGCAAATCCAAAATATTGCA‑G‑‑‑G‑AG‑TTATGTAGGCCTGATAAGCGTAGCGCATCAgaga                        >  2:2876384/3‑97 (MQ=21)
                                                                          aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAATCATGCAAATCCAAAATATTGCA‑G‑‑‑G‑AG‑TTATGTAGGCCTGATAAGCGTAGCGCATCAGGCa                        >  2:2378803/3‑100 (MQ=21)
                                                                          aggCCGGATAAGGCGTTCACGCCGCATCCGGCATACAAAATTGTGTAAATTCAATAAATTGAT‑C‑‑‑G‑AT‑CCACGTAGGCCTGATAAGCGTAGCGCATCAGGc                         >  2:3519958/3‑100 (MQ=21)
                                                                           ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAATTC‑TGCAAATTCAACAAATTGC‑‑G‑‑‑A‑ATCCCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAt                       >  2:3805789/2‑100 (MQ=18)
                                                                           ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCA‑G‑‑‑A‑AA‑CCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCa                        >  2:3286393/2‑100 (MQ=17)
                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCA‑G‑‑‑G‑AG‑CTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCa                        >  1:951293/2‑100 (MQ=17)
                                                                            gCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAATTC‑TGCAAATTCAACAAATTGC‑‑G‑‑‑A‑ATCCCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAtt                      >  2:1277516/1‑100 (MQ=18)
                                                                             ccggataaggcgTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCA‑G‑‑‑T‑TT‑ATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCaat                      >  2:2003096/9‑98 (MQ=25)
                                                                             ccggataaggcgTTCACGCCGCATCCGGCCTACAAAATCATGAAAATTCAATACATTGCA‑A‑‑‑G‑AT‑TTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCacg                      >  1:3401505/9‑98 (MQ=17)
                                                                             ccggataaggcgTTCACGCCGCATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCA‑G‑‑‑G‑AG‑CTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCagt                      >  2:3218106/9‑98 (MQ=17)
                                                                             ccggataaggcgTTCACGCCGCATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCA‑G‑‑‑G‑AG‑CTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCagt                      >  2:2641618/9‑98 (MQ=18)
                                                                             ccggataaggcgTTCACGCCGCATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCA‑G‑‑‑G‑GA‑CTATGTAGGCCTGATAAGCATAGCGCATCAGGCagc                      >  1:2696896/9‑98 (MQ=18)
                                                                             ccggataaggcgTTCACGCCGCATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCA‑G‑‑‑A‑GA‑TTATGTAGGCCTGATAAGCGTAGCGCATCAGGCaat                      >  2:3697232/9‑98 (MQ=17)
                                                                             ccGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑G‑‑‑A‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCAtt                      >  2:1906813/1‑100 (MQ=255)
                                                                             ccGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑G‑‑‑A‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCAtt                      >  1:3546130/1‑100 (MQ=255)
                                                                               ngataaggcgTTCACGCCGCATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGG‑‑‑G‑AC‑TAT‑GTAGGCCTGATAAGCATAGCGCATCAGGCagctt                    >  2:3145180/7‑96 (MQ=18)
                                                                               ggataaggcgTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGCA‑G‑‑‑A‑GT‑TCACGTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                    >  2:3307601/7‑100 (MQ=17)
                                                                               ggataaggcgTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGCA‑G‑‑‑A‑GT‑TCACGTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                    >  2:934387/7‑100 (MQ=17)
                                                                               ggataaggcgTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGCA‑G‑‑‑A‑GT‑TCACGTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                    >  2:1673375/7‑100 (MQ=17)
                                                                               ggataaggcgTTTACGCCGCATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCA‑ATTCA‑AC‑TT‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaatt                     >  2:929691/7‑97 (MQ=21)
                                                                               ggataaggcgTTCACGCCGCATCCGGCCTACAAAATCATGAAAATTCAATACATTGCA‑A‑‑‑G‑AT‑TTTCGTAGGCCTGATAAGCGTAGCGAATCAGGCacgc                     >  2:2779062/7‑96 (MQ=17)
                                                                               ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAATCATGCAAATCCAAAATATTGCA‑G‑‑‑G‑AG‑TTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAtttat                   >  2:2166157/7‑98 (MQ=21)
                                                                               ggataaggcgTTCACGCCGCATCCGACCCACAAAGTCTTGCAAAATCAATAAATTGCA‑G‑‑‑G‑AA‑CGATGTAGGCCTGATAAGCGTAGCGCATCAGGCaacta                    >  1:1013021/7‑96 (MQ=21)
                                                                               ggATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑G‑‑‑A‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCAtttt                    <  2:1220856/100‑1 (MQ=255)
                                                                               ggATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑G‑‑‑A‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCAtttt                    >  1:1435319/1‑100 (MQ=255)
                                                                                gataaggcgTTCACGCCGCATCCGGCCTACAAATCCTTGCGAATTCAATATATTGTATG‑‑‑G‑AA‑ATTC‑‑AGGCCTGATAAGCGTAGCGCATCAGGCttttcac                  >  1:2805802/6‑97 (MQ=37)
                                                                                gataaggcgTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCA‑G‑‑‑T‑TT‑ATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCaatctt                   >  2:637052/6‑95 (MQ=25)
                                                                                gataaggcgTTCACGCCGCATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCA‑G‑‑‑G‑AA‑ACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAgtttt                   >  2:436250/6‑100 (MQ=18)
                                                                                 aTGCGGCGTTAACGCTTTATCCGGCCTACAAAATCGTGCAAATCCAATACATTGCA‑G‑‑‑A‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                  <  2:491520/100‑1 (MQ=255)
                                                                                   gcggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑G‑‑‑A‑AT‑TTTCGTAGGCCTGATACGCGGAGCGCATCAGCCAttttttc                 >  2:663766/1‑97 (MQ=255)
                                                                                    cggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑G‑‑‑A‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCAttttttttt               >  1:6669/1‑96 (MQ=255)
                                                                                    cggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑G‑‑‑A‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct               <  2:2436025/100‑1 (MQ=255)
                                                                                    cggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑G‑‑‑A‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct               >  2:178452/1‑100 (MQ=255)
                                                                                       cgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑G‑‑‑A‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg            <  1:1461450/100‑1 (MQ=255)
                                                                                       cgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑G‑‑‑A‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg            <  1:1056723/100‑1 (MQ=255)
                                                                                        gTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑G‑‑‑A‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgc           <  1:314861/100‑1 (MQ=255)
                                                                                         cttccgatctATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑G‑‑‑A‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgcc          <  1:2592890/91‑1 (MQ=255)
                                                                                          tAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑G‑‑‑A‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCa         >  2:2065443/1‑100 (MQ=255)
                                                                                          tAACGCCTTATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCA‑G‑‑‑A‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCa         <  1:1558964/100‑1 (MQ=255)
                                                                                            aCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑G‑‑‑A‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc       <  2:3528349/100‑1 (MQ=255)
                                                                                             cGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑G‑‑‑A‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGTTTTTTCTATCa      >  1:1495488/1‑100 (MQ=17)
                                                                                             cGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑G‑‑‑A‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc       <  1:515611/99‑1 (MQ=255)
                                                                                              gCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑G‑‑‑A‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAg     >  1:202999/1‑100 (MQ=255)
                                                                                              gCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑G‑‑‑A‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAg     <  1:1743347/100‑1 (MQ=255)
                                                                                               ncTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑G‑‑‑A‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGt    >  1:823598/2‑100 (MQ=255)
                                                                                                  |                                                                                                        
GCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑G‑‑‑A‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCT  >  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: ‑