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

CGACATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAA  >  NC_000913/1824908‑1825202
                                                                                                                                                      |                                                                                                                                                  
ggACATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGATCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTgagaagaga                                                                                                                                                   >  2:534320/2‑150 (MQ=255)
cGACATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTgagaagagaa                                                                                                                                                  <  1:194329/151‑1 (MQ=255)
   cATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGagaaa                                                                                                                                                 >  1:751320/1‑149 (MQ=255)
      tttACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGagaaaagaa                                                                                                                                             <  1:403127/150‑1 (MQ=255)
          cGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAAccg                                                                                                                                        <  1:967428/151‑1 (MQ=255)
              tatCAAAGTATTCCATCCTCGCTCCCGCCCCTTTCTTTTCATTTTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTTTTTGAGAAGAGAAAAGAAAAccgccg                                                                                                                                     <  2:894872/150‑1 (MQ=255)
                tCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATc                                                                                                                                  >  1:574293/1‑151 (MQ=255)
                         tCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCAcc                                                                                                                         >  2:866511/1‑151 (MQ=255)
                          ccATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCg                                                                                                                        >  2:1034786/1‑151 (MQ=255)
                               cTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCAtt                                                                                                                    >  1:60874/1‑150 (MQ=255)
                               cTCGCCCCCGGCCCTTTTTTTTCTTTTTTTTCGTTGAGGTAACGGTGTTTTTTTACCGGCAACTGGGTTGATTTAGTTTAGTGGTGATCAGTTTTTTTTTACAAGGGTTTGAGAAGAGAAAAAAAAACCGCCGATCCTGTCCACCGCa                                                                                                                      <  2:565560/148‑1 (MQ=255)
                                       cGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAAGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAgg                                                                                                           >  2:776438/1‑151 (MQ=255)
                                        ggCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGTTTATCGGCAACTGGGGTGTTTTACGTTAGGGGTGATCGGGTTTTTTTTAAAACTGTTGGAGAAGAGAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAgg                                                                                                           <  1:91684/150‑1 (MQ=255)
                                          cccGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGt                                                                                                          >  2:659238/1‑149 (MQ=255)
                                            ctttCTTTCCTTTGTTATCTTGGAGGTAACGGTGGGTGTTTATGGGCATCTGGGGTGATTTTCTTTTGTGTTAATCCGGTTTTTTTCACAACTTTTGGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGt                                                                                                       <  1:14274/148‑1 (MQ=255)
                                              tttttttCATTGTTTTCGTTCAGGTAACGGTGGTTGTTTACCGGCAACTGGGGTAATTTACGTTAGTGGTGATCGGTTTTTTTTAAAATTTGTTTGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTgg                                                                                                     <  2:163756/147‑1 (MQ=255)
                                              tttttttCATTGTTATCGTTCAGGTACCTGTGGATGTTTACCGGCAACTGGGGTAATTTACGTTAGTGGTGATCAGTTTTTTTTTACAATTTTTGGAAAAGAAAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGt                                                                                                       <  2:162006/145‑1 (MQ=255)
                                               tCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGa                                                                                                    >  2:143796/1‑150 (MQ=255)
                                                 ttttCATTTTTATGGTTAAGGTAACGGGGGTTGTTTTTCGGCAACTCGGGTGATTTCCTTAGGTGGTGACCGGTTTTTTTTTACAACTGTTGGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGAc                                                                                                   <  2:722854/149‑1 (MQ=255)
                                                 ttttCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACa                                                                                                  >  1:723475/1‑150 (MQ=255)
                                                 ttttAATTGTTATCGTTCAGGTACCGGTGGATGGTTTTCGGCACCTGGGGTAATTTACGTTAGTGGTGATCAGGTTTTTTTTACAATTGTTTGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGa                                                                                                    <  2:195862/148‑1 (MQ=255)
                                                     cATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGAcc                                                                                             >  2:1009136/1‑151 (MQ=255)
                                                       ttGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcg                                                                                            <  1:913657/150‑1 (MQ=255)
                                                          ttATCGTTCAGGTAGCCGTGGATGGTTAGCGGCATCTCGGGTTATTTACGTAAGTGTTGATGAGGTTTTTTTTACAACGGTTGGAAAAAAAAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcggc                                                                                          <  2:606379/149‑1 (MQ=255)
                                                             tCGTTCAGGTACCGGTGGATGGTTATCGCCACCTGGGGTGTTTTACGTTAGTGGTGATCAGGTTTTTTTTACAACTGTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTc                                                                                      <  1:57665/150‑1 (MQ=255)
                                                               gTTCAGGTAACGGTGGATGTTTATCGGCAACGGGGGTGTTTTACGTTAGGGGTGATCGGGTTTTTTTTACAACTGTTGGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCtt                                                                                    <  1:866553/150‑1 (MQ=255)
                                                                  cAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCtt                                                                                    >  2:969383/1‑147 (MQ=255)
                                                                           gTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGc                                                                        >  1:897272/1‑150 (MQ=255)
                                                                             ggATGGTTATCGGCACCGGGGGTGATTTACGTTAGTGGTGATCGGGTTTTTTTTACAATTTTTTGAGAAGAAAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTg                                                                      <  2:192409/150‑1 (MQ=255)
                                                                                tGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGaa                                                                    >  1:252202/1‑149 (MQ=255)
                                                                                 gtttACCGGCAACTGGGGTATTTTACTTTAGGGGTAATCGGGTTTTTTTTACAACTGTTGGAAAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAg                                                                  <  1:274000/148‑1 (MQ=255)
                                                                                  tttttCGGCAACTGGGTTATTTTACTTTAGTGGTGATCGGTTTTTTTTTACAATTTTTGGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGa                                                                     <  2:368530/142‑1 (MQ=255)
                                                                                          cAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTcag                                                         >  1:916756/1‑150 (MQ=255)
                                                                                           aaCTGGGGTGATTTACGTTAGTGGTGTTCAGTTTTTTTTTAAAATTGTTGGAAAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTc                                                           <  1:480042/147‑1 (MQ=255)
                                                                                            aCTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTTGCTGAAAGATTATTcagc                                                        >  2:425996/1‑149 (MQ=255)
                                                                                               cGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGGTTTTTGGCTGAAAGATTATTcagcagat                                                    >  2:410393/1‑148 (MQ=255)
                                                                                                gggTTATTTAGTTAGGTGGTAATCAGTTTTTTTTAACAATTTTTTGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGc                                                  <  2:865672/151‑1 (MQ=255)
                                                                                                    gATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAgg                                               >  2:104357/1‑150 (MQ=255)
                                                                                                      tttACTTTAGTGGTAATCAGTTTTTTTTTACAATTTTTTGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAgg                                               <  2:416139/148‑1 (MQ=255)
                                                                                                      tttACGTTAGTGGTGATCAGGTTTTTTTTAAAACTGTTGGAAAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAg                                                <  1:320859/147‑1 (MQ=255)
                                                                                                                  ggtgATCAGTTTTTTTTTACAATTTTTAGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCtttt                                  <  2:417020/150‑1 (MQ=255)
                                                                                                                  ggtgATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCtttt                                  >  1:41527/1‑150 (MQ=255)
                                                                                                                   gtcagcgggttttttttACATCTTTTTAAAAAGAAAAAAAAAACCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCttt                                   <  2:26923/138‑1 (MQ=255)
                                                                                                                       tCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTTcc                               >  2:892460/1‑148 (MQ=255)
                                                                                                                            tttttttaACAATTTTTTGAGAATAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGa                         <  2:651761/148‑1 (MQ=255)
                                                                                                                            cTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGCAGTGGACAAGACCGGCGGTCTTAAGTTTTTTTGCTGAAAGGTTATTCAGCAGTTGCCGGCATTTTACCTTTTTCCCCTGGACg                       >  2:539133/1‑151 (MQ=255)
                                                                                                                              ttttttACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGtt                     >  1:776403/1‑151 (MQ=255)
                                                                                                                                       atcTGTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAg              <  1:985911/148‑1 (MQ=255)
                                                                                                                                       accTGTTGGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCa               <  1:92418/147‑1 (MQ=255)
                                                                                                                                         ttGTTGGAAAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGAc            <  2:400507/149‑1 (MQ=255)
                                                                                                                                         cTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTTCCGCTGGACGTTCTGTCAGACg           >  2:287429/1‑151 (MQ=255)
                                                                                                                                          tttttGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGAc            <  1:835279/147‑1 (MQ=255)
                                                                                                                                          tttttGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGAc            <  1:150102/147‑1 (MQ=255)
                                                                                                                                            tttGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTTAGTTTTTTGGCTGAAAGATTATTCCGCAGTTTCCGGCATTTTACCTTTTTCCCCTGGACGTTTTGTTCGACGCtt        >  2:59057/1‑151 (MQ=255)
                                                                                                                                            ttagaaaagaaaaAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGc          <  1:305834/143‑1 (MQ=255)
                                                                                                                                              tGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTtct      >  1:111609/1‑151 (MQ=255)
                                                                                                                                                agaaaaaaaAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTCCTGTCAGACGCTtct      <  1:825611/149‑1 (MQ=255)
                                                                                                                                                aaaagagaaAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACg           <  1:969426/142‑1 (MQ=255)
                                                                                                                                                 gaagagaaAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGGTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGAGGTTCTGTCAGACGCTTCTCa    >  1:1014521/1‑150 (MQ=255)
                                                                                                                                                  aagaaaaAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCaa   <  2:624571/150‑1 (MQ=255)
                                                                                                                                                   agagaaAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTTCCGCTGGACGTTCTGTCAGACGCTTCTCaaa  >  1:211524/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                  
CGACATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAA  >  NC_000913/1824908‑1825202

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: