Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913811,7210TG38.1% 145.0 / 12.7 87V67G (GTG→GGG) bioCmalonyl‑ACP O‑methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base T (12/40);  new base G (32/0);  total (47/40)
Fisher's exact test for biased strand distribution p-value = 2.13e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AAATGCAGGATATCGACCGTCTGCTGGAGGTGCTGCATGGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCGGTAATTTATCCACGGCACTCCGCGAGCTGTATCGGGTGGTGCGCCCCAAAGGCGTGGTCGCGTTTACCACGCTG  >  NC_000913/811486‑811944
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                   
aaaTGCAGGATATCGCCCTTCTGCTGGAGGTGTTGCAGGGCAACGGTAAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGAc                                                                                                                                                                                                                                 <  2:769242/239‑1 (MQ=255)
 aaTGCAGGATATCGACCGTCTGCTGGAGGTGCTGCATGGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGAc                                                                                                                                                                                                                                 <  1:1989062/238‑1 (MQ=255)
        gATATCGACCGTCTGCTGGGGGTGCTGCATGGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACCCCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGAc                                                                                                                                                                                                                                 <  1:1012610/231‑1 (MQ=255)
         atatCGACCGTCTGCTGGAGGTGCTGCATGGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGGCACGCCCCGGGGACGGCCCTAGa                                                                                                                                                                                                                        >  1:665643/1‑239 (MQ=255)
             cGACCGTCTGCTGGAGGTGCTGCATGGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGCGCCGCGCACGGGGCGGGCTTTGagt                                                                                                                                                                                                                      >  2:1542317/1‑235 (MQ=255)
                ccGTCTGCTGGAGGTGCTGCATGGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGGCACGCGCAGGTGACGGCCTTAGa                                                                                                                                                                                                                        >  1:803293/1‑232 (MQ=255)
                  gTCTGCTGGAGGTGCTGCATGGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGGCACGGGCCGGGGAGGGGCCTTAggcctggtct                                                                                                                                                                                                                >  2:1012610/1‑230 (MQ=255)
                   tCTGCTGGAGGTGCTGCATGGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGCGACGGCCCTAGATCTCTcgcc                                                                                                                                                                                                               >  1:1013664/1‑238 (MQ=255)
                      gctgGAGGTGCTGCATGGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGGACGTCACGCGCGGGGGACGGGCCTAGATCTCTCccctccc                                                                                                                                                                                                           >  2:729128/1‑238 (MQ=255)
                          gAGGTGCTGCATGGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCCGGGGACGGCCCTAGCTCCCTCTCCCCCCAt                                                                                                                                                                                                         >  1:160122/1‑237 (MQ=255)
                             gTGCTGCATGGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAAAGGCACGCGCAGGGGACGGGCCTAGATCTCCCgccgcccaag                                                                                                                                                                                                        >  1:2329420/1‑231 (MQ=255)
                               gctgcATGGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCCCGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAAAGGTTGt                                                                                                                                                                                                   >  1:1564414/1‑238 (MQ=255)
                                    tggacaggggTAAAAAAAAAAGCCGTTGCGGGGGCGTTGGTGCGGGCGGCCGCCCACTATGGGAAAAATGCAAAGCAAAAGCGCCAGAGTGTTGGCGCCTTACTGGAATTGCTTCCACGGCGTAAACACCCCCACGTACGGGCCGCGGGTTGTGAACAGGGGTGGAGGCGCCGCCACGGGCGGGAAGGTCACCCGCAGGTGACGGGCTTAGATCTCTCGCCGCCAATGCTTGTTCAgg                                                                                                                                                                                              <  2:1076270/230‑1 (MQ=255)
                                     tgccAACGTTTAAAAACCAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCa                                                                                                                                                                                            <  1:1542317/236‑1 (MQ=255)
                                       gaaaggtttAATAAACAAGCCATGGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCCGATCTAAAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTCCTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAgg                                                                                                                                                                                              <  2:780006/228‑1 (MQ=255)
                                           cGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCCATGCTTGTTCAGGCACGcccg                                                                                                                                                                                      >  1:1525456/1‑237 (MQ=255)
                                             gTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGGGACGGCCCTAGATCTCTCGCCGCCCAAGCTTGTTCAGGCACGCCAGaa                                                                                                                                                                                    >  1:1706170/1‑239 (MQ=255)
                                                  taaACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGGGACGGCCTTAGATCTCTCGCCGCCCATGCTTGTTCAGGCACGCCCGAAGGAt                                                                                                                                                                                >  1:2246790/1‑238 (MQ=255)
                                                  taaACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGGCACGCGCCGGGGACGGGCTTAGATCTCTCGCCCCCCAAGCCTGTTCAGGCACGCCCGAAGGATg                                                                                                                                                                               >  2:568930/1‑239 (MQ=255)
                                                       aggcaaTTGGAGCGGCATTTGGCCGGGCAGCCGCACCCTAGGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCa                                                                                                                                                                          <  2:2246790/234‑1 (MQ=255)
                                                           cATGGCAGCGGCATTTGGTGGGGAAGCCGCCCACTATGAGCAACATGCAGATCTACAGCGCCAGAGGGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACCCCCACGTACTGGACGGGGGTTGGGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCa                                                                                                                                                                          <  2:1564414/235‑1 (MQ=255)
                                                           aaTTGCAGCGGCTTTTGGGCGGGCAGCCGCACCCTATGAGCAACATGTAGTTCTACGGCGCTAGAGTGCTGCCGCCTTACTGGAAATGCTCCCACAGCGTAAATACCCCCACGTACGGGACGCGGGTTGGGGCCCGGGCTGGGTGGGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGa                                                                                                                                                                        <  2:1013664/236‑1 (MQ=255)
                                                              tGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGGTCAGGCACGCCAGAAGGATGCCGCAGAc                                                                                                                                                                       >  1:246358/1‑235 (MQ=255)
                                                              tGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGGGACGGCCTTAGATCTCTCGCCCCCAATGCTTGTTCAGGCACCCCCGGACGAAGGCGCCGAACAtt                                                                                                                                                                   >  1:1717618/1‑239 (MQ=255)
                                                              tGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGGGACGGCCCTAGATCTCTCGCCGCCCATGCTTGTTCAGGGCCGCCAGAAGGATGcc                                                                                                                                                                             >  1:1595970/1‑229 (MQ=255)
                                                                    ggCATTTGGCCGGGAGGCCGCACACTAGGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCTTTACTGGCAATGCTTCCACAGCGTAAATACCCCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTAt                                                                                                                                                                 <  2:246358/235‑1 (MQ=255)
                                                                     gCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAAAGCCACGCCCCGGGGACGGCCCTAGATCTCCCCCCGCCCATGCTTGTTCAGGGCCGCCACAAGGATGCCGCCGACCATTAttcggc                                                                                                                                                            >  2:739606/1‑234 (MQ=255)
                                                                                      cGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAgg                                                                                                                                                                                  <  1:2321156/200‑1 (MQ=255)
                                                                                           acTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGa                                                                                                                                                                                                                                   <  1:1781230/147‑1 (MQ=255)
                                                                                                 aGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGAc                                                                                                                                                                                                                                 <  2:1803044/142‑1 (MQ=255)
                                                                                                 aGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGAc                                                                                                                                                                                                                                 >  1:1803044/1‑142 (MQ=255)
                                                                                                     aCATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAAAGGCACGCGCAGGGGACGGCCCTAGa                                                                                                                                                                                                                        >  2:1781230/1‑147 (MQ=255)
                                                                                                       aTGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGGACGTCACGCGCAGGTGACGGCCCTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACCCCCGaa                                                                                                                                                                                    >  1:921917/1‑181 (MQ=255)
                                                                                                       aTGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGa                                                                                                                                                                                                                                   >  1:2366120/1‑135 (MQ=255)
                                                                                                       aTGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGa                                                                                                                                                                                                                                   <  2:2366120/135‑1 (MQ=255)
                                                                                                           aGATCTACACCGCCAGAGTGCTGCCGCCTTACTGGCAATGCTCCCACAGCGTAAATACCCCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAgg                                                                                                                                                                                              <  2:1877808/167‑1 (MQ=255)
                                                                                                                tACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAgg                                                                                                                                                                                              <  2:897992/162‑1 (MQ=255)
                                                                                                                tACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAgg                                                                                                                                                                                              >  1:897992/1‑162 (MQ=255)
                                                                                                                    gcgcCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGGCACGCGCGGGGGACGGCCCTAGATCTCCCGCCCCCCATGCCTGTTCCGGCACGCCAGAAGGAAGCCGCAGACCACTATCTGGCGGGGGAGATCGGAACCCCGCCCGTCGCGGCCGCGGCGTTCGATCtt                                                                                                             >  2:752452/1‑238 (MQ=255)
                                                                                                                     cgcCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGGACGTCACGCGCAGGGGGCGGCCCTAGATCTCTCGCCGCCCATGCTTGTTCAGGCACGCCCGGAGGGAGCCGCAGACCATTATCTGGCGGGAGGTATCGAATTCCCGCCGGTTGCGACTGGGAGGCTCgcgcgtg                                                                                                            >  2:814886/1‑232 (MQ=255)
                                                                                                                     cgcCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGGGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACCCCAGaa                                                                                                                                                                                    >  1:1877808/1‑167 (MQ=255)
                                                                                                                             tGCTGACGCCTTGCGGGCAATGCTTCCACGGCGTAAATCCCCCCACGTACTGGAGGCGGGTTGTGGACCTGGCTGGATGGGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAgg                                                                                                                                                                                              <  1:2317662/149‑1 (MQ=255)
                                                                                                                             tGCTGACGCCTTACTGGCAATGCTCCCACACCGTAAATACACCCACGTCCTGGACGCGGGTTGTGGCCCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTcg                                                                                                                                                                                                                 <  2:1791994/130‑1 (MQ=255)
                                                                                                                                       ttACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGGAAGTCCCGCGCAGGGGGCGGCCTTTGAATTCTCGCCGCCAATgg                                                                                                                                                                                                       >  1:1791994/1‑129 (MQ=255)
                                                                                                                                             gCAATGCTTCCACAGCGTAAATACCCCCACGTACGGGACGCGGGTTGTGGCCCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCCCGCGCAGGTGACGGCCTTAGAtctct                                                                                                                                                                                                                   <  1:1397199/112‑1 (MQ=255)
                                                                                                                                                aTGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGa                                                                                                                                                                                                                                   <  1:1969112/94‑1 (MQ=255)
                                                                                                                                                aTGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGGGa                                                                                                                                                                                                                                   >  2:1969112/1‑94 (MQ=255)
                                                                                                                                                       cacgGCATAAATCCCCCCACGTCCTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCGGTAAt                                                                        <  1:2092897/239‑1 (MQ=255)
                                                                                                                                                       cacaGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGGCCCCCGCAGGAGACGGGCTTAGATCTCCCGCCGCCaaa                                                                                                                                                                                                         >  2:1397199/1‑111 (MQ=255)
                                                                                                                                                            cGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTa                                                                                                                                                                  <  1:384525/146‑1 (MQ=255)
                                                                                                                                                             gTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCCCCACTGGCGGGGACGGTACCGGGTGGCGACGGCCTTAGATCTCCCGCCGCCAATg                                                                                                                                                                                                        >  1:167335/1‑107 (MQ=255)
                                                                                                                                                              tAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGGAGGGGGAGGCCCTAGATCTCTCCCCGCCCATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCAGTAGCTGGCGGGAGATATTGAATCCCCGCCGTTAGCGACCGCGACGTTCGATCTTGCATGGGGCGATCCTGCAATGCCGCGGTGCGGGAATTTAtac                                                                  >  2:1144300/1‑236 (MQ=255)
                                                                                                                                                                      cccACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGGAACTCACGCGCCGGGGACGGCCTTAGATCTCTCGCCGCCCAAGCTTGTTCAGGCACGCCCGAAGGATGGCCCAGACCATTATCTCGCGGGGGATATAGCAACCCCGCCCCTTGGGCACTCGGACTTCTATCTTTCCTGGAGGCATCTCTCAGTGCAGTGGGGGGGGAATTTATCTCCGGCACt                                                          >  2:1624898/1‑238 (MQ=255)
                                                                                                                                                                      cccACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGGGACGGCCTTAGATCTCTCCCCCCCCAAGCTTGTTCAGGCCCGCCCGAAGGATGCCGCCGACCATTATCTGGCGGGGGATATCGAATCCCTGCCGTTAGCGAATGCGACGGTCGATCTTGCATGGAGCAATCTCGCCGTGGCGGGGGGGGGGAATTTATTCACCGCACt                                                          >  2:2044895/1‑238 (MQ=255)
                                                                                                                                                                      cccACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGGCACGCGCAGGGGACGGGCTTAGGTCTCTCGCCGCCAATGCTTGTTCCGGCCCGCCAGAAGGATGGCGCCGACCCCTATCTGGGGGGGGGTATCGGATCCCCGCCCGTTGCGACTGCGACGTTCGGACTTGGATGGGGCAATCTCGCAGGGGAGGGGGGGGGGCAATTATCCACgggccc                                                          >  2:904097/1‑234 (MQ=255)
                                                                                                                                                                      cccACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGGCACCCGCAGGGGGCGGCCCTTGATCTCTCGCCGCCCATGCTTGTTTAGGCCCGCCCGGAGGATGCCGCAGACCATTTTCTTGCcggg                                                                                                                                                        >  2:384525/1‑145 (MQ=255)
                                                                                                                                                                      cccACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAAAGGCACGCGCAGGGGACGGCCTTAGATCTCTCGCCGCCAAAGCCTGTTTAGGCACGCCCGAAGGATGCCGCAGACCCTTATCTGGGGGGGGATATCGAAACCCCGCCCTTAGCGGCTGCGCCGTTTGGTCTTGGATTGGGGCAACTCCCCGTGGCGGGGGGCGGGGATTTATccccagccac                                                          >  2:1186114/1‑230 (MQ=255)
                                                                                                                                                                          cGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGGGAACGCCTTAGAACTCTCGCCGCCCATCCTTGTTCAGGCACGCCCGAAGGAAGCCGCAGACCATTATCTGGCGGGGGATATCGCAACCCTGCCCTTTGCGACAGCGACGTTCGATCTTGCATTGGGCAAACCCACAGTGCAGGGGGGCGGGGATTTATTCACGGCACTCcgcc                                                     >  2:570150/1‑238 (MQ=255)
                                                                                                                                                                             tgtggggggggggTGTGGGCCTGGGTGGTAGGGCCCGCCACGGGCGGGAACGTCCCGCGCAGGTGACGGCCTTAGACCCCCCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGCCCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCGGTAATTTATCCACGGCACTCcg                                                       <  2:412802/226‑1 (MQ=255)
                                                                                                                                                                             gcTGGACGCGGTTTGTGGCCCGGGCTGGATGAGCCGCCACTGGCGGGACCGTCACGCGCAGGTGAGGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCCCGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCGGTAATTTATCCACGGCACTCCGCGAg                                                   <  1:904097/237‑1 (MQ=255)
                                                                                                                                                                             agcggggggggggTGGGGCCCTGGGTGGTATGAGCCGCCGTGGGCGGGAGCGTCCCGCGCGGGTGCGGGCTTAAGTTCCTTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGCCCATTATCTGGCGGGAGATGTCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCGGTAATTTATCCACGGCACTCCGCGAg                                                   <  2:1742823/230‑1 (MQ=255)
                                                                                                                                                                             acgggACGGGGGTTGGGGACCGGGCTGGATGGGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCCCGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCGGTAATTTATCCACGGCACTCCGCGAg                                                   <  1:752452/235‑1 (MQ=255)
                                                                                                                                                                              cggggcgcgGGTTGTGGCCCTGGGGGGATGGGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCTTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCGGTAATTTATCCACGGCACTCCGCGAGCt                                                 <  1:729128/234‑1 (MQ=255)
                                                                                                                                                                                 gccgcgGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGa                                                                                                                                                                                                                                   <  1:2363762/59‑1 (MQ=255)
                                                                                                                                                                                 gACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCCGGGGACGGCCCTAGATCTCTCCCCCCCCATGCTTGTTTAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGAGATCGAATCCCTCCCGTTAGCGGCCGCGACGTTCGATCTTGCATGGGGCAATCTCGCCGTGCAGTGGTGCGGGAATTTATCCACCGCCCCCCGCGa                                                    >  2:277697/1‑233 (MQ=255)
                                                                                                                                                                                 gACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGGCCCGCGCCGGGGACCGCCCTTGATCTCTCGCCGCCCATGCTTGTTTAGGCACGCCCCAAGGAGGCCCCAGACCCTTATTTGGCGGGGGAGATCAAAACCCCCCCGCTAGCGCCTGCGCCGTTTGACCTTTCATTGGGCAATCTCGCGGTGCAGTGGGGCGGGGATTTATCCACCGCACTCCCCGAGCCGt                                               >  2:2070945/1‑238 (MQ=255)
                                                                                                                                                                                 gACGCGGGTTGGGGCCCTGGCTGGAGGGGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCGGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCGGTAATTTATCCACGGCACTCCGCGAGCTGt                                               <  1:653371/238‑1 (MQ=255)
                                                                                                                                                                                      gggTTGTGGACCTGGCTGGATGGGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGa                                                                                                                                                                                 <  1:900128/105‑1 (MQ=255)
                                                                                                                                                                                           gtgGACCTGGCTGGATGAGCCGCCACTGGCGGGAACCTCCCCCGCAGGTGACGGCCCTaaa                                                                                                                                                                                                                        >  2:2363762/1‑59 (MQ=255)
                                                                                                                                                                                           ggggcccgggctggAGGGCCGGCCCCGGGCGGGACCGGCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCGGTAATTTATCCACGGCACTCCGCGAGCTGTATCGGgtggtg                                    <  1:1385080/231‑1 (MQ=255)
                                                                                                                                                                                           ggggACCGGGTGGGATGGGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCGGTAATTTATCCACGGCACTCCGCGAGCTGTATCGGgtggtg                                    <  1:2272655/237‑1 (MQ=255)
                                                                                                                                                                                            tgGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTgg                                                                                                                                                             <  1:176764/119‑1 (MQ=255)
                                                                                                                                                                                              gcccgggctggATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCGGTAATTTATCCACGGCACTCCGCGAGCTGTATCGGGTGGTgcg                                  <  1:2104895/233‑1 (MQ=255)
                                                                                                                                                                                                ccTGGCTGGATGAGCCGCCACTGGCGGGAACATCACGCGCAGGTGGCGGCCTTTGATCTCTCGCCGCCAATGCCTGTTCACGCCCGCCAGAAGGAGGCCGCAGAc                                                                                                                                                                       >  2:900128/1‑105 (MQ=255)
                                                                                                                                                                                                  tggctggATGAGCCGCCACTGGCGGGGACGGCACGCGCAGGTGGCGGCCTTTGATCTCTCGCCGCCAATGCTTGTTTAGGGACGCCCGAAGGATGCCGCCGACCCTTATCTGGCGGGGGGTATTGAATCCCTGCCGTTAGCGGACGCGACGGTTGATCTTGCATGGAGGAATCTCgcagtg                                                                                       >  2:1136546/1‑181 (MQ=255)
                                                                                                                                                                                                    gcgggATGGGCCGCCACTGGCGGGAACGTCCCGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGAGGCCGCAGACCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCaa                                                                                                 <  1:1527265/166‑1 (MQ=255)
                                                                                                                                                                                                     ctggATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCa                                                                                                                                                                          <  1:170790/97‑1 (MQ=255)
                                                                                                                                                                                                      tggATGAGCCGCCACTGGCGGGGACCTCCCCCGCAGGGGGACGCCCTTGATCCCTCCCCGCCCACACTTGTTTCGGCACGCCCGAAGGATTCCGCCGACCACTATCTTGCGGGGGatat                                                                                                                                                   >  2:176764/1‑119 (MQ=255)
                                                                                                                                                                                                         gaggggccgccACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAAtc                                                                                               <  1:1993422/162‑1 (MQ=255)
                                                                                                                                                                                                          aTGAGCCGCCACTGGCGGGAACGGCACCCGCAGGGGACGGCCCTTGATCTCTCCCCGCCCATGCTTGTTCAGGCCCGCCAGAAGGATGCCGCAGACCATTATTTGGCGGGGGATATCGAATCCCCGCCGGTAGCGACCGCGACCGTTGGTCTTTCATGGAGCGATCTCCCCGGGCAGGGGGGCGGGAATTTATTCACGGCACCCCGCGCGCTGTATCGGGGGGGGCGCCCCACAggg                       >  2:2104895/1‑236 (MQ=255)
                                                                                                                                                                                                              gccgccACTGGCGGGAACGTCACGCGCAGGGGACGGCCTTAGATCTCTCGCCCCCCATGCTTGTTCAGGCACGCCCGAAGGATGCCGCCGACCCTTATCTGGCGGGGGATATCCAATTCCTGCCGGTTGCGACTGCGACGGTCGATCTTTCATGGAGGCATCTCGCAGGGCCGTGGGGCGGGAATTTTTCCACGGGCCCCCCCCAGCTTGATCGGGGGGGGGGCCCCCAAAGGGGTggc                  >  2:1399266/1‑238 (MQ=255)
                                                                                                                                                                                                               ccgccACTGGCGGGAACGTCACGCGCCGGGGACGGCCTTAGATCTCTCGCCGCCGACGCTTGTTCAGGCACGCCCGGAGGTTGCCCCCGACCCTTATTTGGCGGGAGATATAGAATTCCCGCCCGTGGCGGCTGCGACGTTCGATCTTTCACGGGGCAATCTCgcaggg                                                                                       >  2:1527265/1‑167 (MQ=255)
                                                                                                                                                                                                                   cACTGGCGGGAACCTCACGCGCAGGGGACGGCCTTAGATCTCTCGCCGCCCAAGCTTGTTTAGGCACGCCCGAAGGATGCCGCCGACCCTTATCCGGCGGGGGATATCGAATCCCTGCCGTTTGCGACCGCGACGTTCGATCTTGCATGGAGGAAACCCgcagtgcag                                                                                    >  2:1993422/1‑168 (MQ=255)
                                                                                                                                                                                                                      gggcggcAGCGGCGCGCGAAGGTGCCGGCCTTAGCTCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCGGTAATTTATCCACGGCACTCCGCGAGCTGTATCGGGTGGTGCGCCCCAAAGGCGTGGTCGCGTTTAc         <  2:651513/238‑1 (MQ=255)
                                                                                                                                                                                                                              aCGTCACGCGCAGGTGCCGGCTTTAGATCTCCCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCGGTAATTTATCCACGGCACTCCGCGAGCTGTATCGGGTGGTGCGCCCCAAAGGCGTGGTCGCGTTTACCACGCTg  <  1:814886/238‑1 (MQ=255)
                                                                                                                                                                                                                                gTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCGGTAATTTATCCACGGCACTCCGCGAGCTGTa                                              >  1:1228772/1‑193 (MQ=255)
                                                                                                                                                                                                                                         aGGTGACGGCTTTAGATCTCTCGCCGCCAATGCTTGTTCAGGc                                                                                                                                                                                             <  1:237461/43‑1 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                   
AAATGCAGGATATCGACCGTCTGCTGGAGGTGCTGCATGGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCGGTAATTTATCCACGGCACTCCGCGAGCTGTATCGGGTGGTGCGCCCCAAAGGCGTGGTCGCGTTTACCACGCTG  >  NC_000913/811486‑811944

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