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

CATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAG  >  NC_000913/1824911‑1825207
                                                                                                                                                       |                                                                                                                                                  
cATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGagaaaag                                                                                                                                                    <  1:1123707/151‑1 (MQ=255)
 aTTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGagaaaaga                                                                                                                                                   <  1:87950/151‑1 (MQ=255)
 aTTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGagaaaaga                                                                                                                                                   <  1:352129/151‑1 (MQ=255)
 aTTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGagaaaag                                                                                                                                                    <  2:975443/150‑1 (MQ=255)
   tttACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGagaaaagaa                                                                                                                                                  <  2:478667/150‑1 (MQ=255)
      aCGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAAc                                                                                                                                               <  2:49103/150‑1 (MQ=255)
      aCGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAAc                                                                                                                                               <  1:470115/150‑1 (MQ=255)
      aCGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAAc                                                                                                                                               <  1:1296549/150‑1 (MQ=255)
       cGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAAccg                                                                                                                                             >  1:435015/1‑151 (MQ=255)
        gCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAAc                                                                                                                                               >  2:1221916/1‑148 (MQ=255)
         cATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAAccgcc                                                                                                                                           >  2:831213/1‑151 (MQ=255)
          atatCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAAccgccg                                                                                                                                          >  2:968504/1‑151 (MQ=255)
            ataaaaTTATTCCTTACTCGCTCCCGGCCCTTTCTTTTCATTGTTATCTTTCAGGTAACCGTGGATGGTTATCGGAAACTCGGGTGATTTACGTTAGTGGTGATCAGGTTTTTTTTACAACTGTTGGAGAAGAAAAAAAAAAAccgccg                                                                                                                                          <  2:738278/146‑1 (MQ=255)
            atCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGAt                                                                                                                                        >  1:133695/1‑151 (MQ=255)
              cAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATcc                                                                                                                                      >  2:1389991/1‑151 (MQ=255)
               aaaGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCt                                                                                                                                     >  1:841059/1‑151 (MQ=255)
               aaaGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCt                                                                                                                                     >  2:439548/1‑151 (MQ=255)
                     ttCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCAc                                                                                                                               >  2:871342/1‑151 (MQ=255)
                      tCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCAcc                                                                                                                              >  1:1275770/1‑151 (MQ=255)
                        cATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGc                                                                                                                            >  2:89337/1‑151 (MQ=255)
                           aCTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCAtt                                                                                                                         >  1:1471773/1‑151 (MQ=255)
                           aCTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCAtt                                                                                                                         >  2:1522054/1‑151 (MQ=255)
                             tCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTAc                                                                                                                       >  1:11260/1‑151 (MQ=255)
                                cccccGGCCCGTTTTTTTCTTTGTTATGGTTCAGGTACCTGTGGATGGTTATCGACAACTGGGGTGTTTTACGTTAGTGGTGATCAGTTTTTTTTTACAACTGTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTa                                                                                                                        <  2:291714/145‑1 (MQ=255)
                                cTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCTGCAACTCGGGGGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTAc                                                                                                                       >  1:1088246/1‑148 (MQ=255)
                                  cccGGCCCGTTTTTTTCATTTTTATCGTTCAGGTAACGGTGGATGTTTATCGGCAACTGGGGTGATTTACGTTAGTGGTGATCAGGTTTTTTTTACAATTGTTGGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCaa                                                                                                                  <  1:1252375/151‑1 (MQ=255)
                                   ccGGCCCGTTTTTTTCTTTGTTATCGTTCAGGTACCGGTGGATGTTTATCGGCACCTGGGGTGATTTACTTTAGTGGTGATCGGGTTTTTTTTACAACTGTTTAAGAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAg                                                                                                                 <  1:41287/151‑1 (MQ=255)
                                     ggCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGt                                                                                                               >  1:873631/1‑151 (MQ=255)
                                       cctgTTTTTTCCATTTTTATCGTTGAGGTACCCGGGATGGTTTTTCACCATCTCGGTAAATTTACGTTAGGGGTGATCAGGCTTTTTTCACATCTGTTGAAAAAGAAAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAg                                                                                                                 <  2:247677/144‑1 (MQ=255)
                                       cccGTTTTTTTCTTTTTTTTGTTTCAGGTACCCGTGGATGTTTACCGGAAACTCGGGTAATTTCGGTTAGTGTTGATCGGTTTTTTTTAACAACTTTTGAAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTa                                                                                                              <  2:430515/150‑1 (MQ=255)
                                       cccGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGt                                                                                                               >  2:415968/1‑149 (MQ=255)
                                        ccGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGt                                                                                                            >  2:643209/1‑151 (MQ=255)
                                           tttttttCTTTTTTATCTTTCGGGTACCGGGGGATGTTTATCGGCACCTGGGGTGATTTCCGTTGGTGGTGATCGGGTTTTTTTAACAACTTTTTGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTg                                                                                                           <  1:190072/146‑1 (MQ=255)
                                               ttttATTTTTTTCTTACAGGAACCGGTGGAGGTTTGCCGGCAACGGGGGTATTTTACTTTGGGGGTGATCGGTTTTTTTTTACAACTGTTGGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACaa                                                                                                      <  1:752455/150‑1 (MQ=255)
                                               tttCTTTTTTTTCGTGCAGGTAACCGGGGATGTTTATCGCAAACTGGGGTGATTTACGTTTGTGGTGACCAGGTTTTTTTTACAACTTTTTGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACaa                                                                                                      <  1:98985/150‑1 (MQ=255)
                                                ttCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACa                                                                                                       >  1:1001614/1‑148 (MQ=255)
                                                 tCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGAc                                                                                                   >  2:1517377/1‑151 (MQ=255)
                                                    tttttATCTTTGGGGTACCGGGGGAGGGTTGTCGGAAACTGGGGTATTTTACGTTGGGGGTGATCGGTTTTTTTTTAAATTTTTTGGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcg                                                                                                 <  1:868559/147‑1 (MQ=255)
                                                    ttGTTATCTTTAAGGTAACGGGGGAGGTTTTTCGGAAACGGGGTTGTTTTCTTTAAGTGGTGCTCGTGTTTTTTTTACAGCTGTTAGAGAAGAAAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcg                                                                                                 <  2:798967/150‑1 (MQ=255)
                                                             ttCAGTTACCGGGGGATGGTTAGCGGCAACTGGGGTGATTTACTTTAGTGGTGATAAGTTTTTTTTTACAACTTTTGGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGt                                                                                            <  2:513030/146‑1 (MQ=255)
                                                             ttCAGTTACCCGGGGATGGTTACCGGCAACTGGGGTGATTTACGTTAGTGGTGATCGGGTTTTTTTTACACTTTTTTGAGAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTa                                                                                        <  2:1365877/150‑1 (MQ=255)
                                                                 ggTAACCGTGGATGGTTTTCGGCAACTGGGGTGTTTTACGTTAGGGGTGATCAGGTTTTTTTTACAACTTTTGGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGtt                                                                                    <  1:517969/150‑1 (MQ=255)
                                                                    aacctgggATGGTTACCGGCAACTGGGTTATTTTACTTTAGTGGGAATCGGGTTTTTTTTCAAACTGTTTAAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGtt                                                                                    <  1:92906/141‑1 (MQ=255)
                                                                     aCCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGtttttt                                                                                >  2:155958/1‑150 (MQ=255)
                                                                        gTGGATGTTTATCGGCAACTGGGGTGTTTTACTTTAGTGGTGATCGGGTTTTTTTTACAATTTTTGGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCt                                                                            <  1:1138997/151‑1 (MQ=255)
                                                                         tGGATGTTTCTCGGAATTTCGGATTTTTTACGTTAGTGGTGATTGGTTTTTTTTTACAGCTGTTTGAAAAGAGAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGttttt                                                                                 <  2:1297534/145‑1 (MQ=255)
                                                                           gATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGa                                                                          >  1:954397/1‑150 (MQ=255)
                                                                               gtttATCGAAACCTGGGGTGATTTCCGTTGGTGGTGATCGGGTTTTTTTTACACTTGTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTg                                                                           <  1:47261/144‑1 (MQ=255)
                                                                                tttATCGGAAACGGGGTTAATTTCCGTTGGTGGTGACCAGGTTTTTTTAACAATTGTTGGAAAGGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGaa                                                                         <  2:1084783/146‑1 (MQ=255)
                                                                                 ttATGGGCAACTGGGGTATTTTACGTGGGTGGTGATCGGTTTTTTTTTACATGTGTTGGAGAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGat                                                                     <  2:1565542/150‑1 (MQ=255)
                                                                                 ttATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGatt                                                                    >  1:708694/1‑151 (MQ=255)
                                                                                  tATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTTGCTGAAAGatt                                                                    >  2:1373716/1‑150 (MQ=255)
                                                                                      ggCAACTGGGGTTTTTTCAGTTAGGGGTATTCGGGTTTTTTTAACAATTGTTAAAAAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAg                                                                       <  1:645752/143‑1 (MQ=255)
                                                                                      ggCAACTCGGGTGATTTACGTTAGTGGCGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAAGTTTTTGGCTGAAAGatt                                                                    >  1:1383204/1‑146 (MQ=255)
                                                                                       gatatgtgggTGATTTCTGTTGGGGGTCATCGGATTTTTTTTACATTTTTTGTAAAAAAGAAAAAAAGACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGatta                                                                   <  2:621699/139‑1 (MQ=255)
                                                                                          aCTGGGGTATTTTACGTTGGGGGTGATCAGGTTTTTTTTACACTTGTTTGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTcagcag                                                           <  2:870104/151‑1 (MQ=255)
                                                                                              gggTGTTTTACGTTAGTGGTGATCAGGTTTTTTTTACAACTGTTTGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTg                                                        <  1:1021006/150‑1 (MQ=255)
                                                                                               ggTGATTGCCGTTAGTGTTAATCAGGTTTTTTTTACCATTGTTGAAAAAAAAAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTg                                                        <  2:931822/149‑1 (MQ=255)
                                                                                                    tttACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGGCTTAAGTTTTTTTGCTGAAAGATTATTCAGCCGTTGCAGGCAt                                                 >  2:1080036/1‑151 (MQ=255)
                                                                                                     ttccGTTAGTGGTGCTCGGGTTTTTTTTACATTTTTTGGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCAt                                                 <  2:715793/147‑1 (MQ=255)
                                                                                                     ttACGTTAGTGTTAATCAGGTTTTTTTAACAACTGTTAGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCAt                                                 <  1:147067/150‑1 (MQ=255)
                                                                                                         gTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTAc                                            >  1:1103957/1‑151 (MQ=255)
                                                                                                             gtggtgATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTTAGTTTTTTGGCTGAAAGATTATTCAGCAGGTGCAGGCATTTTACCttt                                        >  2:1217684/1‑151 (MQ=255)
                                                                                                              tggtgATCAGGTTTTTTTTACAATTGTTGGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCt                                          <  2:537471/148‑1 (MQ=255)
                                                                                                              tggtgATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGCCTTTAGTTTTTTGGCTGGAAGATTATTTAGGAGTTGCAGGGATTTTTCCtttt                                       >  1:1170291/1‑151 (MQ=255)
                                                                                                                     cAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGGTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTg                                >  2:190900/1‑151 (MQ=255)
                                                                                                                          tttttttACATCTGTTTGAGAAGAGAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACg                            <  1:1286157/150‑1 (MQ=255)
                                                                                                                          tttttttACACCTTTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGAc                             <  1:467303/149‑1 (MQ=255)
                                                                                                                               ttACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGGTTTTTGGCTGAAAGATTATTCAGCAGTTTCAGGCATTTTACCTTTTTCCGCTGGACGGTCTGt                      >  1:547214/1‑151 (MQ=255)
                                                                                                                                     ttttttAAGAAGAGAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGGCAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAg                   <  2:913424/145‑1 (MQ=255)
                                                                                                                                     cTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACg                >  1:1242390/1‑151 (MQ=255)
                                                                                                                                      ttttggagaagaaaaAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGa                  <  2:1275752/143‑1 (MQ=255)
                                                                                                                                      tGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGCAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGGCGTTCTGTCAGAc                 >  2:1104486/1‑149 (MQ=255)
                                                                                                                                            agaagagaaAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGcgga            <  1:1171454/148‑4 (MQ=255)
                                                                                                                                               agagaaAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTTGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTTCCGCTGGACGTTCTGTCAGACGCTTCTCaaa       >  2:245206/1‑150 (MQ=255)
                                                                                                                                                  aaaaaaaaaaCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTa   <  1:1380233/150‑1 (MQ=255)
                                                                                                                                                   aaaaaaaaaCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAg  <  1:160612/151‑1 (MQ=255)
                                                                                                                                                   aaaaaaaaaCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAg  <  2:874937/151‑1 (MQ=255)
                                                                                                                                                    aaagaaaaCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGAATATTCAGCAGTTGGAGGCATTTTACCTTTTGCCGCTGGACGATCTGTCAGACGCTTCTCacca      >  2:1338584/1‑143 (MQ=255)
                                                                                                                                                       |                                                                                                                                                  
CATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAG  >  NC_000913/1824911‑1825207

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: