Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,027,092:1 +T 5.6% intergenic (‑122/+25) dppA‑II ← / ← dppA‑III dipeptide ABC transporter substrate‑binding protein/dipeptide ABC transporter substrate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,027,0921.T 5.6% 287.5 / 10.3 72intergenic (‑122/+25)dppA‑II/dppA‑IIIdipeptide ABC transporter substrate‑binding protein/dipeptide ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base . (40/28);  new base T (2/2);  total (42/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.24e-01

GCGGAAAGGCAGGTAACGTGGCGCATGGCAAATCCTTTTTCATCTGTTCGCGCCGCCTCCGGGTACGACTGACGGCAACAGTTCGCCGACCCCGTAGCGTATGGCCGCAGGACGACATGAAGGGGTAGCGGGGGTGAAGGGCACCCGCGATTCAGCGGGTGCCTCAGTCAGCCTATTTGTCCACGCTGACGCCGTAGAAGGAATTCAGCGCGAACGGGCTGATCTTGAAGTCCTTCACCTTGGCGCTCATGGGCTGGTAGACAGTGGAGTGGGCGATCGGGGTGATCGGTACCTGC  >  NC_002947/1026945‑1027240
                                                                                                                                                    |                                                                                                                                                     
gcgGAAAGGCAGGTAACGTGGCGCATGGCAAATCCTTTTTCATCTGTTCGCGCCGCCTCCGGGTACGACTGACGGCAACAGTTCGCCGACCCCGTAGCGTATGGCCGCAGGACGACATGAAGGGGTAGCGGGGGTGAAGGGCACCCGCGa                                                                                                                                                     <  1:1674322/150‑1 (MQ=255)
gcgGAAAGGCAGGTAACGTGGCGCATGGCAAATCCTTTTTCATCTGTTCGCGCCGCCTCCGGGTACGACTGACGGCAACAGTTCGCCGACCCCGTAGCGTATGGCCGCAGGACGACATGAAGGGGTAGCGGGGGTGAAGGGCACCCGCGa                                                                                                                                                     <  1:1525021/150‑1 (MQ=255)
  ggAAAGGCAGGTAACGTGGCGCATGGCAAATCCTTTTTCATCTGTTCGCGCCGCCTCCGGGTACGACTGACGGCAACAGTTCGCCGACCCCGTAGCGTATGGCCGCAGGACGACATGAAGGGGTAGCGGGGGTGAAGGGCACCCGCGAtt                                                                                                                                                   >  1:946997/1‑150 (MQ=255)
         cAGGTAACGTGGCGCATGGCAAATCCTTTTTCATCTGCTGACTCCGCCTCCGGGTACGACTGACGGCAACAGTTCGCCGACCCCGTAGCGTATGGCCGCAGGACGACATGAAGGGGTAGCGGGGGTGAAGGGCACCCGCGATTCAGCggg                                                                                                                                            <  1:23059/150‑1 (MQ=255)
          aGGTAACGTGGCGCATGGCAAATCCTTTTTCATCTGTTCGCGCCGCCTCCGGGTACGACTGACGGCAACAGTTCGCCGACCCCGTAGCGTATGGCCGCAGGACGACATGAAGGGGTAGCGGGGGTGAAGGGCACCCGCGATTCAGCGGGt                                                                                                                                           <  1:619451/150‑1 (MQ=255)
          aGGTAACGTGGCGCATGGCAAATCCTTTTTCATCTGTTCGCGCCGCCTCCGGGTACGACTGACGGCAACAGTTCGCCGACCCCGTAGCGTATGGCCGCAGGACGACATGAAGGGGTAGCGGGGGTGAAGGGCACCCGCGATTCAGCGGGt                                                                                                                                           <  2:1444923/150‑1 (MQ=255)
              aaCGTGGCGCATGGCAAATCCTTTTTCATCTGTTCGCGCCGCCTCCGGGTACGACTGACGGCAACAGTTCGCCGACCCCGTAGCGTATGGCCGCAGGACGACATGAAGGGGTAGCGGGGGTGAAGGGCACCCGCGATTCAGCGGGTGCCt                                                                                                                                       >  1:2292553/1‑150 (MQ=255)
                cGTGGCGCATGGCAAATCCTTTTTCATCTGTTCGCGCCGCCTCCGGGTACGACTGACGGCAACAGTTCGCCGACCCCGTAGCGTATGGGCGCAGCACGACCTGCCAGTGCTGCGGGGGTGAAGGGCACCCGCGATTCAGCGGGTGCCtca                                                                                                                                    <  2:1892148/150‑1 (MQ=255)
                   ggCGCATGGCAAATCCTTTTTCATCTGTTCGCGCCGCCTCCGGGTACGACTGACGGCAACAGTTCGCCGACCCCGTAGCGTATGGCCGCAGGACGACATGAAGGGGTACCGGGGGTGAAGGGGACCCGCGATTCAGCGGGTGCCtcagtc                                                                                                                                 <  2:372250/150‑1 (MQ=255)
                      gcATGGCAAATCCTTTTTCATCTGTTCGCGCCGCCTCCGGGTACGACTGACGGCAACAGTTCGCCGACCCCGTAGCGTATGGCCGCAGGACGACATGAAGGGGTAGCGGGGGTGAAGGGCACCCGCGATTCAGCGGGTGCCTCAGTCAGc                                                                                                                              >  2:1417453/1‑150 (MQ=255)
                          ggCAAATCCTTTTTCATCTGTTCGCGCCGCCTCCGGGTACGACTGACGGCAACAGTTCGCCGACCCCGTAGCGTATGGCCGCAGGACGACATGAAGGGGTAGCGGGGGTGAAGGGCACCCGCGATTCAGCGGGTGCCTCAGTCAGCCTAt                                                                                                                          >  2:867090/1‑150 (MQ=255)
                            cAAATCCTTTTTCATCTCTTCCCGCCGCCTCCGGTTACGACTGACGTCAACTGTTCGCCGACCCCGTAGCGTATGGCCGCAGCACGACATGAAGGGGTAGCGGGGGTGAAGGGCACCCGCGATTCAGCGGGTGCCTCAGTCAGCCTAttt                                                                                                                        <  1:1123888/150‑1 (MQ=255)
                               aTCCTTTTTCATCTGTTCGCGCCGTCTCCGGGTACGACTGACGGCAACAGTTCGCCGACCCCGTAGCGTATGGCCGCAGGACGACATGAAGGGGTAGCGGGGGTGAAGGGCACCCGCTGAATCGCGGGTGcc                                                                                                                                        <  1:1939852/132‑1 (MQ=255)
                               aTCCTTTTTCATCTGTTCGCGCCGCCTCCGGGTACGACTGACGGCAACAGTTCGCCGACCCCGTAGCGTATGGCCGCAGGACGACATGAAGGGGTAGCGGGGGTGAAGGGCACCCGCTGAATCGCGGGTGcc                                                                                                                                        >  2:1939852/1‑132 (MQ=255)
                                   tttttCATCTGTTCGCGCCGCCTCCGGGTACGACTGACGGCAACAGTTCGCCGACCCCGTAGCGTATGGCCGCAGGACGACATGAAGGGGTAGCGGGGGTGAAGGGCACCCGCGATTCAGCGGGTGCCTCAGTCAGCCTATTTGTCCACg                                                                                                                 >  1:245178/1‑150 (MQ=255)
                                       tCATCTGTTCGCGCCGCCTCCGGGTACGACTGACGGCAACAGTTCGCCGACCCCGTAGCGTATGGCCGCAGGACGACATGAAGGGGTAGCGGGGGTGAAGGGCACCCGCGATTCAGCGGGTGCCTCAGTCAGCCTATTTGTCCACGCTGa                                                                                                             >  1:341997/1‑150 (MQ=255)
                                           cTGTTCGCGCCGCCTCCGGGTACGACTGACGGCAACAGTTCGCCGACCCCGTAGCGTATGGCCGCAGGACGACATGAAGGGGTAGCGGGGGTGAAGGGCACCCGCGATTCAGCGGGTGCCTCAGTCAGCCTATTTGTCCACGCTGACGcc                                                                                                         >  1:1348157/1‑150 (MQ=255)
                                                 gcgcCGCCACCGGGTACGACTGACGGCAACAGTTCGCCGACCCCGTAGCGTATGGCCGCAGGACGACATGAAGGGGTAGCGGGGGTGAAGGGCACCCGCGATTCAGCGGGTGCCTCAGTCAGCCTATTTGTCCACGCTGACGGCGTAGaa                                                                                                   <  2:225536/150‑1 (MQ=255)
                                                  cgccgcCTCCGGGTACGACTGACGGCAACAGTTCGCCGACCCCGTAGCGTATGGCCGCAGGACGACATGAAGGGGTAGCGGGGGTGAAGGGCACCCGCGATTCAGCGGGTGCCTCAGTCAGCCTATTTGTCCACGCTGACGCCGTAGAAg                                                                                                  <  2:2395388/150‑1 (MQ=255)
                                                                    cTGACGGCAACAGTTCGCCGACCCCGTAGCGTATGGCCGCAGGACGACATGAAGGGGTAGCGGGGGTGAAGGGCACCCGCGATTCAGCGGGTGCCTCAGTCAGCCTATTTGTCCACGCTGACGCCGTAGAAGGAATTCAGCGCGAACggg                                                                                >  1:459434/1‑150 (MQ=255)
                                                                       aCGGCAACAGTTCGCCGACCCCGTAGCGTATGGCCGCAGGACGACATGAAGGGGTAGCGGGGGTGAAGGGCACCCGCGATTCAGCGGGTGCCTCAGTCAGCCTATTTGTCCACGCTGACGCCGTAGAAGGAATTCAGCGCGAACGGGCTg                                                                             >  1:1675408/1‑150 (MQ=255)
                                                                            aaCAGTTCGCCGACCCCGTAGCGTATGGCCGCAGGACGACATGAAGGGGTAGCGGGGGTGAAGGGCACCCGCGATTCAGCGGGTGCCt                                                                                                                                       >  1:1709387/1‑88 (MQ=255)
                                                                            aaCAGTTCGCCGACCCCGTAGCGTATGGCCGCAGGACGACATGAAGGGGTAGCGGGGGTGAAGGGCACCCGCGATTCAGCGGGTGCCt                                                                                                                                       <  2:1709387/88‑1 (MQ=255)
                                                                            aaCAGTTCGCCGACCCCGTAGCGTATGGCCGCAGGACGACATGAAGGGGTAGCGGGGGTGAAGGGCACCCGCGATTCAGCGGGTGCCTCAGTCAGCCTATTTGTCCACGCTGACGCCGTAGAAGGAATTCAGCGCGAACGGGCTGATCtt                                                                        >  1:989204/1‑150 (MQ=255)
                                                                             aCAGTTCGCCGACCCCGTAGCGTATGGCCGCAGGACGACATGAAGGGGTAGCGGGGGTGAAGGGCACCCGCGATTCAGCGGGTGCCTCAGTCAGCCTATTTGTCCACGCTGACGCCGTAGAAGGAATTCAGCGCGAACGGGCTGATCTTg                                                                       >  2:1750588/1‑150 (MQ=255)
                                                                               aGTTCGCCGACCCCGTAGCGTATGGCCGCAGGACGACATGAAGGGGTAGCGGGGGTGAAGGGCACCCGCGATTCAGCGGGTGCCTCAGTCAGCCTATTTGTCCACGCTGACGCCGTAGAAGGAATTCAGCGCGAACGGGCTGATCTTGaa                                                                     <  2:245178/150‑1 (MQ=255)
                                                                                gTTCGCCGACCCCGTAGCGTATGGCCGCAGGACGACATGAAGGGGTAGCGGGGGTGAAGGGCACCCGCGATTCAGCGGGTGCCTCAGTCAGCCTATTTGTCCACGCTGACGCCGTAGAAGGAATTCAGCGCGAACGGGCTGATCTTGAAg                                                                    <  1:1417453/150‑1 (MQ=255)
                                                                                 ttCGCCGACCCCGTAGCGTATGGCCGCAGGACGACATGAAGGGGTAGCGGGGGTGAAGGGCACCCGCGATTCAGCGGGTGCCTCAGTCAGCCTATTTGTCCACGCTGACGCCGTAGAAGGAATTCAGCGCGAACGGGCTGATCTTGAAGt                                                                   >  1:2013425/1‑150 (MQ=255)
                                                                                  tCGCCGACCCCGTAGCGTATGGCCGCAGGACGACATGAAGGGGTAGCGGGGGTGAAGGGCACCCGCGATTCAGCGGGTGCCTCAGTCAGCCTATTTGTCCACGCTGACGCCGTAGAAGGAATTCAGCGCGAACGGGCTGATCTTGAAGTc                                                                  <  2:1348157/150‑1 (MQ=255)
                                                                                    gCCGACCCCGTAGCGTATGGCCGCAGGACGACATGAAGGGGTAGCGGGGGTGAAGGGCACCCGCGATTCAGCGGGTGCCTCAGTCAGCCTATTTGTCCACGCTGACGCCGTAGAAGGAATTCAGCGCGAACGGGCTGATCTTGAAGTCCt                                                                >  1:2078973/1‑150 (MQ=255)
                                                                                     ccGACCCCGTAGCGTATGGCCGCAGGACGACATGAAGGGGTAGCGGGGGTGAAGGGCACCCGCGATTCAGCGGGTGCCTCAGTCAGCCTATTTGTCCACGCTGACGCCGTAGAAGGAATTCAGCGCGAACGGGCTGATCTTGAAGTCCtt                                                               <  2:2290595/150‑1 (MQ=255)
                                                                                          cccGTAGCGTATGGCCGCAGGACGACATGAAGGGGTAGCGGGGGTGAAGGGCACCCGCGATTCAGCGGGTGCCTCAGTCAGCCTATTTGTCCACGCTGACGCCGTAGAAGGAATTCAGCGCGAACGGGCTGa                                                                            <  1:1758775/132‑1 (MQ=255)
                                                                                          cccGTAGCGTATGGCCGCAGGACGACATGAAGGGGTAGCGGGGGTGAAGGGCACCCGCGATTCAGCGGGTGCCTCAGTCAGCCTATTTGTCCACGCTGACGCCGTAGAAGGAATTCAGCGCGAACGGGCTGa                                                                            >  2:1758775/1‑132 (MQ=255)
                                                                                          cccGTAGCGTATGGCCGCAGGACGACATGAAGGCGTAGCGGAGGAGAAGTGCACGCGGGATTCAGCGGGTGCCTAAACTACGCCCTGCTGGGCAGCTGACGCCGTAGAAGGAATTCAGCGCGAACGGGCTGATCTTGAAGTCCTTCACCt                                                          >  1:1749194/1‑150 (MQ=255)
                                                                                           ccGTAGCGTATGGCCGCAGGACGACATGAAGGGGTAGCGGGGGTGAAGGGCACCCGCGATTCAGCGGGTGCCTCAGTCAGCCTATTTGTCCACGCTGACGCCGTAGAAGGAATTCAGCGCGAACGGGCTGATCTTGAAGTCCTTCACCtt                                                         >  2:131379/1‑150 (MQ=255)
                                                                                            cGTAGCGTATGGCCGCAGGACGACATGAAGGGGTAGCGGGGGTGAAGGGCACCCGCGATTCAGCGGGTGCCTCAGTCAGGCTATTTGTCCACGCTGACGCCGTAGAAGGAATTCAGCGCGAACGGACTGATCTTGAAGTCCTTCACCTTg                                                        >  1:2275005/1‑150 (MQ=255)
                                                                                            cGTAGCGTATGGCCGCAGGACGACATGAAGGGGTAGCGGGGGTGAAGGGCACCCGCGATTCAGCGGGTGCCTCAGTCAGCCTATTTGTCCACGCTGACGCCGTAGAAGGAATTCAGCGCGAACGGGCTGATCTTGAAGTCCTTCACCTTg                                                        >  1:1980237/1‑150 (MQ=255)
                                                                                                   tATGGCCGCAGGACGACATGAAGGGGTAGCGGGGGTGAAGGGCACCCGCGATTCAGCGGGTGCCTCAGTCAGCCTATTTGTCCACGCTGACGCCGTAGAAGGAATTCAGCGCGAACGGGCTGATCTTGAAGTCCTTCACCTTGGCGCTCa                                                 >  2:1775583/1‑150 (MQ=255)
                                                                                                       gCCGCAGGACGACATGAAGGGGTAGCGGGGGTGAAGGGCACCCGCGATTCAGCGGGTGCCTCAGTCAGCCTATTTGTCCACGCTGACGCCGTAGAAGGAATTCAGCGCGAACGGGCTGATCTTGAAGTCCTTCACCTTGGCGCTCATggg                                             >  1:2116798/1‑150 (MQ=255)
                                                                                                       gCCGCAGGACGACATGAAGGGGTAGCGGGGGTGAAGGGCACCCGCGATTCAGCGGGTGCCTCAGTCAGCCTATTTGTCCACGCTGACGCCGTAGAAGGAATTCAGCGCGAACGGGCTGATCTTGAAGTCCTTCACCTTGGCGCTCATggg                                             >  1:920916/1‑150 (MQ=255)
                                                                                                         cGCAGGACGACATGAAGGGGTAGCGGGGGTGAAGGGCACCCGCGATTCAGCGGGTGCCTCAGTCAGCCTATTTGTCCACGCTGACGCCGTAGAAGGAATTCAGCGCGAACggg                                                                                <  2:1582036/113‑1 (MQ=255)
                                                                                                         cGCAGGACGACATGAAGGGGTAGCGGGGGTGAAGGGCACCCGCGATTCAGCGGGTGCCTCAGTCAGCCTATTTGTCCACGCTGACGCCGTAGAAGGAATTCAGCGCGAACggg                                                                                >  1:1582036/1‑113 (MQ=255)
                                                                                                         cGCAGGACGACATGAAGGGGTAGCGGGGGTGAAGGGCACCCGCGATTCAGCGGGTGCCTCAGTCAGCCTATTTGTCCACGCTGACGCCGTAGAAGGAATTCAGCGCGAACGGGCTGATCTTGAAGTCCTTCACCTTGGCGCTCATGGGCt                                           >  1:1778711/1‑150 (MQ=255)
                                                                                                          gCAGGACGACATGAAGGGGTAGCGGGGGTGAAGGGCACCCGCGATTCAGCGGGTGCCTCAGTCAGCCTATTTGTCCACGCTGACGCCGTAGAAGGAATTCAGCGCGAACGGGCTGATCTTGAAGTCCTTCACCTTGGCGCTCATGGGCTg                                          >  1:1118359/1‑150 (MQ=255)
                                                                                                          gCAGGACGACATGAAGGGGTAGCGGGGGTGAAGGGCACCCGCGATTCAGCGGGTGCCTCAGTCAGCCTATTTGTCCACGCTGACGCCGTAGAAGGAATTCAGCGCGAACGGGCTGATCTTGAAGTCCTTCACCTTGGCGCTCATGGGCTg                                          <  1:867090/150‑1 (MQ=255)
                                                                                                          gCAGGACGACATGAAGGGGTAGCGGGGGTGAAGGGCACCCGCGATTAAGAAACAGCCTCAGTCAGCCTATTTGTCCACGCTGACGCCGTAGAAGGAATTCAGCGCGAACGGGCTGATCTTGAAGTCCTTCACCTTGGCGCTCATGGGCTg                                          >  1:1376741/1‑150 (MQ=255)
                                                                                                             ggACGACATGAAGGGGTAGCGGGGGTGAAGTGCACCCGCGATTCAGCGGGTGCCTCAGTCAGCCTATTTGTCCACGCTGACGCCGTAGAAGGAATTCAGCGCGAACGGGCTTATCTTGAAGTCCTTCACCTTGGCGCTCATGGGCTGGTa                                       <  2:1698143/150‑1 (MQ=255)
                                                                                                                 gacATGAAGGGGTAGCGGGGGTGAAGGGCACCCGCGATTCAGCGGGTGCCTCAGTCAGCCTATTTGTCCACGCTGACGCCGTAGAAGGAATTCAGCGCGAACGGGCTGATCTTGAAGTCCTTCACCTTGGCGCTCATGGGCTGGTAGACa                                   >  2:2276839/1‑150 (MQ=255)
                                                                                                                 gacATGAAGGGGTAGCGGGGGTGAAGGGCACCCGCGATTCAGCGGGTGCCTCAGTCAGCCTATTTGTCCACGCTGACGCCGTAGAAGGAATTCAGCGCGAACGGGCTGATCTTGAAGTCCTTCACCTTGGCGCTCATGGGCTGGTAGACa                                   >  2:1464092/1‑150 (MQ=255)
                                                                                                                   cATGAAGGGGTAGCGGGGGTGAAGGGCACCCGCGATTCAGCGGGTGCCTCAGTCAGCCTATTTGTCCACGCTGACGCCGTAGAAGGAATTCAGCGCGAACGGGCTGATCTTGAAGTCCTTCACCTTGGCGCTCATGGGCTGGTAGACagt                                 <  1:2276839/150‑1 (MQ=255)
                                                                                                                   cATGAAGGGGTAGCGGGGGTGAAGGGCACCCGCGATTCAGCGGGTGCCTCAGTCAGCCTATTTGTCCACGCTGACGCCGTAGAAGGAATTCAGCGCGAACGGGCTGATCTTGAAGTCCTTCACCTTGGCGCTCATGGGCTGGTAGACagt                                 <  1:1371486/150‑1 (MQ=255)
                                                                                                                    aTGAAGGGGTAGCGGGGGTGAAGGGCACCCGCGATTCAGCGGGTGCCTCAGTCAGCCTATTTGTCCACGCTGACGCCGTAGAAGGAATTCAGCGCGAACGGGCTGATCTTGAAGTCCTTCACCTTGGCGCTCATGGGCTGGTAGACagtg                                <  2:181460/150‑1 (MQ=255)
                                                                                                                     tGAAGGGGTAGCGGGGGTGAAGGGCACCCGCGATTCAGCGGGTGCCTCAGTCAGCCTATTTGTCCACGCTGACGCCGTAGAAGGAATTCAGCGCGAACGGGCTGATCTTGAAGTCCTTCACCTTGGCGCTCATGGGCTGGTAGACagtgg                               >  2:1646404/1‑150 (MQ=255)
                                                                                                                         ggggTAGCGGGGGTGAAGGGCACCCGCGATTCAGCGGGTGCCTCAGTCAGCCTATTTGTCCACGCTGACGCCGTAGAAGGAATTCAgcgc                                                                                       >  1:1606860/1‑90 (MQ=255)
                                                                                                                         ggggTAGCGGGGGTGAAGGGCACCCGCGATTCAGCGGGTGCCTCAGTCAGCCTATTTGTCCACGCTGACGCCGTAGAAGGAATTCAgcgc                                                                                       <  2:1606860/90‑1 (MQ=255)
                                                                                                                          gggTAGCGGGGGTGAAGGGCACCCGCGATTCAGCGGGTGCCTCAGTCAGCCTATTTGTCCACTCTGACGCCGTAGAAGGAATTCAGCGCGAACGGGCTGATCTTGAAGTCCTTCACCTTGGCGCTCATGGGCTGGTAGACagtggagtgg                          >  2:2358163/1‑150 (MQ=255)
                                                                                                                            gTAGCGGGGGTGAAGGGCACCCGCGATTCAGCGGGTGCCTCAGTCAGCCTATTTGTCCACGCTGACGCCGTAGAAGGAATTCAGCGCGAACGGGCTGATCTTGAAGTCCTTCACCTTGGCGCTCATGGGCTGGTAGACAGTGGAGTGGGc                        >  1:1256751/1‑150 (MQ=255)
                                                                                                                            gTAGCGGGGGTGAAGGGCACCCGCGATTCAGCGGGTGCCTCAGTCAGCCTATTTGTCCACGCTGACGCCGGCGAAGGAATTCAGCGCGAACGGGCTGATCTTGAAGTCCTTCACCTTGGCGCTCATGGGCTGGTAGACAGTGGAGTGGGc                        >  1:853197/1‑150 (MQ=255)
                                                                                                                                cGGGGGTGAAGGGCACCCGCGATTCAGCGGGTGCCTCAGTCAGCCTATTTGTCCACGCTGACGCCGTAGAAGGAATTCAGCGCGAACGGGCTGATCTTGAAGTCCTTCACCTTGGCGCTCATGGGCTGGTAGACAGTGGAGTGGGCGATc                    <  2:1694409/150‑1 (MQ=255)
                                                                                                                                cGGGGGTGAAGGGCACCCGCGATTCAGCGGGTGCCTCAGTCAGCCTATTTGTCCACGCTGACGCCGTAGAAGGAATTCAGCGCGAACGGGCTGATCTTGAAGTCCTTCACCTTGGCGCTCATGGGCTGGTAGACAGTGGAGTGGGCGATc                    >  1:1694409/1‑150 (MQ=255)
                                                                                                                                  ggggTGAAGGGCACCCGCGATTCAGCGGGTGCCTCAGTCAGCCTATTTGTCCACGCTGACGCCGTAGAAGGAATTCAGCGCGAACGGGCTGATCTTGAAGTCCTTCACCTTGGCGCTCATGGGCTGGTAGACAGTGGAGTGGGCGATCgg                  <  2:920916/150‑1 (MQ=255)
                                                                                                                                      tGAAGGGCACCCGCGATTCAGCGGGTGCCTCAGTCAGCCTATTTGTCCACGCTGACGCCGTAGAAGGAATTCAGCGCGAACGGGCTGATCTTGAAGTCCTTCACCTTGGCGCTCATGGGCTGGTAGACAGTGGAGTGGGCGATCGGGGTg              >  2:2092919/1‑150 (MQ=255)
                                                                                                                                      actgaggCACCCGCTGAATCGCGGGTGCCTCAGTCAGCCTATTTGTCCACGCTGACGCCGTAGAAGGAATTCAGCGCGAACGGGCTGATCTTGAAGTCCTTCACCTTGGCGCTCATGGGCTGGTAGACAGTGGAGTGGGCGATCggg                 <  2:155805/142‑1 (MQ=255)
                                                                                                                                      actgaggCACCCGCTGAATCGCGGGTGCCTCAGTCAGCCTATTTGTCCACGCTGACGCCGTAGAAGGAATTCAGCGCGAACGGGCTGATCTTGAAGTCCTTCACCTTGGCGCTCATGGGCTGGTAGACAGTGGAGTGGGCGATCggg                 >  1:155805/6‑147 (MQ=255)
                                                                                                                                        aaGGGCACCCGCGATTCAGCGGGTGCCTCAGTCAGCCTATTTGTCCACGCTGACGCCGTAGAAGGAATTCAGCGCGAACGGGCTGATCTTGAAGTCCTTCACCTTGGCGCTCATGGGCTGGTAGACAGTGGAGTGGGCGATCGGGGTGAt            >  1:2004994/1‑150 (MQ=255)
                                                                                                                                         aGGGCACCCGCGATTCAGCGGGTGCCTCAGTCAGCCTATTTGTCCACGCTGACGCCGTAGAAGGAATTCAGCGCGAACGGGCTGATCTTGAAGTCCTTCACCTTGGCGCTCATGGGCTGGTAGACAGTGGAGTGGGCGATCGGGGTGATc           <  1:1464092/150‑1 (MQ=255)
                                                                                                                                             cACCCGCGATTCAGCGGGTGCCTCAGTCAGCCTATTTGTCCACGCTGACGCCGTAGAAGGAATTCAGCGCGAACGGGCTGATCTTGAAGTCCTTCACCTTGACGCTCATGGGCTGGTAGACAGTGGAGTGGGCGATCGGGGTGATCGGTa       >  1:1662290/1‑150 (MQ=255)
                                                                                                                                              aCCCGCGATTCAGCGGGTGCCTCAGTCAGCCTATTTGTCCACGCTGACGCCGTAGAAGGAATTCAGCGCGAACGGGCTGATCTTGAAGTCCTTCACCTTGGCGCTCATGGGCTGGTAGACAGTGGAGTGGGCGATCGGGGTGATCGGTAc      <  2:2004994/150‑1 (MQ=255)
                                                                                                                                                ccGCGATTCAGCGGGTGCCTCAGTCAGCCTATTTGTCCACGCTGACGCCGTAGAAGGAATTCAgcgc                                                                                       >  2:1445821/1‑67 (MQ=255)
                                                                                                                                                ccGCGATTCAGCGGGTGCCTCAGTCAGCCTATTTGTCCACGCTGACGCCGTAGAAGGAATTCAgcgc                                                                                       <  1:1445821/67‑1 (MQ=255)
                                                                                                                                                ccGCGATTCAGCGGGTGCCTCAGTCAGCCTATTTGTCCACGCTGACGCCGTAGAAGGAATTCAGCGCGAACGGGCTGATCTTGAAGTCCTTCACCTTGGCGCTCATGGGCTGGTAGACAGTGGAGTGGGCGATCGGGGTGATCGGTACCt    >  1:2315528/1‑150 (MQ=255)
                                                                                                                                                ccGCGATTCAGCGGGTGCCTCAGTCAGCCTATTTGTCCACGCTGACGCCGCAGAAGGACTTCAGCCCGAACGGGCTGATCTTCAAGTACTTCACCTCGGCGCTCATGGGCCGGTAGACAGTGCAGTGGGCGCTCGGCGCGATCGGCACCt    >  2:81028/1‑150 (MQ=255)
                                                                                                                                                  gcgATTCAGCGGGTGCCTCAGTCAGCCTATTTGTCCACGCTGACGCCGTAGAAGGAATTCAGCGCGAACGGGCTGATCTTGAAGTCCTTCACCTTGGCGCTCATGGGCTGGTAGACAGTGGAGTGGGCGATCGGGGTGATCGGTACCTGc  >  1:102394/1‑150 (MQ=255)
                                                                                                                                                  gcgATTCAGCGGGTGCCTCAGTCAGCCTATTTGTCCACGCTGACGCCGTAGAAGGAATTCAGCGCGAACGGGCTGATCTTGAAGTCCTTCACCTTGGCGCTCATGGGCTGGTAGACAGTGGAGAGGGCGAGCGGGGTGATCGGTACCTGc  >  2:560391/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GCGGAAAGGCAGGTAACGTGGCGCATGGCAAATCCTTTTTCATCTGTTCGCGCCGCCTCCGGGTACGACTGACGGCAACAGTTCGCCGACCCCGTAGCGTATGGCCGCAGGACGACATGAAGGGGTAGCGGGGGTGAAGGGCACCCGCGATTCAGCGGGTGCCTCAGTCAGCCTATTTGTCCACGCTGACGCCGTAGAAGGAATTCAGCGCGAACGGGCTGATCTTGAAGTCCTTCACCTTGGCGCTCATGGGCTGGTAGACAGTGGAGTGGGCGATCGGGGTGATCGGTACCTGC  >  NC_002947/1026945‑1027240

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: