Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,825,0610GA45.5% 76.9 / 41.9 66intergenic (‑124/+79)ves/spycold‑ and stress‑inducible protein/periplasmic ATP‑independent protein refolding chaperone, stress‑induced
Reads supporting (aligned to +/- strand):  ref base G (32/4);  new base A (0/30);  total (32/34)
Fisher's exact test for biased strand distribution p-value = 8.50e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAAC‑CGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAG‑TGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGC  >  NC_000913/1824915‑1825208
                                                                                                                                                    |                                                                                                                                                    
ttACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAAC‑CGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAG‑TGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGagaaaagaaa                                                                                                                                                   <  2:757885/150‑1 (MQ=255)
ttACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAAC‑CGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAG‑TGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGagaaaagaaa                                                                                                                                                   <  2:279572/150‑1 (MQ=255)
  aCGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAAC‑CGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAG‑TGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAA‑cc                                                                                                                                               >  2:315776/1‑151 (MQ=255)
    gCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAAC‑CGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAG‑TGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAA‑ccgc                                                                                                                                             <  1:347518/151‑1 (MQ=255)
     cATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAAC‑CGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAG‑TGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAA‑ccgcc                                                                                                                                            >  1:212596/1‑151 (MQ=255)
        atCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAAC‑CGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAG‑TGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAA‑CCGCCGAt                                                                                                                                         >  1:562799/1‑151 (MQ=255)
          cAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAAC‑CGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAG‑TGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAA‑CCGCCGATcc                                                                                                                                       >  1:797539/1‑151 (MQ=255)
           aaaGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAAC‑CGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAG‑TGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAA‑CCGCCGATCCt                                                                                                                                      >  1:503083/1‑151 (MQ=255)
               tATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAAC‑CGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAG‑TGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAA‑CCGCCGATCCTGTcc                                                                                                                                  >  1:12508/1‑151 (MQ=255)
                aTTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAAC‑CGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAG‑TGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAA‑CCGCCGATCCTGTcc                                                                                                                                  >  2:493941/1‑150 (MQ=255)
                 ttCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAAC‑CGTGGATGGTTATCGGCCACTCGGGTGATTTACGTTAG‑TGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAA‑ACGCCGATCCTGTTCAc                                                                                                                                >  2:679303/1‑151 (MQ=255)
                  tCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAAC‑CGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAG‑TGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAA‑CCGCCGATCCTGTCCAcc                                                                                                                               >  1:16385/1‑151 (MQ=255)
                  tCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAAC‑CGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAG‑TGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAA‑CCGCCGATCCTGTCCAcc                                                                                                                               >  2:254037/1‑151 (MQ=255)
                  cccATCCTCCCTCCCGGCCTTTTTTTTTCATTGTTATCTTTC‑GGTGACGGGGGGATGGTTATCGGCAATGCGGGTTATTTCTGTTGG‑TGGGGTTCAGGTTTTTTTTACAACTTTTTAAGAAGAAAAAAAAAAA‑CCGCCGATCCTGTCCAcc                                                                                                                               <  2:557961/150‑1 (MQ=255)
                    cATACTCGCCCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTACC‑CGGGGATGTTTATCGGCACCTCGGGTTTTTTACGTTAG‑TGGTGATCAGGTTTTTTTTACACCTGTTGGAAAAGAAAAAAAAAAA‑CCGCCGATCCTGTCCACCGc                                                                                                                             <  1:25305/151‑1 (MQ=255)
                            cTCCCGGCCCGTTTTTTTCTTTGTTATCGTTCAGGTAAC‑CGTGGATGGTTATCGCCAACTGGGGTGATTTACGTTAG‑TGGTGATCGGTTTTTTTTTACAACTTTTTGAAAAGAGAAAAAAAAA‑CCGCCGATCCTGTCCACCGCATTACt                                                                                                                       <  1:130870/149‑1 (MQ=255)
                            cTCCCGCCCGTTTCTTTTCATTGTTCTCTTACAGGAAAC‑CGGGGATTTTTCTCTGCAACCGGGGAAATTTACGGTAG‑GTGCAATCAGGTTTTTTTTCCACTTTTGGGAAAAGAAAAAAAAAAG‑CCGCCGATCCTGTCCACCGCATTACt                                                                                                                       <  2:907875/149‑1 (MQ=255)
                              cccGGCTCTTTTTTTTGTTTTTTATCTTTAAGGAAAC‑CGTGGATTGTTTGCGGCAGCTGGGGTAATTTCCGTTGG‑TGGTCAGCGTGTTTTTTTAACAACTGTTGGAGAAAGAAAAAAAAAA‑CCGCCGATCCTGTCCACCGCAt                                                                                                                           <  2:412126/143‑1 (MQ=255)
                               ccGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAAC‑CGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAG‑TGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCa                                                                                                                    >  2:545075/1‑149 (MQ=255)
                                cGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAAC‑CGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAG‑TGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAgg                                                                                                                 >  2:355071/1‑151 (MQ=255)
                                cGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAAC‑CGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAG‑TGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAgg                                                                                                                 >  1:795819/1‑151 (MQ=255)
                                 gcccggTTTTTTTTTTTTTTATAGTTCAGGAAAC‑GGGGGTTGTTTAGCGTAAACTGGGGTATTTTACGTTAG‑TGGTTATCAGGTTTTTTTTACAGCTTTTAGAGAAGAGAAAAAAAAG‑CCGCCGATCCTGTCCACCGCATTACTGCAAgg                                                                                                                 <  2:702771/145‑1 (MQ=255)
                                      gTTCTTTTCATTGTTATCGTTCAGGTAAC‑CGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAG‑TGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTgg                                                                                                           >  1:1004887/1‑151 (MQ=255)
                                          ttttCTTTTTTTTCGTTCAGTTACC‑CGTTGATTGTTACCGCCAACTGGGTTGTTTTACTTTAG‑TGGTACTCAGGTTTTTTTTAAATTTTTATAAAAAAAAAAAAAAAAG‑CCGCCGATCCTGTCCACCGCATTACTGCAATGTAGTGGAc                                                                                                         <  2:697794/149‑1 (MQ=255)
                                               aTTGTTATCGTTCAGGTAAC‑CGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAG‑TGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcg                                                                                                  >  2:162334/1‑151 (MQ=255)
                                                tttttATCGTCCAGGGAAC‑GGTGTATGTTTATCGGCAACTCGGGTGATTTACGTTGG‑TGGTGATCAGGTTTTTTTCAAAATTGTTTGAAAAAAGAAAAGAAAC‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcg                                                                                                  <  2:1076058/147‑1 (MQ=255)
                                                  gTTATCGTTCAGGTAAC‑CGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAG‑TGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcggcg                                                                                               >  2:1167952/1‑151 (MQ=255)
                                                   ttatcgtcggggcagc‑cGTAGTTGGTTATCGTCAACTTGGTTTATTTACTTGAG‑TGGTGTTAATGTTTTTTTACCATCTATTTGAAAAAAAAAAAAAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcggcg                                                                                               <  1:736249/135‑1 (MQ=255)
                                                    tATCGTTCAGGTAAC‑CGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAG‑TGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGt                                                                                             >  1:467748/1‑151 (MQ=255)
                                                     aTCGTTCAGGTAAC‑CGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAG‑TGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTc                                                                                            >  1:57405/1‑151 (MQ=255)
                                                          tCAGGTAAC‑CGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAG‑TGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCtt                                                                                          >  2:1116331/1‑148 (MQ=255)
                                                             gttacc‑cGTGGATGTTTATCGGCACCTCGGTTGATTTATTTTAG‑TGGTATTCAGTTTTTTTTTACAATTGTTTGAGAAAAAAAAAAAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTaa                                                                                        <  2:412965/142‑1 (MQ=255)
                                                             ggTACC‑CGTGGATGTTTATCGGCAACGGGGGTATTTTACTTTAG‑GGGTGATCAGTTTTTTTTTACAACTTTTTGAGAAGAAAAAAAAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGttt                                                                                    <  1:944344/151‑1 (MQ=255)
                                                                  c‑cGTGGATGTTTATCGGCATCTCGGGTGTTTTACTTTAG‑TGGTGATCAGTTTTTTTTTACAATTGTTTGAAAAGAAAAAAAAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTgg                                                                               <  1:76043/151‑1 (MQ=255)
                                                                   cggTGGATGTTTATCGGAAACGGGGGTGATTTACTTTAG‑TGGTGATCGGGTTTTTTTTACAATTTTTGGAAAAGAAAAAAAAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTg                                                                                <  1:1271287/148‑1 (MQ=255)
                                                                       ggATGGTTTCCGGCATCTGGGGTTATTTTGTTTAG‑TGGTTATCAGTTTTTTTTTCCAACTTTTAGAGAAGAGAAAAAAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTg                                                                            <  2:731369/150‑1 (MQ=255)
                                                                       ggATGGTTCCCGGAAGCTGGGGTATTTTGCTTGAG‑TGGTGATCGGCTTTTTTTTACATATGTTTGAGAAAAAAAAAAAAAC‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGttttt                                                                                  <  1:517150/144‑1 (MQ=255)
                                                                             ttATCGGCAACTCGGGTGTTTTACGTTAG‑TGGTGATCAGTTTTTTTTTACATCTGTTGGAGAAGAAAAAAAAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGa                                                                       <  1:185490/149‑1 (MQ=255)
                                                                               atcgccaccggggggGATTTA‑GTTAGTTGGTGACCAGTTTTTTTTTACAACTGTTGGGAAAGAAAAAAAAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGatta                                                                    <  2:513914/139‑1 (MQ=255)
                                                                                  ggCACCTGGGGTGATTTACGTTAG‑TGGTGATCAGGTTTTTTTTACAACTGTTTAAGAAGAAAAAAAAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTc                                                                 <  1:466218/150‑1 (MQ=255)
                                                                                  gacaatttgggTGATTTACGTATG‑TGGTGATCGGGTTTTTTTTACATATTTTAGAGAAAAAAAAAGAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGattat                                                                   <  2:706031/140‑1 (MQ=255)
                                                                                  gacAACTGGGGTGATGTAGGTTAG‑GGGTAGTCGTGTTTTTTTAACAGTTTTTTGGAAAGA‑AAAAGAAAAGCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGatta                                                                    <  2:845721/145‑1 (MQ=255)
                                                                                    cAACTCGGGTGATTTACGTTAG‑TGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTcagc                                                              >  1:1092414/1‑151 (MQ=255)
                                                                                     acctggggTGTTTTACTTTAG‑TGGTGAGCAGGTTTTTTTTACAACTTTTGGAAAAGAAAAAAAAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTcagca                                                             <  1:975006/146‑1 (MQ=255)
                                                                                            gTGATTTACGTTAG‑TGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGc                                                        >  1:824079/1‑149 (MQ=255)
                                                                                             tGATTTACGTTAG‑GGGTGATCAGGCTTTTTTTACAACTGGTTGAGAAGAGAAAAGAACA‑CCGCCGATCCTGTCAACCGCATTACTGCAAGGTAGTGGACACGACCGGCGGTCTTTAGATTTTTGGCGGAAAGGTTATTCAGCCGTTGCAgg                                                     >  2:602144/1‑151 (MQ=255)
                                                                                                tttagtttAG‑GGGTGATCAGTTTTTTTTTACATTTGTTGAAAAAGAAAAAAAAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCa                                                       <  1:233587/140‑1 (MQ=255)
                                                                                                tttACTTTAG‑TGGTGATCAGGTTTTTTTTACAACTTTTGGAGAAGAAAAAAAAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAg                                                      <  1:1203296/147‑1 (MQ=255)
                                                                                                tttACGTTAG‑TGGTGATCAGGTTTTTTTTACAACTGTTTGAGAAGAAAAAAAAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAgg                                                     <  1:636270/148‑1 (MQ=255)
                                                                                                     gTTAG‑TGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTTTTCAGCAGTTGCAGGCAtttt                                               >  2:1031293/1‑149 (MQ=255)
                                                                                                      ttAG‑TGGTGATCGGGTTTTTTTAAAAACTGTTTGAAAAGAAAAAAAAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCAtt                                                 <  1:723852/146‑1 (MQ=255)
                                                                                                       tAG‑TGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCt                                           >  1:359856/1‑151 (MQ=255)
                                                                                                            ggtgATCGTGCTTTTTGTCCATTTGTTAGAAAAGAAAAAAAAAAG‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCtttt                                        <  2:7397/150‑1 (MQ=255)
                                                                                                                    ggttttttttACAACTTTTTGAAAAGAAAAAAAAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTg                                 <  1:355228/146‑1 (MQ=255)
                                                                                                                          ttttACAACTGTTTGAGAAGAGAAAAGAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGGTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCt                         >  1:690436/1‑151 (MQ=255)
                                                                                                                          ttttACAACTGTTTGAGAAGAGAAAAGAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGAAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTTAGGAGTTGCAGGCATTTTACCTTTTGCCGCCGGACGt                            >  1:845409/1‑148 (MQ=255)
                                                                                                                              aCAACTGTTTGAGAAGAGAAAAGAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGATTTTTGGCTGAAAGATTATTTAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGt                       >  1:266314/1‑149 (MQ=255)
                                                                                                                                  cTGTTTGAGAAGAGAAAAGAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTTCCGCTGGACGTTCTGTCAGAc                  >  2:1166765/1‑150 (MQ=255)
                                                                                                                                  cTGTTTGAGAAGAGAAAAGAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGAc                  >  1:59909/1‑150 (MQ=255)
                                                                                                                                   tttttGAAAAGAAAAAAAAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGAc                  <  1:954569/147‑1 (MQ=255)
                                                                                                                                     tttGAAAAGAAAAAAAAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACg                 <  1:800006/148‑1 (MQ=255)
                                                                                                                                       tGAGAAGAGAAAAGAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTTTTCAGCAGTTGCCGGACATTTACCTTTTGCCGCTGGAAGTTCTGTAAAACGCTtct            >  2:1168304/1‑151 (MQ=255)
                                                                                                                                        gagaagagaaAAGAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTtct            >  1:155427/1‑150 (MQ=255)
                                                                                                                                        gaaaaaaaaaaaaaaaa‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCt               <  1:800104/139‑1 (MQ=255)
                                                                                                                                         aaaagaaaaAAAAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTtc             <  1:1144118/146‑1 (MQ=255)
                                                                                                                                           aagagaaAAGAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCaa         >  1:425879/1‑150 (MQ=255)
                                                                                                                                           aagaaaaaagaaaa‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTtctc           <  2:269262/148‑1 (MQ=255)
                                                                                                                                           aagaaaaAAAAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCa          <  1:1204659/149‑1 (MQ=255)
                                                                                                                                            agaaaaAAAAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCaaaa       <  2:170718/151‑1 (MQ=255)
                                                                                                                                            agaaaaAAAAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCaaa        <  1:1194935/150‑1 (MQ=255)
                                                                                                                                             aaaaaaaaaaac‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCaaaa       <  1:59051/147‑1 (MQ=255)
                                                                                                                                              aaaaaaaaaaa‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAAtt     <  1:494361/149‑1 (MQ=255)
                                                                                                                                               gaaaagaaaa‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGACTTTAGGTTTTTGGCTGAAAGATTATTCAGGAGGTTCAGGCATTTTACCTTTTGGCGCTGGACGTTCTGTTAGACGCTTCTCAAAAtt     >  1:1284396/1‑150 (MQ=255)
                                                                                                                                               aaaaaaaaaa‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCaaaa       <  2:662388/147‑1 (MQ=255)
                                                                                                                                               aaaaaaaaaa‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTa    <  1:315926/150‑1 (MQ=255)
                                                                                                                                                aaaaaaaaa‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAg   <  1:254176/151‑1 (MQ=255)
                                                                                                                                                 aaagaaaa‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGc  >  1:510778/1‑151 (MQ=255)
                                                                                                                                                 aaaaaaaa‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGc  <  2:525742/151‑1 (MQ=255)
                                                                                                                                                  aagaaaa‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGc  >  1:14752/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAAC‑CGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAG‑TGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGC  >  NC_000913/1824915‑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: ‑