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

ATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGC  >  NC_000913/1824912‑1825208
                                                                                                                                                        |                                                                                                                                                   
aTTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGagaaaaga                                                                                                                                                    <  2:586859/151‑1 (MQ=255)
  tttACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGagaaaagaa                                                                                                                                                   >  2:172000/1‑150 (MQ=255)
  tttACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGagaaaaga                                                                                                                                                    >  2:677240/1‑149 (MQ=255)
   ttACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGagaaaagaaa                                                                                                                                                  <  1:913542/150‑1 (MQ=255)
      cGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAAccg                                                                                                                                              >  1:1163194/1‑151 (MQ=255)
        cATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAAccgcc                                                                                                                                            <  1:67813/151‑1 (MQ=255)
        cATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAAccgcc                                                                                                                                            >  2:1047459/1‑151 (MQ=255)
         atatCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAAccgccg                                                                                                                                           >  1:1183615/1‑151 (MQ=255)
           atCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGAt                                                                                                                                         >  2:1177137/1‑151 (MQ=255)
              aaaGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCt                                                                                                                                      >  1:681943/1‑151 (MQ=255)
              aaaGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCt                                                                                                                                      >  1:1248222/1‑151 (MQ=255)
              aaaGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCt                                                                                                                                      >  2:115793/1‑151 (MQ=255)
               aaGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTg                                                                                                                                     >  2:909401/1‑151 (MQ=255)
               aaGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTg                                                                                                                                     >  1:1453271/1‑151 (MQ=255)
                aGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGt                                                                                                                                    >  2:270191/1‑151 (MQ=255)
                   aTTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGCTTGAGAAGAGAAAAGAAGACCGCCGATCCTGTcc                                                                                                                                  >  1:1136781/1‑150 (MQ=255)
                    ttCCATACCCCCCCCCGGCCCTTTTTTTTCTTTGTTATCGTTCAGGTACCCGTGGATGGTTATCGGCAACTGGGGTGTTTTACTTTAGTGGTGATCAGTTTTTTTTTACAACTTTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTcc                                                                                                                                  <  1:677396/149‑1 (MQ=255)
                          aCTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCAtt                                                                                                                          >  2:817850/1‑151 (MQ=255)
                           cccGCTCCCGGCCTGTTCTTTTCATTGTTATCGTTCAGGTACCGGTGGTTGTTTACGGGAAACTGGGGTGTTTTGCTTTGGTGGTGACGGGGTTTTTTTTACACTTTTTGGAGAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCAt                                                                                                                           <  1:720588/147‑1 (MQ=255)
                           cccGCCCCCGCCCTTTTTTTTTTATTTTTATCGTTCAGGTAACGGTGGATGTTTATCGGCAACTGGGGTTATTTACTTTGGTGGTGACCGGTTTTTTTTTACAACTGTTGGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGc                                                                                                                             <  2:924603/145‑1 (MQ=255)
                            tCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTAc                                                                                                                        >  1:597862/1‑151 (MQ=255)
                            tCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTAc                                                                                                                        >  1:264857/1‑151 (MQ=255)
                             cccccccGGCCCTTTTTTTTTTTTTTTATCGTTCAGGAAACGGTGGATGTTTATCGGCACCTGGGGTGTTTTACTTGAGTGGTGCTCAGTTTTTTTTTCCACCTGTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGc                                                                                                                             <  2:245611/141‑1 (MQ=255)
                                tCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGc                                                                                                                     >  1:578335/1‑150 (MQ=255)
                                 cccGGCCTTTTTTTTTAATTTTTATCGTTCAGGAACCGGTGGATGGTTATCGGAACCTGGGGTGATTTAGTTTGGTGGTGATCAGGTTTTTTTTACAACTTTTGGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACt                                                                                                                       <  1:118532/147‑1 (MQ=255)
                                         gTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTgg                                                                                                           >  1:128794/1‑151 (MQ=255)
                                          tttttttCTTTGTTATCGTTCAGGTACCGGTGGTGGTTTATCGGCAACTGGGTTTATTTACTTTGGTGGGATCCAGGTTTTTTTAACATTTTTTTGAAAAAAAAAAAAAAACCCGCCGATCCTGTCCACCGCATTACt                                                                                                                       <  2:1067134/135‑1 (MQ=255)
                                           tCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGAc                                                                                                         >  1:708804/1‑151 (MQ=255)
                                            cTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGAc                                                                                                         >  1:556108/1‑150 (MQ=255)
                                             ttttCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACa                                                                                                        >  1:123326/1‑150 (MQ=255)
                                                      ttATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcggcgg                                                                                              >  2:1052840/1‑151 (MQ=255)
                                                        aTCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTc                                                                                            >  1:1018405/1‑151 (MQ=255)
                                                           gTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCtt                                                                                          >  2:811832/1‑150 (MQ=255)
                                                            ttCAGGTAACGGGGGTTGGTTATCGCCAACTCGGGTGATTTACTTTGGTGGTGATCAGTTTTTTTTTACACTTTTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCtt                                                                                          <  1:1355159/149‑1 (MQ=255)
                                                             tCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGGCtt                                                                                          >  1:889460/1‑148 (MQ=255)
                                                               aGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGtt                                                                                     >  2:977355/1‑151 (MQ=255)
                                                                gttAACGGTGGATGGTTATCGGCAACTGGGGTGATTTACTTTAGTGGTGATCAGGTTTTTTTTACAACTGTTTGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGtt                                                                                     <  1:188094/148‑1 (MQ=255)
                                                                 gaaccgggggATGTTTATCGGCACCTGGGGTTATTTACGTTAGTGGGGACCAGTTTTTTTTTAAAAGTTTTTGAAAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGt                                                                                      <  2:1392133/140‑1 (MQ=255)
                                                                 gTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGtttt                                                                                   >  1:384993/1‑151 (MQ=255)
                                                                   aaCGGTGGAGGTTTACCGGCATCTGGGGTGTTTTACGTTAGTGGAGACCAGGTTTTTTTTACATCTTTTTAAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGtttt                                                                                   <  1:1306215/149‑1 (MQ=255)
                                                                   aaCCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGtttttt                                                                                 >  1:107311/1‑151 (MQ=255)
                                                                    aCCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGtttttt                                                                                 >  2:607019/1‑150 (MQ=255)
                                                                     cggTGGATGTTTTTCGCCACCGGGGGTTATTTACTTGTGTGGGGATCGGTTTTTTTTAACACTTGTTGAAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCtt                                                                                          <  2:1456555/138‑1 (MQ=255)
                                                                      cGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGGTTTTTgg                                                                               >  1:783349/1‑150 (MQ=255)
                                                                       gTGGATGTTTATCGGCACCTCGGGTGATTTACGTTAGTGGTGATCAGGTTTTTTTTACACCTGTTTGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGtt                                                                                     <  1:1052888/143‑1 (MQ=255)
                                                                       gTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTTAGGTTTTTGGCt                                                                             >  2:1309841/1‑151 (MQ=255)
                                                                       gTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTgg                                                                               >  1:962655/1‑149 (MQ=255)
                                                                       gTGGATGGTTATCGGCAACTCGGGCGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGCAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGACTTGAGgttttt                                                                                 >  1:1351676/1‑147 (MQ=255)
                                                                        tGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTTAGtttttttggt                                                                             >  2:937813/1‑146 (MQ=255)
                                                                             ggTTATCGGCAACGGGGGTGTTTTACTTTGGTGGTGATCGGTTTTTTTTTACAATTTTTAAAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGaa                                                                          <  2:1298683/148‑1 (MQ=255)
                                                                                           ggggTGTTTTATGTTGGGGGTAAGCAGTTTTTTTTAACATCTTTTGGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTcagc                                                              <  2:662028/145‑1 (MQ=255)
                                                                                           ggggTGATTTACTTGAGGGGTGACGGGTTTTTTTTAAAAACTGTTGGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTcagcag                                                            <  2:1086462/147‑1 (MQ=255)
                                                                                            gggTGTTTTACGTTAGTGGTGATCAGTTTTTTTTTACACCTTTTGGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGtt                                                          <  1:312423/149‑1 (MQ=255)
                                                                                            gggTAATTTACTTTGGTGTTAATCAGTTTTTTTTTAAAATTTTTAAAAAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGtt                                                          <  2:779292/149‑1 (MQ=255)
                                                                                             gttgttttattttAGTGGTAATCAGTTTTTTTTTCAATTTTTTGGAAAAAAAAAAAAAAACCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCa                                                       <  2:728981/140‑1 (MQ=255)
                                                                                              gTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTTGCTGAAAGATTATTCAGAAGtt                                                          >  1:172094/1‑147 (MQ=255)
                                                                                                  tttcttttAGTGGTAATAAGTTTTTTTTTAAACTTGTTGAAAAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCa                                                   <  1:681870/143‑1 (MQ=255)
                                                                                                  tttATGTTGGTGGTGACCAGGTTTTTTTTACAACTGTTGGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCa                                                   <  1:1162479/150‑1 (MQ=255)
                                                                                                  tttACTTTGGTTGTGATCGGCTTTTTTTACCACTTTTTAAAAAAGAAAAAAAAAATCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCa                                                   <  1:727231/150‑1 (MQ=255)
                                                                                                  tttACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGCCTTACGATTTTTTGTTCGAAGAATTTTCCGCAGGTGccg                                                      >  2:637566/1‑145 (MQ=255)
                                                                                                     aCGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGCTTTAAGTTTTTTGGCTGAAAGATTATTCCGCAGGTGGAGGGAtttt                                               >  2:602888/1‑151 (MQ=255)
                                                                                                      cGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTTAGCAGTTGCAGGCAtttt                                               >  2:1091271/1‑150 (MQ=255)
                                                                                                       actttggtggatatcggCTTTTTTTCATAATTTTTTGAAAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCAttt                                                <  1:458827/135‑1 (MQ=255)
                                                                                                            ggggataccctgttttttttACAACTTTTGGAGAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGTACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCtt                                          <  2:1375743/137‑1 (MQ=255)
                                                                                                             ggtgATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGCCTTTATTTTTTTTGTTGAAAGGTTATTAACACGTTGCCGGCCTTTTTCCtttt                                        >  2:972927/1‑150 (MQ=255)
                                                                                                              ggtgtgcggtttttttttACATTTGTTGGAAAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCtttt                                        <  1:719386/139‑1 (MQ=255)
                                                                                                                 aTCAGGTTTTTTTTACATTTTTTGGAAAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTg                                       <  2:1183557/147‑1 (MQ=255)
                                                                                                                    cAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGGTGCAGGCATTTTACCTTTTGCCGc                                   >  1:33981/1‑149 (MQ=255)
                                                                                                                    cAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAAATTATTCAGCAGTTGCAGGCATTTTTCCTTTTGCCGCTg                                 >  1:974956/1‑151 (MQ=255)
                                                                                                                     aGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTgg                                >  1:1141634/1‑151 (MQ=255)
                                                                                                                      gtttttttttACAACTTTTGGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCt                                  <  1:1346133/145‑1 (MQ=255)
                                                                                                                      gttttttttaACACCTTTTGGAAAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGa                               <  1:766440/148‑1 (MQ=255)
                                                                                                                      ggttttttttACATTTGTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCtt                                          <  2:488096/137‑1 (MQ=255)
                                                                                                                          ttttttACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGtt                           >  1:572059/1‑151 (MQ=255)
                                                                                                                           tttttACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGGCGGCATTTTACCTTTTTCCGCTGGACg                             >  2:1001663/1‑148 (MQ=255)
                                                                                                                            ttttACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGAAGTGGACAAGACCGGCGGTCTTAAGCTTTTTCGCTGAAGCGTTATTTAACAATTGCAGGCATTTTTACTTTTGCCGCTGGGCggt                           >  2:1450343/1‑147 (MQ=255)
                                                                                                                                 aCAGCTTTTGGAGAAAAGAAAGGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCa                     <  1:1446964/151‑1 (MQ=255)
                                                                                                                                 aCAAGTGTTAGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTc                      <  1:33319/150‑1 (MQ=255)
                                                                                                                                  cAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGACTTTAGGTTTTTGGGTGAAAGAATTTTTAGCAGTTGCAGGCATTTTACCTTTTGGCCCCGGACGTTccttc                      >  2:306894/1‑145 (MQ=255)
                                                                                                                                   aaCTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTTGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTc                      >  2:197532/1‑148 (MQ=255)
                                                                                                                                      cttttagaaaagaaaaAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTc                      <  2:1398159/137‑1 (MQ=255)
                                                                                                                                       gTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGc                >  2:139635/1‑150 (MQ=255)
                                                                                                                                         ttagaaaagaaaaAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGAc                  <  1:811853/141‑1 (MQ=255)
                                                                                                                                         ttGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTTAGGTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCtt              >  1:820531/1‑150 (MQ=255)
                                                                                                                                             agaagagaaAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTTGCTGAAAGATTATTCAGCAGTTGCAGGGCTTTTAACTTTTGTCGCTGGACGTTCTGTCAGACGCTtc             >  2:1119862/1‑148 (MQ=255)
                                                                                                                                             agaagagaaAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTTCCGCTGGACGTTCTGTCAGACGCTtctc           >  1:793067/1‑150 (MQ=255)
                                                                                                                                             agaagaaaaAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTtctc           <  1:172002/150‑1 (MQ=255)
                                                                                                                                             aaaaaaaaaaaaaaaaCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTtct            <  2:112441/142‑1 (MQ=255)
                                                                                                                                               aagagaaAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCCTTTTACCTTTTTCCGCTGGACGTTCTGTCAGACGCTTCTCaa         >  2:447291/1‑150 (MQ=255)
                                                                                                                                               aagaaaaAAAAAACCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCaaa        <  1:1253039/151‑1 (MQ=255)
                                                                                                                                               aagaaaaAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTGCCTTTTGCCGCTGGACGTTCTGTCAGACGCTtc             <  2:407691/146‑1 (MQ=255)
                                                                                                                                                agagaaAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCaaa        >  1:1138340/1‑150 (MQ=255)
                                                                                                                                                  agaaaagaaaACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAAtt     >  1:327451/1‑151 (MQ=255)
                                                                                                                                                    aaaagaaaaCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTa    <  2:1468035/150‑1 (MQ=255)
                                                                                                                                                    aaaaaaaaaCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTa    <  2:1224919/150‑1 (MQ=255)
                                                                                                                                                    aaaaaaaaaCCGCCGATCCTGTCCACCGCAATACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAg   <  1:819122/151‑1 (MQ=255)
                                                                                                                                                     aaagaaaaCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGc  <  1:630610/151‑1 (MQ=255)
                                                                                                                                                        |                                                                                                                                                   
ATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGC  >  NC_000913/1824912‑1825208

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