Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009135,1590GA38.3% 75.2 / 25.3 60intergenic (+139/‑75)thrC/yaaXL‑threonine synthase/DUF2502 family putative periplasmic protein
Reads supporting (aligned to +/- strand):  ref base G (29/8);  new base A (0/23);  total (29/31)
Fisher's exact test for biased strand distribution p-value = 3.58e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCG  >  NC_000913/5013‑5305
                                                                                                                                                      |                                                                                                                                                   
atcaGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATa                                                                                                                                                 <  1:418360/151‑1 (MQ=255)
   aGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATaaca                                                                                                                                             <  2:1051707/151‑1 (MQ=255)
      aaaaTCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATaacaac                                                                                                                                           <  1:486009/150‑1 (MQ=255)
      aaaaTCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAAcc                                                                                                                                          >  2:478870/1‑151 (MQ=255)
                      ctcAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGt                                                                                                                          >  2:469360/1‑151 (MQ=255)
                          aTCAGGCGGGTTTTGTTTTTAGGCACCCGGGTTTTTTTATGAAAAAATTAGGGAAAAAAATGACAGGAAAAAAGGAAAATTTCTAAAAAAAGGCGTTAACTTAGAGATTAGGATTGGGGAGAAAAACAACCGCCGTTCTCATCGAGTAAtc                                                                                                                      <  2:832030/151‑1 (MQ=255)
                              ggCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCg                                                                                                                  >  2:901554/1‑151 (MQ=255)
                                cgggttttGTTTTTAGGCAGCCCGGTTTTTTAATGAAAAAATTATGGAGAAAAATGAAGGGGAAAAAGGAAAATTTCTAAAAAAAGGGGTTAACTTAAAATTGGGGATTGCGAAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCgg                                                                                                                 <  1:295785/145‑1 (MQ=255)
                                   agTTTGCTTTTGTGCAGGCCGGCTTTTTTATGAAGAAATTATGGTGAAAAATGACAGGTAAAAATGAGTAATTCTCAATAAATGCGGTAACTTAGAGAGTAGGAGTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGatat                                                                                                             >  1:481999/2‑151 (MQ=255)
                                     ttttCTTTAATGCAGCCGGTTTTTTTAAGAAAAAAATTAGGGAAAAAAAGAAGGGGAAAAAAGGAAAAATTCACAAAAAATGCGTTCTCATAGAGTTGAGGTGGGCGAAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGata                                                                                                              <  1:442570/148‑1 (MQ=255)
                                     tttGTTTTTAAGCAGCCCGTTTTTTTTATAAAAAAATTATGGAAAAAAAAAACAGGAAAAAAGGAGAAATTCTAAATATATGGGTTAACTTAGAAATAAGGATTGCGGAAAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATc                                                                                                            <  2:580127/150‑1 (MQ=255)
                                         cTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGAccc                                                                                                       >  1:169440/1‑151 (MQ=255)
                                           tttaggcgtcccGTCTTTTTTAAGAAAAAATTTTGAAAAAAATGGAAAGGAAAAAAAGAGAAATCATCAAAAACGGGGTTATTTTAGGGATTATGGTGGCGAAAAACAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCa                                                                                                      <  2:532310/141‑1 (MQ=255)
                                             tagcaaGCCGGGTTTTTTTAAAAAAAAATTATGGAAAAAAAAGACAGGAAAAAAGGAAAAATTCTAAATAAATGCGTTAACTAAGAGATTAGGATGGGGGAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATa                                                                                                    <  2:10977/145‑1 (MQ=255)
                                             tATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATa                                                                                                    >  2:795496/1‑150 (MQ=255)
                                              aTGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATaa                                                                                                   >  1:950117/1‑150 (MQ=255)
                                               gacaccctgttttttttAGAAAAAAATTAGGAAAAAAAAGAAAGGAAAAAAAGAAAAACTTCTCAATAATTGCGTTAAATTAAAATTTAGATTGGCGAAAAAAAACACCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACg                                                                                                 <  1:493377/140‑1 (MQ=255)
                                                  aGCCCGTTTTGTTTATAAAAAAATTATGGAAAAAAAGGAAAGGAAAAAAGGAAAAATTATCAATAAGTGCGGTAAATAAGAGATTAGGAGGGCGAAGAATACCAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAAc                                                                                                  <  1:940768/147‑1 (MQ=255)
                                                  aGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGc                                                                                              >  2:684153/1‑151 (MQ=255)
                                                    ccgggtttttttAAAAAGAAATTAGGAAAAAAAATGACGGGAAAAAAGAAAAAATTCTAAATAAACGCGGTAACTAAGAGAGAAGGATGGGGGAAAAAAACAGCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACg                                                                                                 <  1:302050/140‑1 (MQ=255)
                                                         tttttttAAAAAAAAATTAGGAAGAAAAAGGACGGAAAAAAAGGAAAAATTCTAAAAAAATGGGTTAACTAGGAGATTAGGATTGCGAAAAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGa                                                                                         <  1:1090574/148‑1 (MQ=255)
                                                            ttttATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATaaa                                                                                     >  2:886386/1‑150 (MQ=255)
                                                            ttttAGAAAAAAATTATGGAAAAAAATGACAGGAAAAAAGGAAAATTTCAAAAAAAATGCGTTAATTTAGAGATTAGGATTGGGAAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATaaa                                                                                     <  2:693075/150‑1 (MQ=255)
                                                                 tGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGt                                                                               >  1:40159/1‑151 (MQ=255)
                                                                    aaaaaaTTAGGAAAAAAAATGAAAGGAAAAAAGGAAAAATTCTAAAAAAAGGCGGTAACTTAGAGATTAGGATGGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGt                                                                               <  2:602446/146‑1 (MQ=255)
                                                                       aaaTTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCt                                                                          >  1:93680/1‑151 (MQ=255)
                                                                        aaTGATGGAAAAAAAGGACGGAAAAAAGGGAGAAATTCCAAATAAGTGCGTTAATTAAAAATTGAGTAGGGGTGAGAAAAAAAACCGCCGTTCTCAGCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCtg                                                                         <  1:301904/151‑1 (MQ=255)
                                                                          ttATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCtgtg                                                                       >  2:1079367/1‑151 (MQ=255)
                                                                          ttATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCtgtg                                                                       >  1:449245/1‑151 (MQ=255)
                                                                            aTGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGa                                                                      >  2:625068/1‑150 (MQ=255)
                                                                             ggaagaAAAAGCAAAGGAAAAAAAGAAAACTTCTAAAAAAATGCGGAATTTTAGATATGAGATTGGGGAAAAAAACCACCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGTCAATGATAAAAGGAGTAACCtggg                                                                       <  1:230151/145‑3 (MQ=255)
                                                                                 aaaaaaGGAGAGGAAAAAAGGAAAAATCCTCAAAAAAGGCGTAACCTAAGAAATTAGAATTGAGGAAAAAAAAACCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGaaaa                                                                   <  1:451991/147‑1 (MQ=255)
                                                                                  aaaaaTGAGGGAAAAAAGGGAAAAATTCTAAAAAAAGGGGGTAATTAAGAAATTAGGTTGGGGAAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGa                                                                <  2:697503/150‑1 (MQ=255)
                                                                                  aaaaaTGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGaaaa                                                                   >  2:827380/1‑147 (MQ=255)
                                                                                  aaaaaTGACAGGAAAAAAGGATAAATTCACAATAAAGGCGGTACCTTAGAGATTAGGATTGCGGAGAAAAACACCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGaaa                                                                    <  2:405676/146‑1 (MQ=255)
                                                                                     aaagaaggggAAAAAGGAGAAATTTTAAATAAGGGCGTTATTGAGAAGATTGGAATGGCGAAAAACAACAGCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCa                                                            <  2:1127577/144‑1 (MQ=255)
                                                                                     aaTGAAAGGAAAAAAGGAAAAATTCTAAAAAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCa                                                            <  2:800321/151‑1 (MQ=255)
                                                                                       tGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAAt                                                          >  1:988559/1‑151 (MQ=255)
                                                                                             ggAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCTAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTAt                                                      >  2:722908/1‑149 (MQ=255)
                                                                                               aaaaaGGAAAATTTCAAAAAAACTCCGCTAATTTGAAGATTAGGATGGTGAAAAATAAGAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTa                                                       <  1:769566/146‑1 (MQ=255)
                                                                                                aaaaGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTAc                                                 >  2:161404/1‑151 (MQ=255)
                                                                                                aaaaGGAAAAATTCACAATAAAGGCGTTATCTAAGAGTTTAGGTTGGCGAAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTAc                                                 <  1:280967/151‑1 (MQ=255)
                                                                                                       gaAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACt                                          >  1:451993/1‑151 (MQ=255)
                                                                                                                 aaaaaaTGCGGTAATTTAGAGATTAGGATTGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTgg                                  <  2:987450/147‑1 (MQ=255)
                                                                                                                 aaTAAAGGCGTTAATTGAGATATTGGGTTTGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGt                                 <  2:1030871/151‑1 (MQ=255)
                                                                                                                 aaTAAAGGCGGAAACTTAGAGATTAGGATTGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTg                                   <  2:16211/149‑1 (MQ=255)
                                                                                                                    aaaTGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCt                              >  2:384198/1‑151 (MQ=255)
                                                                                                                       tGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATAGAGAAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTgg                            >  1:998845/1‑150 (MQ=255)
                                                                                                                             aaCTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTcc                     >  2:812321/1‑151 (MQ=255)
                                                                                                                              aCTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTccc                    >  1:142384/1‑151 (MQ=255)
                                                                                                                                 ttaaagttttgggATTGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCg                         <  2:1108855/134‑1 (MQ=255)
                                                                                                                                 tAGAGATTAGGATTGCGGAGAATAAGAACCGCCGTAATCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATACAAGGAGTAACCAGTGAAAAAGATGCAATCTATCGTACTCCCACTTTCCCGGGATCTGGTGGATCCCa                   >  1:49866/1‑149 (MQ=255)
                                                                                                                                 tAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATg                 >  1:561123/1‑151 (MQ=255)
                                                                                                                                   gagaTTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGgc               >  1:643703/1‑151 (MQ=255)
                                                                                                                                     agaTTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGgc               >  1:659468/1‑150 (MQ=255)
                                                                                                                                      agtgattgtatGCGAAAAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCAta                 <  1:41897/138‑2 (MQ=255)
                                                                                                                                           aggtgggaggAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCAt                  <  1:309834/134‑1 (MQ=255)
                                                                                                                                            ggATTGCGGAGAATAACAACCGCCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGcac          >  1:120651/1‑149 (MQ=255)
                                                                                                                                               tgggaaagaAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTgg                            <  1:459955/122‑1 (MQ=255)
                                                                                                                                               tggCGGAAAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAgg       <  2:216224/147‑1 (MQ=255)
                                                                                                                                                  cGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCg  >  1:731091/1‑151 (MQ=255)
                                                                                                                                                      |                                                                                                                                                   
ATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCG  >  NC_000913/5013‑5305

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