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

CCGACATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAA  >  NC_000913/1824907‑1825202
                                                                                                                                                      |                                                                                                                                                  
ccGACATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTgagaagaga                                                                                                                                                   <  1:990618/151‑1 (MQ=255)
   aCATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGGGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGagaaa                                                                                                                                                 >  2:1091148/1‑150 (MQ=255)
    cATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGagaaa                                                                                                                                                 >  1:128947/1‑149 (MQ=255)
     aTTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGagaaa                                                                                                                                                 >  2:1063395/1‑148 (MQ=255)
          aCGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAAcc                                                                                                                                         <  1:941989/151‑1 (MQ=255)
          aCGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAAcc                                                                                                                                         <  2:1168519/151‑1 (MQ=255)
            gCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAAccgc                                                                                                                                       <  1:47502/151‑1 (MQ=255)
                atCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAAAAAACCGCCGa                                                                                                                                    <  2:768136/150‑1 (MQ=255)
                    aaGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTg                                                                                                                               >  2:738891/1‑151 (MQ=255)
                     aGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGt                                                                                                                              >  1:330175/1‑151 (MQ=255)
                           caatccTCGCTCCCGGCCGGTTTTTTTCTTTTTTATCTTTCAGGTACCCGTGTTTGTTTTCCGCCACCTGGTTTTTTTTACTTTAGTGGTGATCGTTTTTTTTTTCCACCTTTTTGAAAAAAAAAAAAAAAACCGCCGATCCTGTCCAcc                                                                                                                         <  2:917884/145‑1 (MQ=255)
                            cATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGc                                                                                                                       >  2:371499/1‑151 (MQ=255)
                            cATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGc                                                                                                                       >  1:991269/1‑151 (MQ=255)
                               aCTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCAtt                                                                                                                    >  1:1001272/1‑151 (MQ=255)
                                cTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGGGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCAtt                                                                                                                    >  1:361383/1‑150 (MQ=255)
                                    cccccGGCCCTTTTTTTTCATTTTTATCGTTAAGGTAACCGGGGATGGTTATCGGCAACTGGGGTGATTTACGTTAGTGGTGATCAGGTTTTTTTTACAACTGTTTGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTg                                                                                                                <  1:473727/148‑1 (MQ=255)
                                     tCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGc                                                                                                               >  2:317763/1‑150 (MQ=255)
                                      cccGGCCCTTTTTTTTCATTTTTATCGTTCAGGTACCCGTGGATGTTTATCGGCAACTCGGGTAATTTACTTTAGTGGTGATCAGGTTTTTTTTACAACTTTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTa                                                                                                                   <  2:471042/145‑1 (MQ=255)
                                       ccGGCCCGTTCTTTTCATTTTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTGGGGTGATTTACGTTAGGGGTGATCAGGTTTTTTTTACACCTGTTGGAAAAGAGAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAg                                                                                                            <  1:881010/151‑1 (MQ=255)
                                           cccGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAgg                                                                                                           >  2:306406/1‑148 (MQ=255)
                                            ccGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAg                                                                                                        >  2:1094703/1‑150 (MQ=255)
                                              gTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTgg                                                                                                     >  1:31664/1‑151 (MQ=255)
                                                  ttttCTTTTTTATCTTTCAGGAACCCGTGGTGGTTTATCGGCACCTCGGTTATTTTCCGTTAGTGGTGATCGGGTTTTTTTTCAACTTTTTGAAAAAGAGAAAAAAAACCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACa                                                                                                  <  2:956800/150‑1 (MQ=255)
                                                  ttttCATTTTTATCTTTCAGGTACCCGTGGATGGTTATCGGCACCTGGGGTGATTTACGTTAGTGGTGATCAGGTTTTTTTTACAACTTTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACa                                                                                                  <  1:1042405/150‑1 (MQ=255)
                                                  ttttCATTTTTATCTTTCAGGTACCCGTGGATGGTTATCGGCAACTCGGGTGTTTTACGTTAGTGGTGATCGGGTTTTTTTAACAACTGTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGAc                                                                                                   <  1:190141/149‑1 (MQ=255)
                                                  ttttCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACa                                                                                                  >  2:1177439/1‑150 (MQ=255)
                                                        ttGTTTTCTTTCAGGTAACGGGGGATGGTTATCGGCAACTGGGGTGATTTACGTTGGTGGTGATCAGGCTTTTTTTAAATTTGTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGAcc                                                                                             <  1:1159610/149‑1 (MQ=255)
                                                         ttttATCTTTCAGGTACCGGGGGTTGTTTATCGCCACCTGGGGTGTTTTAGGTTATTGGTATTCGGTTTTTTTTAACAATTTTTTGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcgg                                                                                           <  1:474452/148‑1 (MQ=255)
                                                          gTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcggcg                                                                                         >  1:1090595/1‑151 (MQ=255)
                                                           ttatttttCAGGTAACGGTGTTTGTTTATGGGCAACGGGGGTATTTTACGTTAGGGGTGATCAGTTTTTTTTTACACCTGTTGGAAAAGAAAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcggcg                                                                                         <  2:106634/144‑1 (MQ=255)
                                                           ttATCTTTCGGTTACCGGGGGATGTTTATCGAAAACTGGGGTGTTTTACTTTAGTGGTGATCGGGTTTTTTTTACAACTGTTTGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcggc                                                                                          <  1:1092259/149‑1 (MQ=255)
                                                           ttATCGTTCAGGTAACGGGGGTTGTTTATCGGCAACGGGGGTGTTTTATTTTAGTGGGTATCAGTTTTTTTTAACACTTGTTGGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcggcg                                                                                         <  2:905324/150‑1 (MQ=255)
                                                            tatgtttCAGGTACCGGGGGATGTTTATCGGCAACTGGGGTGTTTTACGTTAGGGGTGATCGGTTTTTTTTTACACTTGTTGGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcggcgg                                                                                        <  1:799067/145‑1 (MQ=255)
                                                            tATTGTTCGGGTAACGGTGGTGGGTTAGCGGAAACTGGGGTGATTTACGTTAGTGGTGATCAGCTTTTTTTTACAACTGTTAGAAAAGAGAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcggcgg                                                                                        <  1:353281/150‑1 (MQ=255)
                                                               cgttaggttaacgttGGATGTTTTTCGTCACCTTGGGTAATTTGCGTATGGGGTGATCAGGTTTTTTTTAAAACTTTTTGAAAAAAGAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCtt                                                                                    <  1:483925/137‑1 (MQ=255)
                                                               cGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCtt                                                                                    >  1:249238/1‑151 (MQ=255)
                                                                 ttCAGGTACCCGTGGATGGTTACCGGCAACTGGGTTGATTTACGTTAGTGGTGATCAGTTTTTTTTTACACCTTTTGGAGAAGAAAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTa                                                                                   <  1:993478/150‑1 (MQ=255)
                                                                   caggaaccggTGGATGTTTATCGCCAACTCGGTTGATTTCGGTTAGTGGTGACCGGGTTTTTTTTCCAACTGTTTGAGAAGAAAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAg                                                                                 <  2:1094333/141‑1 (MQ=255)
                                                                     ggaaacgggggATGTTTAGCGGCAACGGGGGTGATTAACGTTGGTGGGGAGCAGTTTTTTTTAACAATTGTTGGAGAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGttt                                                                              <  2:915495/142‑1 (MQ=255)
                                                                           cGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTTAGtttttt                                                                           >  1:6734/1‑148 (MQ=255)
                                                                              ggATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCt                                                                       >  2:795477/1‑149 (MQ=255)
                                                                              ggATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTg                                                                      >  1:284105/1‑150 (MQ=255)
                                                                               gATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTg                                                                      >  1:187286/1‑149 (MQ=255)
                                                                                aTGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGGTTTTTGGCTGaa                                                                    >  2:110549/1‑150 (MQ=255)
                                                                                 tGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGaa                                                                    >  1:8848/1‑149 (MQ=255)
                                                                                  ggtttggggCAACGCGGGTGATTTACGTTAGTGGTGTTCAGTCTTTTTTTACAATTGTTTGAGAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGa                                                                 <  2:963346/144‑1 (MQ=255)
                                                                                    ttttCGCCAACTCGGGTTATTTACGTTAGGGGGGATCAGGCTTTTTTTAAAACTTTTTGAGAAGAAAAAAAAAATCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGat                                                                <  1:537406/147‑1 (MQ=255)
                                                                                    ttACCGGCACCTCGGGTGATTTACGTTAGTGGTGATCAGTTTTTTTTTACATCTTTTTGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGat                                                                <  1:1145196/150‑1 (MQ=255)
                                                                                    ttACCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGTTTTTTTTACACCTGTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGaa                                                                    <  1:161673/146‑1 (MQ=255)
                                                                                         ggaaaCTGGGGTGATTTACTTTAGTGGTAATCAGGTTTTTTTTACAACTGTTGGAGAAAAAAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTca                                                          <  1:1063422/148‑1 (MQ=255)
                                                                                         ggCAACTGGGGTGTTTTAGGTTAGTGGTGTTCGGGTTTTTTTTACAATTGTTGGAGAAGAGAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTc                                                           <  2:1014339/150‑1 (MQ=255)
                                                                                          gCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGattatt                                                            >  2:490916/1‑148 (MQ=255)
                                                                                           cAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTcagc                                                        >  1:476589/1‑151 (MQ=255)
                                                                                           cAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTcag                                                         >  2:353179/1‑150 (MQ=255)
                                                                                            acctgggttaaTTTACGTTAGTGGTGATCGGGTTTTTTTTACAATTGTTTGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTcag                                                         <  2:567763/139‑1 (MQ=255)
                                                                                            aagtCGGGTGATTTACGTTAGGGGTGATCAGGTTTTTTTAACAGCTGTTAGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTca                                                          <  1:682792/145‑1 (MQ=255)
                                                                                            aaCTGGGGTAATTTACTTTAGTGGTGATCAGGTTTTTTTAACAACTTTTTGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTcagca                                                       <  1:31781/151‑1 (MQ=255)
                                                                                                 gggTGATTTCGGTGAGGGGGGACCGGTTTTTTTTTCCAAGTGTTGGAAAGAAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGt                                                     <  2:804921/148‑1 (MQ=255)
                                                                                                      ttttACGTTGGTGGTAATCAGGTTTTTTTTAAAATTTTTGGAAAAGAAAAAAAAAACCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAgg                                               <  2:121799/148‑1 (MQ=255)
                                                                                                      aTTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGc                                              >  1:495306/1‑150 (MQ=255)
                                                                                                       tttagtttAGTGGTGATCAGTTTTTTTTTACACCTGTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAg                                                <  1:110700/141‑1 (MQ=255)
                                                                                                       tttATGTTAGTGTTGATCGGGTTTTTTTTAAAATTTTTGGAAAAGAGAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGt                                                     <  2:979850/142‑1 (MQ=255)
                                                                                                        ttATGTTAGTGGTGATCTGTTTTTTTTTACAACTTTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCAt                                            <  2:935136/150‑1 (MQ=255)
                                                                                                           cGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGATGAAAGATTATTCAGCAGTTGCAGGCAtttt                                         >  2:542881/1‑150 (MQ=255)
                                                                                                             ttggtggtaATCAGTTTTTTTTTACAATTTTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTa                                        <  2:37879/146‑1 (MQ=255)
                                                                                                             ttAGTGTTGATCAGGTTTTTTTAACAACTGTTTGAGAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTa                                        <  1:852963/149‑1 (MQ=255)
                                                                                                             ttAGTGGTGATAAGGTTTTTTTTACAACTTTTTGAAAAGAGAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTa                                        <  1:203471/149‑1 (MQ=255)
                                                                                                                  ggtgATCAGGTTTTTTTTACAACTGTTTGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCtttt                                  <  1:662990/150‑1 (MQ=255)
                                                                                                                  ggtaatcggtttttttttACACCTGTTTGAGAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCt                                     <  2:422387/136‑1 (MQ=255)
                                                                                                                      aTCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGcc                               >  2:256632/1‑149 (MQ=255)
                                                                                                                                 tttACAATTGTTTAAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCt                   <  1:68930/150‑1 (MQ=255)
                                                                                                                                  ttACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGt                 >  2:446371/1‑151 (MQ=255)
                                                                                                                                  ttACAACTGTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTg                  <  2:338200/150‑1 (MQ=255)
                                                                                                                                   tACAACCGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGGTCTGt                 >  1:969700/1‑150 (MQ=255)
                                                                                                                                    aCAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTg                  >  1:731968/1‑148 (MQ=255)
                                                                                                                                     cAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCCGGACGTTCTGt                 >  2:791296/1‑148 (MQ=255)
                                                                                                                                     cAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGATTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTc                >  2:650464/1‑149 (MQ=255)
                                                                                                                                      agctttttGAGAAGAAAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAg              <  2:31662/145‑1 (MQ=255)
                                                                                                                                         tttttGAAAAGAAAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCa               <  1:258141/144‑1 (MQ=255)
                                                                                                                                         tGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGc          >  1:579556/1‑151 (MQ=255)
                                                                                                                                                agaagagaaAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTtctc     >  2:141838/1‑150 (MQ=255)
                                                                                                                                                agaagagaaAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGGTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTtc       >  1:886327/1‑148 (MQ=255)
                                                                                                                                                aaaagaaaaAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTtctc     <  2:256095/148‑1 (MQ=255)
                                                                                                                                                   agagaaAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGGTGCAGGCATTTTACCTTTTGCCGCTGGACGTTTTGTCAGAAGCTTCTCaaa  >  2:950273/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                  
CCGACATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAA  >  NC_000913/1824907‑1825202

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