Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009135,1590GA45.2% 63.8 / 37.0 62intergenic (+139/‑75)thrC/yaaXL‑threonine synthase/DUF2502 family putative periplasmic protein
Reads supporting (aligned to +/- strand):  ref base G (29/5);  new base A (0/28);  total (29/33)
Fisher's exact test for biased strand distribution p-value = 7.17e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTG  >  NC_000913/5013‑5303
                                                                                                                                                  |                                                                                                                                                
atcaGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATa                                                                                                                                              <  2:13359/151‑1 (MQ=255)
  caGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATaac                                                                                                                                            <  2:297081/151‑1 (MQ=255)
  caGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATaa                                                                                                                                             >  1:249761/1‑150 (MQ=255)
      aaaaTCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATaacaac                                                                                                                                         <  2:879989/150‑1 (MQ=255)
      aaaaTCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAAcc                                                                                                                                        >  2:8770/1‑151 (MQ=255)
        aaTCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAAccgc                                                                                                                                      >  1:129079/1‑151 (MQ=255)
         aTCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAAccgcc                                                                                                                                     >  1:658452/1‑151 (MQ=255)
              ttCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTtct                                                                                                                                >  2:346479/1‑151 (MQ=255)
              ttCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTtct                                                                                                                                >  2:165619/1‑151 (MQ=255)
                         aaTCAGGCGGGGTTTGTTTTAACGCAGCCCGGTTTTTTTATGAAAAAATTAGGGAAAAAAAGGACGGGGAAAAAGGAAAAATTCTCAAAAAAGGCGGTACTTTAGAGATTAGGATGGCGGAAAATAACAACCGCCGTTCTCATCGAGTaa                                                                                                                      <  2:119686/150‑1 (MQ=255)
                         aaTCAGGCCGGGTTTTTTTTTAGGCCGCCGGGTTTTTTTATAAAAAAATTAGGGAAAAAAATGACAGGAAAAAAGGAAAAATTCTAAAAAAAGGCGTTAATTAAGAAATTAGGATTGCGGAAAAAAACAACCGCCGTTCTCATCg                                                                                                                           <  1:916524/145‑1 (MQ=255)
                            cAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCTGAGAATAACAACCGCCGTTCTCATCGAGTAAtctc                                                                                                                  >  2:794855/1‑151 (MQ=255)
                              ggcgggttttttttttATGCACCCCGGTTTTTTTATAAAAAAATTATGAAAAAAAAGGACGGGAAAAAAGGAAAATTTCTCAAAAAAGGGGGTATCTTAGAGATAAGTTTTGGGGAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTcc                                                                                                                 <  2:599380/138‑1 (MQ=255)
                              ggCCGTGTTTGTTTTTATACACCCCGGTTTTTTTATAAAAAAATTATGGAAAAAAAGGACAGGAAAAAAGGAGAATTTCTCAATAAATGCGGTATTTAAGAGATGAGGATGGCGGAAAATAACAACCGCCGTTCTCATCGAGTAATCTCCg                                                                                                                <  1:421176/151‑1 (MQ=255)
                              ggCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCg                                                                                                                >  1:201536/1‑151 (MQ=255)
                                 cGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGat                                                                                                             >  1:173949/1‑151 (MQ=255)
                                     tttGTTTTTATGCAGCCCGGTTTTTTTAGAAAAAAATTATGGAAAAAAAGGAGAGGAAAAAAGGAAAAATTCTAAAAAATTGGGGTATTTTAGAGATTAGGTTGGGGAAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGatat                                                                                                           <  1:610359/149‑1 (MQ=255)
                                        gCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGAcc                                                                                                      >  2:62282/1‑151 (MQ=255)
                                         cTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGCAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGAATGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGAccc                                                                                                     >  1:716706/1‑151 (MQ=255)
                                          ttttATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGAccc                                                                                                     >  1:1001500/1‑150 (MQ=255)
                                           tttATCCACCCCGTTTTTTTAATAAAAAAATTATGAAAAAAAAGGACAGGAAAAAAGGAAAAATTCTAAAAAAAGGCGGTAACTAAGAGATTAGGATTGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCa                                                                                                    <  2:96176/150‑1 (MQ=255)
                                            ttATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCAt                                                                                                   >  1:844501/1‑150 (MQ=255)
                                              acgccccccgttttttttAAGAAAAAATTGGGAAGAAAAAGAAGAGGAAAAAAGAAAAAATTCTAAATAAATGCATTACCTTAGAATTTAGGATTGCGGAAAAAACAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAAc                                                                                                <  2:710021/139‑1 (MQ=255)
                                                gCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACgg                                                                                              >  2:804815/1‑151 (MQ=255)
                                                         cTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGAt                                                                                      >  2:133955/1‑150 (MQ=255)
                                                          ttttttATGAAAAAATTAGGGAAAAAAAGGACAGGGAAAAAGAAAAAATCCTCAAAAAATGGGGTATTTGAGAGATTGGGGTGGCGGAGAAAAACAGCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATa                                                                                     <  1:990723/150‑1 (MQ=255)
                                                             tttATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATaaa                                                                                   >  2:531390/1‑149 (MQ=255)
                                                               tATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAgg                                                                                >  2:383136/1‑150 (MQ=255)
                                                                aTGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAgg                                                                                >  1:61793/1‑149 (MQ=255)
                                                                   aagaaATTATGGAAAAAAATGACAGGAAAAAAGGAAAAATTCTAAAAAAATGCGTTAATTAAGAGTTTAGGATTGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGt                                                                             <  1:472759/149‑1 (MQ=255)
                                                                          tATGGAAAAAAATGACAGGAAAAAAGGAAAATTTATAAAAAAAGGCGGTAATTTAGAGATTAGGATTGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCtgtg                                                                     <  2:206056/150‑1 (MQ=255)
                                                                             gaaaaaaaaTGAAAGGAAAAAAGGAAAAATTCTAAAAAAATGGGTTAACTTAGAGATTAGGATTGGGGAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGa                                                                    <  2:701436/144‑1 (MQ=255)
                                                                               aaaaaaaTGAAGGGAAAAAAGGAAAAATTCTAAATAAAGGCGGTAACTTAGAGATTAGGATGGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCtgtg                                                                     <  2:239186/143‑1 (MQ=255)
                                                                               aaaaaaaTGAAGGGAAAAAAGAAAAAATTCTCAAAAAAGGCGTTAACTTAGAAATTAGATTTGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAg                                                               <  2:941699/149‑1 (MQ=255)
                                                                                 aaaaaTAAGAGGAAAAAAGAAAAAATTCAAAAAAAAGGCGGTATTTAAGAGATTAGGATGGGGGAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGAt                                                             <  2:77265/151‑1 (MQ=255)
                                                                                  aaaaTGACAGGAAAAAAGGAAAAATTCAAAAAAAATGCGGTAACTTAGAGATTAGGATGGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATg                                                            <  2:249821/151‑1 (MQ=255)
                                                                                  aaaaGGACAGGAAAAAAGGAAAAATTCTAAATAAAGGCGGTAATTTAGAGATGAGGATGGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATg                                                            <  2:525406/151‑1 (MQ=255)
                                                                                  aaaaGGACAGGAAAAAAGGAAAAATTCTAAAAAAAGGCGGTAACTAAGAGATGAGGTTTGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATg                                                            <  2:853486/151‑1 (MQ=255)
                                                                                       gACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATc                                                       >  2:421341/1‑151 (MQ=255)
                                                                                           ggaaaaaaGGAAAAATTCTAAAAAAATGCGGTAACTAAGAGATTAGGATGGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCt                                                      <  2:456257/145‑1 (MQ=255)
                                                                                              aaaaaGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGt                                                 >  2:425170/1‑150 (MQ=255)
                                                                                              aaaaaGGAAAAATTCTAAAAAAATGCGTTAACTTAGAAATTAGGATTGGGAAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATc                                                   <  1:795988/148‑1 (MQ=255)
                                                                                              aaaaaGGAAAAATTCTAAAAAAAGGCGGTAATTTAGAGATAAGGATGGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCg                                                  <  2:196373/149‑1 (MQ=255)
                                                                                                aaaGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACt                                              >  2:88247/1‑151 (MQ=255)
                                                                                                 aaGGAGAAATTCTCAAAAAATGCGATAACTTAGAGATGAGGATCGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTc                                             <  2:276254/151‑1 (MQ=255)
                                                                                                 aaGGAAAAATTCTCAAAAAAGGCGGAACTTTAGAAATGAGGATTGCGAAAAATAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTc                                             <  1:547432/151‑1 (MQ=255)
                                                                                                  aGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCg                                            >  2:401417/1‑151 (MQ=255)
                                                                                                   ggAAAAATTCTCAATAAAGGCGTTAACTTAGAGATGAGGATGGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCg                                            <  1:353147/150‑1 (MQ=255)
                                                                                                       aaaTTCTAAAAAAGTGGGGTAACTAAGAGATTAGTATTGGGGAAAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACtt                                       <  2:124913/151‑1 (MQ=255)
                                                                                                       aaaTTCAAAAAAAATGAGATAATTAAGATATTGGGATTGGGAAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACtt                                       <  1:165567/151‑1 (MQ=255)
                                                                                                         aTTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTc                                     >  2:895611/1‑151 (MQ=255)
                                                                                                          ttCTAAAAAAATGGGGTAATTTAGAGATTAGGATGGGGGAAAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTc                                     <  2:129127/150‑1 (MQ=255)
                                                                                                                aTAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGtt                              >  2:451458/1‑151 (MQ=255)
                                                                                                                  aaaTGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCt                            >  1:782685/1‑151 (MQ=255)
                                                                                                                   aaTGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTg                           >  2:970374/1‑151 (MQ=255)
                                                                                                                   aaTGCGGTAACTTAGAGATTAGGATTGCGGAAAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTg                           <  2:201438/151‑1 (MQ=255)
                                                                                                                    aTGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTgg                          >  2:146826/1‑151 (MQ=255)
                                                                                                                           aattaagagATTAGGATTGTGAAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCt                     <  1:794739/144‑1 (MQ=255)
                                                                                                                            aCTTAGAGATTAGGATTGCGGAAAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTccc                  <  2:113080/151‑1 (MQ=255)
                                                                                                                               tAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATg               >  1:422933/1‑151 (MQ=255)
                                                                                                                                        ggATTGCGAAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGcaca       <  1:682700/150‑1 (MQ=255)
                                                                                                                                            tGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGCCAGCACAGGCTg  >  1:491065/1‑151 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
ATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTG  >  NC_000913/5013‑5303

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