Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,825,0610GA46.2% 87.0 / 52.2 80intergenic (‑124/+79)ves/spycold‑ and stress‑inducible protein/periplasmic ATP‑independent protein refolding chaperone, stress‑induced
Reads supporting (aligned to +/- strand):  ref base G (34/9);  new base A (0/37);  total (34/46)
Fisher's exact test for biased strand distribution p-value = 2.18e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGC  >  NC_000913/1824917‑1825208
                                                                                                                                                 |                                                                                                                                                    
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)
                                                                                                                                           aaaaaaaaaaaCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGAGTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCGCAAAAt      <  2:840747/148‑1 (MQ=255)
                                                                                                                                             aaaagaaaaCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTTGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGAAGCTTCTCAAAAtt     >  2:607140/1‑149 (MQ=255)
                                                                                                                                             aaaagaaaaCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTTAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAAtt     >  1:646323/1‑149 (MQ=255)
                                                                                                                                             aaaagaaaaCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCACAAtt     >  2:163919/1‑149 (MQ=255)
                                                                                                                                             aaaaaaaaaCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTa    <  2:438492/150‑1 (MQ=255)
                                                                                                                                              aaagaaaacCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCAATTTACCTTTTGCCGCTGGGCGTTCTGTCAGACACTTCTCAAAATTa    <  2:610805/150‑1 (MQ=255)
                                                                                                                                              aaagaaaaCCGCCGATCCTGTCCACCGCATTGCTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGc  >  1:386245/1‑151 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
ACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGC  >  NC_000913/1824917‑1825208

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: