Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,693,3900TG44.0% 36.3 / 146.6 91K616T (AAA→ACA) bcsAcellulose synthase, catalytic subunit
Reads supporting (aligned to +/- strand):  ref base T (28/23);  new base G (25/15);  total (53/38)
Fisher's exact test for biased strand distribution p-value = 5.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.32e-01

AGCATCTCGGTTGGCGGGCCATAGAAGTAGCGCCAGATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGCCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGATATACCACGCCAGCACCGTTTCGTAGATTTCACTCCAGAAAGAGTGGCGATATTTGCCCTGGATCTTGGAGTTGGTCAGGCTGGCATGGATCATATGCGGC  >  NC_000913/3693240‑3693536
                                                                                                                                                      |                                                                                                                                                  
aGCATCTCGGTTGGCGGGCCATAGAAGTAGCGCCAGATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGCCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAAtt                                                                                                                                                    <  1:1003616/151‑1 (MQ=255)
aGCATCTCGGTTGGCGGGCCATAGAAGTAGCGCCAGATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGCCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAAtt                                                                                                                                                    <  2:675720/151‑1 (MQ=255)
       cGGTTGGCGGGCCATAGAAGTAGCGCCAGATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGCCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGcc                                                                                                                                                <  1:717583/148‑1 (MQ=255)
       cGGTTGGCGGGCCATAGAAGTAGCGCCAGATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGCCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGc                                                                                                                                                 <  1:213741/147‑1 (MQ=255)
       cGGTTGGCGGGCCATAGAAGTAGCGCCAGATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGCCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGACCACCTTTGGCGGTGACGTTAAATTTGCCttt                                                                                                                                             <  1:681359/151‑1 (MQ=255)
             gCGGGCCATAGAAGTAGCGCCAGATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGCCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATGTGCCTTTGTGCgg                                                                                                                                       <  1:499075/151‑1 (MQ=255)
              cGGGCCATAGAAGTAGCGCCAGATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGCCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCggg                                                                                                                                      >  2:467672/1‑151 (MQ=255)
                              cgcCAGATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGCCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGccacca                                                                                                                      <  2:1011429/151‑1 (MQ=255)
                                   gATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGCCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATGTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTc                                                                                                                 >  1:459927/1‑151 (MQ=255)
                                     tGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGCCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCgg                                                                                                               >  2:1084929/1‑151 (MQ=255)
                                      gCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGCCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCggt                                                                                                              <  2:145722/151‑1 (MQ=255)
                                       ccTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGCCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCggt                                                                                                              >  2:372976/1‑150 (MQ=255)
                                            cGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGCCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCg                                                                                                        <  2:31419/151‑1 (MQ=255)
                                            cGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGCCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCg                                                                                                        <  2:674742/151‑1 (MQ=255)
                                            cGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGCCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATGTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCg                                                                                                        >  1:1084734/1‑151 (MQ=255)
                                             gCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGCCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCg                                                                                                        >  1:1125277/1‑150 (MQ=255)
                                              cAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGCCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGa                                                                                                       <  1:1108825/150‑1 (MQ=255)
                                               aaCGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGCCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATGTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGata                                                                                                     <  2:645501/151‑1 (MQ=255)
                                               aaCGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGCCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATGTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGat                                                                                                      >  2:730584/1‑150 (MQ=255)
                                                aCGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGCCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGatat                                                                                                    <  2:159596/151‑1 (MQ=255)
                                                   cccACCAGGTTGAGCAGGACAAGGAAGATGTAGGGCCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATGTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGATATAcc                                                                                                 <  2:632987/151‑1 (MQ=255)
                                                    ccaccaGGTTGAGCAGGACAAGGAAGATGTAGGGCCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATGTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGATATAcc                                                                                                 >  1:1104507/1‑150 (MQ=255)
                                                     caccaGGTTGAGCAGGACAAGGAAGATGTAGGGCCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATGTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGATATACCAc                                                                                               <  2:1061337/151‑1 (MQ=255)
                                                       ccaGGTTGAGCAGGACAAGGAAGATGTAGGGCCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATGTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGATATACCACGc                                                                                             >  1:211594/1‑151 (MQ=255)
                                                        caGGTTGAGCAGGACAAGGAAGATGTAGGGCCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATGTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGATATACCACGcc                                                                                            >  1:828058/1‑151 (MQ=255)
                                                          ggTTGAGCAGGACAAGGAAGATGTAGGGCCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATGTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGATATACCACGCCa                                                                                           <  1:744754/150‑1 (MQ=255)
                                                              gAGCAGGACAAGGAAGATGTAGGGCCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGCGGTGCGATATACCACGCCAGCAcc                                                                                      >  2:922062/1‑151 (MQ=255)
                                                               aGCAGGACAAGGAAGATGTAGGGCCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGATATACCACGCCAGCACCg                                                                                     >  1:320460/1‑151 (MQ=255)
                                                               aGCAGGACAAGGAAGATGTAGGGCCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGATATACCACGCCAGCACCg                                                                                     >  1:97580/1‑151 (MQ=255)
                                                                gCAGGACAAGGAAGATGTAGGGCCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATGTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGATATACCACGCCAGCACCGt                                                                                    >  2:250592/1‑151 (MQ=255)
                                                                  aGGACAAGGAAGATGTAGCGCCGCGAGATCCCCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGATATACCACGCCAGCACCGttt                                                                                  <  2:224159/151‑1 (MQ=255)
                                                                   ggACAAGGAAGATGTAGGGCCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGATATACCACGCCAGCACCGTTTc                                                                                 >  1:1083688/1‑151 (MQ=255)
                                                                    gACAAGGAAGATGTAGGGCCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATGTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGATATACCACGCCAGCACCGTTTCg                                                                                >  2:303230/1‑151 (MQ=255)
                                                                      cAAGGAAGATGTAGGGCCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATGTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGATATACCACGCCAGCACCGTTTCGt                                                                               >  2:98404/1‑150 (MQ=255)
                                                                       aaGGAAGATGTAGGGCCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGATATACCACGCCAGCACCGTTTCGTAg                                                                             <  1:609800/151‑1 (MQ=255)
                                                                       aaGGAAGATGTAGGGCCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGATATACCACGCCAGCACCGTTTCGTAg                                                                             <  2:798154/151‑1 (MQ=255)
                                                                       aaGGAAGATGTAGGGCCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATGTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGATATACCACGCCAGCACCGTTTCGTa                                                                              <  2:273496/150‑1 (MQ=255)
                                                                          gAAGATGTAGGGCCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATGTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGATATACCACGCCAGCACCGTTTCGTAGAtt                                                                          >  2:694967/1‑151 (MQ=255)
                                                                           aaGATGTAGGGCCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGATATACCACGCCAGCACCGTTTCGTAGAttt                                                                         >  1:252592/1‑151 (MQ=255)
                                                                            aGATGTAGGGCCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGTCGGTGACGTTAATTGTGCCTTTGGGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGTTATACCACGCCAGCACCGTTTCGTAGATTTc                                                                        >  2:888439/1‑151 (MQ=255)
                                                                              aTGTAGGGCCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGATATACCACGCCAGCACCGTTTCGTAGATTTCAc                                                                      <  1:493546/151‑1 (MQ=255)
                                                                                  aGGGCCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATGTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGATATACCACGCCAGCACCGTTTCGTAGATTTCACTcc                                                                   >  2:1109569/1‑150 (MQ=255)
                                                                                  aGGGCCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATGTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGATATACCACGCCAGCACCGTTTCGTAGATTTCACTCCa                                                                  <  1:518900/151‑1 (MQ=255)
                                                                                    ggCCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATGTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGATATACCACGCCAGCACCGTTTCGTAGATTTCACTCCAGa                                                                <  2:420690/151‑1 (MQ=255)
                                                                                     gCCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGATATACCACGCCAGCACCGTTTCGTAGATTTCACTCCAGaa                                                               <  1:149753/151‑1 (MQ=255)
                                                                                       cgcgAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGATATACCACGCCAGCACCGTTTCGTAGATTTCACTCCAGaa                                                               >  2:89986/1‑149 (MQ=255)
                                                                                        gcgAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGATATACCACGCCAGCACCGTTTCGTAGATTTCACTCCAGaa                                                               >  2:588068/1‑148 (MQ=255)
                                                                                         cgAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGATATACCACGCCAGCACCGTTTCGTAGATTTCACTCCAGAAag                                                             >  2:153773/1‑149 (MQ=255)
                                                                                           agaTCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGATATACCACGCCAGCACCGTTTCGTAGATTTCACTCCAGAAAGAGTg                                                         >  2:613451/1‑151 (MQ=255)
                                                                                           agaTCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATGTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGATATACCACGCCAGCACCGTTTCGTAGATTTCACTCCAGAAAGAGTg                                                         >  1:344207/1‑151 (MQ=255)
                                                                                            gaTCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATGTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGATATACCACGCCAGCACCGTTTCGTAGATTTCACTCCAGAAAGAGTgg                                                        >  2:467341/1‑151 (MQ=255)
                                                                                            gaTCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATGTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGATATACCACGCCAGCACCGTTTCGTAGATTTCACTCCAGAAAGAGTgg                                                        >  2:765938/1‑151 (MQ=255)
                                                                                             aTCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGATATACCACGCCAGCACCGTTTCGTAGATTTCACTCCAGAAAGAGTGGc                                                       >  2:147707/1‑151 (MQ=255)
                                                                                                aCCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATGTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGATATACCACGCCAGCACCGTTTCGTAGATTTCACTCCAGAAAGAGTGGCGat                                                    >  2:85828/1‑151 (MQ=255)
                                                                                                 cccAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATGTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGATATACCACGCCAGCACCGTTTCGTAGATTTCACTCCAGAAAGAGTGGCGat                                                    <  2:424104/150‑1 (MQ=255)
                                                                                                  ccAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGATATACCACGCCAGCACCGTTTCGTAGATTTCACTCCAGAAAGAGTGGCGatat                                                  <  1:205191/151‑1 (MQ=255)
                                                                                                   cAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGGTATACCACGCCAGCACCGTTTCGTAGATTGCACTCCAGAAAGAGTGGCGatat                                                  >  2:40412/1‑150 (MQ=255)
                                                                                                    aGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGATATACCACGCCAGCACCGTTTCGTAGATTTCACTCCAGAAAGAGTGGCGATAttt                                                >  1:957171/1‑151 (MQ=255)
                                                                                                    aGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGATATACCACGCCAGCACCGTTTCGTAGATTTCACTCCAGAAAGAGGGGCGATAttt                                                >  1:582580/1‑151 (MQ=255)
                                                                                                      tCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGATATACCACGCCAGCACCGTTTCGTAGATTTCACTCCAGAAAGAGTGGCGATATTTGc                                              >  2:903146/1‑151 (MQ=255)
                                                                                                          cGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGATATACCACGCCAGCACCGTTTCGTAGATTTCACTCCAGAAAGAGTGGCGATATTTGCCCTg                                          >  2:981069/1‑151 (MQ=255)
                                                                                                          cGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATGTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGATATACCACGCCAGCACCGTTTCGTAGATTTCACTCCAGAAAGAGTGGCGATATTTGCCCt                                           <  1:1101347/150‑1 (MQ=255)
                                                                                                           gTACTCTTCTTCCACCAGTCCGCCTTTGGCGGTGACGTTAAATGTGCCTTTGTGCGGGTTAATCAGCGCCTCCAGCGTCGGTGGTGCGATATACCACGCCAGCACCGTTTCGTAGATTTCACTCCAGAAAGAGTGGCGATATTTGCCCTg                                          <  2:386764/150‑1 (MQ=255)
                                                                                                           gTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGATATACCACGCCAGCACCGTTTCGTAGATTTCACTCCAGAAAGAGTGGCGATATTTGCCCTgg                                         >  1:769183/1‑151 (MQ=255)
                                                                                                           gTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATGTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGATATACCACGCCAGCACCGTTTCGTAGATTTCACTCCAGAAAGAGTGGCGATATTTGCCCTgg                                         >  1:306025/1‑151 (MQ=255)
                                                                                                           gTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATGTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGATATACCACGCCAGCACCGTTTCGTAGATTTCACTCCAGAAAGAGTGGCGATATTTGCCCTgg                                         <  1:700414/151‑1 (MQ=255)
                                                                                                              ctctTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATGTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGATATACCACGCCAGCACCGTTTCGTAGATTTCACTCCAGAAAGAGTGGCGATATTTGCCCTGGATc                                      >  2:473336/1‑151 (MQ=255)
                                                                                                                  tcttcCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGATATACCACGCCAGCACCGTTTCGTAGATTTCACTCCAGAAAGAGTGGCGATATTTGCCCTGGATCTTgg                                  >  1:631511/1‑151 (MQ=255)
                                                                                                                     tcCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGATATACCACGCCAGCACCGTTTCGTAGATTTCACTCCAGAAAGAGTGGCGATATTTGCCCTGGATCTTgg                                  <  1:1113536/148‑1 (MQ=255)
                                                                                                                     tcCACCAGTCCACCTTTGGCGGTGACGTTAAATGTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGAGCGATATACCACGCCAGCACCGTTTCGTAGATTTCACTCCAGAAAGAGTGGCGATATTTGCCCTGGAGCTTGGAGt                               >  2:1130639/1‑151 (MQ=255)
                                                                                                                      ccaccaGTCCACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGATATACCACGCCAGCACCGTTTCGTAGATTTCACTCCAGAAAGAGTGGCGATATTTGCCCTGGATCTTGGAGtt                              <  2:631457/151‑1 (MQ=255)
                                                                                                                      ccaccaGTCCACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGATATACCACGCCAGCACCGTTTCGTAGATTTCACTCCAGAAAGAGTGGCGATATTTGCCCTGGATCTTGGAGtt                              <  2:522166/151‑1 (MQ=255)
                                                                                                                      ccaccaGTCCACCTTTGGCGGTGACGTTAAATGTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGATATACCACGCCAGCACCGTTTCGTAGATTTCACTCCAGAAAGAGTGGCGATATTTGCCCTGGATCTTGGAGtt                              >  1:1025926/1‑151 (MQ=255)
                                                                                                                         ccaGTCCACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGATATACCACGCCAGCACCGTTTCGTAGATTTCACTCCAGAAAGAGTGGCGATATTTGCCCTGGATCTTGGAGTTGGt                           >  1:376982/1‑151 (MQ=255)
                                                                                                                          caGTCCACCTTTGGCGGTGACGTTAAATGTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGATATACCACGCCAGCACCGTTTCGTAGATTTCACTCCAGAAAGAGTGGCGATATTTGCCCTGGATCTTGGAGTTGGTc                          >  2:83390/1‑151 (MQ=255)
                                                                                                                           aGTCCACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGATATACCACGCCAGCACCGTTTCGTAGATTTCACTCCAGAAAGAGTGGCGATATTTGCCCTGGATCTTGGAGTTGGTCa                         <  2:432198/151‑1 (MQ=255)
                                                                                                                            gTCCACCTTTGGCGGTGACGTTAAATGTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGATATACCACGCCAGCACCGTTTCGTAGATTTCACTCCAGAAAGAGTGGCGATATTTGCCCTGGATCTTGGAGTTGGTc                          >  2:386670/1‑149 (MQ=255)
                                                                                                                            gTCCACCTTTGGCGGTGACGTTAAATGTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGATATACCACGCCAGCACCGTTTCGTAGATTTCACTCCAGAAAGAGTGGCGATATTTGCCCTGGATCTTGGAGTTGGTCAg                        <  1:951406/151‑1 (MQ=255)
                                                                                                                             tCCACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGATATACCACGCCAGCACCGTTTCGTAGATTTCACTCCAGAAAGAGTGGCGATATTTGCCCTGGATCTTGGAGTTGGTCAgg                       >  2:1115704/1‑151 (MQ=255)
                                                                                                                                   tttGGCGGTGACGTTAAATGTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGATATACCACGCCAGCACCGTTTCGTAGATTTCACTCCAGAAAGAGTGGCGATATTTGCCCTGGATCTTGGAGTTGGTCAGGCTGGc                  >  2:265671/1‑150 (MQ=255)
                                                                                                                                     tGGCGGTGACGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGATATACCACGCCAGCACCGTTTCGTAGATTTCACTCCAGAAAGAGTGGCGATATTTGCCCTGGATCTTGGAGTTGGTCAGGCTGGCAt                <  1:976687/150‑1 (MQ=255)
                                                                                                                                     tGGCGGTGACGTTAAATGTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGATATACCACGCCAGCACCGTTTCGTAGATTTCACTCCAGAAAGAGTGGCGATATTTGCCCTGGATCTTGGAGTTGGTCAGGCTGGCAt                <  1:587718/150‑1 (MQ=255)
                                                                                                                                      ggCGGTGACGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGATATACCACGCCAGCACCGTTTCGTAGATTTCACTCCAGAAAGAGTGGCGATATTTGCCCTGGATCTTGGAGTTGGTCAGGCTGGCATg               <  1:1026157/150‑1 (MQ=255)
                                                                                                                                      ggCGGTGACGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGATATACCACGCCAGCACCGTTTCGTAGATTTCACTCCAGAAAGAGTGGCGATATTTGCCCTGGATCTTGGAGTTGGTCAGGCTGGCATg               <  1:921663/150‑1 (MQ=255)
                                                                                                                                       gCGGTGACGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGATATACCACGCCAGCACCGTTTCGTAGATTTCACTCCAGAAAGAGTGGCGATATTTGCCCTGGATCTTGGAGTTGGTCAGGCTGGCATgg              >  2:963920/1‑150 (MQ=255)
                                                                                                                                         ggTGACGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGATATACCACGCCAGCACCGTTTCGTAGATTTCACTCCAGAAAGAGTGGCGATATTTGCCCTGGATCTTGGAGTTGGTCAGGCTGGCATGGATc           >  2:548256/1‑151 (MQ=255)
                                                                                                                                              cGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGATATACCACGCCAGCACCGTTTCGTAGATTTCACTCCAGAAAGAGTGGCGATATTTGCCCTGGATCTTGGAGTTGGTCAGGCTGGCATGGATCATATg      >  1:948228/1‑151 (MQ=255)
                                                                                                                                              cGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGATATACCACGCCAGCACCGTTTCGTAGATTTCACTCCAGAAAGAGTGGCGATATTTGCCCTGGATCTTGGAGTTGGTCAGGCTGGCATGGATCATATg      <  2:1040357/151‑1 (MQ=255)
                                                                                                                                               gTTAAATTTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGATATACCACGCCAGCACCGTTTCGTAGATTTCACTCCAGAAAGAGTGGCGATATTTGCCCTGGATCTTGGAGTTGGTCAGGCTGGCATGGATCATATGc     >  2:962784/1‑151 (MQ=255)
                                                                                                                                                 tAAATGTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGATATACCACGCCAGCACCGTTTCGTAGATTTCACTCCAGAAAGAGTGGCGATATTTGCCCTGGATCTTGGAGTTGGTCAGGCTGGCATGGATCATATGCgg   >  2:581938/1‑151 (MQ=255)
                                                                                                                                                  aaaTGTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGATATACCACGCCAGCACCGTTTCGTAGATTTCACTCCAGAAAGAGTGGCGATATTTGCCCTGGATCTTGGAGTTGGTCAGGCTGGCATGGATCATATGCGgc  >  2:752409/1‑151 (MQ=255)
                                                                                                                                                  aaaTGTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGATATACCACGCCAGCACCGTTTCGTAGATTTCACTCCAGAAAGAGTGGCGATATTTGCCCTGGATCTTGGAGTTGGTCAGGCTGGCATGGATCATATGCGgc  <  1:154564/151‑1 (MQ=255)
                                                                                                                                                   aaTTTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGATATACCACGCCAGCACCGTTTCGTAGATTTCACTCCAGAAAGAGTGGCGATATTTGCCCTGGATCTTGGAGTTGGTCAGGCTGGCATGGATCATATGCGgc  >  1:277070/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                  
AGCATCTCGGTTGGCGGGCCATAGAAGTAGCGCCAGATGCCTACCGCAACGCCCACCAGGTTGAGCAGGACAAGGAAGATGTAGGGCCGCGAGATCACCCAGTCGACGTACTCTTCTTCCACCAGTCCACCTTTGGCGGTGACGTTAAATTTGCCTTTGTGCGGGTTAATCAGCGCCACCAGCGTCGGTGGTGCGATATACCACGCCAGCACCGTTTCGTAGATTTCACTCCAGAAAGAGTGGCGATATTTGCCCTGGATCTTGGAGTTGGTCAGGCTGGCATGGATCATATGCGGC  >  NC_000913/3693240‑3693536

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