Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,825,0610GA45.6% 65.6 / 28.2 57intergenic (‑124/+79)ves/spycold‑ and stress‑inducible protein/periplasmic ATP‑independent protein refolding chaperone, stress‑induced
Reads supporting (aligned to +/- strand):  ref base G (25/6);  new base A (0/26);  total (25/32)
Fisher's exact test for biased strand distribution p-value = 7.52e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTA‑GTGGTGATCAGGCTTTTTTTACAA‑CTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAA  >  NC_000913/1824914‑1825202
                                                                                                                                                     |                                                                                                                                             
tttACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTA‑GTGGTGATCAGGCTTTTTTTACAA‑CTGTTTGAGAAGagaaaagaa                                                                                                                                             <  1:403127/150‑1 (MQ=255)
    cGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTA‑GTGGTGATCAGGCTTTTTTTACAA‑CTGTTTGAGAAGAGAAAAGAAAAccg                                                                                                                                        <  1:967428/151‑1 (MQ=255)
        tatCAAAGTATTCCATCCTCGCTCCCGCCCCTTTCTTTTCATTTTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTA‑GTGGTGATCAGGCTTTTTTTACAA‑CTTTTTGAGAAGAGAAAAGAAAAccgccg                                                                                                                                     <  2:894872/150‑1 (MQ=255)
          tCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTA‑GTGGTGATCAGGCTTTTTTTACAA‑CTGTTTGAGAAGAGAAAAGAAAACCGCCGATc                                                                                                                                  >  1:574293/1‑151 (MQ=255)
                   tCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTA‑GTGGTGATCAGGCTTTTTTTACAA‑CTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCAcc                                                                                                                         >  2:866511/1‑151 (MQ=255)
                    ccATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTA‑GTGGTGATCAGGCTTTTTTTACAA‑CTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCg                                                                                                                        >  2:1034786/1‑151 (MQ=255)
                         cTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTA‑GTGGTGATCAGGCTTTTTTTACAA‑CTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCAtt                                                                                                                    >  1:60874/1‑150 (MQ=255)
                         cTCGCCCCCGGCCCTTTTTTTTCTTTTTTTTCGTTGAGGTAACGGTGTTTTTTTACCGGCAACTGGGTTGATTTAGTTTA‑GTGGTGATCAGTTTTTTTTTACAA‑GGGTTTGAGAAGAGAAAAAAAAACCGCCGATCCTGTCCACCGCa                                                                                                                      <  2:565560/148‑1 (MQ=255)
                                 cGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTA‑GTGGTGATCAAGCTTTTTTTACAA‑CTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAgg                                                                                                           >  2:776438/1‑151 (MQ=255)
                                  ggCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGTTTATCGGCAACTGGGGTGTTTTACGTTA‑GGGGTGATCGGGTTTTTTTTAAAA‑CTGTTGGAGAAGAGAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAgg                                                                                                           <  1:91684/150‑1 (MQ=255)
                                    cccGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTA‑GTGGTGATCAGGCTTTTTTTACAA‑CTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGt                                                                                                          >  2:659238/1‑149 (MQ=255)
                                      ctttCTTTCCTTTGTTATCTTGGAGGTAACGGTGGGTGTTTATGGGCATCTGGGGTGATTTTCTTTT‑GTGTTAATCCGGTTTTTTTCACAA‑CTTTTGGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGt                                                                                                       <  1:14274/148‑1 (MQ=255)
                                        tttttttCATTGTTTTCGTTCAGGTAACGGTGGTTGTTTACCGGCAACTGGGGTAATTTACGTTA‑GTGGTGATC‑GGTTTTTTTTAAAATTTGTTTGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTgg                                                                                                     <  2:163756/147‑1 (MQ=255)
                                        tttttttCATTGTTATCGTTCAGGTACCTGTGGATGTTTACCGGCAACTGGGGTAATTTACGTTA‑GTGGTGATCAGTTTTTTTTTACAA‑TTTTTGGAAAAGAAAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGt                                                                                                       <  2:162006/145‑1 (MQ=255)
                                         tCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTA‑GTGGTGATCAGGCTTTTTTTACAA‑CTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGa                                                                                                    >  2:143796/1‑150 (MQ=255)
                                           ttttCATTTTTATGGTTAAGGTAACGGGGGTTGTTTTTCGGCAACTCGGGTGATTTCCTTAG‑GTGGTGACCGGTTTTTTTTTACAA‑CTGTTGGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGAc                                                                                                   <  2:722854/149‑1 (MQ=255)
                                           ttttCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTA‑GTGGTGATCAGGCTTTTTTTACAA‑CTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACa                                                                                                  >  1:723475/1‑150 (MQ=255)
                                           ttttAATTGTTATCGTTCAGGTACCGGTGGATGGTTTTCGGCACCTGGGGTAATTTACGTTA‑GTGGTGATCAGGTTTTTTTTACAA‑TTGTTTGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGa                                                                                                    <  2:195862/148‑1 (MQ=255)
                                               cATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTA‑GTGGTGATCAGGCTTTTTTTACAA‑CTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGAcc                                                                                             >  2:1009136/1‑151 (MQ=255)
                                                 ttGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTA‑GTGGTGATCAGGCTTTTTTTACAA‑CTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcg                                                                                            <  1:913657/150‑1 (MQ=255)
                                                    ttATCGTTCAGGTAGCCGTGGATGGTTAGCGGCATCTCGGGTTATTTACGTAA‑GTGTTGATGAGGTTTTTTTTACAA‑CGGTTGGAAAAAAAAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcggc                                                                                          <  2:606379/149‑1 (MQ=255)
                                                       tCGTTCAGGTACCGGTGGATGGTTATCGCCACCTGGGGTGTTTTACGTTA‑GTGGTGATCAGGTTTTTTTTACAA‑CTGTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTc                                                                                      <  1:57665/150‑1 (MQ=255)
                                                         gTTCAGGTAACGGTGGATGTTTATCGGCAACGGGGGTGTTTTACGTTA‑GGGGTGATCGGGTTTTTTTTACAA‑CTGTTGGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCtt                                                                                    <  1:866553/150‑1 (MQ=255)
                                                            cAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTA‑GTGGTGATCAGGCTTTTTTTACAA‑CTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCtt                                                                                    >  2:969383/1‑147 (MQ=255)
                                                                     gTGGATGGTTATCGGCAACTCGGGTGATTTACGTTA‑GTGGTGATCAGGCTTTTTTTACAA‑CTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGc                                                                        >  1:897272/1‑150 (MQ=255)
                                                                       ggATGGTTATCGGCACCGGGGGTGATTTACGTTA‑GTGGTGATCGGGTTTTTTTTACAA‑TTTTTTGAGAAGAAAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTg                                                                      <  2:192409/150‑1 (MQ=255)
                                                                          tGGTTATCGGCAACTCGGGTGATTTACGTTA‑GTGGTGATCAGGCTTTTTTTACAA‑CTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGaa                                                                    >  1:252202/1‑149 (MQ=255)
                                                                           gtttACCGGCAACTGGGGTATTTTACTTTA‑GGGGTAATCGGGTTTTTTTTACAA‑CTGTTGGAAAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAg                                                                  <  1:274000/148‑1 (MQ=255)
                                                                            tttttCGGCAACTGGGTTATTTTACTTTA‑GTGGTGATCGGTTTTTTTTTACAA‑TTTTTGGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGa                                                                     <  2:368530/142‑1 (MQ=255)
                                                                                    cAACTCGGGTGATTTACGTTA‑GTGGTGATCAGGCTTTTTTTACAA‑CTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTcag                                                         >  1:916756/1‑150 (MQ=255)
                                                                                     aaCTGGGGTGATTTACGTTA‑GTGGTGTTCAGTTTTTTTTTAAAA‑TTGTTGGAAAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTc                                                           <  1:480042/147‑1 (MQ=255)
                                                                                      aCTCGGGTGATTTACGTTA‑GTGGTGATCAGGCTTTTTTTACAA‑CTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTTGCTGAAAGATTATTcagc                                                        >  2:425996/1‑149 (MQ=255)
                                                                                         cGGGTGATTTACGTTA‑GTGGTGATCAGGCTTTTTTTACAA‑CTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGGTTTTTGGCTGAAAGATTATTcagcagat                                                    >  2:410393/1‑148 (MQ=255)
                                                                                          gggTTATTTA‑GTTAGGTGGTAATCAGTTTTTTTTAACAA‑TTTTTTGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGc                                                  <  2:865672/151‑1 (MQ=255)
                                                                                              gATTTACGTTA‑GTGGTGATCAGGCTTTTTTTACAA‑CTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAgg                                               >  2:104357/1‑150 (MQ=255)
                                                                                                tttACTTTA‑GTGGTAATCAGTTTTTTTTTACAA‑TTTTTTGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAgg                                               <  2:416139/148‑1 (MQ=255)
                                                                                                tttACGTTA‑GTGGTGATCAGGTTTTTTTTAAAA‑CTGTTGGAAAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAg                                                <  1:320859/147‑1 (MQ=255)
                                                                                                            ggtgATCAGTTTTTTTTTACAA‑TTTTTAGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCtttt                                  <  2:417020/150‑1 (MQ=255)
                                                                                                            ggtgATCAGGCTTTTTTTACAA‑CTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCtttt                                  >  1:41527/1‑150 (MQ=255)
                                                                                                             gtcagcgggttttttttACAT‑CTTTTTAAAAAGAAAAAAAAAACCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCttt                                   <  2:26923/138‑1 (MQ=255)
                                                                                                                 tCAGGCTTTTTTTACAA‑CTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTTcc                               >  2:892460/1‑148 (MQ=255)
                                                                                                                      tttttttaACAA‑TTTTTTGAGAATAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGa                         <  2:651761/148‑1 (MQ=255)
                                                                                                                      cTTTTTTTACAA‑CTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGCAGTGGACAAGACCGGCGGTCTTAAGTTTTTTTGCTGAAAGGTTATTCAGCAGTTGCCGGCATTTTACCTTTTTCCCCTGGACg                       >  2:539133/1‑151 (MQ=255)
                                                                                                                        ttttttACAA‑CTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGtt                     >  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)
                                                                                                                                    tttttaaaaaaaaaaaaaaaacCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGAc            <  2:164578/136‑1 (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)
                                                                                                                                                     |                                                                                                                                             
TTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTA‑GTGGTGATCAGGCTTTTTTTACAA‑CTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAA  >  NC_000913/1824914‑1825202

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: ‑