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

TTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGC  >  NC_000913/1824915‑1825208
                                                                                                                                                    |                                                                                                                                                    
ttACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGagaaaagaaa                                                                                                                                                   <  2:757885/150‑1 (MQ=255)
ttACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGagaaaagaaa                                                                                                                                                   <  2:279572/150‑1 (MQ=255)
  aCGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAAcc                                                                                                                                               >  2:315776/1‑151 (MQ=255)
    gCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAAccgc                                                                                                                                             <  1:347518/151‑1 (MQ=255)
     cATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAAccgcc                                                                                                                                            >  1:212596/1‑151 (MQ=255)
        atCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGAt                                                                                                                                         >  1:562799/1‑151 (MQ=255)
          cAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATcc                                                                                                                                       >  1:797539/1‑151 (MQ=255)
           aaaGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCt                                                                                                                                      >  1:503083/1‑151 (MQ=255)
               tATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTcc                                                                                                                                  >  1:12508/1‑151 (MQ=255)
                aTTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTcc                                                                                                                                  >  2:493941/1‑150 (MQ=255)
                 ttCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCCACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAAACGCCGATCCTGTTCAc                                                                                                                                >  2:679303/1‑151 (MQ=255)
                  tCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCAcc                                                                                                                               >  1:16385/1‑151 (MQ=255)
                  tCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCAcc                                                                                                                               >  2:254037/1‑151 (MQ=255)
                  cccATCCTCCCTCCCGGCCTTTTTTTTTCATTGTTATCTTTCGGTGACGGGGGGATGGTTATCGGCAATGCGGGTTATTTCTGTTGGTGGGGTTCAGGTTTTTTTTACAACTTTTTAAGAAGAAAAAAAAAAACCGCCGATCCTGTCCAcc                                                                                                                               <  2:557961/150‑1 (MQ=255)
                    cATACTCGCCCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTACCCGGGGATGTTTATCGGCACCTCGGGTTTTTTACGTTAGTGGTGATCAGGTTTTTTTTACACCTGTTGGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGc                                                                                                                             <  1:25305/151‑1 (MQ=255)
                            cTCCCGGCCCGTTTTTTTCTTTGTTATCGTTCAGGTAACCGTGGATGGTTATCGCCAACTGGGGTGATTTACGTTAGTGGTGATCGGTTTTTTTTTACAACTTTTTGAAAAGAGAAAAAAAAACCGCCGATCCTGTCCACCGCATTACt                                                                                                                       <  1:130870/149‑1 (MQ=255)
                            cTCCCGCCCGTTTCTTTTCATTGTTCTCTTACAGGAAACCGGGGATTTTTCTCTGCAACCGGGGAAATTTACGGTAGGTGCAATCAGGTTTTTTTTCCACTTTTGGGAAAAGAAAAAAAAAAGCCGCCGATCCTGTCCACCGCATTACt                                                                                                                       <  2:907875/149‑1 (MQ=255)
                              cccGGCTCTTTTTTTTGTTTTTTATCTTTAAGGAAACCGTGGATTGTTTGCGGCAGCTGGGGTAATTTCCGTTGGTGGTCAGCGTGTTTTTTTAACAACTGTTGGAGAAAGAAAAAAAAAACCGCCGATCCTGTCCACCGCAt                                                                                                                           <  2:412126/143‑1 (MQ=255)
                               ccGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCa                                                                                                                    >  2:545075/1‑149 (MQ=255)
                                cGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAgg                                                                                                                 >  2:355071/1‑151 (MQ=255)
                                cGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAgg                                                                                                                 >  1:795819/1‑151 (MQ=255)
                                 gcccggTTTTTTTTTTTTTTATAGTTCAGGAAACGGGGGTTGTTTAGCGTAAACTGGGGTATTTTACGTTAGTGGTTATCAGGTTTTTTTTACAGCTTTTAGAGAAGAGAAAAAAAAGCCGCCGATCCTGTCCACCGCATTACTGCAAgg                                                                                                                 <  2:702771/145‑1 (MQ=255)
                                      gTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTgg                                                                                                           >  1:1004887/1‑151 (MQ=255)
                                          ttttCTTTTTTTTCGTTCAGTTACCCGTTGATTGTTACCGCCAACTGGGTTGTTTTACTTTAGTGGTACTCAGGTTTTTTTTAAATTTTTATAAAAAAAAAAAAAAAAGCCGCCGATCCTGTCCACCGCATTACTGCAATGTAGTGGAc                                                                                                         <  2:697794/149‑1 (MQ=255)
                                               aTTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcg                                                                                                  >  2:162334/1‑151 (MQ=255)
                                                tttttATCGTCCAGGGAACGGTGTATGTTTATCGGCAACTCGGGTGATTTACGTTGGTGGTGATCAGGTTTTTTTCAAAATTGTTTGAAAAAAGAAAAGAAACCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcg                                                                                                  <  2:1076058/147‑1 (MQ=255)
                                                  gTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcggcg                                                                                               >  2:1167952/1‑151 (MQ=255)
                                                   ttatcgtcggggcagccGTAGTTGGTTATCGTCAACTTGGTTTATTTACTTGAGTGGTGTTAATGTTTTTTTACCATCTATTTGAAAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcggcg                                                                                               <  1:736249/135‑1 (MQ=255)
                                                    tATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGt                                                                                             >  1:467748/1‑151 (MQ=255)
                                                     aTCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTc                                                                                            >  1:57405/1‑151 (MQ=255)
                                                          tCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCtt                                                                                          >  2:1116331/1‑148 (MQ=255)
                                                             gttacccGTGGATGTTTATCGGCACCTCGGTTGATTTATTTTAGTGGTATTCAGTTTTTTTTTACAATTGTTTGAGAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTaa                                                                                        <  2:412965/142‑1 (MQ=255)
                                                             ggTACCCGTGGATGTTTATCGGCAACGGGGGTATTTTACTTTAGGGGTGATCAGTTTTTTTTTACAACTTTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGttt                                                                                    <  1:944344/151‑1 (MQ=255)
                                                                  ccGTGGATGTTTATCGGCATCTCGGGTGTTTTACTTTAGTGGTGATCAGTTTTTTTTTACAATTGTTTGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTgg                                                                               <  1:76043/151‑1 (MQ=255)
                                                                   cggTGGATGTTTATCGGAAACGGGGGTGATTTACTTTAGTGGTGATCGGGTTTTTTTTACAATTTTTGGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTg                                                                                <  1:1271287/148‑1 (MQ=255)
                                                                       ggATGGTTTCCGGCATCTGGGGTTATTTTGTTTAGTGGTTATCAGTTTTTTTTTCCAACTTTTAGAGAAGAGAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTg                                                                            <  2:731369/150‑1 (MQ=255)
                                                                       ggATGGTTCCCGGAAGCTGGGGTATTTTGCTTGAGTGGTGATCGGCTTTTTTTTACATATGTTTGAGAAAAAAAAAAAAACCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGttttt                                                                                  <  1:517150/144‑1 (MQ=255)
                                                                             ttATCGGCAACTCGGGTGTTTTACGTTAGTGGTGATCAGTTTTTTTTTACATCTGTTGGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGa                                                                       <  1:185490/149‑1 (MQ=255)
                                                                               atcgccaccggggggGATTTAGTTAGTTGGTGACCAGTTTTTTTTTACAACTGTTGGGAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGatta                                                                    <  2:513914/139‑1 (MQ=255)
                                                                                  ggCACCTGGGGTGATTTACGTTAGTGGTGATCAGGTTTTTTTTACAACTGTTTAAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTc                                                                 <  1:466218/150‑1 (MQ=255)
                                                                                  gacaatttgggTGATTTACGTATGTGGTGATCGGGTTTTTTTTACATATTTTAGAGAAAAAAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGattat                                                                   <  2:706031/140‑1 (MQ=255)
                                                                                  gacAACTGGGGTGATGTAGGTTAGGGGTAGTCGTGTTTTTTTAACAGTTTTTTGGAAAGAAAAAGAAAAGCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGatta                                                                    <  2:845721/145‑1 (MQ=255)
                                                                                    cAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTcagc                                                              >  1:1092414/1‑151 (MQ=255)
                                                                                     acctggggTGTTTTACTTTAGTGGTGAGCAGGTTTTTTTTACAACTTTTGGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTcagca                                                             <  1:975006/146‑1 (MQ=255)
                                                                                            gTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGc                                                        >  1:824079/1‑149 (MQ=255)
                                                                                             tGATTTACGTTAGGGGTGATCAGGCTTTTTTTACAACTGGTTGAGAAGAGAAAAGAACACCGCCGATCCTGTCAACCGCATTACTGCAAGGTAGTGGACACGACCGGCGGTCTTTAGATTTTTGGCGGAAAGGTTATTCAGCCGTTGCAgg                                                     >  2:602144/1‑151 (MQ=255)
                                                                                                tttagtttAGGGGTGATCAGTTTTTTTTTACATTTGTTGAAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCa                                                       <  1:233587/140‑1 (MQ=255)
                                                                                                tttACTTTAGTGGTGATCAGGTTTTTTTTACAACTTTTGGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAg                                                      <  1:1203296/147‑1 (MQ=255)
                                                                                                tttACGTTAGTGGTGATCAGGTTTTTTTTACAACTGTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAgg                                                     <  1:636270/148‑1 (MQ=255)
                                                                                                     gTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTTTTCAGCAGTTGCAGGCAtttt                                               >  2:1031293/1‑149 (MQ=255)
                                                                                                      ttAGTGGTGATCGGGTTTTTTTAAAAACTGTTTGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCAtt                                                 <  1:723852/146‑1 (MQ=255)
                                                                                                       tAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCt                                           >  1:359856/1‑151 (MQ=255)
                                                                                                            ggtgATCGTGCTTTTTGTCCATTTGTTAGAAAAGAAAAAAAAAAGCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCtttt                                        <  2:7397/150‑1 (MQ=255)
                                                                                                                    ggttttttttACAACTTTTTGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTg                                 <  1:355228/146‑1 (MQ=255)
                                                                                                                          ttttACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGGTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCt                         >  1:690436/1‑151 (MQ=255)
                                                                                                                          ttttACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGAAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTTAGGAGTTGCAGGCATTTTACCTTTTGCCGCCGGACGt                            >  1:845409/1‑148 (MQ=255)
                                                                                                                              aCAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGATTTTTGGCTGAAAGATTATTTAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGt                       >  1:266314/1‑149 (MQ=255)
                                                                                                                                  cTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTTCCGCTGGACGTTCTGTCAGAc                  >  2:1166765/1‑150 (MQ=255)
                                                                                                                                  cTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGAc                  >  1:59909/1‑150 (MQ=255)
                                                                                                                                   tttttGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGAc                  <  1:954569/147‑1 (MQ=255)
                                                                                                                                     tttGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACg                 <  1:800006/148‑1 (MQ=255)
                                                                                                                                       tGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTTTTCAGCAGTTGCCGGACATTTACCTTTTGCCGCTGGAAGTTCTGTAAAACGCTtct            >  2:1168304/1‑151 (MQ=255)
                                                                                                                                        gagaagagaaAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTtct            >  1:155427/1‑150 (MQ=255)
                                                                                                                                        gaaaaaaaaaaaaaaaaCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCt               <  1:800104/139‑1 (MQ=255)
                                                                                                                                         aaaagaaaaAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTtc             <  1:1144118/146‑1 (MQ=255)
                                                                                                                                           aagagaaAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCaa         >  1:425879/1‑150 (MQ=255)
                                                                                                                                           aagaaaaaagaaaaCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTtctc           <  2:269262/148‑1 (MQ=255)
                                                                                                                                           aagaaaaAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCa          <  1:1204659/149‑1 (MQ=255)
                                                                                                                                            agaaaaAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCaaaa       <  2:170718/151‑1 (MQ=255)
                                                                                                                                            agaaaaAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCaaa        <  1:1194935/150‑1 (MQ=255)
                                                                                                                                             aaaaaaaaaaacCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCaaaa       <  1:59051/147‑1 (MQ=255)
                                                                                                                                              aaaaaaaaaaaCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAAtt     <  1:494361/149‑1 (MQ=255)
                                                                                                                                               gaaaagaaaaCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGACTTTAGGTTTTTGGCTGAAAGATTATTCAGGAGGTTCAGGCATTTTACCTTTTGGCGCTGGACGTTCTGTTAGACGCTTCTCAAAAtt     >  1:1284396/1‑150 (MQ=255)
                                                                                                                                               aaaaaaaaaaCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCaaaa       <  2:662388/147‑1 (MQ=255)
                                                                                                                                               aaaaaaaaaaCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTa    <  1:315926/150‑1 (MQ=255)
                                                                                                                                                aaaaaaaaaCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAg   <  1:254176/151‑1 (MQ=255)
                                                                                                                                                 aaagaaaaCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGc  >  1:510778/1‑151 (MQ=255)
                                                                                                                                                 aaaaaaaaCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGc  <  2:525742/151‑1 (MQ=255)
                                                                                                                                                  aagaaaaCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGc  >  1:14752/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGC  >  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: