Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,326,4500TG43.8% 104.2 / 16.2 65G114G (GGT→GGGatoBacetyl‑CoA acetyltransferase
Reads supporting (aligned to +/- strand):  ref base T (7/29);  new base G (28/0);  total (36/29)
Fisher's exact test for biased strand distribution p-value = 7.97e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCGCGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGATTACCGCCGAAAAC  >  NC_000913/2326303‑2326600
                                                                                                                                                   |                                                                                                                                                      
ggCGCGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATg                                                                                                                                                     <  2:1048366/151‑1 (MQ=255)
    cgTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCGCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTAGGGaaa                                                                                                                                                 <  2:664719/151‑1 (MQ=255)
           cACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAg                                                                                                                                          <  1:707361/151‑1 (MQ=255)
            aCTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAAAATGAg                                                                                                                                          >  2:938960/1‑150 (MQ=255)
               gTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGtt                                                                                                                                        <  1:94002/149‑1 (MQ=255)
                    aaaGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGcccc                                                                                                                                 <  1:127546/151‑1 (MQ=255)
                     aaGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGccccc                                                                                                                                <  2:27886/151‑1 (MQ=255)
                     aaGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGccccc                                                                                                                                <  2:292591/151‑1 (MQ=255)
                      aGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCt                                                                                                                               <  2:92050/151‑1 (MQ=255)
                          ggCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCtactt                                                                                                                           <  1:1064571/151‑1 (MQ=255)
                             tGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTGTGGGAAAAATGAGGTTAGCCCCCtactt                                                                                                                           >  1:898457/1‑148 (MQ=255)
                             tGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAATATGAGATTAGCCCCCtactta                                                                                                                          >  2:800961/1‑149 (MQ=255)
                                 aGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAATATGAGGTTAGCCCCCtcccttccc                                                                                                                       >  1:982681/1‑140 (MQ=255)
                                  gAAACGGTGTTCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATg                                                                                                                   <  2:371063/151‑1 (MQ=255)
                                     aCGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAATATGAGTTTAGCCCCCCACTTTCTCGATGCa                                                                                                                 >  2:357551/1‑150 (MQ=255)
                                      cGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCaaa                                                                                                               <  1:940055/151‑1 (MQ=255)
                                       ggTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTTTGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCaaaa                                                                                                              <  2:1044628/151‑1 (MQ=255)
                                       ggTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAATATGAGTTTAGCCCCCTACTTACTCGATGCaaa                                                                                                               >  2:885521/1‑150 (MQ=255)
                                         tgtgCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGTGGGGAAGGGAAATATTAGGTTAGACCCCTACTTACTCAATgacaaa                                                                                                              >  2:798313/1‑144 (MQ=255)
                                             cGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGc                                                                                                        <  1:210284/151‑1 (MQ=255)
                                                aTTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGGAAATATGAGTTTTACCCCCCACTTACTCGAAGCAAAAGCACGctc                                                                                                      >  2:957093/1‑150 (MQ=255)
                                                aTTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGctct                                                                                                     <  2:378879/151‑1 (MQ=255)
                                                aTTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAATATGAGTTTAGCCCCCCACTTACTCGATGCAAAAGCACGccc                                                                                                      >  2:210814/1‑148 (MQ=255)
                                                   cACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAAAATTGGGTTAGGCCCCATCTTACTCGATGGAAAAACCCGc                                                                                                        >  2:888728/1‑145 (MQ=255)
                                                     cGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTAGGGAGAATAGGAGGTTTGCCCCCCACTTTCTCGATTCAAAAAGACGctcttggt                                                                                                 >  1:754516/1‑146 (MQ=255)
                                                           aTAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAATATGAGGTTAGCCCCCCACTTTCTCGATGCCAAAGCAAGCCCTTGTTATCGTCt                                                                                          >  1:168425/1‑151 (MQ=255)
                                                             aaaGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTg                                                                                        <  1:1024611/151‑1 (MQ=255)
                                                             aaaGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAATATTAGTTTAGGCCCCTACTTTATCGCTGCAAAAACACGCTCTGGGTTTcttctt                                                                                         >  1:578480/1‑150 (MQ=255)
                                                             aaaGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAAAATGAGGTTAGCCCCCTACTTACTCGATGCAAAAGCACCCCCCTGTTATCg                                                                                             >  2:921548/1‑146 (MQ=255)
                                                               aGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGGTGGAAAATATGAGTTTAGCCCCCTACTTTCTCGATGCAAAAGCACGCTCTGGTTATCGTCTTg                                                                                        >  1:433724/1‑149 (MQ=255)
                                                                  aTGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAATATGAGTTTTGCCCCCTACTTTCTCGATGCAAAAGAAAGCTCTGGTTATTGTTTTGgag                                                                                     >  2:34654/1‑149 (MQ=255)
                                                                    gtgGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACg                                                                                  <  1:834655/150‑1 (MQ=255)
                                                                     tgGGTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACgg                                                                                 <  1:177288/150‑1 (MQ=255)
                                                                      ggTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGAc                                                                               <  2:267988/151‑1 (MQ=255)
                                                                         tCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAg                                                                             <  1:613326/150‑1 (MQ=255)
                                                                              tCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTAt                                                                       >  2:254365/1‑151 (MQ=255)
                                                                                ttAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAATATGGGTTTTGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGGTATCGTCTTTGAGACGGACAGGTTTAt                                                                       >  1:446543/1‑149 (MQ=255)
                                                                                  aaaaGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAAAATGAGGTTAGCCCCCCACTTTCTCGGTGGAAAAGCAAGCCCTTGGGATCGTCTTGGGGAAGGGCAGGtgt                                                                         >  1:426637/1‑143 (MQ=255)
                                                                                  aaaaGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGGCTTGGAGACGGACAGGTTTATGAc                                                                    >  2:278839/1‑150 (MQ=255)
                                                                                  aaaaGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGAGTGGAAAAAATAAGGTTATCCCCCTACTTTCTTGATGCAAAAGCACGCTCTGGTTATCGTTTTTGAGACGGACAGGATTAAGACg                                                                   >  2:529861/1‑151 (MQ=255)
                                                                                     aGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAAATTAAGGTTAGGCCCCTACTTACCCGATGCAAAAGCACGCTCTGGGTTTCGTCTTGGAGACGGGCAGGgtt                                                                         >  2:834873/1‑142 (MQ=255)
                                                                                          ggCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGGGACGGACAGGGTTATGAAGTAATACt                                                            >  2:641509/1‑150 (MQ=255)
                                                                                                tGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAACACGCTCTGGTTATCGTCTTGGAGACGGACAGGGTTATGACGTAATcc                                                             >  2:22159/1‑143 (MQ=255)
                                                                                                 gccgccCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGAt                                                     <  2:965221/150‑1 (MQ=255)
                                                                                                  ccgccCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATgg                                                   <  2:20876/151‑1 (MQ=255)
                                                                                                  ccgccCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAATATGAAGTTAGCCCCCCACCTACTCGGTGCGAAAAGACGCACTGGGTTTCCGATTGGGGAGGGGAAGGGTTATTACGTAATccccggcgacg                                                    >  1:100681/1‑141 (MQ=255)
                                                                                                     cccAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCt                                                <  1:460409/151‑1 (MQ=255)
                                                                                                     cccAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCt                                                <  2:417195/151‑1 (MQ=255)
                                                                                                       cAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTTCTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGa                                              <  1:943882/151‑1 (MQ=255)
                                                                                                            cATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAAAATGAGGTTAGCCCCCTACTTTCTCGATGCAAAAGCACGCTCTGGGTATTGTTTTGGAGAAGGACAGGGTTATGAAGTAATCCTGCGCGATgggc                                                 >  1:1030310/1‑141 (MQ=255)
                                                                                                              ttCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGGTGGAAAATATGGGGTTAGCCCCCCACTTAATCGATGCAAAAGCACGCTCTGGGTATCGGCTTTGAGACGGGCAGGGTTATGACGTAATCCTGCGCGGTGGGCTGAtgtggg                                        >  1:1071230/1‑148 (MQ=255)
                                                                                                                  ggcaggTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAAAATGAGGTTAGCCCCCTACTTACTCGATGCAAAAACACGCTCTGGGTATCGGCTTGGAGACGGACAGGGTTTTGACGGAATCCTGCGCGATGGCCTGATGTGCGCCAcc                                   >  1:558483/1‑151 (MQ=255)
                                                                                                                     aggTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCAt                                <  1:1115067/151‑1 (MQ=255)
                                                                                                                       gTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGGGACGGACAGGTTTATGACGTAATCCTGCGCGATGGGCTGATGTGGGCCACCCATgg                              >  2:270812/1‑151 (MQ=255)
                                                                                                                             cgcAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCa                        <  2:282372/151‑1 (MQ=255)
                                                                                                                              gcAGAGCATTGTGGCGGGGGGGATGGAAAATATGAGTTTAGCCCCCCACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGGTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGGTATCa                        >  2:448096/1‑150 (MQ=255)
                                                                                                                              gcAGAGCATTGTGGCGGGGGGGATGGAAAAAATGAGGTTAGCCCCCCACTTTCTCGATGCAAAAACACGCTCTGGTTATCGTCTTGGAGACGGACAGGGTTATGACGTAATCCCGCGCGATGGCCTGATGTGCGCCACCCCTGGTTAt                          >  1:14074/1‑148 (MQ=255)
                                                                                                                               cAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCata                      <  1:1105256/151‑1 (MQ=255)
                                                                                                                               cAGAGCATTGTGGCGGGGGGGATGGAAAATATGAGTTTTGACCCCTACTTACTCGAAGCAAAAGCACGCTCTGGGTTTCGGCTTGGGGGGGGACAAGTTTATGACGTAATTCTGCGGGATGGCCCGTTGGGCGCCAACCATGGtttt                          >  2:1068447/1‑145 (MQ=255)
                                                                                                                                  agCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATgg                   <  2:258319/151‑1 (MQ=255)
                                                                                                                                      ttGTGGCGGGGGGGGTGGGAAAAATGAGTTTAGCCCCCTACTTACTCGATGCGAAAGCACGCTCTGGGTTTTGTTTTTGGGAAGGGGAGGGTTTTGACGTAATCCTGCGCGGTGGCCCGTTGTGCGCCCCCCATGGGTATCATATtgggg                >  1:241469/1‑149 (MQ=255)
                                                                                                                                            cGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGATTAccgc         >  1:427182/1‑151 (MQ=255)
                                                                                                                                             ggggggTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGATTAccgcc        >  1:479915/1‑151 (MQ=255)
                                                                                                                                                 ggTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGATTACCGCCGaa     >  1:565/1‑150 (MQ=255)
                                                                                                                                                  gTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGATTACCGCCGaaaa   <  2:465658/151‑1 (MQ=255)
                                                                                                                                                   tATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGATTACCGCCGAAAAc  >  2:252216/1‑151 (MQ=255)
                                                                                                                                                   |                                                                                                                                                      
GGCGCGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGATTACCGCCGAAAAC  >  NC_000913/2326303‑2326600

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