Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,151,6860GT42.1% 43.4 / 18.4 57intergenic (+40/+25)yegJ/yegKuncharacterized protein/ser/thr phosphatase‑related protein
Reads supporting (aligned to +/- strand):  ref base G (20/13);  new base T (8/16);  total (28/29)
Fisher's exact test for biased strand distribution p-value = 6.10e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.03e-01

GATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCG  >  NC_000913/2151538‑2151835
                                                                                                                                                         |                                                                                                                                                       
gATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAAttgtt                                                                                                                                                       <  1:311389/151‑1 (MQ=255)
    tataCGCGCGGGGATAAAATTTACGGCGGTTATAACAATTATTCTTTTTTTTTTACCTATCCAAAAAAAGAAGACGACGAGCTAAAAGCAAAAATAGTTTGCTTATCAGCCGCTATTTTTTTTTTTTTAAATGGGGGGTTGAAttttt                                                                                                                                                       >  2:1048446/1‑148 (MQ=255)
     ataCGCGCGGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCtt                                                                                                                                                  <  2:992500/151‑1 (MQ=255)
            cgGGGATAAAATTTACGGCGGTTTTTCCATTTATCCCTTTTTTTTCGCCCATCCCAAAAAAGAAAGCGACGAGCTACAAGCAAAACTTGTTTTCTCATCAGCCGATATTTTTTTTTTTTAAAATGGGGGGTGGAATTTTGtttttctaccc                                                                                                                                           >  1:1133451/1‑144 (MQ=255)
            cgGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGc                                                                                                                                           <  1:778075/151‑1 (MQ=255)
                     aaTTTACGGCGTTTATACCATTGATCCTTTATTGGTCCCCTATCAAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTAAATTGTTGTCTTCAACCGCCATGAGaaa                                                                                                                                  <  2:161064/151‑1 (MQ=255)
                           cggcggTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGGTGTAATTTTTTTTTTTAACccccctaggaa                                                                                                                                   >  1:867005/1‑138 (MQ=255)
                                gTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGcc                                                                                                                      >  1:745436/1‑151 (MQ=255)
                                 ttATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTGGGAATTTTTTCTTTAACCGACATGAGaaaatatatcttccc                                                                                                                      >  2:354279/1‑139 (MQ=255)
                                      ccATTGTTCTTTTATTGGTCCCTTATCAAAAAAAAGAAGCCGACGAGCTAAGAGAAAAATTAGTTCGCTAATCAGCCTTTATATTTTTTTTTTAAACGGGGGGTTGAATTGTTGTTTTCAACCGCCATGAGAAACTCACTCTGCCcaca                                                                                                                  <  1:1689/149‑1 (MQ=255)
                                      caaTTGATCTTTTATTGGTCCCCTATCAAAAAGAAGAAGCCGACGAGTTAAGAGAAAAATTATTTCGTTAATCAGCCTTTATTTTTATTGTTAAAACTGGCGTTTAAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCcaca                                                                                                                  <  1:810563/147‑1 (MQ=255)
                                         ttGATCTTTTATTGGTCCCCTACCAAAAAAAAGAAGCCGACGAGCAAAGAGAAAAATTAGTCCGCTAATCAGCCTTTATATTTATTTTTTAACCGGGCGGTTGAATTTTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAgc                                                                                                                <  2:418414/148‑1 (MQ=255)
                                         ttGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTTTTCAACCGCCATGAGAAAATCACTCTGCCCACACCGcc                                                                                                             >  1:463126/1‑151 (MQ=255)
                                              ccTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTTAATTTTTTTTTTCAAATGCCATGAAAAAATCACTCCTCCCcacacgac                                                                                                              >  2:352487/1‑142 (MQ=255)
                                              ccTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGATGAATTGTTGTCTTCAACCGCCCTGAGAAAATCAATCTGACCCCAACGcc                                                                                                             >  2:56044/1‑146 (MQ=255)
                                                tttATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTTTTTTCTTTAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGcc                                                                                                       >  2:237572/1‑150 (MQ=255)
                                                 tttttgtcccccTATCAAAAAAAAAAAGCCGAGAACAAAAAAAAAAATTTAGTTCCTAAACCACCCTTTTTTTTTTTTGTTAAAACGGGCGGTTAAATTGTTGTCTCAAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCa                                                                                                      <  2:448919/140‑1 (MQ=255)
                                                        tCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTTTTCAAACGCCATGAGGAAATCACTCTGCCCACAGCGCCAGTGGCCAAGtttt                                                                                               >  2:527434/1‑150 (MQ=255)
                                                              accaaaaaGAAGAAGCCGACGGGCTAAGAGAAAAATTATTTCGCTAATCAGCCTTTATTTTTTTTGTTTAAACGGGGGTTTAATTTTTTTTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTa                                                                                              <  2:157232/141‑1 (MQ=255)
                                                                  aaaaGAAGACGCCGACGAGATAAGAGAAAAATTAGTTCGCTAATCAGCTTTTATATTTATTGTTTAATCGGGGGTTTAAATTTTTGTCTCAATCCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTatcat                                                                                          <  1:1064280/145‑1 (MQ=255)
                                                                    aagaagaagCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGt                                                                                   >  1:741982/1‑150 (MQ=255)
                                                                     agaagaagCCGACGAGCAAAGAGAAAAATTAGTCCGCTAATCAGCCTTTATTTTTTTTGTTAAAACTGGCGGTTAAATTTTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTAc                                                                                 <  1:390174/151‑1 (MQ=255)
                                                                      gaagaagCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTTTCTTCAACCCCCATGAGAAACTCACTCTGCCCCCAGCGCCAGTGCCAACGTTTTTTCATCATCCGTAc                                                                                 >  2:1123228/1‑150 (MQ=255)
                                                                       aaaaagCGGACAACCAAAGAAAAAAATTAGTTCGTTAATCGGCCTTTTTTTTTTTTTTTTAACCGGGCGTTTGATTTGTTGTTTCAAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACg                                                                                <  1:1202647/147‑1 (MQ=255)
                                                                          aagCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAAAGTTTTATCATCATCCGTACGCTCg                                                                            >  1:1177057/1‑151 (MQ=255)
                                                                              cgacgaGCAAAAAGAAAATTTAGTCCGATAACCACCCTTTATTTTTATTTTTTAACCGGGGGGTTAATTTGTTGTCTAAAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGttt                                                                         <  2:163288/150‑1 (MQ=255)
                                                                               gacgaGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTTAACCGCCCTGAGAAACTCACTCTTCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGgtt                                                                         >  2:602563/1‑149 (MQ=255)
                                                                                   aGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTTTTGTCTTCAACCCCCATTAGAAACTCACTCTGCCCACAGCGCCAGTGCCCAAGTTTTTTTATCATCCCTACGCTCGTTTaacccc                                                                   >  2:784681/1‑146 (MQ=255)
                                                                                       aagaaaaaaaTTATTCCGTTAATCGGCCTTTATTTTTATTTTTAAACCGGGGGTTTAAATTTTTGTCTTAAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGc                                                              <  2:1218371/145‑1 (MQ=255)
                                                                                       aaaagaaaaaTTATTCCGTTAACCAGCCTTTATTTTTTTTTTTTAAACGGGGGGTTAATTTTTTTTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCgg                                                                 <  1:411395/143‑1 (MQ=255)
                                                                                             aaaaTTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGAc                                                         >  1:90207/1‑150 (MQ=255)
                                                                                                  ttAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTTAATTGTTGTCTTCAACCGCCATGAGAAAATCACTCTGCCCCCAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGGCaaaa                                                     >  1:875294/1‑150 (MQ=255)
                                                                                                   tAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTTAATTTTTGTCTTCAACCGCCCTGAGAAACTCACTCTTCCCACAGCGCCAGTGCCCACGTTTTTTCATCCTCCGTACCCCCGTTTACCGCCGGGCTGGACaaaaa                                                    >  2:503750/1‑150 (MQ=255)
                                                                                                   tAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAAACGCCATGAGAAAATCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACaaaaa                                                    >  1:634580/1‑150 (MQ=255)
                                                                                                        cGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAAATGTTGTCTTCAACCGCCCTGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGGTGGACAAAAACCGCtt                                              >  1:663603/1‑151 (MQ=255)
                                                                                                           tAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTTTTTTCTTCAACCGCCATTAGAAACTCACTCTTCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCCTACGCTCGTTTACCGCCGGGCTGGACAAAAAATGCttt                                             >  2:214835/1‑149 (MQ=255)
                                                                                                           tAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCttt                                             >  1:542179/1‑149 (MQ=255)
                                                                                                            aaTCAGCCTTTATATTTTTTGTTTAAACGGGCGTTTGATTTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTaa                                           <  1:719035/150‑1 (MQ=255)
                                                                                                             accacccGTTATATTTTTTTTTAATCCGGGGGGTTAAATTTTTCTTTACAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACaaa                                                      <  2:826752/133‑1 (MQ=255)
                                                                                                                 gctgTTATTTTTTTTGTTAAACCTGGCGTTTGATTTTTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACa                                         <  2:147269/144‑1 (MQ=255)
                                                                                                                  ctgTTATTTTTATTTTTAATACGGGTGTTTGACTGTTTGTCTAAACCCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAAtt                                      <  1:501729/147‑1 (MQ=255)
                                                                                                                  ccGTTATATTTATTGTTTAAACTGGCGGTTTGATTGTTGTCTTCAACCGCCCTGAGAAACTCACTCTTCCCACAGCGCCAGTGCCAACGGTTTTTCATCATCCGTACGCCCGTTTACCCCCGGGCTGGACAAAAACTGCTTTAACaaat                                      >  2:799239/1‑147 (MQ=255)
                                                                                                                   cGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAAACGCCATGAGAAACTCACTTTGCCCCCAGCGCCAGTGCCAACGATTTATCATCATCCCTACGCTCCATTAACGCCGGGCTGGACAAAAACTGCTTTAACAATTc                                     >  1:37877/1‑149 (MQ=255)
                                                                                                                      tataTTTTTTTTTAAAACGGGGGTTTGATTTTTTTTCTTCACCCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGt                                 <  1:1098967/150‑1 (MQ=255)
                                                                                                                       ttatttttttttttAAACCGGGGGTTTAATTTTTTTTCTCCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAgg                                  <  2:1131376/146‑1 (MQ=255)
                                                                                                                         tatttattGTTTAAACTGGCGGCTTAAATGTTGTCTTCAAACGCCATGAGAAACTCACTCTTCCCCCAGCGCCCGTGCCAACGGTTTATTATCATCCCTACCCCCGTTTACCCCCGGGCCGGGACAAAAATGCTTTTACAAATTAGGtaa                               >  2:688300/1‑150 (MQ=255)
                                                                                                                           tttatgTTTTAAACGGGGGTTTAAATTTTTGTCTCCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGt                                 <  1:567110/146‑1 (MQ=255)
                                                                                                                            ttttttgttAAACCGGGGGTTTAATTTTTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGtaat                              <  1:979304/142‑1 (MQ=255)
                                                                                                                            ttattGTTTAAACTGGCGGTTTAATTGTTGTCTTTAACCGCCATGAGAAAATCACCCTGCCCACAGCGCCAGTGCCAACGGTTTATCATCATTCGTACGCTCGGTTACCGCCGGGCTGGACAAAAAATGCTTTAACAAATCAGGgaataa                            >  2:1191398/1‑150 (MQ=255)
                                                                                                                               ttGTTTAAACTGGCGGTTTAATTGTTGTCTTCAACCCCCCTGAGGAACTCACTCTGCCCACAGCGCCAGTGGCAACGTTTTATCATCATCCGTACGCCCGTTTACCGCCGGGCTGGACAAAAAATGCTTTAAAAATTTCGGTAATAAATc                         >  2:501629/1‑150 (MQ=255)
                                                                                                                               ttGTTAAACCTGGCGGTTGATTTTTTGTCTCCACCCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAAt                          <  2:875288/149‑1 (MQ=255)
                                                                                                                                  ttaagccgggCTGTTGTTTTTTTGTCTACAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAg                      <  2:336412/142‑1 (MQ=255)
                                                                                                                                     aaccggggtgTTAAATTGTTGTCTACAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAgtt                    <  2:942964/140‑1 (MQ=255)
                                                                                                                                     aaaCTGGCGGTTGATTTTTTGTCTTCACCCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAg                      <  2:951208/147‑1 (MQ=255)
                                                                                                                                          gggggTTGAATTGTTGTCTCCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAgttgtt                 <  1:987052/144‑1 (MQ=255)
                                                                                                                                          gggggTTAAGTTGTGTTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTg             <  2:1133473/148‑1 (MQ=255)
                                                                                                                                             ggTTGATTTGTTTTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGt          <  1:976749/151‑1 (MQ=255)
                                                                                                                                                  ttaatttgtttTCTTAACCCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGtt         <  2:349079/145‑1 (MQ=255)
                                                                                                                                                  gAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTtgct      >  1:344754/1‑151 (MQ=255)
                                                                                                                                                     atttgtttTTTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGCCAAAAACTGCTTTAACAATTCAGGTAATAAATCGATTTGTTCCTGCGTtgct      <  2:1204948/148‑1 (MQ=255)
                                                                                                                                                     aTTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTtgctgc    >  1:453174/1‑151 (MQ=255)
                                                                                                                                                        tgttgtCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTgccg  >  1:1170847/1‑151 (MQ=255)
                                                                                                                                                        tgttgtCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTgccg  >  1:749620/1‑151 (MQ=255)
                                                                                                                                                         gttttCTTATACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCTTtgct      <  2:666277/146‑1 (MQ=255)
                                                                                                                                                         |                                                                                                                                                       
GATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCG  >  NC_000913/2151538‑2151835

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