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

ACATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATT  >  NC_000913/1824910‑1825205
                                                                                                                                                    |                                                                                                                                                      
aCATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGagaaa                                                                                                                                                     >  2:1143455/1‑150 (MQ=255)
     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                                                                                                              >  1:795819/1‑151 (MQ=255)
                                     cGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAgg                                                                                                              >  2:355071/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)
                                                                                                                                              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)
                                                                                                                                                    gaaaagaaaaCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGACTTTAGGTTTTTGGCTGAAAGATTATTCAGGAGGTTCAGGCATTTTACCTTTTGGCGCTGGACGTTCTGTTAGACGCTTCTCAAAAtt  >  1:1284396/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
ACATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATT  >  NC_000913/1824910‑1825205

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: