Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,326,4500TG40.9% 122.2 / 13.2 66G114G (GGT→GGGatoBacetyl‑CoA acetyltransferase
Reads supporting (aligned to +/- strand):  ref base T (5/34);  new base G (27/0);  total (32/34)
Fisher's exact test for biased strand distribution p-value = 8.22e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGATTACCGCCGAAAA  >  NC_000913/2326306‑2326599
                                                                                                                                                |                                                                                                                                                      
gcgTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGaa                                                                                                                                                 <  2:206798/151‑1 (MQ=255)
   tCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCCGGCGAAGAGCATTGTGGCGGGGGGTAAGGgaaa                                                                                                                                               >  2:924458/1‑150 (MQ=255)
   tCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGaaa                                                                                                                                                <  2:727465/149‑1 (MQ=255)
    cAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGgaaa                                                                                                                                               >  2:1125368/1‑149 (MQ=255)
     aGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAatat                                                                                                                                            <  1:541422/151‑1 (MQ=255)
     aGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGgaaaaa                                                                                                                                             >  1:776718/1‑148 (MQ=255)
            gTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGttt                                                                                                                                      <  1:803432/150‑1 (MQ=255)
            gTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGttt                                                                                                                                      <  1:1070422/150‑1 (MQ=255)
            gTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTa                                                                                                                                     <  2:992476/151‑1 (MQ=255)
            gTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGGTGGAAAATATGAGttt                                                                                                                                      >  2:1151441/1‑150 (MQ=255)
             ttAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAATATGAGttt                                                                                                                                      >  1:551706/1‑149 (MQ=255)
             ttAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAATATGAGgtt                                                                                                                                      >  1:329223/1‑149 (MQ=255)
               aaaaaGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGc                                                                                                                                   <  1:1244984/150‑1 (MQ=255)
                        gCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCtactta                                                                                                                         <  2:1304579/151‑1 (MQ=255)
                        gCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGGTGGAAAAAATGAGTTTAGCCCCCtac                                                                                                                            >  2:529640/1‑148 (MQ=255)
                             cAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGAAGGGAAAAATGAGGTTAGCCCCCcactt                                                                                                                          >  2:272859/1‑145 (MQ=255)
                                  aCGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGGTGGAAAAAATGGGAGTAGGCCCCCTCTTACTCg                                                                                                                     >  2:190802/1‑145 (MQ=255)
                                   cGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCaaa                                                                                                              <  1:1303032/151‑1 (MQ=255)
                                     gtgtGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAAAATGAGTTTAGCCCCCTACTTACTCGATGCaaaa                                                                                                             >  2:1257826/1‑150 (MQ=255)
                                          cGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGGTGGGAAAAATGGGGTTAACCCCCTCCTTACTCGGTGCaaa                                                                                                              >  1:320074/1‑144 (MQ=255)
                                                 aCGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGAATTGTGGCGGGGGGGATGGAAAAGATGAGGTTAGGCCCCCACCTACTAGATGGACAAGCTCGCGCTGGtt                                                                                                >  2:761187/1‑151 (MQ=255)
                                                     tCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAATATGAGGTTTGCCCCCTACTTATTCCGTGCCAAAGGACGCCCTGGTTATCg                                                                                            >  1:113773/1‑151 (MQ=255)
                                                      cAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGt                                                                                                 <  1:916234/145‑1 (MQ=255)
                                                      cAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGGTGGAAAAAATGAGTTTAGCCCCCTACTTACTTCATGCAAAAAAACGCACTGGGTATCg                                                                                            >  2:804219/1‑150 (MQ=255)
                                                      cAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAATATGAGTTTAGCCCCCCACTTACTCGATGCAAAAGCACGCCCTgggttt                                                                                              >  1:199150/1‑144 (MQ=255)
                                                      cAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAATATGAGGTTAGACCCCCACCTTCTCGATGCAAAAGCACGCTCTGGt                                                                                                 >  2:916670/1‑145 (MQ=255)
                                                       aaTAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAAGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGGAAATATGAGGTTTGCCCCCCCCTTACTCGCTTCAAAAGCAAGCGCTGGTTAttg                                                                                            >  2:328912/1‑147 (MQ=255)
                                                        aTAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCt                                                                                         <  1:370099/151‑1 (MQ=255)
                                                             gTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGgag                                                                                    <  1:308761/151‑1 (MQ=255)
                                                             gTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGgag                                                                                    <  2:76060/151‑1 (MQ=255)
                                                              tATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAAAATTAGGTTAGCCCCCTACTTTCTCGATGCAAAAGCACGCCCTGGTTATCGTCTTg                                                                                       >  2:83227/1‑147 (MQ=255)
                                                                  tgGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAAAATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGGGACgg                                                                                >  1:436654/1‑150 (MQ=255)
                                                                     ttCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAACACGCTCTGGTTATCGTCTTGGAGACgg                                                                                >  1:548009/1‑147 (MQ=255)
                                                                      tCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAg                                                                            <  1:803435/150‑1 (MQ=255)
                                                                          gTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAAGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTCGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTa                                                                       <  2:1281696/151‑1 (MQ=255)
                                                                             ttAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAACACGCTCTGGTTATCGTCTTGGAGACGGACAGGGTTAt                                                                      >  1:85681/1‑149 (MQ=255)
                                                                               aaaaGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAAAAATGGTTTAGCCCCCCACGGTATAGGTGCAAAAACACGCACTGGTTGTCGTTTTTGGGACGGACAGGGTTATGACg                                                                  >  2:702113/1‑151 (MQ=255)
                                                                                aaaGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAAAATGAGGTTTACCCCCCACTTTCTCGCTGCAAAAACACGCCCTGGGTATTGTCTTTGAGACGGACAGGGTTTTTACg                                                                  >  2:529054/1‑150 (MQ=255)
                                                                                 aaGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTa                                                                <  1:1217180/151‑1 (MQ=255)
                                                                                   gtgtGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAAt                                                              <  1:495406/151‑1 (MQ=255)
                                                                                    tgtgGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAACACGCCCTGGGTATCGTCTTGGGGACGGACAGGTTTATGACGTaa                                                               >  1:283052/1‑149 (MQ=255)
                                                                                       ggCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCt                                                           <  2:1308781/150‑1 (MQ=255)
                                                                                          gcTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTgcgc                                                       <  2:777081/151‑1 (MQ=255)
                                                                                              gccgccCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGAt                                                    <  1:1133240/150‑1 (MQ=255)
                                                                                              gccgccCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGAt                                                    <  1:997812/150‑1 (MQ=255)
                                                                                                    cAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTGTGGAAAATATGAGTTTCGCCCCCCACATACTTGCTGCAAAAGCACGCTCTAGTTATTGTCTTGGAGACGGACAGGTTTATGACGTAATCCTTCGCGCTgg                                                  >  2:1152037/1‑146 (MQ=255)
                                                                                                    cAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAATATGAGTTTAGCCCCCTACTTACTCGCTGCAAAAACACGCCCCTGGTATAGTCTTGGAGACGGGAAGGATTATGACGGAATCCTGCGCGATGGcc                                                >  2:566819/1‑148 (MQ=255)
                                                                                                     aGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGAt                                            <  1:1272710/151‑1 (MQ=255)
                                                                                                     aGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTAAGGAAAATATGAGTTTATCCCCCTACTTACTCGAAGCAAAAGCACGCTCTGGGTATCGTCTTGGAGACGGAAAGGGTTTTGACGGAATCATGCGCGATGGGCTGAt                                            >  1:171642/1‑151 (MQ=255)
                                                                                                        ccATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGAtgtg                                         <  2:465553/151‑1 (MQ=255)
                                                                                                          aTTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGTTGGGAAATAAGAGGTTTGCCCCCTCCTTACTTCATGTAAAACCACGCTCTGGGTATCGGCTTGGAGACGGGAAGGTTTTTGACGTAAACCTGCGTGATTGCCTgttctg                                         >  2:547011/1‑144 (MQ=255)
                                                                                                              aggcaggTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCAc                                   <  2:783859/151‑1 (MQ=255)
                                                                                                                  aggTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCAt                               <  1:62965/151‑1 (MQ=255)
                                                                                                                  aggTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAATATGAGTTTTGCCCCCTACTTACTCGGTGCAAAAAGACGCGCTGGTTATCGTCTTGGGGACGGACAGGGTTATTAAGAAATCCTGCGCGATGGGCTGATGTGCGCCAccc                                 >  2:1272723/1‑149 (MQ=255)
                                                                                                                    gTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAATATGAGGTTAGCCCCCCACTTACTCGATGCAAAAAGACGCCCTGGTTATCGTCTTGGAGACGGACAGGTTTATTACGTAATCCTGCGCGGTGGCCTGCTGTGCGCCACCCATgg                             >  1:1279332/1‑151 (MQ=255)
                                                                                                                      cAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGtt                           <  1:1257125/151‑1 (MQ=255)
                                                                                                                      cAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGt                            <  1:252706/150‑1 (MQ=255)
                                                                                                                             agagCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCatat                    <  2:853539/151‑1 (MQ=255)
                                                                                                                              gagCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATg                   <  2:348003/151‑1 (MQ=255)
                                                                                                                                  aTTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGa               <  2:1146364/151‑1 (MQ=255)
                                                                                                                                  aTTGTGGCGGGGGGGGTGGGAAAAATTAGTTTAGCCCCCCACTTAATCGATGCAAAAGCACGCTCTGGGTATCGTCTTGGAGACGGGCAGGGTTTTGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGGTATCATATggggg               >  2:299250/1‑150 (MQ=255)
                                                                                                                                    tgtgGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGAt              <  1:1150787/150‑1 (MQ=255)
                                                                                                                                      tgGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGATTAc           >  1:1201532/1‑151 (MQ=255)
                                                                                                                                          ggggggTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGATTAccgcc       <  1:924010/151‑1 (MQ=255)
                                                                                                                                              ggTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGATTACCGCCGaaa   <  1:1124737/151‑1 (MQ=255)
                                                                                                                                               gTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGATTACCGCCGaaaa  <  1:538438/151‑1 (MQ=255)
                                                                                                                                                tATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGCTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGATTACCGCCGaaaa  >  1:358680/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                      
GCGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGATTACCGCCGAAAA  >  NC_000913/2326306‑2326599

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