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

GAAATGCAGGATATCGACCGTCTGCTGGAGGTGCTGCATGGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCGGTAATTTATCCACGGCACTCCGCGAGCTGTATCGGGTGGTGCGCCCCAAAGGCGTGGTCGCGTTTACCACGC  >  NC_000913/811485‑811942
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                    
gAAATGCAGGAGAGCGGCCGTCTGCGGGAGGTGCGGAAGGGCACCGTTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACCCCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGGGCCGCCACTGGCGGGAACGTCACGCGCAGGTGa                                                                                                                                                                                                                                   <  2:1072104/239‑1 (MQ=255)
  aaTGCAGGATATCGACCGTCTGCTGGAGGTGCTGCATGGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGAc                                                                                                                                                                                                                                  <  1:1919428/238‑1 (MQ=255)
      agggttaccGACCGTCGGCTGGAGGTGCGGCAGGGAAACGTTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGa                                                                                                                                                                                                                                   <  2:884558/225‑1 (MQ=255)
                   gTCTGCTGGAGGTGCTGCATGGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGGCCCGGCCCGGTGACGGCCTTAGAttctacgc                                                                                                                                                                                                                 >  1:34852/1‑231 (MQ=255)
                   gTCTGCTGGAGGTGCTGCATGGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGGCCCGCGCAGGGGACGGCCTTAGATATCTggcc                                                                                                                                                                                                                >  1:568417/1‑239 (MQ=255)
                   gTCTGCTGGAGGTGCTGCATGGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAAAGTCAAGCGCCGGGGGACGCCTTAGGTCTCTcgc                                                                                                                                                                                                                 >  1:1241547/1‑238 (MQ=255)
                       gctgGAGGTGCTGCATGGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGGCACGCGCGCGGGAAGGCCTGCAGTCTCTCgccgccg                                                                                                                                                                                                            >  2:673422/1‑238 (MQ=255)
                               tgctgcATGGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTATGGACCTGGCTGGATGAGCCGCCACTGGCGGGAAAGGCACGCGCAGGGGACGGCCTTAGATCTCTCGCCGCCAATGCgtgt                                                                                                                                                                                                  >  1:1764998/1‑239 (MQ=255)
                                gctgcATGGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCCCGCGCAGGGGACGGCCCTAGATCTCTCGGCGCCCAAGCCTGt                                                                                                                                                                                                  >  1:469145/1‑238 (MQ=255)
                                     tgggCAACGGTAAATAAACAAGCCTTTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATCCCCCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGAGGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAgg                                                                                                                                                                                             <  1:171378/236‑1 (MQ=255)
                                     aTGGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAAAGGCACGCGCAGGGGACGGCCTTAGATCTCTCGCCCCCCAAGCTTGTTTAgg                                                                                                                                                                                             >  1:74998/1‑238 (MQ=255)
                                      tGGCACCGGTAAATAAACAAGCCATTGCGGGGGCATTTGGTGGGGCTGCCGCACACTATGAGCAACAGGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGTTTCCACAGCGTAAATACACCCACGGACGGGACGCGGGTTGTGGCCCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCGTAGATCTCTCGCCGCCAATGCTTGTTCAGGCa                                                                                                                                                                                           <  2:323537/239‑1 (MQ=255)
                                       ggCAACGGTAAATAACAAAGCATTTGCAGCGGCATTGGGTCGGGGAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGGGCTGACGCCTTCCTGGCAATGCTTCCACAGCGTAAATACCCCCACGTACTGGACGCGGGTTGGGGACCGGGCGGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCa                                                                                                                                                                                           <  2:2223379/238‑1 (MQ=255)
                                         cAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGCACGTCACGCGACGGTGGCGGCCTTAGGTCTCTCGCCCCCCATGGTTGGTAACGCACg                                                                                                                                                                                         >  1:1005790/1‑238 (MQ=255)
                                          aaCGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGc                                                                                                                                                                                        <  1:1699149/238‑1 (MQ=255)
                                              gTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGGCACGCGCAGGGGACGGCCTTCGGTCTCTCCCCGCCCATGCTTGTTCCgggcagcccgaa                                                                                                                                                                                   >  1:1554575/1‑229 (MQ=255)
                                                 aataaACAAGCCAGGGCAGCGACATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTCCCACAGCGTAAATACCCCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGAc                                                                                                                                                                                                                                  <  1:1171739/191‑1 (MQ=255)
                                                        aaGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGAAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCa                                                                                                                                                                         <  1:915384/239‑1 (MQ=255)
                                                         aGCCATTGCAGCGGCATTTGGTCTGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGGCCCACGCAGGGGCCGGGCTTAGATCTCTCTCTCCCCAGCCCTGCGCACGCACGCCAGACGGGTGCCgtca                                                                                                                                                                        >  2:744485/1‑236 (MQ=255)
                                                         aGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGGACGGTAAGCGGCGGGGACGGCCTTAGGTCTCTCGCCCCCCATGCCTTGTCTAGGCCCGCCGGaggaggcggcc                                                                                                                                                                          >  2:719035/1‑228 (MQ=255)
                                                         aGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGGAAGGCCCCCGCCGGGGAAGGGCCTCGAATCTCTCCCGCCAATGCCTTGTCCAGCCGCGCCCAAGGACGCCGCa                                                                                                                                                                         >  2:418489/1‑238 (MQ=255)
                                                         aGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAAAGGCACGCGCCGGGGGAGGCCCTAGATCTCTCCCCGCCCATGCCTGTTCAGGCACCGCCGAAGGATGgcgccg                                                                                                                                                                        >  2:1362532/1‑237 (MQ=255)
                                                          gCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGCCGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAgg                                                                                                                                                                                             <  1:396418/217‑1 (MQ=255)
                                                           ccATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGGCACGGCCAGGGGGAGGCCCTAGaa                                                                                                                                                                                                                        >  2:1171739/1‑190 (MQ=255)
                                                            cgttGCAGGTGCTTTTGGTCGGGCGGCCCGCCCCTATGGGCAACGGTCTGTGCCACAGCGCCAGGGTGCGGACGCCTTACTGGCAGTGCTTCCAGAGGGTAAATCCACCCCGGTACGGGCCGCGGGTTGGGGACCTGGGGGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGGGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGa                                                                                                                                                                       <  2:2111036/235‑1 (MQ=255)
                                                            cATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGa                                                                                                                                                                       <  1:161720/237‑1 (MQ=255)
                                                                  aGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGa                                                                                                                                                                                                                                   >  2:1857840/1‑173 (MQ=255)
                                                                  aGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGa                                                                                                                                                                                                                                   <  1:1857840/173‑1 (MQ=255)
                                                                          ttGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACGGCGCCGGAGTGCTGACGCCTTACGGGCAATGCTTCCACAGCGTAAATACCCCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGga                                                                                                                                                      <  1:304887/239‑1 (MQ=255)
                                                                             gTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTTGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGa                                                                                                                                                                                                                                   <  1:364441/162‑1 (MQ=255)
                                                                                gggCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAAGACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGatat                                                                                                                                                 <  1:2160805/238‑1 (MQ=255)
                                                                                        gaacacTAGGAGAAACATGCAGATCTACAGCGCCAGAGTGCTGACGCTTTACGGGCAATGCTTCCACAGCGTAAATCCCCCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGa                                                                                                                                                                       <  2:1683183/207‑1 (MQ=255)
                                                                                        gCACCCTATGAGCAACATGCAGATCTACAGCCCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACCCCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCa                                                                                                                                                                         <  1:1123631/207‑1 (MQ=255)
                                                                                            acTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACCCCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAgg                                                                                                                                                                                             <  1:1893136/183‑1 (MQ=255)
                                                                                              tATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACCCCCCCGTACTGGCCGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGga                                                                                                                                                      <  1:642483/219‑1 (MQ=255)
                                                                                              tATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGa                                                                                                                                                                                                                                   <  1:890006/145‑1 (MQ=255)
                                                                                                tGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGGACGGCACGCGCCTGGGAAGGCCCTAGATCTCTCGCCGCCCACGCTTGTTCTAGCCCGCCAGAAAGGTGCCTCCGACCCTTGTCTGGCGGGGGGTATCGAATCCCAGCCGTCAg                                                                                                                                >  2:1713459/1‑239 (MQ=255)
                                                                                                tGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGCCACCCGCAGGGGACGGCCCTTGGACTCTCGCCGCCCAAGCTTGTTCGTGGACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGGGATATCCAATCCCTGTCGttc                                                                                                                                 >  2:1956017/1‑237 (MQ=255)
                                                                                                  aGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGc                                                                                                                                                                          <  1:408051/196‑1 (MQ=255)
                                                                                                  aGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACCCGCAGGGGACGGCCTTAGATCTCTCTCCGCCAATGCTTGTTCAGGCACGCCCGAAGGATGCCGGAGCCCCTTATccg                                                                                                                                                             >  1:1683183/1‑207 (MQ=255)
                                                                                                  aGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGGCACGCGCAGGGGAGGGGCTTAGGTCTCTCGCCGCCCATGCTTGTTCAGGCACGCCAGAAGGGTGCCGCAGACCCTTATc                                                                                                                                                               >  2:1123631/1‑207 (MQ=255)
                                                                                                        aTGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGGACGTCCCGCGCAGGGGACGGGCTTAGGTCTCCCGCCGCCCATGCTTTGTCAGGCACGCCAGACGGATGCCGCAGACCACTATCTGGCGGGAGAGATCGCATCCCTGCCGTTAGCGACcgcg                                                                                                                        >  2:555279/1‑239 (MQ=255)
                                                                                                        aTGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAAAGGCCCGCGCCGGGGGCGGCCCTAGa                                                                                                                                                                                                                         >  2:890006/1‑145 (MQ=255)
                                                                                                        aTGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAAAAGCACGCGGCGGGGACGGCCTTTGATCTCTCGCCGCCAAAGCTTGTTCAGGCCCCCCAGAAGGCTGCCGCAGACCCTTTTc                                                                                                                                                               >  1:1970783/1‑201 (MQ=255)
                                                                                                        aTGCAAATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGCACGTCACGCCCCGGGGACGGCCTTTGATCCCTCCCCGCCCATGCTTGTTCCGGCACCGCAGACGGGAGGCGCCGACCCTTTTCCGGCGGGGGATTTTGCATACCCGCCCCTAGCGAccccg                                                                                                                        >  2:1557663/1‑235 (MQ=255)
                                                                                                         tGCAGATCTACAGCGCCAGAGTGCTGACGCCTTCCTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGa                                                                                                                                                                                                                                   <  1:2060280/134‑1 (MQ=255)
                                                                                                            aGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACCTCCCGCGCAGGGGGCGGGCTTAGAACTCTCGCCCCCCATGCTTTGTCAGGCCCGGCCGGAGGACGGCGCCGACCaattt                                                                                                                                                                >  2:408051/1‑192 (MQ=255)
                                                                                                                   cAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGGCCCGCGCAGGGGACGGCCTTAGa                                                                                                                                                                                                                         >  2:2060280/1‑134 (MQ=255)
                                                                                                                         cAGAGTGCTGACGCCTTACTGGCAAGGCTTCCACGGCGTAAATACACCCACGTGCTGGACGCGGGTTGTGGACCTGGCTGTATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGc                                                                                                                                                                          <  1:1431731/173‑1 (MQ=255)
                                                                                                                                    cGCCTTACTGGCAATGCTCCCACAGCGTAAATACCCCCCCGTGCGGGGCGCGGGTTGTGGACCTGGCGGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCg                                                                                           <  1:826637/239‑1 (MQ=255)
                                                                                                                                     gCCTTACTGGCAATGCTTCCACAGCGTAAATACCCCCACGTACTGGACGCGGGTTGTGAACCTGGCGGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAAt                                                                                                                                                                                                         <  1:823870/131‑1 (MQ=255)
                                                                                                                                      ccTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGa                                                                                                                                                                                                                                   <  1:1730433/105‑1 (MQ=255)
                                                                                                                                      ccTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGa                                                                                                                                                                                                                                   >  2:1730433/1‑105 (MQ=255)
                                                                                                                                               cAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGGGGCGGCCTTAGATCCCTCGCCGCCACTGCTTGTTCAg                                                                                                                                                                                              >  2:823870/1‑131 (MQ=255)
                                                                                                                                                     ttCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAgg                                                                                                                                                                                             <  1:2079344/126‑1 (MQ=255)
                                                                                                                                                         acaGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGa                                                                                                                                                                       <  1:876975/144‑1 (MQ=255)
                                                                                                                                                             cGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGATGGATGAGCCGCCACTGGCGGGAACGTCCCGCGCAGGTGACGGCCTTAGATCTCTCgccgcc                                                                                                                                                                                                             >  1:1440014/1‑104 (MQ=255)
                                                                                                                                                             cGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGATGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCgccgcc                                                                                                                                                                                                             <  2:1440014/104‑1 (MQ=255)
                                                                                                                                                               tAAATCCACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGGGCGGGGGACGGCCTTTGATCTCTCGCCGCCCCAGCTTGTTCAAGCCCCCCCGAAGGATGCCGCAGACCATTATCTGGGCGGGGGAATTCAACCCCCTCCCGTAGCGACCGCGTCCTTCGGTCCTTCATGGAGCACTCCCCCACGGCAGGGGTGGCGGCATTTATcc                                                                 >  2:1606494/1‑238 (MQ=255)
                                                                                                                                                               tAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAAAGGCACGCGCAGGGGACGGCCTTAGATCTCTCGCCCCCCATGCTTGTTCAGGCCCGCCCGaa                                                                                                                                                                                   >  2:2079344/1‑126 (MQ=255)
                                                                                                                                                                aaaTACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCa                                                                                                                                                                         <  1:998302/135‑1 (MQ=255)
                                                                                                                                                                   tACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGGAAGGCACGCGCAGGGGACGGCCCTTGATCTCTCCCCGCCACTGCTTGTTCAGGGCCGCCCGAAGGATGGCGCAGACCATTATCTg                                                                                                                                                            >  2:876975/1‑144 (MQ=255)
                                                                                                                                                                       cccACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACCGCACGCGCCGGGGGCGGCCTTAGATCTCTCCCCGCCCATGCTTGTTCAGGCCCCCCCGAAGGATGCCGCAGACCCTTATCTGGCGGGGGGTATCGGATCCCCGCCGGTTGCGCCTGCGAAGTTCGATATTGCATGGGGCAATCCTCCAGTGCAGAGGGGCGGGAAGTTATTCACCGGACt                                                         >  2:1471429/1‑238 (MQ=255)
                                                                                                                                                                         cccGTGCTGGACGGGGGTTGTGGACCTGGCTGGAGGGGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCTTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCGGTAATTTATCCACGGCACTCcg                                                      <  1:673422/237‑1 (MQ=255)
                                                                                                                                                                         cACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACCCGCAGGGGACGGCCCTAGATCTCTCGCCGCCCATGCTTGTTCAGGCACGCCCGAAGGGGGCCGCCGACCATTATCTGGCGGGGGGTATCGAATCCCCGCCGGTTGCGGACGCGGAGGTCGACTTTTGCTGGGGCACATTCGCAGTGCGTGGGGGGGGGAATTTTTCCACGGcgcacc                                                       >  2:910962/1‑233 (MQ=255)
                                                                                                                                                                          aCGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTcgc                                                                                                                                                                                                                 <  2:356689/87‑1 (MQ=255)
                                                                                                                                                                          aCGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTcgc                                                                                                                                                                                                                 >  1:356689/1‑87 (MQ=255)
                                                                                                                                                                          aCGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACCTCACGCGCCGGTGGCGGCCTTAGATCTCTCCCCGCCCATGCTTTTTCAAGCACGCCAGAAGGAAGCCGCCGACCCTTa                                                                                                                                                                 >  2:998302/1‑133 (MQ=255)
                                                                                                                                                                            gTACGGGACGCGGGTTGGGGCCCGGGCTGGATGGGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCGGTAATTTATCCACGGCACTCCGCGa                                                   <  1:1362532/239‑1 (MQ=255)
                                                                                                                                                                             tACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGc                                                                                                                                                                          <  1:1647366/121‑1 (MQ=255)
                                                                                                                                                                              cgtGGACGGGGTTGTGGGGCTGGGGTGGCGGACCCGCCACGGGAGGGACGGTCGCGGGAGGGGGAGGGCCTTCGATCTCCCGCCGCATAGGCTTGGTAAGGCGCGCCAGAAGGAGGCCGCAGACAATTATCTGGCGGGAGATTTCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCGGTAATTTATCCACGGCACTCCGCGAg                                                  <  2:71388/236‑1 (MQ=255)
                                                                                                                                                                               ccgggccggggTTGGTGGCCTTGGCGGGAGGGCCCGCCACGGGCGGGAACGTCCCGCGAAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCGGTAATTTATCCACGGCACTCCGCGAg                                                  <  1:1713459/230‑1 (MQ=255)
                                                                                                                                                                                  gACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGa                                                                                                                                                                       <  1:2035565/119‑1 (MQ=255)
                                                                                                                                                                                  gACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGa                                                                                                                                                                       >  2:2035565/1‑119 (MQ=255)
                                                                                                                                                                                       gggTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCCCGCGCCGGGGGCGGGCCTAGATCTCTTGCCGCCCAAGATTGTTTAGCCACGCCAGAAGCGGGCCGCAGACCacttt                                                                                                                                                                >  2:1647366/1‑117 (MQ=255)
                                                                                                                                                                                          ttGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCa                                                                                                                                                                         <  1:1333218/109‑1 (MQ=255)
                                                                                                                                                                                                      ggctggATGAGCCGCCACTGGCGGGAACGGCACGCGCAGGGGACGGCCTTAGATCTCTCGCCCCCCATGCTTGTTCAGGCACGCCCGAAGGATGGCGCAGACCCTTATc                                                                                                                                                               >  2:1333218/1‑109 (MQ=255)
                                                                                                                                                                                                       cggggaggggccgccACTGGCGGGACCGTCACGCGCAGTTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCCCGCCAGAAGAGTGCCGCAGACCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCGGTAATTTATCCACGGCACTCCGCGAGCTGTATCGGGTGGTgcg                               <  1:582223/223‑1 (MQ=255)
                                                                                                                                                                                                           gATGGGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCgccgcc                                                                                                                                                                                                             <  1:623737/60‑1 (MQ=255)
                                                                                                                                                                                                            tagacccgccCCGGGCGGGACAGTCCCGCCCGGGTGACGGCTTAAGATCTCTCCCGCCTAAGGCTTGTGCAGCCCCCCCAGAAGTATGCGGCATCCCATTATCTGGGGGACAATTTGGAATACGTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCGGTAATTTATCCACGGCACTCCGCGAGCTGTATCGGGTGGTGCGCCCCAAATGCg                    <  2:568417/233‑1 (MQ=255)
                                                                                                                                                                                                             tGAGCCGCCACTGGCGGGAACGTCCCGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCGGTAATTTATCCACGGCACTCCGCGAGCTGTATCGGGTGGTGCGCCCCAAAGGCGTg                  <  1:1471429/239‑1 (MQ=255)
                                                                                                                                                                                                                gccgccCCTGGGGGGACCGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCGGTAATTTATCCACGGCACTCCGCGAGCTGTATCGGGTGGTGCGCCCCAAAGGCGTGGt                <  1:1452466/238‑1 (MQ=255)
                                                                                                                                                                                                                acacccACGGGCGGGAACTTCCCGCGCAGGTGACGGCTTTAGATCTTCCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGCCCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCGGTAATTTATCCACGGCACTCCGCGAGCTGTATCGGGTGGTGCGCCCCAAAGGCGTGGt                <  1:706241/234‑1 (MQ=255)
                                                                                                                                                                                                                 ccgccACTGGCGGGAACGTCACGCGCCGGGGACGGCCTTAGATCTCTCGCCGCCCATGCTTGTTCAGGCACGCCCGAAGGATGCCGCAGACCATTATCTGGCGGGGGATATCGAATCCCCGCCGTTTTCGACAGCGACGTTTGATCTTGCATGGGGCAATCTCGCCGTGGAGGGGGGCGGGAAATTATTCTCCGCACTCCCCCAGCGTGAATGGGGGGGGGGCCCCCAAGGCGCgggtc              >  2:1403698/1‑236 (MQ=255)
                                                                                                                                                                                                                     cACTGGCGGGAACGTCACCCGCAGGTGACGGCCCTAGATCTCTCGCCGCCAATGCTTGTTCAGGCGCGCCCGAAGGATGCCGCAGACCATTATCTGGCGGGGGATATCGAATCCCCGCCGTTAGCGGCTGCGACGTTCGATCTTTCATGGAGCAATACCGCAGTGCAGTGGGGCGGTAATTTTTCCACCGCCCTCCGCCGCCTGGTTCGGGGGGTGGGGCCCCAAAgggggggcgtc             >  2:1350480/1‑228 (MQ=255)
                                                                                                                                                                                                                     cACTGGCGGGAAAGTCACGCGCAGGTGACGGCCTTAGAATTCCCGCCGCCAATGCTTGtt                                                                                                                                                                                                 >  2:623737/1‑60 (MQ=255)
                                                                                                                                                                                                                              gAACGTCACGCGCAGGTGCCGGCTTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGCCCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCGGTAATTTATCCACGGCACTCCGCGAGCTGTATCGGGTGGTGCGCCCCAAAGGCGTGGTCGCGTTTACCACGc  <  1:2056199/238‑1 (MQ=255)
                                                                                                                                                                                                                              gAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTAAGGCCCGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCGGTAATTTATCCACGGCACTCCGCGAGCTGTATCGGGTGGTGCGCCCCAAAGGCGTGGTCGCGTTTACCACGc  <  1:160715/238‑1 (MQ=255)
                                                                                                                                                                                                                               aaCGTCACGCGCAGGGGACGGCCTTAGATCTCTCGCCGCCCATGCTTTTTCAGGCACGCCCGAAGGATGCCGCAGACCATTATCTGGCGGGGGGTATCGAATCCCTGCCGTTAGCGACCGCGACGGTCGGTCTTGCATGGAGCAATCTCCCCGTGCCGTGGGGCGGGAATTTTTTCCCGGCACCCCCCGCGCCTGATCGGGTGGTGGGCCCCCAAGGACGGGTcgccgtttcc      >  2:706241/1‑226 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                    
GAAATGCAGGATATCGACCGTCTGCTGGAGGTGCTGCATGGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCGGTAATTTATCCACGGCACTCCGCGAGCTGTATCGGGTGGTGCGCCCCAAAGGCGTGGTCGCGTTTACCACGC  >  NC_000913/811485‑811942

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: