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

CGACATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGC‑AACTCGGGTGATTTAC‑GTTAGTGGTGATCAGGCTTTTTTTA‑CAACTGTTT‑GAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAA  >  NC_000913/1824908‑1825203
                                                                                                                                                        |                                                                                                                                                   
cGACATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGC‑AACTCGGGTGATTTAC‑GTTAGTGGTGATCAGGCTTTTTTTA‑CAACTGTTT‑gagaagagaa                                                                                                                                                   <  2:232214/151‑1 (MQ=255)
  aCATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGC‑AACTCGGGTGATTTAC‑GTTAGTGGTGATCAGGCTTTTTTTA‑CAACTGTTT‑GAGAAGagaaa                                                                                                                                                  >  2:682498/1‑150 (MQ=255)
    aTTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGC‑AACTCGGGTGATTTAC‑GTTAGTGGTGATCAGGCTTTTTTTA‑CAACTGTTT‑GAGAAGagaaaaga                                                                                                                                               <  1:400190/151‑1 (MQ=255)
     ttttACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGC‑AACTCGGGTGATTTAC‑GTTAGTGGTGATCAGGC‑TTTTTTA‑CAACTGTTT‑GAGAAGagaaaagaa                                                                                                                                              >  1:1140408/1‑150 (MQ=255)
       ttACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGC‑AACTCGGGTGATTTAC‑GTTAGTGGTGATCAGGCTTTTTTTA‑CAACTGTTT‑GAGAAGagaaaagaaa                                                                                                                                             >  1:959987/1‑150 (MQ=255)
        taagCATATCAAAGTATTCCATACCCGCTCCCCGACCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGC‑AACTCGGGTGATTTAC‑GTTAGTGGTTGTCAGGCTTTTTTTA‑CAACTGTTT‑GAGAAGagaaaag                                                                                                                                                <  1:787642/143‑1 (MQ=255)
        tACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGC‑AACTCGGGTGATTTAC‑GTTAGTGGTGATCAGGCTTTTTTTA‑CAACTGTTT‑GAGAAGAGAAAAGAAAAc                                                                                                                                           >  2:35995/1‑151 (MQ=255)
        tACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGC‑AACTCGGGTGATTTAC‑GTTAGTGGTGATCAGGCTTTTTTTA‑CAACTGTTT‑GAGAAGAGAAAAGAAAAc                                                                                                                                           >  2:1128030/1‑151 (MQ=255)
           gCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGC‑AACTCGGGTGATTTAC‑GTTAGTGGTGATCAGGCTTTTTTTA‑CAACTGTTT‑GAGAAGAGAAAAGAAAAccgc                                                                                                                                        >  1:314993/1‑151 (MQ=255)
             atatCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGC‑AACTCGGGTGATTTAC‑GTTAGTGGTGATCAGGCTTTTTTTA‑CAACTGTTT‑GAGAAGAGAAAAGAAAAccgccg                                                                                                                                      >  2:547642/1‑151 (MQ=255)
              tatCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGC‑AACTCGGGTGATTTAC‑GTTAGTGGTGATCAGGCTTTTTTTA‑CAACTGTTT‑GAGAAGAGAAAAGAAAAccgccg                                                                                                                                      >  1:861499/1‑150 (MQ=255)
                 cAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGC‑AACTCGGGTGATTTAC‑GTTAGTGGTGATCAGGCTTTTTTTA‑CAACTGTTT‑gagaagagaa                                                                                                                                                   <  1:887908/134‑1 (MQ=255)
                 cAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGC‑AACTCGGGTGATTTAC‑GTTAGTGGTGATCAGGCTTTTTTTA‑CAACTGTTT‑gagaagagaa                                                                                                                                                   >  2:887912/1‑134 (MQ=255)
                       aTTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGC‑AACTCGGGTGATTTAC‑GTTAGTGGTGATCAGGCTTTTTTTA‑CAACTGTTT‑GAGAAGAGAAAAGAAAACCGCCGATCCTGTcc                                                                                                                             >  1:1262884/1‑150 (MQ=255)
                        ttCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGC‑AACTCGGGTGATTTAC‑GTTAGTGGTGATCAGGCTTTTTTTA‑CAACTGTTT‑GAGAAGAGAAAAGAAAACCGCCGATCCTGTCCAc                                                                                                                           >  1:1041277/1‑151 (MQ=255)
                        ttCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGC‑AACTCGGGTGATTTAC‑GTTAGTGGTGATCAGGCTTTTTTTA‑CAACTGTTT‑GAGAAGAGAAAAGAAAACCGCCGATCCTGTCCAc                                                                                                                           >  1:416158/1‑151 (MQ=255)
                          ccatccccccccccGCCCCTTTTTTTTCTTTGTTATCTTTCAGGAACCCGTGGATGGTTATCGGC‑AACTGGGGTGTTTTAC‑TTTGGTGGTGATCAGGTTTTTTTTA‑CAATTTTTT‑GAAAAAAAAAAAAAAAACCGCCGATCCTGTCCAcc                                                                                                                          <  2:664738/139‑1 (MQ=255)
                            aTACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGC‑AACTCGGGTGATTTAC‑GTTAGTGGTGATCAGGCTTTTTTTA‑CAACTGTTT‑GAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGc                                                                                                                        >  2:52065/1‑150 (MQ=255)
                             tACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGC‑AACTCGGGTGATTTAC‑GTTAGTGGTGATCAGGCTTTTTTTA‑CAACTGTTT‑GAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCAt                                                                                                                      >  2:151416/1‑151 (MQ=255)
                               cTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGC‑AACTCGGGTGATTTAC‑GTTAGTGGTGATCAGGCTTTTTTTA‑CAACTGTTT‑GAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCAtt                                                                                                                     >  1:804352/1‑150 (MQ=255)
                                      ccGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGC‑AACTCGGGTGATTTAC‑GTTAGTGGTGATCAGGCTTTTTTTA‑CAACTGTTT‑GAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCa                                                                                                               >  2:384248/1‑149 (MQ=255)
                                       cGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGC‑AACTCGGGTGATTTAC‑GTTAGTGGTGATCAGGCTTTTTTTA‑CAACTGTTT‑GAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAgg                                                                                                            >  1:421003/1‑151 (MQ=255)
                                        ggccttttttttttatttttaTGTTTGAGGTACCCGGGGAGGTTTAGCGGC‑AACTGGGGTGATTTGC‑GTTAGTGGTGATCAGGTTTTTTTTA‑CAGCTTTTT‑GAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCAt                                                                                                                      <  1:1231595/131‑1 (MQ=255)
                                          cccGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGC‑AACTCGGGTGATTTAC‑GTTAGTGGTGATCAGGCTTTTTTTA‑CAACTGTTT‑GAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGt                                                                                                           >  1:1225720/1‑149 (MQ=255)
                                              ttCTTTTTATTGTTTTCGTTCAGGTACCGGTGGATGTTTATCGGC‑AACTGGGGTATTTTAC‑GTTAGTGGTGATCGGGTTTTTTTTA‑CATTTGTTT‑GAAAAAAAAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGt                                                                                                        <  2:590432/148‑1 (MQ=255)
                                              ttCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGC‑AACTCGGGTGATTTAC‑GTTAGTGGTGATCAGGCTTTTTTTA‑CAACTGTTT‑GAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTgg                                                                                                      >  1:1197916/1‑150 (MQ=255)
                                                cTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGC‑AACTCGGGTGATTTAC‑GTTAGTGGTGATCAGGCTTTTTTTA‑CAACTGTTT‑GAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGAc                                                                                                    >  2:414519/1‑150 (MQ=255)
                                                       tttttttCGTTAAGGAAACGTTGGTTGGTTAGCGGC‑‑ACTGGGGGTGTTTACTGTTGGGGGGGAGC‑GGTTTTTTTTACCAGCTGTTT‑GAAAAGAAAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGAcc                                                                                              <  1:995155/143‑1 (MQ=255)
                                                          ttttCGTTCGGGTACCCGGGGTTGGTTATCGCC‑AACTGGGGTGATTTAC‑GTTAGTGGTGATCAGTTTTTTTTAA‑CAATTTTTG‑GAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcggc                                                                                           <  1:252559/146‑1 (MQ=255)
                                                             tCGTTCAGGTAACCGTGGATGGTTATCGGC‑AACTCGGGTGATTTAC‑GTTAGTGGTGATCAGGCTTTTTTTA‑CAACTGGTT‑GAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcggcgggc                                                                                       >  1:904132/1‑148 (MQ=255)
                                                                ttCAGGTAACCGTGGATGGTTATCGGC‑AACTCGGGTGATTTAC‑GTTAGTGGTGATCAGGCTTTTTTTA‑CAACTGTTT‑GAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCtt                                                                                     >  2:863353/1‑149 (MQ=255)
                                                                ttCAGGTAACCGGGGTTGTTTACCGGC‑AACTGGGGTGATTTAC‑TTTAGGGGTAATCGGGTTTTTTTTA‑CACCTTTTG‑GAAAAGAAAAAAAAAACCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGt                                                                                        <  1:667371/146‑1 (MQ=255)
                                                                tgcgggTAACGGGGGTTGTTTATCGGC‑ACCTGGGGTTTTTTAC‑GTTAGGGGTGATCAGTTTTTTTTTA‑CATTTTTTG‑GAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGt                                                                                        <  2:1055768/142‑1 (MQ=255)
                                                                 tCAGGTAACCGTGGATGGTTATCGGC‑AACTCGGGTGATTTAC‑GTTAGTGGTGATCAGGCTTTTTTTA‑CAACTGTTT‑GAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCtt                                                                                     >  1:140278/1‑148 (MQ=255)
                                                                    ggTACCGGTGGATGTTTATCGGC‑ACCTGGGGTATTTTAC‑TTAGGTGGTGATC‑GGTTTTTTTTA‑AAAACGTTTAAAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTaa                                                                                   <  1:810763/147‑1 (MQ=255)
                                                                       aaCCGGGGATTGTTTGCGGC‑AACTGGGGTGATTTAC‑GTTAGGGGTGGTGGGTTTTTTTTTA‑CAACTTTTT‑AAAAAGAGAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGtttttt                                                                            <  1:429589/151‑1 (MQ=255)
                                                                        cgcGTGGATGTTTACCGGA‑AACTGGGATATTTTAC‑GTTGGCGGTAATCGATTTTTTTTTA‑CCCTGGTTG‑GAGAAAAAAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAg                                                                                  <  1:228055/142‑1 (MQ=255)
                                                                         cggTGGTTGGTTACCGCC‑ACCTGGGGTGTTTTAC‑TTTGGTGGTGATCGGTTTTTTTTTA‑AAACTGTTG‑GAAAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTgg                                                                          <  1:1106667/149‑1 (MQ=255)
                                                                                tgtttATCGGA‑AACTGGGGTGTTTTAC‑TTTAGTGGTGATCAGTTTTTTTTAA‑CAATTGTTT‑GAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGaaa                                                                    <  1:960742/147‑1 (MQ=255)
                                                                                tGGTTATCGGC‑AACTCGGGTGATTTAC‑GTTAGTGGTGATCAGGCTTTTTTTA‑CAACTGTTT‑GAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTTAGTTTTTTGGCTGaa                                                                     >  1:165923/1‑149 (MQ=255)
                                                                                 gtttATGGGC‑ATCTGGGGTGTTTTAT‑TTTGGTGGTGTTCCGTTTTTTTTTA‑AATGTGTTT‑GAAAAAAAAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTg                                                                           <  2:857788/140‑1 (MQ=255)
                                                                                  gTTATCGGC‑AACTCGGGTGATTTAC‑GTTAGTGGTGATCAGGCTTTTTTTA‑CAACTGTTT‑GAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTTGCTGaa                                                                     >  1:311810/1‑147 (MQ=255)
                                                                                  gTTATCGGC‑AACTCGGGTGATTTAC‑GTTAGTGGTGATCAGGCTTTTTTTA‑CAACTGTTT‑GAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTTGCTGaa                                                                     >  2:644662/1‑147 (MQ=255)
                                                                                     ttCGGCAAACTGGGGTGTTTTAC‑TTTGGTGGTGATCAGGTTTTTTTTA‑CAATTTTTA‑GAAAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGatt                                                                <  1:966334/149‑1 (MQ=255)
                                                                                      tCGGC‑AACTCGGGTGATTTAC‑GTTAGTGGTGATCAGGCTTTTTTTA‑CAACTGTTT‑GAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGGAAGattatt                                                             >  1:111799/1‑151 (MQ=255)
                                                                                       cGGC‑AACTCGGGTGATTTAC‑GTTAGTGGTGATCAGGCTTTTTTTA‑CAACTGTTT‑GAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGattatt                                                             >  1:555459/1‑150 (MQ=255)
                                                                                        ggC‑AACTGGGGTGATTTAC‑GTTAGTGGTGATCAGGTTTTTTTTA‑CAACTGTTA‑GAGAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTca                                                           <  1:890833/151‑1 (MQ=255)
                                                                                         ggaaaCTGGGGTGATTTAC‑GTTAGTGGTGATCAGTTTTTTTTAA‑AAACTGTTG‑GAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGattatt                                                             <  1:719901/146‑1 (MQ=255)
                                                                                         gC‑AACTCGGGTGATTTAC‑GTTAGTGGTGATCAGGCTTTTTTTA‑CAACTGTTT‑GAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGattatt                                                             >  1:394213/1‑148 (MQ=255)
                                                                                            aaCTGGGGTGATTTAC‑GTTGGTGGTTATCAGTTTTTTTTTA‑AAACTTTTG‑AAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTcag                                                          <  1:404629/149‑1 (MQ=255)
                                                                                             cctggggTGTTTTAC‑TTTGGTGGTGATCAGGTTTTTTTTA‑AAATTGTTG‑GAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTcagca                                                        <  2:1075832/146‑1 (MQ=255)
                                                                                             aCTCGGGTGATTTAC‑GTTAGTGGTGATCAGGCTTTTTTTA‑CAACTGTTT‑GAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTcagc                                                         >  2:6691/1‑149 (MQ=255)
                                                                                               tCGGGTGTTTTAC‑TTTGGTGGTAATCAGTTTTTTTTTA‑CAATTGTTG‑GAGAAGAAAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTcagcag                                                       <  1:1231592/149‑1 (MQ=255)
                                                                                                    tGATTTAC‑GTTAGTGGTGATCAGGCTTTTTTCA‑CAACTGTTT‑GAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAgg                                                >  2:625006/1‑151 (MQ=255)
                                                                                                       tttAC‑TTTAGTGGTAATCGGGTTTTTTTTA‑CAATTTTTT‑GAAAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAg                                                 <  1:112176/147‑1 (MQ=255)
                                                                                                          aC‑GTTAGTGGTGATCGGGTTTTTTTTA‑CAACTGTTT‑GAAAAAAAAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCAttt                                           <  2:654444/150‑1 (MQ=255)
                                                                                                          aC‑GTTAGTGGTGATCAGGCTTTTTTTA‑CAACTGTTT‑GAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCAtttt                                          >  2:252420/1‑151 (MQ=255)
                                                                                                                    gtgATCTGTTTTTTTTTA‑CACTTTTTG‑GAGAAGAAAAAAAAAACCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCttt                                    <  1:1221498/148‑1 (MQ=255)
                                                                                                                     tgATCAGGCTTTTTTTA‑CAACTGTTT‑GAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGGCCTAAGATTTTTTGCTTAAAGAATTTTCAGGAGTTGTAGGCATTTTaatttgt                                   >  2:586412/1‑142 (MQ=255)
                                                                                                                        tCAGGTTTTTTTTA‑CACCTTTTG‑AAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGc                              <  2:839729/150‑1 (MQ=255)
                                                                                                                        tCAGGCTTTTTTTA‑CAACTGTTT‑GAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCCGGCATTTTACCTTTTGCCGc                              >  2:878486/1‑150 (MQ=255)
                                                                                                                         cAGGCTTTTTTTA‑CAACTGTTT‑GAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTTGCTTAAAAATTATTCAGCAGTTGCAGGGATTTTACCTTTTGCCGCTg                            >  1:1170973/1‑151 (MQ=255)
                                                                                                                          aGGCTTTTTTTA‑CAACTGTTT‑GAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTgg                           >  2:89759/1‑151 (MQ=255)
                                                                                                                            gCTTTTTTTA‑CAACTGTTT‑GAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTTCCGCTgg                           >  1:644022/1‑149 (MQ=255)
                                                                                                                            gCTTTTTTTA‑CAACTGTTT‑GAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGCTTTAAGGTTTTTTGCTTCAAGCTTATTCAGCAGTTGAAGGCATTTTACCATTTGcc                                >  1:94287/1‑144 (MQ=255)
                                                                                                                              tttttttA‑CAACGTTGT‑GAAAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTGCCTTTTGCCGCTGGACg                        <  1:785852/150‑1 (MQ=255)
                                                                                                                              ttttttaA‑CAACTTTTG‑AAAAAGAGAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGc                              <  1:161745/144‑1 (MQ=255)
                                                                                                                               ttttttA‑CAACTGTTT‑GAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGtt                      >  2:363741/1‑151 (MQ=255)
                                                                                                                               ttttttA‑CAACTGTTT‑GAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGt                       <  1:599092/150‑1 (MQ=255)
                                                                                                                                tttttC‑CAACTTTTT‑GAAAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGtt                      <  2:1178703/150‑1 (MQ=255)
                                                                                                                                tttttA‑CAACTGTTT‑GAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGtt                      >  2:942611/1‑150 (MQ=255)
                                                                                                                                   ttA‑CAACTGTTT‑GAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGt                  >  2:1235310/1‑151 (MQ=255)
                                                                                                                                    aa‑CAATTGTTG‑GAAAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGt                  <  1:46253/149‑1 (MQ=255)
                                                                                                                                     a‑CAACTGTTT‑GAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTc                 <  2:602020/150‑1 (MQ=255)
                                                                                                                                      aaaaCTGTTA‑GAGAAGAAAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGt                  <  2:39647/147‑1 (MQ=255)
                                                                                                                                        aattGTTA‑GAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTc                 <  1:557695/145‑1 (MQ=255)
                                                                                                                                           tGTTT‑GAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTTCCGCTGGACGTTCTGTCAGACg            >  1:7764/1‑150 (MQ=255)
                                                                                                                                            gTTT‑GAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGc           >  1:1172729/1‑150 (MQ=255)
                                                                                                                                              ttggagaagaaaaAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACg            <  1:1052036/145‑1 (MQ=255)
                                                                                                                                               t‑GAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTtct       >  2:212679/1‑151 (MQ=255)
                                                                                                                                                  agaagagaaAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTTGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTTCCGCTGGACGTTCTGTCAGACGCTtctc      >  1:574704/1‑150 (MQ=255)
                                                                                                                                                  aaaagaaaaAAAAAACCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGAc             <  1:160811/141‑1 (MQ=255)
                                                                                                                                                    aagaaaaaagaaaaCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCaaa   <  1:169117/151‑1 (MQ=255)
                                                                                                                                                      gagaaAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTTGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCaaaa  >  2:341059/1‑150 (MQ=255)
                                                                                                                                                        |                                                                                                                                                   
CGACATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGC‑AACTCGGGTGATTTAC‑GTTAGTGGTGATCAGGCTTTTTTTA‑CAACTGTTT‑GAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAA  >  NC_000913/1824908‑1825203

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