Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,825,0610GA40.7% 77.3 / 31.7 60intergenic (‑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/3);  new base A (0/24);  total (32/28)
Fisher's exact test for biased strand distribution p-value = 1.35e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGC  >  NC_000913/1824911‑1825208
                                                                                                                                                        |                                                                                                                                                    
cATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGagaaaag                                                                                                                                                      <  1:518684/151‑1 (MQ=255)
 aTTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGagaaaaga                                                                                                                                                     <  1:34222/151‑1 (MQ=255)
     tatgcttatCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGAAACCCTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGCGGCGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAAc                                                                                                                                                >  2:749348/7‑151 (MQ=255)
     tACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAgaaaagaaaa                                                                                                                                                  <  2:216326/150‑1 (MQ=255)
     tACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAAc                                                                                                                                                >  1:220512/1‑151 (MQ=255)
        gCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAAccgc                                                                                                                                             >  2:670013/1‑151 (MQ=255)
          atatCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAAccgccg                                                                                                                                           >  1:160230/1‑151 (MQ=255)
            ttaaaaGTTTTCCATGCTCCCTCCCGGCCCGTTCTTTTCTTTTTTTTCTTTCAGGTACCCGTGGATGGTTATCGGCACCTCGGGTAATTAACGTAAGTGGTGATCAGGCTTTTTTTACACCTGTTTGAGAAGAAAAAAAAAAAccgccg                                                                                                                                           <  2:826223/146‑1 (MQ=255)
              cAAAGTATTCCATACTCGCTCCCCGCCCGCTCTTTTCATTGGCATCGTTCCGGTAACCGTGGATGCTTATCGGCAACTCGGGTGATTTACGTTTGTGGTGATTATGCTTTTTTTACAACTGTTTGAGAAGAGAAACGAAAACCGCCGATcc                                                                                                                                       >  2:642832/1‑151 (MQ=255)
                     ttCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCAc                                                                                                                                >  1:613564/1‑151 (MQ=255)
                            cTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCAtt                                                                                                                          >  2:337105/1‑150 (MQ=255)
                             tCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTAc                                                                                                                        >  1:773967/1‑151 (MQ=255)
                                ccccccccccTTTTTTTTTATTTTTATCTTGCAGTGACCCGTGGTTGTTTAGCGCAAACTGGGGTTATTTCTGTTAGTGGGGATCTGTTTTTTTTTACATGTGTTTGGGAGAAAAAAATAACACCGCCGATCCTGTCCACCGCATTACTGc                                                                                                                     <  2:417935/144‑1 (MQ=255)
                                  cccgcccgGTTTTTTTTATTTTTATCGTTCGGTTAGCGTTGGTTGGTTATGGGCACCTGGGGTGTTTTAGGTGAGTGTTGTTCGGGTTTTTTTAACAATTGTTTGAGAAGAAAAAAGAAAACCGCCGATCCTGTCCACCGCAtt                                                                                                                          <  1:126050/144‑1 (MQ=255)
                                    cGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAgg                                                                                                                 >  2:760642/1‑151 (MQ=255)
                                    cGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAgg                                                                                                                 >  1:429406/1‑151 (MQ=255)
                                      gCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGt                                                                                                                >  1:924684/1‑150 (MQ=255)
                                       cccTTTCTTTTCTTTTTTATGTTGCAGAAAACGGGGGTTGGTTGTCGGCACCGGGGTGGATTCACGTTAGTGGTGCTCAGTTTTTTTTAAAATTTTTTGAAAAAAAAAAAAAAAATCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAg                                                                                                              <  1:646474/151‑1 (MQ=255)
                                       cccGTTTTTTTATTTTTTATCTTCCAGGAACCCGGGGATGGTTACCGGCAACTGGGGTGATTTACGTTAGGGGTGATCGGTTTTTTTTAACAACTGTTGGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAgg                                                                                                                 <  1:1143335/148‑1 (MQ=255)
                                       cccGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGTTTACCGGCAACCGGGGTGATTTACGTTGGTGGTGATCGGTTTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAg                                                                                                              <  1:990968/151‑1 (MQ=255)
                                         cGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTg                                                                                                            >  1:1044983/1‑151 (MQ=255)
                                              ttttCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACa                                                                                                        >  1:386667/1‑150 (MQ=255)
                                               tttCATTGTTATCTTTCGGTAACCGGGGGTTGTTTATGGGCACCTGGGGTGTTTTATTTTGGTGGTGATCGGTTTTTTTTTACACTTGTTTGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACaa                                                                                                       <  2:868364/150‑1 (MQ=255)
                                                    tttttttCGTTCAGGTAACGGGGGTTGTTTATCGGCACCTGGGTTTTTTTACTTTAGGGGTGATCGGGTTTTTTTTACATCTGTTAGAGAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcg                                                                                                  <  2:737952/144‑1 (MQ=255)
                                                    tttttATCTTTAAGGTAGCCGTGGATGGTTACCGGCAATTGGGGTGATTTACGTTAGGGTTGATCAGTTTTTTTTTCCATCTGTTTGAGAAGAGAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGAcc                                                                                                   <  2:337831/146‑1 (MQ=255)
                                                    ttGTTTTCGTTGGGATAACCGTGGAGTTTTACCGGCAACTGGGGTGATTAGCTTTATGGGTGAGCTGCTTTTTTTTACACTTGTTTGAAAAGAGAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGAcc                                                                                                   <  1:1213046/149‑1 (MQ=255)
                                                     tGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcgg                                                                                                 >  2:814979/1‑150 (MQ=255)
                                                        tATCGTTCAGGTAACCGTGGATGTTTAGCGGAACCTGGGGTGATTTACGTTGGTGGTGATCAGGCTTTTTTTACAATTTTTTGAAAAGAAAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcggcgg                                                                                              <  2:160228/150‑1 (MQ=255)
                                                            gTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCtt                                                                                          >  1:60293/1‑150 (MQ=255)
                                                            gTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGGGATTTACGTTAGCGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGGCtt                                                                                          >  2:931044/1‑150 (MQ=255)
                                                              tCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGt                                                                                             >  2:994932/1‑145 (MQ=255)
                                                              tCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTa                                                                                         >  2:705753/1‑149 (MQ=255)
                                                                 ggTACCGGTGGATGGTTACCGGCACCTGGGGTGTTTTGCGTTGGTGGTGATCAGGTTTTTTTTACACCTGTTTGAGAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGtt                                                                                     <  1:582802/150‑1 (MQ=255)
                                                                   gtaccggTGGTTGTTTATCGGAACCTGGGTTGTTTTACTTTAGGGGTGATCAGTTTTTTTTTACAACTTTTGGAAAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGtttt                                                                                   <  1:350011/145‑1 (MQ=255)
                                                                         gTGGATGTTTATCGGCACCTGGGTTGTTTTACTTTGGTGGTGATCAGTTTTTTTTTACATTTTTTGGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTg                                                                                <  2:145480/148‑1 (MQ=255)
                                                                           ggATGTTTATCGGCAACTGGGGTGATTTACGTTAGGGGTGACCAGGTTTTTTTTACAATTGTTGGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGa                                                                           <  2:59760/151‑1 (MQ=255)
                                                                            gATGGTTATCGGCAACGGGGGTTTTTTACGTTGGTGGTGATCAGTTTTTTTTAACACTTTTTGGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTg                                                                            <  2:65067/149‑1 (MQ=255)
                                                                                   aTCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTCAGCTTTTTTGCTGAAAGattat                                                                   >  2:248030/1‑151 (MQ=255)
                                                                                     cGGCAACTCGGGCGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCAGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGattatt                                                                  >  2:507171/1‑150 (MQ=255)
                                                                                              ggttgttttAGGTTAGTGGTGACAAGTTTTTTTTTACAATTTTTAGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCGGCAGtt                                                          <  2:98537/143‑1 (MQ=255)
                                                                                              gggTGTTTTACTTTAGTGGTGACCGGTTTTTTTTTACACCGTTTTGAGAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGtt                                                          <  2:867966/149‑1 (MQ=255)
                                                                                              gggTGTTTTACGTTGGTGGTGATCAGGTTTTTTTTACAACTGTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGtt                                                          <  1:435959/149‑1 (MQ=255)
                                                                                                    tttcggTTGGTGGTGCACGGGTTTTTTTTACAATTTTTTGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGc                                                    <  2:82768/144‑1 (MQ=255)
                                                                                                    tttACGTGGGTGGTGATCAGTTTTTTTTAAAAATTGTTGGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGc                                                    <  2:741876/149‑1 (MQ=255)
                                                                                                      tACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCCGCCGTTGCAGGCAttt                                                >  1:1197691/1‑151 (MQ=255)
                                                                                                         gTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCAtttt                                               >  2:530223/1‑149 (MQ=255)
                                                                                                         gTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGGTTTTTGGCTGAAAGATTATTCCGCAGTTTCAGGCAtt                                                 >  1:1016515/1‑147 (MQ=255)
                                                                                                             gtggtgATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTTGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCttt                                         >  1:596420/1‑151 (MQ=255)
                                                                                                             gtggtgATCAGTTTTTTTTAAAACTGGTGTAAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCttt                                         <  2:1071737/151‑1 (MQ=255)
                                                                                                              tggtgctcggattttttttACAATTTTTGGAAAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCtt                                          <  2:1021275/137‑1 (MQ=255)
                                                                                                              tggggagcagtttttttttACACTTTTTTGAGAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCttt                                         <  2:679992/138‑1 (MQ=255)
                                                                                                               gttattcggtttttttttACACTGGTTAAAGAGAAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCAtt                                                 <  1:1050963/130‑1 (MQ=255)
                                                                                                               ggtgATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCtttt                                        >  1:445277/1‑150 (MQ=255)
                                                                                                                 tgATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGcc                                     >  1:105381/1‑151 (MQ=255)
                                                                                                                    tCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGc                                   >  1:326075/1‑150 (MQ=255)
                                                                                                                          tttttttACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTTCCGCTGGACg                             >  1:407821/1‑150 (MQ=255)
                                                                                                                            ttttaACACCTGTTGGAAAAAAGAAAAAAAAGCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGtt                           <  2:1005539/150‑1 (MQ=255)
                                                                                                                             ttttACAACTGTTAGAAAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTc                          <  1:505102/150‑1 (MQ=255)
                                                                                                                              tttACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGGTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCCGGACGTTCt                         >  1:855156/1‑150 (MQ=255)
                                                                                                                                  cAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTc                      >  1:916638/1‑149 (MQ=255)
                                                                                                                                    aactttataagaaaaaaaAAAAAAGCCGCCGATCCTGTCCCCCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGt                       <  2:970081/138‑1 (MQ=255)
                                                                                                                                      tGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGGTGAAAGATTATTCAGCAGTTGCCGGCCTTTTACCTTTTGCCGCTGGACGGTCTGTCAGACGc                >  1:695380/1‑151 (MQ=255)
                                                                                                                                        gTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGc                >  1:46339/1‑150 (MQ=255)
                                                                                                                                                 aagaaaagaaaACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCGTTTGCCGCTGGACTTTCTGTCAGAc                  <  2:382666/138‑1 (MQ=255)
                                                                                                                                                  agaaaaaaaaACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTTCCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCaaa        <  2:554641/148‑1 (MQ=255)
                                                                                                                                                  agaaaaaaaaACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATCATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTtc             <  2:373409/143‑1 (MQ=255)
                                                                                                                                                   aaaaagaaaaCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTa    <  1:670022/150‑1 (MQ=255)
                                                                                                                                                   aaaaaaaaaaCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAAtt     <  2:616180/149‑1 (MQ=255)
                                                                                                                                                     aaagaaaaCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGc  <  2:1044965/151‑1 (MQ=255)
                                                                                                                                                     aaaaaaaaCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGc  <  1:222949/151‑1 (MQ=255)
                                                                                                                                                        |                                                                                                                                                    
CATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGC  >  NC_000913/1824911‑1825208

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