Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_bme_pgi 4,235,408 Δ1 bp intergenic (‑2175/‑499) lysC ← / → yjbE lysine‑sensitive aspartokinase 3/extracellular polysaccharide production threonine‑rich protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_bme_pgi4,235,4070A.97.1% 280.5 / ‑0.7 68intergenic (‑2174/‑500)lysC/yjbElysine‑sensitive aspartokinase 3/extracellular polysaccharide production threonine‑rich protein
Reads supporting (aligned to +/- strand):  ref base A (2/0);  new base . (30/36);  total (32/36)
Fisher's exact test for biased strand distribution p-value = 2.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCAATGCCTTCGATCAATGGGGCGTGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCG‑GCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATT  >  NC_000913_3_bme_pgi/4235259‑4235550
                                                                                                                                                    |                                                                                                                                                
tCAATGCCTTCGATCCATGGGGCGGGGTACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGa                                                                                                                                                  <  1:8626/149‑1 (MQ=255)
    tGCCTTCGATCAATGGGGCGTGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGGGCGAAGGAAAGCTTGatcat                                                                                                                                              >  2:164448/1‑145 (MQ=255)
     gCCTTCGATCAATGGGGCGTGGAACTGGGTAAGGAACTGGGGCAGGCTTTTTTCCCCGTTCTACCTGGCAAGGAAGGGGCGAGTGGGCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑atcatc                                                                                                                                            <  2:717865/149‑1 (MQ=255)
        ttCGATCAATGGGGCGTGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGACGGAAAGCTTG‑ATCAtcgtc                                                                                                                                         >  2:65683/1‑149 (MQ=255)
        ttCGATCAATGGGGCGTGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCAtcgtc                                                                                                                                         >  1:172478/1‑149 (MQ=255)
          cGATCAATGGGGCGTGGAACGGGGTAAGGAACGGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGa                                                                                                                                       <  2:158021/149‑1 (MQ=255)
                    ggCGTGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAGGGGGCGAGTGGTAGCGACGCTTCGACGCGAGGGCTTGTCGCGCAGCTGCATGCGCGAAGGGAAGCTTG‑ATCATCGTCGATATGTAGGcc                                                                                                                             >  2:397634/1‑149 (MQ=255)
                          gAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATaa                                                                                                                       <  1:1001513/149‑1 (MQ=255)
                          gAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATaa                                                                                                                       <  2:83640/149‑1 (MQ=255)
                           aaCTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGGGCGAAGGCAAGCTTG‑ATCATCGTCGATATGCAGGCCGGATAAg                                                                                                                      >  2:306825/1‑149 (MQ=255)
                           aaCTGGGAAAGGCACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAg                                                                                                                      <  1:423853/149‑1 (MQ=255)
                             cTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCCTGCGCGAAGGAAAGCTTG‑ACAATCAGCGATATCTAGGCCGGATAATGc                                                                                                                    >  2:415990/1‑149 (MQ=255)
                              tGGGGAAGGAACTGGCGACGGAACTGATGCGGGTTGTATCGGGCGAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCg                                                                                                                   <  1:521439/149‑1 (MQ=255)
                               gggTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGt                                                                                                                  >  2:407739/1‑149 (MQ=255)
                                ggTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGtt                                                                                                                 >  2:549748/1‑149 (MQ=255)
                                     ggAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGTCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTGGGCCGGATAAGGCGTTCAcgc                                                                                                            >  1:909207/1‑149 (MQ=255)
                                       acgtGGCGACGGAACTTGTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCAcgccg                                                                                                          <  2:438411/146‑1 (MQ=255)
                                        acgggCGACGGAAGTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCAcgccgc                                                                                                         <  2:194982/146‑1 (MQ=255)
                                           ggCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATc                                                                                                      <  1:481921/149‑1 (MQ=255)
                                             cGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCg                                                                                                    <  1:694410/149‑1 (MQ=255)
                                              gACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGAAACCg‑g                                                                                                  >  1:565767/1‑149 (MQ=255)
                                                 ggAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGAATAAGGCGTTAAAGCCGCATCCG‑GCaa                                                                                               >  1:781198/1‑149 (MQ=255)
                                                   aaCTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCG‑GCAAcc                                                                                             >  1:538308/1‑149 (MQ=255)
                                                      ttttGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCG‑GCAACCGAt                                                                                          >  2:646234/1‑149 (MQ=255)
                                                      ttttGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTGGTCGCGCATCTGCATGCGCGAAGGAAAGCTT‑AATCATCGTCGATATGTAGGCCGGTTAAGCCGTTCACGCCGCATACG‑GGACCCGAt                                                                                          >  2:300335/1‑149 (MQ=255)
                                                          gCCGGTTGTATCGGGCAAGGAAGGGGCCAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCG‑GCAACCGATGCCt                                                                                      <  2:312934/149‑1 (MQ=255)
                                                            cGGTTGTATCGGGCAAGGAAGGGGCGAGTGGGCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCG‑GCAACCGATGCCTGa                                                                                    <  2:72426/149‑1 (MQ=255)
                                                             ggtttgaTCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCTTTCACGCCGCATCCG‑GCAACCGATGCCTGAt                                                                                   <  2:303208/143‑1 (MQ=255)
                                                             ggTTGTATCGGGCAAGGAAGGGGTGGGTGGTCGCGACGCTTCCATGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGAGATGTAGGCCGGATAAGGCGTTCACGCCGCATCCG‑GCAACCGATGCCTGAt                                                                                   <  1:49580/149‑1 (MQ=39)
                                                                     cGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCG‑GCAACCGATGCCTGATGCGAcgcg                                                                           <  2:509066/149‑1 (MQ=255)
                                                                      gggCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATCTGTAGGCCGGCTAAGGCGATCACGCCGCATCCG‑GCAACCGATACCTGATGCGACCCgg                                                                          >  1:961242/1‑149 (MQ=255)
                                                                      gggCAAGGAAGGGGCGAGTGGTCGCGACGCTTCAACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCG‑GCAACCGATGCCTGATGCGACGCgg                                                                          >  2:677461/1‑149 (MQ=255)
                                                                       ggCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCG‑GCAACCGATGCCTGATGCGACGCGGt                                                                         >  1:647738/1‑149 (MQ=255)
                                                                             gaaggGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCAcgccg                                                                                                          >  1:483354/1‑111 (MQ=255)
                                                                             gaaggGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCAcgccg                                                                                                          <  2:483354/111‑1 (MQ=255)
                                                                               aggGGCGAGTGTGTGCGGCGCTTCCCCGCAAGTGCTCGGCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGTATCCG‑GCAACCGATGCCTGATGCGACGCGGTCGCGTCtt                                                                 <  2:699395/149‑1 (MQ=25)
                                                                                 gggCGAGTGGTCGGGACGCGTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCCTCCG‑GCAACCGATGCCTGATAAGACGCGACCGCGTcgcat                                                               <  2:190149/149‑5 (MQ=38)
                                                                                    cGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCG‑GCAACCGATGCCTGATAAGACGCGACCGCGTcgcatcag                                                            <  1:407739/149‑8 (MQ=38)
                                                                                     gAGTGGTCGCGACGCTTCCACGCAAGGGGTTGTCGCGCATCTGGATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCG‑GCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAgg                                                           <  1:562064/149‑1 (MQ=255)
                                                                                        tGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCCCGCCGCATCCG‑GCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCt                                                        >  1:641178/1‑149 (MQ=255)
                                                                                        tGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCG‑GCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCt                                                        >  1:913560/1‑149 (MQ=255)
                                                                                        tGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCG‑GCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCt                                                        >  2:216621/1‑149 (MQ=255)
                                                                                         ggTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCG‑GCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTa                                                       <  2:594423/149‑1 (MQ=255)
                                                                                           tCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCGGGATAAGGCGTTCACGCCGCATCCG‑GCCACAGATGCCTGATGTGACGCGGTGGCGTCTTATCATGCCTatc                                                     >  1:969089/1‑147 (MQ=38)
                                                                                            cgcgACGCTTCCACGCAAGGGCTGGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCG‑GCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAg                                                    <  1:250799/149‑1 (MQ=255)
                                                                                             gcgACGCTTCCACGCAAGGGCTTGTGGCGCATCTGCATAAGCGAAGGAAAGCTT‑AATCATCGTCGATATGGAGGCCGGAGAAGGCCGTCTCGCCGCAACCGCGGACCCGGT‑ACTGATGGCACGCGGTCGCGTCTTCACAGGCCTCCatg                                                   >  2:112187/1‑147 (MQ=255)
                                                                                             gcgACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCG‑GCAACCGATGCCTGCTGCGACGCGGTCGCGTCTTATCAGGCCTACAgg                                                   >  2:653118/1‑149 (MQ=255)
                                                                                                aCGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCG‑GCAACCGATGCCTGATGCGACGCGGGCGCGTCTTATCAGGCCTACAGGTCg                                                >  2:989199/1‑149 (MQ=255)
                                                                                                    tgccACGAGAGGGCGTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCG‑GCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGc                                            <  2:19652/147‑1 (MQ=25)
                                                                                                        aCGCAAGGGCTTGTAGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCG‑GCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCACGCCTCCAGGTAGATTACGat                                        >  1:130960/1‑149 (MQ=25)
                                                                                                         cGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCG‑GCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGTTCGATTCCGgta                                       >  2:23327/1‑149 (MQ=255)
                                                                                                         cGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCG‑GCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGata                                       <  2:909207/149‑1 (MQ=255)
                                                                                                          gCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCG‑GCAACCGATGCCTGATGCGACGCGGTcgc                                                                      >  2:198969/1‑117 (MQ=255)
                                                                                                          gCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCG‑GCAACCGATGCCTGATGCGACGCGGTcgc                                                                      <  1:198969/117‑1 (MQ=255)
                                                                                                          gCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCG‑GCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGt                                                  <  2:441266/137‑1 (MQ=255)
                                                                                                          gCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCG‑GCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGt                                                  >  1:441266/1‑137 (MQ=255)
                                                                                                          gCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCG‑GCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGatat                                      <  1:677461/149‑1 (MQ=255)
                                                                                                                  ttGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCAGTGGGGATATGGAGGCCCGACAAAGAGTAACGGCCCCTTGCT‑GTTACCGATGCCTGAGGCGAGCACGTCCCGTCTTATCGGGCCTACAGGTGGCGTCCGATACGTACAGCg                              >  2:162967/1‑149 (MQ=255)
                                                                                                                     tCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCG‑GCAACCGATACCTGATGCGACGCGGTc                                                                        <  2:592870/104‑1 (MQ=255)
                                                                                                                     tCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCG‑GCAACCGATACCTGATGCGACGCGGTc                                                                        >  1:592870/1‑104 (MQ=255)
                                                                                                                        cgcATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCG‑GCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATACAGATCTGTACATCGTATTCg                        <  2:497146/149‑1 (MQ=255)
                                                                                                                        cgcATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCG‑GCAACCGATGACTGACAGTGCTCGTTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCg                        >  1:488286/1‑149 (MQ=37)
                                                                                                                            tcggCAGGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCG‑GCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCaa                    <  2:635626/146‑1 (MQ=255)
                                                                                                                              tGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATCTGTAGGCCGGATAAGGCGTTCACGCCGCATCCG‑GCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAAtt                  <  1:415990/149‑1 (MQ=255)
                                                                                                                              ggCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCG‑GCAACCGATGCCTGATGCGACGCGGTCGAGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAAtt                  <  2:259550/148‑1 (MQ=255)
                                                                                                                                 aTGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCG‑GCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAAt               <  1:925532/149‑1 (MQ=255)
                                                                                                                                 aTGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCG‑GCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAAt               <  1:4258/149‑1 (MQ=255)
                                                                                                                                  tGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCG‑GCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATCCGGCAATTAATa              <  1:403203/149‑1 (MQ=255)
                                                                                                                                       gAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCG‑GCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATCTGTACATCGTATTCGGCAATTAATACATAg         >  2:858008/1‑149 (MQ=255)
                                                                                                                                            gaaGCTTG‑ATCATCGTCGATATGTACGCCGGATAAGGCGTTCACCCCGCATCCG‑GCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACga    <  1:917238/148‑1 (MQ=255)
                                                                                                                                              aGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCG‑GCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACgatt  <  1:164448/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
TCAATGCCTTCGATCAATGGGGCGTGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCG‑GCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATT  >  NC_000913_3_bme_pgi/4235259‑4235550

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