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

CGGCGCGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGATTACCGCCGA  >  NC_000913/2326302‑2326596
                                                                                                                                                    |                                                                                                                                                  
cGGCGCGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTa                                                                                                                                                   >  1:513165/1‑150 (MQ=255)
   cgcgTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGa                                                                                                                                               <  2:546112/151‑1 (MQ=255)
    gcgTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGaa                                                                                                                                              <  3:328974/151‑1 (MQ=255)
    gcgTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGaa                                                                                                                                              <  2:822758/151‑1 (MQ=255)
       tCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGaaaa                                                                                                                                            >  1:1049942/1‑150 (MQ=255)
         aGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAatat                                                                                                                                         <  3:368590/151‑1 (MQ=255)
           gCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAatat                                                                                                                                         >  4:450898/1‑149 (MQ=255)
            cACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGaaaaaa                                                                                                                                          >  3:390302/1‑145 (MQ=255)
            cACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAata                                                                                                                                          >  2:939611/1‑147 (MQ=255)
                gTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGttt                                                                                                                                   <  1:1231924/150‑1 (MQ=255)
                  tAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAg                                                                                                                                 <  2:1279369/150‑1 (MQ=255)
                   aaaaaGCGGGCTGGCGGAACCGGTGTGGGGTTTCACGGTCAATAAAGAAGGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCGGCCCAGGCAATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGTATGGAAAATATGAGTTTAGccc                                                                                                                              <  3:350734/151‑1 (MQ=255)
                      aaGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGccccc                                                                                                                            <  2:1131297/151‑1 (MQ=255)
                        gCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGGTGGGAAATATGAGGTTAGgcccc                                                                                                                            >  3:172114/1‑149 (MQ=255)
                         cGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCtac                                                                                                                         <  3:133885/151‑1 (MQ=255)
                           ggCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCtactt                                                                                                                       <  1:699095/151‑1 (MQ=255)
                            gCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCtactta                                                                                                                      <  2:1267349/151‑1 (MQ=255)
                                  aGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAATATGGGGTTAGCCCCCCACTTACTCGAt                                                                                                                >  1:852625/1‑151 (MQ=255)
                                     aaCGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAATATGAGGTTAGCCCCCCACTTACTCGa                                                                                                                 >  1:1137670/1‑147 (MQ=255)
                                        ggTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCaaaa                                                                                                          <  2:1151120/151‑1 (MQ=255)
                                        ggTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGGTGGGAAATAAGAATTTAGCCCCCTCCTTACTAGTGGCaaa                                                                                                           >  2:1074727/1‑150 (MQ=255)
                                           gtgCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGGTGGGAAATATGGGGTTAACCCCCTTCCTACTCgcagccaaa                                                                                                          >  4:197233/1‑140 (MQ=255)
                                             gCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTATTTACTCGATGCAAAAGCACg                                                                                                     <  2:906535/151‑1 (MQ=255)
                                               ggATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCa                                                                                                       >  4:148424/1‑147 (MQ=255)
                                                  ttCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGGTGGGAAAAATGAGTTTTGCCCCCCACTTACTCGATGCAAAAGAACGctct                                                                                                 >  2:1055347/1‑150 (MQ=255)
                                                    cACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGGTGGAAAATATGAGTTTAGCCCCCTCCTTACTCGATGCAAAAACACGCTCTgg                                                                                               >  4:453682/1‑150 (MQ=255)
                                                        gTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATc                                                                                          <  2:626398/151‑1 (MQ=255)
                                                        gTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAAAATGAGTTTAGCCCCCTACTTACTCGATGCAAAAACACGCTCTGGTTAt                                                                                           >  4:269674/1‑150 (MQ=255)
                                                          cAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGt                                                                                        <  2:1059210/151‑1 (MQ=255)
                                                           aaTAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAATATGAGTTTAGCCCCCCACTTACTCGATGCAAAAGAACGCCCTGGTTATCg                                                                                         >  2:1203095/1‑149 (MQ=255)
                                                            aTAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCt                                                                                      <  1:759685/151‑1 (MQ=255)
                                                            aTAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAATATGAGTTTAGGCCCCTACTTACTCGATGCAAAAGCACGCTCTGGGTATCGt                                                                                        >  1:518763/1‑149 (MQ=255)
                                                            aTAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAATATGAGTTTAGCCCCCCACTTACTCGATGCCAAAACACGCTCTGGTTATCg                                                                                         >  2:975864/1‑148 (MQ=255)
                                                              aaaGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAACACGCGCTGGGTATTGTCtt                                                                                     >  3:254291/1‑150 (MQ=255)
                                                               aaGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTgg                                                                                   <  2:888329/151‑1 (MQ=255)
                                                               aaGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAAAATGAGTTTAGCCCCCCACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTgg                                                                                   >  3:46760/1‑151 (MQ=255)
                                                                aGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGga                                                                                  <  3:68760/151‑1 (MQ=255)
                                                                 gTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGgag                                                                                 <  4:441674/151‑1 (MQ=255)
                                                                  tATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAATATGAGGTTAGCCCCCTCCTTACTCGATGCAAAAGCACGCCCTGGGTATCGGCTTgggg                                                                                 >  3:129251/1‑148 (MQ=255)
                                                                   aTGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGAc                                                                               <  1:1065703/151‑1 (MQ=255)
                                                                    tgtgGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTTCTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGgaga                                                                                <  1:484280/149‑1 (MQ=255)
                                                                    tgtgGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAAAATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGgaga                                                                                >  1:573036/1‑149 (MQ=255)
                                                                    tgtgGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAAAATGAGTTTAGCCCCCCACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCtt                                                                                     >  3:257261/1‑144 (MQ=255)
                                                                           cGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAg                                                                         <  1:715860/149‑1 (MQ=255)
                                                                            gggTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAACACGCCCTGGGTATCGGCTTGGGGACGGACAgg                                                                        >  1:595750/1‑149 (MQ=255)
                                                                             ggTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAATATGAGTTTAGCCCCCCACTTACTCGATGCAAAAACACGCTCTGGTTATCGTCTTTGAGACGGACAGGttt                                                                     >  1:699076/1‑151 (MQ=255)
                                                                               tCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAATATGAGTTTAGCCCCCTACTTTCTCGATGCAAAAGCACGCCCTGGGTATCGTCTTGGAGACGGGCAGGGTTAt                                                                   >  3:212039/1‑151 (MQ=255)
                                                                                cTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAATATGAGTTTAGCCCCCCACTTACTCGGTGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGgtttt                                                                   >  2:149993/1‑148 (MQ=255)
                                                                                  tAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATg                                                                  <  2:491343/149‑1 (MQ=255)
                                                                                  tAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAAAATGAGGTTAGCCCCCCACTTACTCGATGCAAAAGCACGCGCTGGGTATCGTCTTGGAGGCGGaacgggtttt                                                                   >  2:285818/1‑139 (MQ=255)
                                                                                    aaaGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAAAATGAGGTTAGCCCCCTACTTACTCGGTGTAAAAAGACGCCATGGGTTTCGTCTTTGGGGCGGACAGGgtt                                                                     >  3:192651/1‑144 (MQ=255)
                                                                                     aaGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTa                                                             <  1:488532/151‑1 (MQ=255)
                                                                                        tgtgGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAATATGAGGTTAGCCCCCCACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTTGAGACGGACAGGTTTATGACGTAAt                                                           >  1:312504/1‑150 (MQ=255)
                                                                                         gtgGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATcc                                                         <  2:1260988/151‑1 (MQ=255)
                                                                                         gtgGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATcc                                                         <  1:533139/151‑1 (MQ=255)
                                                                                           ggCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTg                                                       <  1:125664/151‑1 (MQ=255)
                                                                                              gcTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGGAAATATGAGTTTAGCCCCCTACTTACTCCATGCAAAAGCACGCTCTGGGTTTCGTCTTGGAGACGGGCAGGGTTATTACGGAATCCTgcg                                                     >  1:1065457/1‑150 (MQ=255)
                                                                                                    cgccCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAAAATGAGGTTTGCCCCCTACTTACTCGGTGCAAAAGCACGCTCTGGTTATCGGCTTTGGGACGGACAGGGTTATGACGTTATCCTGCGCGATgg                                               >  1:667827/1‑150 (MQ=255)
                                                                                                           gCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCCCTGGTTATCGTCTTGGGGACGGACCGGGTTATGACGTAATCCTGCGCGATGGCCTGAtgt                                       >  4:169152/1‑151 (MQ=255)
                                                                                                              aTTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTGTGGAAAAAATGAGTTTAGCCCCCCACTTACTCGATGCAAAAGAACGCCCTGGTTATCGTCTTTGGGACGGACAGGGTTATGACGTAATCCTGCGCGGTGGCCCGATGTgcg                                    >  2:1118351/1‑151 (MQ=255)
                                                                                                                  aggcaggTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAAAATGAGGTTAGCCCCCCTCTTACTCGATGCAAAAACAAGCTCTGGGTATCGTTTTTGAGACGGACAGGGTTATGACGTAATCCTGCGCGATGGCCTGGTGTGCGCCAc                                >  1:595192/1‑151 (MQ=255)
                                                                                                                  aggcaggTCAGGCGCAGAGCATTGTAGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACGCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGTGCCa                                 <  3:169150/150‑1 (MQ=255)
                                                                                                                   ggcaggTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCAcc                               <  1:1118557/151‑1 (MQ=255)
                                                                                                                   ggcaggTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAACACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCAcc                               >  3:315189/1‑151 (MQ=255)
                                                                                                                       ggTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAACACGCTCTGGGTATCGTCTTGGAGACGGACAGGGTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCAt                            >  1:407651/1‑150 (MQ=255)
                                                                                                                       ggTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAATATGAGTTTAGCCCCCCACTTTCTCGATGCAAAAACACGCTCTGGTTATTGTCTTGGGGACGGACAGGTTTATGAAGTAATCCTGCGCGATGGCCTGATGGGCGCCACCCATg                           >  3:403704/1‑151 (MQ=255)
                                                                                                                           aGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTa                       <  3:213705/151‑1 (MQ=255)
                                                                                                                               gcAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCat                   <  2:518670/151‑1 (MQ=255)
                                                                                                                                cAGAGCATTGTGGCGGGGGGGGTGGGAAAAATGAGGTTAGACCCCCACTTTTTCGATGCAAAAAAACGCCCTGGTTATTGTGTTGGGGAAGGACCGGGTTTTGACGCAAAATTGCGCGATGGTCTTGTGTGCGCCCACCATGGGTa                       >  2:847586/1‑146 (MQ=255)
                                                                                                                                        tgtgGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGAt           <  1:149770/150‑1 (MQ=255)
                                                                                                                                         gtgGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGATTa         <  3:455493/151‑1 (MQ=255)
                                                                                                                                               gggggTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGATTAccgccg   >  2:433002/1‑151 (MQ=255)
                                                                                                                                                ggggTATGGAAAATATGAGTTTAGCCCCCTGCTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAAGCCTGCGCGATGGCCTGATGTGCGCCGCCCATGGTTATCATATGGGGATTACCGCCGa  <  3:243584/151‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CGGCGCGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGATTACCGCCGA  >  NC_000913/2326302‑2326596

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