Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,235,478 T→C 5.2% intergenic (+71/‑428) 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,4780TC 5.2% 568.2 / 13.7 174intergenic (+71/‑428)pgi/yjbEglucosephosphate isomerase/extracellular polysaccharide production threonine‑rich protein
Reads supporting (aligned to +/- strand):  ref base T (77/88);  new base C (5/4);  total (82/92)
Fisher's exact test for biased strand distribution p-value = 7.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.10e-01

Only 100 of 178 total aligned reads displayed.

AAAGATGATAAAGAAATCAGCAGCCACGATAGCTCGACCAATGGTCTGATTAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTA  >  NC_000913/4235330‑4235626
                                                                                                                                                    |                                                                                                                                                    
aaaGATGATAAAGAAATCAGCAGCCACGATAGCTCGACCAATGGTCTGATTAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGt                                                                                                                                                      >  2:503703/1‑149 (MQ=255)
aaaGATGATAAAGAAATCAGCAGCCACGATAGCTCGACCAATGGTCTGATTAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGt                                                                                                                                                      >  2:928481/1‑149 (MQ=255)
   gatgatAAAGAAATCAGCAGCCACGATAGCTCGACCAATGGTCTGATTAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTcgc                                                                                                                                                   <  3:1042839/149‑1 (MQ=255)
   gatgatAAAGAAATCAGCAGCCACGATAGCTCGACCAATGGTCTGATTAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTcgc                                                                                                                                                   <  4:194324/149‑1 (MQ=255)
        tAAAGAAATCAGCAGCCACGATAGCTCGACCAATGGTCTGATTAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCtt                                                                                                                                              <  4:1156582/149‑1 (MQ=255)
        tAAAGAAATCAGCAGCCACGATAGCTCGACCAATGGTCTGATTAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATACGGCAACCGATGCCTGATGCGACGCGGTCGCGTCtt                                                                                                                                              >  3:1314512/1‑149 (MQ=255)
           aGAAATCAGCAGCCACGATAGCTCGACCAATGGTCTGATTAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATc                                                                                                                                           <  1:1045775/149‑1 (MQ=255)
           aGAAATCAGCAGCCACGATAGCTCGACCAATGGTCTGATTAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATc                                                                                                                                           >  1:38965/1‑149 (MQ=255)
            gAAATCAGCAGCCACGATAGCTCGACCAATGGTCTGATTAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCa                                                                                                                                          >  3:1085894/1‑149 (MQ=255)
            gAAATCAGCAGCCACGATAGCTCGACCAATGGTCTGATTAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCa                                                                                                                                          <  2:139899/149‑1 (MQ=255)
             aaaTCAGCAGCCACGATAGCTCGACCAATGGTCTGATTAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATAAGACGCGACCGCGTcgcatcag                                                                                                                                         <  4:794773/149‑8 (MQ=38)
             aaaTCAGCAGCCACGATAGCTCGACCAATAGTCTGATTAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAg                                                                                                                                         <  4:518005/149‑1 (MQ=255)
               accagcagCCACGATAGCTCGACCAATGGTCTGATTAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGGTATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGc                                                                                                                                       <  2:604033/147‑1 (MQ=255)
               aTCAGCAGCCACGATAGCTCGACCAATGGTCTGATTAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGc                                                                                                                                       <  4:611740/149‑1 (MQ=255)
                tCAGCAGCCACGATAGCTCGACCAATGGTCTGATTAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGcc                                                                                                                                      >  1:873767/1‑149 (MQ=255)
                   gcagcCACGATAGCTCGACCAATGGTCTGATTAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTAc                                                                                                                                   >  1:892807/1‑149 (MQ=255)
                      gcCACGATAGCTCGACCAATGGTCTGATTAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAgg                                                                                                                                >  2:272190/1‑149 (MQ=255)
                       ccACGATAGCTCGACCAATGGTCTGATTAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGt                                                                                                                               <  1:559097/149‑1 (MQ=255)
                       ccACGATAGCTCGACCAATGGTCTGATTAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGt                                                                                                                               <  4:1085894/149‑1 (MQ=255)
                         aCGATAGCTCGACCAATGGTCTGATTAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCg                                                                                                                             <  1:1357995/149‑1 (MQ=255)
                         aCGATAGCTCGACCAATGGTCTGATTAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATAAGACGCGACCGCGTCGCATCAGGcatcgg                                                                                                                                  >  3:1432974/1‑139 (MQ=37)
                          cGATAGCTCGACCAATGGTCTGATTAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTAc                                                                                                                                   >  4:1419615/1‑142 (MQ=255)
                          cGATAGCTCGACCAATGGTCTGATTAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTAc                                                                                                                                   <  3:1419615/142‑1 (MQ=255)
                           gATAGCTCGACCAATGGTCTGATTAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGAt                                                                                                                           >  2:1246353/1‑149 (MQ=255)
                              aGCTCGACCAATGGTCTGATTAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTAc                                                                                                                                   >  1:263894/1‑138 (MQ=255)
                                 tCGACCAATGGTCTGATTAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGat                                                                                                                     >  4:736934/1‑149 (MQ=255)
                                 tCGACCAATGGTCTGATTAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATAAGACGCGACCGCGTcgcatc                                                                                                                                           <  3:262784/127‑6 (MQ=37)
                                 tCGACCAATGGTCTGATTAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATAAGACGCGACCGCGTcgcatc                                                                                                                                           >  4:262784/1‑122 (MQ=37)
                                  cGACCAATGGTCTGATTAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGata                                                                                                                    <  3:802186/149‑1 (MQ=255)
                                   gACCAATGGTCTGATTAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGatat                                                                                                                   <  3:45789/149‑1 (MQ=255)
                                     ccATTGGTCTGATTAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGt                                                                                                                 <  3:736934/149‑1 (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:316314/1‑149 (MQ=255)
                                              tGATTAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTAt                                                                                                        >  2:171195/1‑149 (MQ=255)
                                              tGATTAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTAt                                                                                                        >  2:1177596/1‑149 (MQ=255)
                                               gATTAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGata                                                                                                                    >  4:1197412/1‑136 (MQ=255)
                                               gATTAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGata                                                                                                                    <  3:1197412/136‑1 (MQ=255)
                                                aTTAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTc                                                                                                      >  4:401548/1‑149 (MQ=255)
                                                aTTAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTc                                                                                                      >  1:1134829/1‑149 (MQ=255)
                                                aTTAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTc                                                                                                      >  1:1199921/1‑149 (MQ=255)
                                                  tAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCgg                                                                                                    <  3:1409616/149‑1 (MQ=255)
                                                    aCCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCa                                                                                                  >  3:73708/1‑149 (MQ=255)
                                                    aCCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCa                                                                                                  >  4:338463/1‑149 (MQ=255)
                                                     ccGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCaa                                                                                                 >  2:104990/1‑149 (MQ=255)
                                                     ccGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCaa                                                                                                 <  1:928481/149‑1 (MQ=255)
                                                       gCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAAtt                                                                                               >  2:938940/1‑149 (MQ=255)
                                                       gCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAAtt                                                                                               <  3:1460189/149‑1 (MQ=255)
                                                         tataAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTaa                                                                                             <  2:798394/149‑1 (MQ=255)
                                                          ataAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAAt                                                                                            <  2:107550/149‑1 (MQ=255)
                                                           taAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATa                                                                                           >  3:883703/1‑149 (MQ=255)
                                                            aaaGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCATTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATa                                                                                           <  4:134397/148‑1 (MQ=255)
                                                              aGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACAt                                                                                        <  4:1440566/149‑1 (MQ=255)
                                                               gCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATa                                                                                       <  3:1364573/149‑1 (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)
                                                                           tAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACgattgatt                                                                           <  2:1134829/149‑1 (MQ=255)
                                                                              tcatcGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTaaa                                                                        >  3:994726/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)
                                                                                 tcgtcgATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCg                                                                                                     <  2:286376/117‑1 (MQ=255)
                                                                                   gtcgATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAAcc                                                                   <  4:1181419/149‑1 (MQ=255)
                                                                                   gtcgATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAAcc                                                                   <  2:253371/149‑1 (MQ=255)
                                                                                   gtcgATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAAcc                                                                   <  1:171195/149‑1 (MQ=255)
                                                                                       atatGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTtaa                                                               >  4:989513/1‑149 (MQ=255)
                                                                                        tatGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTtaat                                                              <  3:1041313/149‑1 (MQ=255)
                                                                                               gCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAAc                                                                    >  4:1428288/1‑136 (MQ=255)
                                                                                                  ggATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCg                                                    >  4:545430/1‑149 (MQ=255)
                                                                                                   gATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGa                                                   <  4:1029961/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)
                                                                                                                     cATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCAttttt                                 >  4:733574/1‑149 (MQ=255)
                                                                                                                      aTCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCg                                         >  1:212142/1‑140 (MQ=255)
                                                                                                                      aTCAGGCATCGGTTGCCGGATGCGGCGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCAtttttt                                >  3:181360/1‑149 (MQ=255)
                                                                                                                         cGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTAt                             <  1:663286/149‑1 (MQ=255)
                                                                                                                          ggCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATa                                                                                       <  3:1054059/90‑1 (MQ=255)
                                                                                                                          ggCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATa                                                                                       >  4:1054059/1‑90 (MQ=255)
                                                                                                                          ggCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCAttttttatca                            <  2:873767/149‑5 (MQ=255)
                                                                                                                          ggCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATc                            >  2:1338819/1‑149 (MQ=255)
                                                                                                                          ggCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACATAACGTCGGCAttgttattaa                            <  2:609198/149‑6 (MQ=39)
                                                                                                                                cGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATaa                      >  3:373050/1‑149 (MQ=255)
                                                                                                                                   tGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTtaat                                                              <  1:1396238/106‑1 (MQ=255)
                                                                                                                                    gCCTGATAAGACGCGACCGCGTCGCATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATcc                  >  4:469108/1‑149 (MQ=255)
                                                                                                                                       tGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGAGTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCt               >  4:836792/1‑149 (MQ=255)
                                                                                                                                        gATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCtt              >  3:1086557/1‑149 (MQ=255)
                                                                                                                                        gATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCtt              <  2:231845/149‑1 (MQ=255)
                                                                                                                                             gACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTg         >  1:525526/1‑149 (MQ=255)
                                                                                                                                               cgcgGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTa       <  2:19252/149‑1 (MQ=255)
                                                                                                                                                gcgGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTaa      <  2:445745/149‑1 (MQ=255)
                                                                                                                                                gcgGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTaa      <  3:251053/149‑1 (MQ=255)
                                                                                                                                                    tCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTa  >  3:519254/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AAAGATGATAAAGAAATCAGCAGCCACGATAGCTCGACCAATGGTCTGATTAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTA  >  NC_000913/4235330‑4235626

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