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.

AAATGCAGGATATCGACCGTCTGCTGGAGGTGCTGCATGGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGG‑ATGAGCCGCCACTGGCGGGAACGTCACGCGCA‑GGTGA‑CGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCG‑AATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCA‑TGGAGCAATCTCGCAGTGCAGTGGTGCGGTAATTTATCCACGGCACTCCGCGAGCTGTATCGGGTGGTGCG‑CCCCAAAGGCGTGGTCGCGTTTACCACGCTG  >  NC_000913/811486‑811944
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                   
aaaTGCAGGATATCGCCCTTCTGCTGGAGGTGTTGCAGGGCAACGGTAAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGG‑ATGAGCCGCCACTGGCGGGAACGTCACGCGCA‑GGTGA‑c                                                                                                                                                                                                                                 <  2:769242/239‑1 (MQ=255)
 aaTGCAGGATATCGACCGTCTGCTGGAGGTGCTGCATGGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGG‑ATGAGCCGCCACTGGCGGGAACGTCACGCGCA‑GGTGA‑c                                                                                                                                                                                                                                 <  1:1989062/238‑1 (MQ=255)
        gATATCGACCGTCTGCTGGGGGTGCTGCATGGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACCCCCACGTACTGGACGCGGGTTGTGGACCTGGCTGG‑ATGAGCCGCCACTGGCGGGAACGTCACGCGCA‑GGTGA‑c                                                                                                                                                                                                                                 <  1:1012610/231‑1 (MQ=255)
         atatCGACCGTCTGCTGGAGGTGCTGCATGGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGG‑ATGAGCCGCCACTGGCGGGAACGGCACGCCCC‑GGGGA‑CGGCCCTAGa                                                                                                                                                                                                                        >  1:665643/1‑239 (MQ=255)
             cGACCGTCTGCTGGAGGTGCTGCATGGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGG‑ATGAGCCGCCACTGGCGGGAACGCGCCGCGCA‑CGGGG‑CGGGCTTTGagt                                                                                                                                                                                                                      >  2:1542317/1‑235 (MQ=255)
                ccGTCTGCTGGAGGTGCTGCATGGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGG‑ATGAGCCGCCACTGGCGGGAACGGCACGCGCA‑GGTGA‑CGGCCTTAGa                                                                                                                                                                                                                        >  1:803293/1‑232 (MQ=255)
                  gTCTGCTGGAGGTGCTGCATGGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGG‑ATGAGCCGCCACTGGCGGGAACGGCACGGGCC‑GGGGAGGGGCCTTAggcctggtct                                                                                                                                                                                                                >  2:1012610/1‑230 (MQ=255)
                   tCTGCTGGAGGTGCTGCATGGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGG‑ATGAGCCGCCACTGGCGGGAACGTCACGCGCA‑GGCGA‑CGGCCCTAGATCTCTcgcc                                                                                                                                                                                                               >  1:1013664/1‑238 (MQ=255)
                      gctgGAGGTGCTGCATGGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGG‑ATGAGCCGCCACTGGCGGGGACGTCACGCGCG‑GGGGA‑CGGGCCTAGATCTCTCccctccc                                                                                                                                                                                                           >  2:729128/1‑238 (MQ=255)
                          gAGGTGCTGCATGGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGG‑ATGAGCCGCCACTGGCGGGAACGTCACGCGCC‑GGGGA‑CGGCCCTAGCTCCCTCTCCCCCCAt                                                                                                                                                                                                         >  1:160122/1‑237 (MQ=255)
                             gTGCTGCATGGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGG‑ATGAGCCGCCACTGGCGGGAAAGGCACGCGCA‑GGGGA‑CGGGCCTAGATCTCCCgccgcccaag                                                                                                                                                                                                        >  1:2329420/1‑231 (MQ=255)
                               gctgcATGGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGG‑ATGAGCCGCCACTGGCGGGAACGTCCCGCGCA‑GGTGA‑CGGCCTTAGATCTCTCGCCGCCAAAGGTTGt                                                                                                                                                                                                   >  1:1564414/1‑238 (MQ=255)
                                    tggacaggggTAAAAAAAAAAGCCGTTGCGGGGGCGTTGGTGCGGGCGGCCGCCCACTATGGGAAAAATGCAAAGCAAAAGCGCCAGAGTGTTGGCGCCTTACTGGAATTGCTTCCACGGCGTAAACACCCCCACGTACGGGCCGCGGGTTGTGAACAGGGGTGG‑AGGCGCCGCCACGGGCGGGAAGGTCACCCGCA‑GGTGA‑CGGGCTTAGATCTCTCGCCGCCAATGCTTGTTCAgg                                                                                                                                                                                              <  2:1076270/230‑1 (MQ=255)
                                     tgccAACGTTTAAAAACCAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGG‑ATGAGCCGCCACTGGCGGGAACGTCACGCGCA‑GGTGA‑CGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCa                                                                                                                                                                                            <  1:1542317/236‑1 (MQ=255)
                                       gaaaggtttAATAAACAAGCCATGGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCCGATCTAAAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTCCTGGACGCGGGTTGTGGACCTGGCTGG‑ATGAGCCGCCACTGGCGGGAACGTCACGCGCA‑GGTGA‑CGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAgg                                                                                                                                                                                              <  2:780006/228‑1 (MQ=255)
                                           cGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGG‑ATGAGCCGCCACTGGCGGGAACGTCACGCGCA‑GGTGA‑CGGCCTTAGATCTCTCGCCGCCCATGCTTGTTCAGGCACGcccg                                                                                                                                                                                      >  1:1525456/1‑237 (MQ=255)
                                             gTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGG‑ATGAGCCGCCACTGGCGGGAACGTCACGCGCA‑GGGGA‑CGGCCCTAGATCTCTCGCCGCCCAAGCTTGTTCAGGCACGCCAGaa                                                                                                                                                                                    >  1:1706170/1‑239 (MQ=255)
                                                  taaACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGG‑ATGAGCCGCCACTGGCGGGAACGTCACGCGCA‑GGGGA‑CGGCCTTAGATCTCTCGCCGCCCATGCTTGTTCAGGCACGCCCGAAGGAt                                                                                                                                                                                >  1:2246790/1‑238 (MQ=255)
                                                  taaACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGG‑ATGAGCCGCCACTGGCGGGAACGGCACGCGCC‑GGGGA‑CGGGCTTAGATCTCTCGCCCCCCAAGCCTGTTCAGGCACGCCCGAAGGATg                                                                                                                                                                               >  2:568930/1‑239 (MQ=255)
                                                       aggcaaTTGGAGCGGCATTTGGCCGGGCAGCCGCACCCTAGGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGG‑ATGAGCCGCCACTGGCGGGAACGTCACGCGCA‑GGTGA‑CGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCa                                                                                                                                                                          <  2:2246790/234‑1 (MQ=255)
                                                           cATGGCAGCGGCATTTGGTGGGGAAGCCGCCCACTATGAGCAACATGCAGATCTACAGCGCCAGAGGGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACCCCCACGTACTGGACGGGGGTTGGGGACCTGGCTGG‑ATGAGCCGCCACTGGCGGGAACGTCACGCGCA‑GGTGA‑CGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCa                                                                                                                                                                          <  2:1564414/235‑1 (MQ=255)
                                                           aaTTGCAGCGGCTTTTGGGCGGGCAGCCGCACCCTATGAGCAACATGTAGTTCTACGGCGCTAGAGTGCTGCCGCCTTACTGGAAATGCTCCCACAGCGTAAATACCCCCACGTACGGGACGCGGGTTGGGGCCCGGGCTGG‑GTGGGCCGCCACTGGCGGGAACGTCACGCGCA‑GGTGA‑CGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGa                                                                                                                                                                        <  2:1013664/236‑1 (MQ=255)
                                                              tGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGG‑ATGAGCCGCCACTGGCGGGAACGTCACGCGCA‑GGTGA‑CGGCCTTAGATCTCTCGCCGCCAATGCTTGGTCAGGCACGCCAGAAGGATGCCGCAGAc                                                                                                                                                                       >  1:246358/1‑235 (MQ=255)
                                                              tGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGG‑ATGAGCCGCCACTGGCGGGAACGTCACGCGCA‑GGGGA‑CGGCCTTAGATCTCTCGCCCCCAATGCTTGTTCAGGCACCCCCGGACGAAGGCGCCGAACAtt                                                                                                                                                                   >  1:1717618/1‑239 (MQ=255)
                                                              tGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGG‑ATGAGCCGCCACTGGCGGGAACGTCACGCGCA‑GGGGA‑CGGCCCTAGATCTCTCGCCGCCCATGCTTGTTCAGGGCCGCCAGAAGGATGcc                                                                                                                                                                             >  1:1595970/1‑229 (MQ=255)
                                                                    ggCATTTGGCCGGGAGGCCGCACACTAGGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCTTTACTGGCAATGCTTCCACAGCGTAAATACCCCCACGTACTGGACGCGGGTTGTGGACCTGGCTGG‑ATGAGCCGCCACTGGCGGGAACGTCACGCGCA‑GGTGA‑CGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTAt                                                                                                                                                                 <  2:246358/235‑1 (MQ=255)
                                                                     gCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGG‑ATGAGCCGCCACTGGCGGGAAAGCCACGCCCC‑GGGGA‑CGGCCCTAGATCTCCCCCCGCCCATGCTTGTTCAGGGCCGCCACAAGGATGCCGCCGACCATTAttcggc                                                                                                                                                            >  2:739606/1‑234 (MQ=255)
                                                                                      cGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGG‑ATGAGCCGCCACTGGCGGGAACGTCACGCGCA‑GGTGA‑CGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAgg                                                                                                                                                                                  <  1:2321156/200‑1 (MQ=255)
                                                                                           acTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGG‑ATGAGCCGCCACTGGCGGGAACGTCACGCGCA‑GGTGa                                                                                                                                                                                                                                   <  1:1781230/147‑1 (MQ=255)
                                                                                                 aGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGG‑ATGAGCCGCCACTGGCGGGAACGTCACGCGCA‑GGTGA‑c                                                                                                                                                                                                                                 <  2:1803044/142‑1 (MQ=255)
                                                                                                 aGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGG‑ATGAGCCGCCACTGGCGGGAACGTCACGCGCA‑GGTGA‑c                                                                                                                                                                                                                                 >  1:1803044/1‑142 (MQ=255)
                                                                                                     aCATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGG‑ATGAGCCGCCACTGGCGGGAAAGGCACGCGCA‑GGGGA‑CGGCCCTAGa                                                                                                                                                                                                                        >  2:1781230/1‑147 (MQ=255)
                                                                                                       aTGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGG‑ATGAGCCGCCACTGGCGGGGACGTCACGCGCA‑GGTGA‑CGGCCCTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACCCCCGaa                                                                                                                                                                                    >  1:921917/1‑181 (MQ=255)
                                                                                                       aTGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGG‑ATGAGCCGCCACTGGCGGGAACGTCACGCGCA‑GGTGa                                                                                                                                                                                                                                   >  1:2366120/1‑135 (MQ=255)
                                                                                                       aTGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGG‑ATGAGCCGCCACTGGCGGGAACGTCACGCGCA‑GGTGa                                                                                                                                                                                                                                   <  2:2366120/135‑1 (MQ=255)
                                                                                                           aGATCTACACCGCCAGAGTGCTGCCGCCTTACTGGCAATGCTCCCACAGCGTAAATACCCCCACGTACTGGACGCGGGTTGTGGACCTGGCTGG‑ATGAGCCGCCACTGGCGGGAACGTCACGCGCA‑GGTGA‑CGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAgg                                                                                                                                                                                              <  2:1877808/167‑1 (MQ=255)
                                                                                                                tACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGG‑ATGAGCCGCCACTGGCGGGAACGTCACGCGCA‑GGTGA‑CGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAgg                                                                                                                                                                                              <  2:897992/162‑1 (MQ=255)
                                                                                                                tACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGG‑ATGAGCCGCCACTGGCGGGAACGTCACGCGCA‑GGTGA‑CGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAgg                                                                                                                                                                                              >  1:897992/1‑162 (MQ=255)
                                                                                                                    gcgcCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGG‑ATGAGCCGCCACTGGCGGGAACGGCACGCGCG‑GGGGA‑CGGCCCTAGATCTCCCGCCCCCCATGCCTGTTCCGGCACGCCAGAAGGAAGCCGCAGACCACTATCTGGCGGGGGAGATCG‑GAACCCCGCCCGTCGCGGCCGCGGCGTTCGATCtt                                                                                                             >  2:752452/1‑238 (MQ=255)
                                                                                                                     cgcCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGG‑ATGAGCCGCCACTGGCGGGGACGTCACGCGCA‑GGGGG‑CGGCCCTAGATCTCTCGCCGCCCATGCTTGTTCAGGCACGCCCGGAGGGAGCCGCAGACCATTATCTGGCGGGAGGTATCG‑AATTCCCGCCGGTTGCGACTGGGAGGCTCgcgcgtg                                                                                                            >  2:814886/1‑232 (MQ=255)
                                                                                                                     cgcCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGG‑ATGAGCCGCCACTGGCGGGAACGTCACGCGCA‑GGGGA‑CGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACCCCAGaa                                                                                                                                                                                    >  1:1877808/1‑167 (MQ=255)
                                                                                                                             tGCTGACGCCTTGCGGGCAATGCTTCCACGGCGTAAATCCCCCCACGTACTGGAGGCGGGTTGTGGACCTGGCTGG‑ATGGGCCGCCACTGGCGGGAACGTCACGCGCA‑GGTGA‑CGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAgg                                                                                                                                                                                              <  1:2317662/149‑1 (MQ=255)
                                                                                                                             tGCTGACGCCTTACTGGCAATGCTCCCACACCGTAAATACACCCACGTCCTGGACGCGGGTTGTGGCCCTGGCTGG‑ATGAGCCGCCACTGGCGGGAACGTCACGCGCA‑GGTGA‑CGGCCTTAGATCTCTcg                                                                                                                                                                                                                 <  2:1791994/130‑1 (MQ=255)
                                                                                                                                       ttACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGG‑ATGAGCCGCCACTGGCGGGGAAGTCCCGCGCA‑GGGGG‑CGGCCTTTGAATTCTCGCCGCCAATgg                                                                                                                                                                                                       >  1:1791994/1‑129 (MQ=255)
                                                                                                                                             gCAATGCTTCCACAGCGTAAATACCCCCACGTACGGGACGCGGGTTGTGGCCCTGGCTGG‑ATGAGCCGCCACTGGCGGGAACGTCCCGCGCA‑GGTGA‑CGGCCTTAGAtctct                                                                                                                                                                                                                   <  1:1397199/112‑1 (MQ=255)
                                                                                                                                                aTGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGG‑ATGAGCCGCCACTGGCGGGAACGTCACGCGCA‑GGTGa                                                                                                                                                                                                                                   <  1:1969112/94‑1 (MQ=255)
                                                                                                                                                aTGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGG‑ATGAGCCGCCACTGGCGGGAACGTCACGCGCA‑GGGGa                                                                                                                                                                                                                                   >  2:1969112/1‑94 (MQ=255)
                                                                                                                                                       cacgGCATAAATCCCCCCACGTCCTGGACGCGGGTTGTGGACCTGGCTGG‑ATGAGCCGCCACTGGCGGGAACGTCACGCGCA‑GGTGA‑CGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCG‑AATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCA‑TGGAGCAATCTCGCAGTGCAGTGGTGCGGTAAt                                                                        <  1:2092897/239‑1 (MQ=255)
                                                                                                                                                       cacaGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGG‑ATGAGCCGCCACTGGCGGGAACGGCCCCCGCA‑GGAGA‑CGGGCTTAGATCTCCCGCCGCCaaa                                                                                                                                                                                                         >  2:1397199/1‑111 (MQ=255)
                                                                                                                                                            cGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGG‑ATGAGCCGCCACTGGCGGGAACGTCACGCGCA‑GGTGA‑CGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTa                                                                                                                                                                  <  1:384525/146‑1 (MQ=255)
                                                                                                                                                             gTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGG‑ATGAGCCCCCACTGGCGGGGACGGTACCGGGT‑GGCGA‑CGGCCTTAGATCTCCCGCCGCCAATg                                                                                                                                                                                                        >  1:167335/1‑107 (MQ=255)
                                                                                                                                                              tAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGG‑ATGAGCCGCCACTGGCGGGAACGTCACGCGGA‑GGGGG‑AGGCCCTAGATCTCTCCCCGCCCATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCAGTAGCTGGCGGGAGATATTG‑AATCCCCGCCGTTAGCGACCGCGACGTTCGATCTTGCA‑TGGGGCGATCCTGCAATGCCGCGGTGCGGGAATTTAtac                                                                  >  2:1144300/1‑236 (MQ=255)
                                                                                                                                                                      cccACGTACTGGACGCGGGTTGTGGACCTGGCTGG‑ATGAGCCGCCACTGGCGGGGAACTCACGCGCC‑GGGGA‑CGGCCTTAGATCTCTCGCCGCCCAAGCTTGTTCAGGCACGCCCGAAGGATGGCCCAGACCATTATCTCGCGGGGGATATAGCAACCCCGCCCCTTGGGCACTCGGAC‑TTCTATCTTTCC‑TGGAGGCATCTCTCAGTGCAGTGGGGGGGGAATTTATCTCCGGCACt                                                          >  2:1624898/1‑238 (MQ=255)
                                                                                                                                                                      cccACGTACTGGACGCGGGTTGTGGACCTGGCTGG‑ATGAGCCGCCACTGGCGGGAACGTCACGCGCA‑GGGGA‑CGGCCTTAGATCTCTCCCCCCCCAAGCTTGTTCAGGCCCGCCCGAAGGATGCCGCCGACCATTATCTGGCGGGGGATATCG‑AATCCCTGCCGTTAGCGAATGCGACGGTCGATCTTGCA‑TGGAGCAATCTCGCCGTGGCGGGGGGGGGGAATTTATTCACCGCACt                                                          >  2:2044895/1‑238 (MQ=255)
                                                                                                                                                                      cccACGTACTGGACGCGGGTTGTGGACCTGGCTGG‑ATGAGCCGCCACTGGCGGGAACGGCACGCGCA‑GGGGA‑CGGGCTTAGGTCTCTCGCCGCCAATGCTTGTTCCGGCCCGCCAGAAGGATGGCGCCGACCCCTATCTGGGGGGGGGTATCG‑GATCCCCGCCCGTTGCGACTGCGACGTTCGGACTTGGA‑TGGGGCAATCTCGCAGGGGAGGGGGGGGGGCAATTATCCACgggccc                                                          >  2:904097/1‑234 (MQ=255)
                                                                                                                                                                      cccACGTACTGGACGCGGGTTGTGGACCTGGCTGG‑ATGAGCCGCCACTGGCGGGAACGGCACCCGCA‑GGGGG‑CGGCCCTTGATCTCTCGCCGCCCATGCTTGTTTAGGCCCGCCCGGAGGATGCCGCAGACCATTTTCTTGCcggg                                                                                                                                                        >  2:384525/1‑145 (MQ=255)
                                                                                                                                                                      cccACGTACTGGACGCGGGTTGTGGACCTGGCTGG‑ATGAGCCGCCACTGGCGGGAAAGGCACGCGCA‑GGGGA‑CGGCCTTAGATCTCTCGCCGCCAAAGCCTGTTTAGGCACGCCCGAAGGATGCCGCAGACCCTTATCTGGGGGGGGATATCG‑AAACCCCGCCCTTAGCGGCTGCGCCGTTTGGTCTTGGATTGGGGCAA‑CTCCCCGTGGCGGGGGGCGGGGATTTATccccagccac                                                          >  2:1186114/1‑230 (MQ=255)
                                                                                                                                                                          cGTACTGGACGCGGGTTGTGGACCTGGCTGG‑ATGAGCCGCCACTGGCGGGAACGTCACGCGCA‑GGGGA‑ACGCCTTAGAACTCTCGCCGCCCATCCTTGTTCAGGCACGCCCGAAGGAAGCCGCAGACCATTATCTGGCGGGGGATATCG‑CAACCCTGCCCTTTGCGACAGCGACGTTCGATCTTGCA‑TTGGGCAAACCCACAGTGCAGGGGGGCGGGGATTTATTCACGGCACTCcgcc                                                     >  2:570150/1‑238 (MQ=255)
                                                                                                                                                                             tgtggggggggggTGTGGGCCTGGGTGGTAGGGCCCGCCACGGGCGGGAACGTCCCGCGCA‑GGTGA‑CGGCCTTAGACCCCCCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGCCCATTATCTGGCGGGAGATATCG‑AATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCA‑TGGAGCAATCTCGCAGTGCAGTGGTGCGGTAATTTATCCACGGCACTCcg                                                       <  2:412802/226‑1 (MQ=255)
                                                                                                                                                                             gcTGGACGCGGTTTGTGGCCCGGGCTGG‑ATGAGCCGCCACTGGCGGGACCGTCACGCGCA‑GGTGA‑GGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCCCGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCG‑AATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCA‑TGGAGCAATCTCGCAGTGCAGTGGTGCGGTAATTTATCCACGGCACTCCGCGAg                                                   <  1:904097/237‑1 (MQ=255)
                                                                                                                                                                             agcggggggggggTGGGGCCCTGGGTGGTATGAGCCGCCGTGGGCGGGAGCGTCCCGCGCG‑GGTGC‑GGGCTTAAGTTCCTTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGCCCATTATCTGGCGGGAGATGTCG‑AATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCA‑TGGAGCAATCTCGCAGTGCAGTGGTGCGGTAATTTATCCACGGCACTCCGCGAg                                                   <  2:1742823/230‑1 (MQ=255)
                                                                                                                                                                             acgggACGGGGGTTGGGGACCGGGCTGG‑ATGGGCCGCCACTGGCGGGAACGTCACGCGCA‑GGTGA‑CGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCCCGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCG‑AATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCA‑TGGAGCAATCTCGCAGTGCAGTGGTGCGGTAATTTATCCACGGCACTCCGCGAg                                                   <  1:752452/235‑1 (MQ=255)
                                                                                                                                                                              cggggcgcgGGTTGTGGCCCTGGGGGG‑ATGGGCCGCCACTGGCGGGAACGTCACGCGCA‑GGTGA‑CGGCTTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCG‑AATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCA‑TGGAGCAATCTCGCAGTGCAGTGGTGCGGTAATTTATCCACGGCACTCCGCGAGCt                                                 <  1:729128/234‑1 (MQ=255)
                                                                                                                                                                                 gccgcgGGTTGTGGACCTGGCTGG‑ATGAGCCGCCACTGGCGGGAACGTCACGCGCA‑GGTGa                                                                                                                                                                                                                                   <  1:2363762/59‑1 (MQ=255)
                                                                                                                                                                                 gACGCGGGTTGTGGACCTGGCTGG‑ATGAGCCGCCACTGGCGGGAACGTCACGCGCC‑GGGGA‑CGGCCCTAGATCTCTCCCCCCCCATGCTTGTTTAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGAGATCG‑AATCCCTCCCGTTAGCGGCCGCGACGTTCGATCTTGCA‑TGGGGCAATCTCGCCGTGCAGTGGTGCGGGAATTTATCCACCGCCCCCCGCGa                                                    >  2:277697/1‑233 (MQ=255)
                                                                                                                                                                                 gACGCGGGTTGTGGACCTGGCTGG‑ATGAGCCGCCACTGGCGGGAACGGCCCGCGCC‑GGGGA‑CCGCCCTTGATCTCTCGCCGCCCATGCTTGTTTAGGCACGCCCCAAGGAGGCCCCAGACCCTTATTTGGCGGGGGAGATCA‑AAACCCCCCCGCTAGCGCCTGCGCCGTTTGACCTTTCA‑TTGGGCAATCTCGCGGTGCAGTGGGGCGGGGATTTATCCACCGCACTCCCCGAGCCGt                                               >  2:2070945/1‑238 (MQ=255)
                                                                                                                                                                                 gACGCGGGTTGGGGCCCTGGCTGG‑AGGGGCCGCCACTGGCGGGAACGTCACGCGCA‑GGTGA‑CGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCGGGCGGGAGATATCG‑AATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCA‑TGGAGCAATCTCGCAGTGCAGTGGTGCGGTAATTTATCCACGGCACTCCGCGAGCTGt                                               <  1:653371/238‑1 (MQ=255)
                                                                                                                                                                                      gggTTGTGGACCTGGCTGG‑ATGGGCCGCCACTGGCGGGAACGTCACGCGCA‑GGTGA‑CGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGa                                                                                                                                                                                 <  1:900128/105‑1 (MQ=255)
                                                                                                                                                                                           gtgGACCTGGCTGG‑ATGAGCCGCCACTGGCGGGAACCTCCCCCGCA‑GGTGA‑CGGCCCTaaa                                                                                                                                                                                                                        >  2:2363762/1‑59 (MQ=255)
                                                                                                                                                                                           ggggcccgggctgg‑AGGGCCGGCCCCGGGCGGGACCGGCACGCGCA‑GGTGA‑CGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCG‑AATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCA‑TGGAGCAATCTCGCAGTGCAGTGGTGCGGTAATTTATCCACGGCACTCCGCGAGCTGTATCGGgtggtg                                    <  1:1385080/231‑1 (MQ=255)
                                                                                                                                                                                           ggggACCGGGTGGG‑ATGGGCCGCCACTGGCGGGAACGTCACGCGCA‑GGTGA‑CGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCG‑AATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCA‑TGGAGCAATCTCGCAGTGCAGTGGTGCGGTAATTTATCCACGGCACTCCGCGAGCTGTATCGGgtggtg                                    <  1:2272655/237‑1 (MQ=255)
                                                                                                                                                                                            tgGACCTGGCTGG‑ATGAGCCGCCACTGGCGGGAACGTCACGCGCA‑GGTGA‑CGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTgg                                                                                                                                                             <  1:176764/119‑1 (MQ=255)
                                                                                                                                                                                              gcccgggctgg‑ATGAGCCGCCACTGGCGGGAACGTCACGCGCA‑GGTGA‑CGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCG‑AATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCA‑TGGAGCAATCTCGCAGTGCAGTGGTGCGGTAATTTATCCACGGCACTCCGCGAGCTGTATCGGGTGGTgcg                                  <  1:2104895/233‑1 (MQ=255)
                                                                                                                                                                                                ccTGGCTGG‑ATGAGCCGCCACTGGCGGGAACATCACGCGCA‑GGTGG‑CGGCCTTTGATCTCTCGCCGCCAATGCCTGTTCACGCCCGCCAGAAGGAGGCCGCAGAc                                                                                                                                                                       >  2:900128/1‑105 (MQ=255)
                                                                                                                                                                                                  tggctgg‑ATGAGCCGCCACTGGCGGGGACGGCACGCGCA‑GGTGG‑CGGCCTTTGATCTCTCGCCGCCAATGCTTGTTTAGGGACGCCCGAAGGATGCCGCCGACCCTTATCTGGCGGGGGGTATTG‑AATCCCTGCCGTTAGCGGACGCGACGGTTGATCTTGCA‑TGGAGGAATCTCgcagtg                                                                                       >  2:1136546/1‑181 (MQ=255)
                                                                                                                                                                                                    gcggg‑ATGGGCCGCCACTGGCGGGAACGTCCCGCGCA‑GGTGA‑CGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGAGGCCGCAGACCATTATCTGGCGGGAGATATCG‑AATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCA‑TGGAGCaa                                                                                                 <  1:1527265/166‑1 (MQ=255)
                                                                                                                                                                                                     ctgg‑ATGAGCCGCCACTGGCGGGAACGTCACGCGCA‑GGTGA‑CGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCa                                                                                                                                                                          <  1:170790/97‑1 (MQ=255)
                                                                                                                                                                                                      tgg‑ATGAGCCGCCACTGGCGGGGACCTCCCCCGCAGGGGGA‑CGCCCTT‑GATCCCTCCCCGCCCACACTTGTTTCGGCACGCCCGAAGGATTCCGCCGACCACTATCTTGCGGGGGatat                                                                                                                                                   >  2:176764/1‑119 (MQ=255)
                                                                                                                                                                                                         gaggggccgccACTGGCGGGAACGTCACGCGCA‑GGTGA‑CGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCG‑AATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCA‑TGGAGCAAtc                                                                                               <  1:1993422/162‑1 (MQ=255)
                                                                                                                                                                                                          aTGAGCCGCCACTGGCGGGAACGGCACCCGCA‑GGGGA‑CGGCCCTTGATCTCTCCCCGCCCATGCTTGTTCAGGCCCGCCAGAAGGATGCCGCAGACCATTATTTGGCGGGGGATATCG‑AATCCCCGCCGGTAGCGACCGCGACCGTTGGTCTTTCA‑TGGAGCGATCTCCCCGGGCAGGGGGGCGGGAATTTATTCACGGCACCCCGCGCGCTGTATCGGGGGGGGCG‑CCCCACAggg                       >  2:2104895/1‑236 (MQ=255)
                                                                                                                                                                                                              gccgccACTGGCGGGAACGTCACGCGCA‑GGGGA‑CGGCCTTAGATCTCTCGCCCCCCATGCTTGTTCAGGCACGCCCGAAGGATGCCGCCGACCCTTATCTGGCGGGGGATATCC‑AATTCCTGCCGGTTGCGACTGCGACGGTCGATCTTTCA‑TGGAGGCATCTCGCAGGGCCGTGGGGCGGGAATTTTTCCACGGGCCCCCCCCAGCTTGATCGGGGGGGGGGCCCCCAAAGGGGTggc                  >  2:1399266/1‑238 (MQ=255)
                                                                                                                                                                                                               ccgccACTGGCGGGAACGTCACGCGCC‑GGGGA‑CGGCCTTAGATCTCTCGCCGCCGACGCTTGTTCAGGCACGCCCGGAGGTTGCCCCCGACCCTTATTTGGCGGGAGATATAG‑AATTCCCGCCCGTGGCGGCTGCGACGTTCGATCTTTCA‑CGGGGCAATCTCgcaggg                                                                                       >  2:1527265/1‑167 (MQ=255)
                                                                                                                                                                                                                   cACTGGCGGGAACCTCACGCGCA‑GGGGA‑CGGCCTTAGATCTCTCGCCGCCCAAGCTTGTTTAGGCACGCCCGAAGGATGCCGCCGACCCTTATCCGGCGGGGGATATCG‑AATCCCTGCCGTTTGCGACCGCGACGTTCGATCTTGCA‑TGGAGGAAACCCgcagtgcag                                                                                    >  2:1993422/1‑168 (MQ=255)
                                                                                                                                                                                                                      gggcggcAGCGGCGCGCGAA‑GGTGC‑CGGCCTTAGCTCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCG‑AATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCA‑TGGAGCAATCTCGCAGTGCAGTGGTGCGGTAATTTATCCACGGCACTCCGCGAGCTGTATCGGGTGGTGCG‑CCCCAAAGGCGTGGTCGCGTTTAc         <  2:651513/238‑1 (MQ=255)
                                                                                                                                                                                                                              aCGTCACGCGCA‑GGTGC‑CGGCTTTAGATCTCCCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCG‑AATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCA‑TGGAGCAATCTCGCAGTGCAGTGGTGCGGTAATTTATCCACGGCACTCCGCGAGCTGTATCGGGTGGTGCG‑CCCCAAAGGCGTGGTCGCGTTTACCACGCTg  <  1:814886/238‑1 (MQ=255)
                                                                                                                                                                                                                                gTCACGCGCA‑GGTGA‑CGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCG‑AATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCA‑TGGAGCAATCTCGCAGTGCAGTGGTGCGGTAATTTATCCACGGCACTCCGCGAGCTGTa                                              >  1:1228772/1‑193 (MQ=255)
                                                                                                                                                                                                                                         a‑GGTGA‑CGGCTTTAGATCTCTCGCCGCCAATGCTTGTTCAGGc                                                                                                                                                                                             <  1:237461/43‑1 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                   
AAATGCAGGATATCGACCGTCTGCTGGAGGTGCTGCATGGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGG‑ATGAGCCGCCACTGGCGGGAACGTCACGCGCA‑GGTGA‑CGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCG‑AATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCA‑TGGAGCAATCTCGCAGTGCAGTGGTGCGGTAATTTATCCACGGCACTCCGCGAGCTGTATCGGGTGGTGCG‑CCCCAAAGGCGTGGTCGCGTTTACCACGCTG  >  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: ‑