Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,235,477 G→A 5.1% intergenic (+70/‑429) pgi → / → yjbE glucosephosphate isomerase/extracellular polysaccharide production threonine‑rich protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,235,4770GA 5.1% 581.6 / 13.3 176intergenic (+70/‑429)pgi/yjbEglucosephosphate isomerase/extracellular polysaccharide production threonine‑rich protein
Reads supporting (aligned to +/- strand):  ref base G (79/88);  new base A (5/4);  total (84/92)
Fisher's exact test for biased strand distribution p-value = 7.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.86e-01

Only 100 of 177 total aligned reads displayed.

AAAGATGATAAAGAAATCAGCAGCCACGATAGCTCGACCAATGGTCTGATTAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTT  >  NC_000913/4235330‑4235625
                                                                                                                                                   |                                                                                                                                                    
aaaGATGATAAAGAAATCAGCAGCCACGATAGCTCGACCAATGGTCTGATTAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGt                                                                                                                                                     >  3:754468/1‑149 (MQ=255)
aaaGATGATAAAGAAATCAGCAGCCACGATAGCTCGACCAATGGTCTGATTAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGt                                                                                                                                                     >  2:503703/1‑149 (MQ=255)
  aGATGATAAAGAAATCAGCAGCCACGATAGCTCGACCAATGGTCTGATTAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTcg                                                                                                                                                   >  3:1440566/1‑149 (MQ=255)
   gatgatAAAGAAATCAGCAGCCACGATAGCTCGACCAATGGTCTGATTAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTcgc                                                                                                                                                  <  4:194324/149‑1 (MQ=255)
       atAAAGAAATCAGCAGCCACGATAGCTCGACCAATGGTCTGATTAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCt                                                                                                                                              <  3:238416/149‑1 (MQ=255)
        tAAAGAAATCAGCAGCCACGATAGCTCGACCAATGGTCTGATTAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCtt                                                                                                                                             <  3:1053374/149‑1 (MQ=255)
        tAAAGAAATCAGCAGCCACGATAGCTCGACCAATGGTCTGATTAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCtt                                                                                                                                             <  4:1156582/149‑1 (MQ=255)
         aaaGAAATCAGCAGCCACGATAGCTCGACCAATGGTCTGATTAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTa                                                                                                                                            <  4:526814/149‑1 (MQ=255)
           aGAAATCAGCAGCCACGATAGCTCGACCAATGGTCTGATTAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATc                                                                                                                                          >  1:38965/1‑149 (MQ=255)
           aGAAATCAGCAGCCACGATAGCTCGACCAATGGTCTGATTAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATc                                                                                                                                          <  1:1045775/149‑1 (MQ=255)
            gAAATCAGCAGCCACGATAGCTCGACCAATGGTCTGATTAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCa                                                                                                                                         <  2:139899/149‑1 (MQ=255)
             aaaTCAGCAGCCACGATAGCTCGACCAATGGTCTGATTAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATAAGACGCGACCGCGTcgcatcag                                                                                                                                        <  4:794773/149‑8 (MQ=38)
               aTCAGCAGCCACGATAGCTCGACCAATGGTCTGATTAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGc                                                                                                                                      <  4:611740/149‑1 (MQ=255)
                tCAGCAGCCACGATAGCTCGACCAATGGTCTGATTAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGcc                                                                                                                                     >  1:873767/1‑149 (MQ=255)
                   gcagcCACGATAGCTCGACCAATGGTCTGATTAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTAc                                                                                                                                  <  2:862903/149‑1 (MQ=255)
                       ccACGATAGCTCGACCAATGGTCTGATTAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAg                                                                                                                                        >  4:1195594/1‑139 (MQ=255)
                       ccACGATAGCTCGACCAATGGTCTGATTAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGt                                                                                                                              <  4:1085894/149‑1 (MQ=255)
                       ccACGATAGCTCGACCAATGGTCTGATTAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGt                                                                                                                              <  1:559097/149‑1 (MQ=255)
                         aCGATAGCTCGACCAATGGTCTGATTAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCg                                                                                                                            <  1:1357995/149‑1 (MQ=255)
                         aCGATAGCTCGACCAATGGTCTGATTAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATAAGACGCGACCGCGTCGCATCAGGcatcgg                                                                                                                                 <  4:1432974/144‑6 (MQ=38)
                          cGATAGCTCGACCAATGGTCTGATTAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTAc                                                                                                                                  <  3:1419615/142‑1 (MQ=255)
                          cGATAGCTCGACCAATGGTCTGATTAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTAc                                                                                                                                  >  4:1419615/1‑142 (MQ=255)
                           gATAGCTCGACCAATGGTCTGATTAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGAt                                                                                                                          >  2:1246353/1‑149 (MQ=255)
                              aGCTCGACCAATGGTCTGATTAACCGCTATAAAGCGTGGCGCGTTTAAGCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTAc                                                                                                                                  <  2:263894/138‑1 (MQ=255)
                                 tCGACCAATGGTCTGATTAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGat                                                                                                                    >  4:736934/1‑149 (MQ=255)
                                 tCGACCAATGGTCTGATTAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATAAGACGCGACCGCGTcgcatc                                                                                                                                          >  4:262784/1‑122 (MQ=37)
                                  cGACCAATGGTCTGATTAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGata                                                                                                                   <  3:802186/149‑1 (MQ=255)
                                   gACCAATGGTCTGATTAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGatat                                                                                                                  <  3:45789/149‑1 (MQ=255)
                                   gACCAATGGTCTGATTAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGatat                                                                                                                  >  3:1042635/1‑149 (MQ=255)
                                     ccATTGGTCTGATTAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGt                                                                                                                <  3:736934/149‑1 (MQ=255)
                                      cAATGGTCTGATTAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTa                                                                                                               >  2:54426/1‑149 (MQ=255)
                                        aTGGTCTGATTAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACa                                                                                                             >  4:1479243/1‑149 (MQ=255)
                                          ggTCTGATTAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCATGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATc                                                                                                           <  3:917938/149‑1 (MQ=255)
                                            tCTGATTAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAg                                                                                                                                        <  2:1363329/118‑1 (MQ=255)
                                              tGATTAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTAt                                                                                                       >  2:1177596/1‑149 (MQ=255)
                                              tGATTAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTAt                                                                                                       >  2:316314/1‑149 (MQ=255)
                                               gATTAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGata                                                                                                                   >  4:1197412/1‑136 (MQ=255)
                                               gATTAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGata                                                                                                                   <  3:1197412/136‑1 (MQ=255)
                                                aTTAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTc                                                                                                     >  1:1134829/1‑149 (MQ=255)
                                                aTTAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTc                                                                                                     >  4:401548/1‑149 (MQ=255)
                                                aTTAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTc                                                                                                     >  1:1199921/1‑149 (MQ=255)
                                                  tAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCgg                                                                                                   <  3:1409616/149‑1 (MQ=255)
                                                    aCCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCa                                                                                                 >  4:338463/1‑149 (MQ=255)
                                                     ccGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCaa                                                                                                <  1:928481/149‑1 (MQ=255)
                                                      cGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATAAGACGCGACCGCGTCGCATCa                                                                                                                                         <  2:282168/107‑1 (MQ=25)
                                                       gCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAAtt                                                                                              >  2:938940/1‑149 (MQ=255)
                                                         tataAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTaa                                                                                            <  2:798394/149‑1 (MQ=255)
                                                          ataAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAAt                                                                                           <  2:107550/149‑1 (MQ=255)
                                                           taAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATa                                                                                          >  3:883703/1‑149 (MQ=255)
                                                            aaaGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATa                                                                                          >  3:134397/1‑148 (MQ=255)
                                                            aaaGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCATTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATa                                                                                          <  4:134397/148‑1 (MQ=255)
                                                              aGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACAt                                                                                       <  4:1440566/149‑1 (MQ=255)
                                                               gCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACaca                                                                                      <  4:1017013/149‑3 (MQ=255)
                                                               gCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATa                                                                                      >  3:152016/1‑149 (MQ=255)
                                                               gCGTGGCGCGGTTAACCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATa                                                                                      >  2:1014159/1‑149 (MQ=255)
                                                                 gTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAAtt                                                                                              >  4:673110/1‑139 (MQ=255)
                                                                 gTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAAtt                                                                                              <  3:673110/139‑1 (MQ=255)
                                                                         gTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACgattga                                                                            >  1:31862/1‑149 (MQ=255)
                                                                          ttAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACgattgat                                                                           <  2:1030069/149‑1 (MQ=255)
                                                                           tAATTACCGTCGTTATGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAACCGATGCCGGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACTTAGCACggtttatt                                                                          <  2:1472844/149‑1 (MQ=25)
                                                                           tAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACgattgatt                                                                          <  2:1134829/149‑1 (MQ=255)
                                                                              tcatcGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTaaa                                                                       <  3:278126/149‑1 (MQ=255)
                                                                                atcGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATa                                                                     >  1:231845/1‑149 (MQ=255)
                                                                                 tcgtcgATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTAt                                                                                                       >  2:1419503/1‑114 (MQ=255)
                                                                                 tcgtcgATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTAt                                                                                                       <  1:1419503/114‑1 (MQ=255)
                                                                                 tcgtcgATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCg                                                                                                    >  1:286376/1‑117 (MQ=255)
                                                                                 tcgtcgATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATaa                                                                    >  3:458132/1‑149 (MQ=255)
                                                                                   gtcgATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAAcc                                                                  <  1:171195/149‑1 (MQ=255)
                                                                                   gtcgATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAAcc                                                                  <  2:253371/149‑1 (MQ=255)
                                                                                       atatGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTtaa                                                              >  4:989513/1‑149 (MQ=255)
                                                                                        tatGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTtaat                                                             <  3:1041313/149‑1 (MQ=255)
                                                                                               gCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAAc                                                                   <  3:1428288/136‑1 (MQ=255)
                                                                                                  ggATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCg                                                   >  4:545430/1‑149 (MQ=255)
                                                                                                   gATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGa                                                  <  4:1029961/149‑1 (MQ=255)
                                                                                                      aaGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCAGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGt                                               <  2:125225/149‑1 (MQ=255)
                                                                                                        ggCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTa                                                                                                        <  3:229928/90‑1 (MQ=255)
                                                                                                        ggCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTa                                             >  2:749517/1‑149 (MQ=255)
                                                                                                        ggCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTa                                             <  1:503703/149‑1 (MQ=255)
                                                                                                         gCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTAt                                            >  4:251053/1‑149 (MQ=255)
                                                                                                           gTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTtaat                                                             <  3:4004/130‑1 (MQ=255)
                                                                                                            ttCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTc                                         <  1:1435409/149‑1 (MQ=255)
                                                                                                               aCGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGc                                      >  2:103086/1‑149 (MQ=255)
                                                                                                                  ccgcATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCATATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACATAACGTCGGCAtt                                   <  3:1066971/149‑1 (MQ=38)
                                                                                                                     cATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCAttttt                                >  3:261962/1‑149 (MQ=255)
                                                                                                                     cATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCAttttt                                >  4:733574/1‑149 (MQ=255)
                                                                                                                      aTCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCg                                        >  1:212142/1‑140 (MQ=255)
                                                                                                                         cGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTAt                            <  1:663286/149‑1 (MQ=255)
                                                                                                                          ggCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATa                                                                                      >  4:1054059/1‑90 (MQ=255)
                                                                                                                          ggCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACATAACGTCGGCAttgttattaa                           <  2:609198/149‑6 (MQ=39)
                                                                                                                                   tGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTtaat                                                             <  1:1396238/106‑1 (MQ=255)
                                                                                                                                    gCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATcc                 <  2:892807/149‑1 (MQ=255)
                                                                                                                                    gCCTGATAAGACGCGACCGCGTCGCATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATcc                 >  4:469108/1‑149 (MQ=255)
                                                                                                                                       tGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGAGTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCt              >  4:836792/1‑149 (MQ=255)
                                                                                                                                        gATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCtt             <  2:231845/149‑1 (MQ=255)
                                                                                                                                        gATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCtt             >  3:1086557/1‑149 (MQ=255)
                                                                                                                                             gACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTg        >  1:525526/1‑149 (MQ=255)
                                                                                                                                               cgcgGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTa      <  2:19252/149‑1 (MQ=255)
                                                                                                                                                gcgGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTaa     >  1:890407/1‑149 (MQ=255)
                                                                                                                                                gcgGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTaa     <  2:445745/149‑1 (MQ=255)
                                                                                                                                                   gTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAAttt  >  3:1199270/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
AAAGATGATAAAGAAATCAGCAGCCACGATAGCTCGACCAATGGTCTGATTAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTT  >  NC_000913/4235330‑4235625

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: