Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,825,0560GA43.9% 78.4 / 34.8 66intergenic (‑119/+84)ves/spycold‑ and stress‑inducible protein/periplasmic ATP‑independent protein refolding chaperone, stress‑induced
Reads supporting (aligned to +/- strand):  ref base G (29/8);  new base A (0/29);  total (29/37)
Fisher's exact test for biased strand distribution p-value = 1.50e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAA  >  NC_000913/1824909‑1825201
                                                                                                                                                   |                                                                                                                                                 
gACATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTgagaagagaa                                                                                                                                                 >  2:750685/1‑150 (MQ=255)
gACATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGagaaa                                                                                                                                                <  2:857836/151‑1 (MQ=255)
  cATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGGTGGTTATCGGCAACTCGGGTGATTTACGTTTGTGGCGATCAGGCTTTTTTTACAACTGTTTgagaagag                                                                                                                                                   >  1:431074/1‑146 (MQ=255)
   aTTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGagaaa                                                                                                                                                >  1:290502/1‑148 (MQ=255)
    ttttACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGagaaaaga                                                                                                                                             >  2:920234/1‑150 (MQ=255)
      ttACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGagaaaagaaa                                                                                                                                           >  2:856412/1‑150 (MQ=255)
       tACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAgaaaagaaaa                                                                                                                                          <  2:378729/150‑1 (MQ=255)
         cGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAAccg                                                                                                                                       <  1:1004792/151‑1 (MQ=255)
         cGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAAccg                                                                                                                                       <  2:339522/151‑1 (MQ=255)
           aatatAAAAGTATTCCATACTAGCTCCCGGACCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAAccgc                                                                                                                                      >  2:293360/2‑150 (MQ=255)
               taaaaGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACTGTGGATGGTTATCGGCAACTCGGGTGATTAACGTTAGTGGTGATCAGGTTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAAccgccg                                                                                                                                    <  2:227154/146‑1 (MQ=255)
                 aaaGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCt                                                                                                                               >  2:50215/1‑151 (MQ=255)
                        tCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCAcc                                                                                                                        >  1:264046/1‑151 (MQ=255)
                           aTACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGc                                                                                                                      >  1:879700/1‑150 (MQ=255)
                           aTACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGc                                                                                                                      >  1:932400/1‑150 (MQ=255)
                           aTACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGc                                                                                                                      >  2:810850/1‑150 (MQ=255)
                            tccTCGCTCCCGGCCTGTTTTTTTCTTTTTTATCTTTCAGTTACCCGGGGTTGTTTATCGGCATCGGGGTTGTTTTACTTTGTTGGTGACCAGTTTTTTTTTCCAACTGTTTAAGAAAAAAAAAGAAAACCGCCGATCCTGTCCACCg                                                                                                                       <  2:327851/146‑1 (MQ=255)
                              cTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCAtt                                                                                                                   >  1:572376/1‑150 (MQ=255)
                               tCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTa                                                                                                                  >  2:512983/1‑150 (MQ=255)
                                 gCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTg                                                                                                               >  2:709898/1‑151 (MQ=255)
                                   tCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGc                                                                                                              >  2:298809/1‑150 (MQ=255)
                                   tCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGc                                                                                                              >  2:930145/1‑150 (MQ=255)
                                    cccGGCCCTTTTTTTTCATTGTTATCTTTCAGGTACCGGTGGAGGGTTATCGGCAACTCGGGTGTTTTACTTTGGGGTGGATCAGGTTTTTTTTACAACTTTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCaa                                                                                                            <  1:628310/151‑1 (MQ=255)
                                     ccGCCCCTTTTTTTTTATTGTTTTCGTTCAGATACCCGTGGATGTTTACCGGCACCTGGGGTATTTTACGTTAGTGGTGATCAGTCTTTTTTTACACTTGTTGGGGAAGAAAAAAAAAATCCGCCGATCCTGTCCACCGCAtt                                                                                                                   <  2:891137/143‑1 (MQ=255)
                                       ggCCCGTTCTTTTCTTTTTTATCGTTCAGGTACCCGTGGATGGTTATCGGCACCTGGGGTGATTTACGTTAGTGGTGATCAGGTTTTTTTTACAACTGTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGt                                                                                                         <  1:958963/151‑1 (MQ=255)
                                            gTTTTTTTCATTGTTATCGTTCAGGTACCGGGGGTTGTTTATCGGCAACTGGGGTGATTTACGTTAGTGGTGATCAGGTTTTTTTTACAACTTTTGGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGt                                                                                                      <  1:437610/149‑1 (MQ=255)
                                             tttttttCATTGTTTTCGTTAAGGTACCCGTGGTTGTTTATCGGCACCTCGGTTGATTTATGTGAGTGGTAATCATGTTTTTTTTACAACTGTTTGAGAAGAAAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGt                                                                                                      <  1:748964/145‑1 (MQ=255)
                                              tCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGAc                                                                                                  >  2:183320/1‑151 (MQ=255)
                                                ttttCTTTTTTATCGTTCAGGAACCGGTGGATGGTTATCGGCAACTGGGGTGATTTACGTTAGTGGTGATCAGGTTTTTTTTACAACTTTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACa                                                                                                 <  1:431640/150‑1 (MQ=255)
                                                 tttCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACa                                                                                                 >  2:202535/1‑149 (MQ=255)
                                                 tttCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACa                                                                                                 >  1:15324/1‑149 (MQ=255)
                                                      tttttATCTTTCAGGTAACCGTGGATGGTTATCGGCAACGGGGTTGATTTACGTTAGGGGTGATCGGTCTTTTTTTACAATTTTTTGAAAAAAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGAcc                                                                                            <  2:166952/146‑1 (MQ=255)
                                                      ttgtttgctttCAGGTAACCGGGGATGGTTATCGGCACCTGGGTTGATTTACGTTAGTGGTGATCAGGTTTTTTTTACACCTGTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGAcc                                                                                            <  2:442015/140‑1 (MQ=255)
                                                         ttATCTTTCAGTTAACCGTGGATGGTTATCGGCAACTGGGGTGATTTGCGTTAGTGGTGATCAGGTTTTTTTTACAACTGTTTGAGAAGAAAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcggc                                                                                         <  1:610354/149‑1 (MQ=255)
                                                         ttATCGTTCAGGTAACCGTGGATGGTTATCGGCAACGCGGGTGATTTACGTTAGTGGTGATCAGGTTTTTTTTACAACTGTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcggcg                                                                                        <  1:380377/150‑1 (MQ=255)
                                                          tATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGt                                                                                      >  1:65193/1‑151 (MQ=255)
                                                               ttCAGGTAACCGTGGATGGTTATCGGCAACTGGGGTGATTTACGTTAGTGGTGATCAGGTTTTTTTTACAACTTTTTGAGAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCt                                                                                    <  1:889813/148‑1 (MQ=255)
                                                               ttCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCtt                                                                                   >  2:894240/1‑149 (MQ=255)
                                                                          gggtgtgtttgtgggCAACGCGGGTTATTGTCGTGAGTGGTGATCAGGTTTTTTTTACAATTTTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGttttt                                                                           <  1:1005594/133‑1 (MQ=255)
                                                                          gTGGATGGTTAGCGAAAACTCGGTTGATTTATGTTAGTGGTGATCAGGTTTTTTTTAAATCTGTTGGAAAAGAAAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGc                                                                       <  2:278237/150‑1 (MQ=255)
                                                                            ggATGGTTATCGGCAACTCGGGTGATTTACGTTAGGGGTGATCAGTCTTTTTTTACAACTGTTGGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGa                                                                    <  1:798198/151‑1 (MQ=255)
                                                                                  ttATCGGCAACTGGGTTAATTTACGTTGGTGGTGATCAGGTTTTTTTTACAACTGTTTGAAAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGat                                                               <  2:370213/150‑1 (MQ=255)
                                                                                     tCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGattatt                                                           >  1:950761/1‑151 (MQ=255)
                                                                                     tCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGattatt                                                           >  2:532919/1‑151 (MQ=255)
                                                                                       ggCAACTCGGGTGATTTACGTTAGTGGTGATCAGGTTTTTTTTACAACTTTTTGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGatta                                                             <  1:526014/147‑1 (MQ=255)
                                                                                          aaCTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTca                                                         >  2:113886/1‑148 (MQ=255)
                                                                                             tCGGGTGATTTACTTTAGTGGTGATCGGGCTTTTTTTACAACTGTTTAAAAAGAGAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTcagca                                                      <  1:416112/148‑1 (MQ=255)
                                                                                               gggTGATTTGCGTTATTGGTATTCGGCCTTTTTTAACAACTTTTTGGAAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGc                                                 <  2:121844/151‑1 (MQ=255)
                                                                                               gggTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGc                                                 >  1:426357/1‑151 (MQ=255)
                                                                                                     tttACTTTGTTGGTGATCGGTTTTTTTTCAAAATTGTTTGAGAAAAAAAAAAAAAGCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAg                                               <  1:750753/147‑1 (MQ=255)
                                                                                                      ttACTTTAGTGGTAATCAGGTTTTTTTAACAACTGTTGGAGAAGAGAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCAt                                           <  1:561194/150‑1 (MQ=255)
                                                                                                      ttACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGc                                             >  1:197629/1‑148 (MQ=255)
                                                                                                         ggTTAGTGGTGATCGGGTTTTTTTTACACTTGTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTa                                       <  2:641509/150‑1 (MQ=255)
                                                                                                           ttAGGGGTGATCAGTTTTTTTTAACAACTGTTGGAGAAGAAAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTAc                                      <  1:856494/150‑1 (MQ=255)
                                                                                                              gtggtgCTCCGGTTTTTTTTACAGCTTTTGGAGAAGAAAAAGGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCttt                                  <  1:252412/151‑1 (MQ=255)
                                                                                                              gtggtgATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGGTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCttt                                  >  2:824589/1‑151 (MQ=255)
                                                                                                              gtggtgATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGCCTTAAGTTTTTTGCCCGAAAGATTATTCAGCAGTTTCCGGCATTTTACCttt                                  >  2:695590/1‑151 (MQ=255)
                                                                                                                ggtgATCAGTTTTTTTTTACAACTGTTTGAAAAGAAAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCtt                                   <  1:810926/148‑1 (MQ=255)
                                                                                                                ggtgATCAGGTTTTTTTTACAACTGTGTGAAAAGAGAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTg                                <  1:833112/151‑1 (MQ=255)
                                                                                                                      cAGGTTTTTTTCAAATTTGTTTGAGAAGAAAAAAAAAAGCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCt                           <  1:426821/150‑1 (MQ=255)
                                                                                                                      cAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGAAACTGCAGGCATTTTTCCTTTTGCCGCt                           >  2:110633/1‑150 (MQ=255)
                                                                                                                      cAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGGTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCt                           >  2:14643/1‑150 (MQ=255)
                                                                                                                       aGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTgg                         >  1:603710/1‑151 (MQ=255)
                                                                                                                        ggttttttttACAACTTTTTGAGAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTgg                         <  2:912506/147‑1 (MQ=255)
                                                                                                                         gCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGAc                       >  2:603492/1‑151 (MQ=255)
                                                                                                                           tttttttACAACTGTTTGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACg                      <  1:124502/150‑1 (MQ=255)
                                                                                                                           tttttttACAACTGTTGAAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGa                        <  1:470616/148‑1 (MQ=255)
                                                                                                                             tttttACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGtt                    >  2:478553/1‑150 (MQ=255)
                                                                                                                               tttACACCTTTTTGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCt                  <  2:284129/150‑1 (MQ=255)
                                                                                                                                  aCAACTGTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCa              <  1:866881/151‑1 (MQ=255)
                                                                                                                                   cAACTGTTTGAGAAGAAAAAAGAAACCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTc               <  1:614922/149‑1 (MQ=255)
                                                                                                                                    acttGTTTGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAg             <  2:314708/147‑1 (MQ=255)
                                                                                                                                    aaCTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTc               >  1:790621/1‑148 (MQ=255)
                                                                                                                                               aagaaaaAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCaa  <  1:486773/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GACATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAA  >  NC_000913/1824909‑1825201

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