Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913811,7210TG32.9% 100.0 / 10.6 74V67G (GTG→GGG) bioCmalonyl‑acyl carrier protein methyltransferase
Reads supporting (aligned to +/- strand):  ref base T (7/42);  new base G (20/4);  total (28/46)
Fisher's exact test for biased strand distribution p-value = 1.26e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CTGCGCATGAAATGCAGGATATCGACCGTCTGCTGGAGGTGCTGCATGGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCGGTAATTTATCCACGGCACTCCGCGAGCTGTATCGGGTGGTGCGCCCCAAAGGCGTGGTCGCGTTTACCACGCTGGTGCAGGGATCGTTACCCGAACT  >  NC_000913/811477‑811967
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                        
cTGCGCATGAAATGCAGGATATCGACCGTCTGCTGGAGGTGCTGCATGGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGGACGTCACGCGCAGGGGACGGc                                                                                                                                                                                                                                                    >  1:73807/1‑251 (MQ=255)
   cgcATGAAAGGCAGGATATCGACCGTCTGCTGGAGGTGCTGCATGGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCtt                                                                                                                                                                                                                                                 <  1:131620/251‑1 (MQ=255)
          aaTGCAGGATATCGACCGTCTGCTGGAGGTGCTGCATGGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACCCTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGGCACGCGGAGGGGACGGGCCTAGAtccc                                                                                                                                                                                                                                          >  1:417755/1‑249 (MQ=255)
          aaTGCAGGATATCGACCGTCTGCTGGAGGTGCTGCATGGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGAtctc                                                                                                                                                                                                                                          <  1:99323/251‑1 (MQ=255)
          aaTGCAGGATATCGACCGTCTGCTGGAGGTGCTGCATGGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGGGACGGCCTTAGAtctc                                                                                                                                                                                                                                          >  1:193806/1‑251 (MQ=255)
          aaTGCAGGATATCGACCGTCTGCTGGAGGTGCTGCATGGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACGGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGAtctc                                                                                                                                                                                                                                          <  1:644465/251‑1 (MQ=255)
           aTGCAGGATATCGACCGTCTGCTGGAGGTGCTGCATGGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGAtctct                                                                                                                                                                                                                                        <  1:286497/251‑1 (MQ=255)
             gCAGGTTATCGACCGTCTGCTGGAGGTGCTGCATGGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTcg                                                                                                                                                                                                                                      <  1:25778/251‑1 (MQ=255)
               aGGATATCGCCCGTCTGCTGGAGGTGCTGCATGGCAACGGTTAATAAACAAGCCAGTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTcgcc                                                                                                                                                                                                                                    <  1:506479/251‑1 (MQ=255)
               aGGATACCGACCGTCTGCGGGAGGGTCTGCAGGGAAACGGTGAATAAATAAGCCATTGCAGCGGTATTTGGTCGGGCAGCCGCCCACGATGAGCAACATGCAGATCTACAGCGCCAGAGGGTTGACGCCTTACTGGCTATGCTTCCACATCGTAAATACCCCCCTGTTCTGTATGCGTGGTGTGGACCGGGCTGGAGGAGCCGCCACTGTCGTGATCTTCACGCTCAGGTGACGGCCTTGGATCTCTcgcc                                                                                                                                                                                                                                    <  2:596891/251‑1 (MQ=255)
                 gATATCGACCGTCTGCTGGAGGTGCTGCATGGCAACGGTTAAAAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTcgccgc                                                                                                                                                                                                                                  <  1:580358/251‑1 (MQ=255)
                    atCGACCGTCTGCTGGAGGTGCTGCATGGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGGACGTCCCGCGCAGGGGGCGGGCCTAGATCTCTCgccgccc                                                                                                                                                                                                                                >  1:270183/1‑249 (MQ=255)
                    atCGACCGTCTGCTGGAGGTGCGGCATGGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCGACAGCGCCAGAGTGCGTACGCCTTACTGGCAATGCTTCCACAGCGTAAATCCCCCCACGTACTGGACGCGGGTTGTGGACCTGGCTGTATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCaa                                                                                                                                                                                                                               <  1:633023/251‑1 (MQ=255)
                       gACCGTCTGCTGGAGGTGCTGCATGGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTATATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGc                                                                                                                                                                                                                            <  1:458400/251‑1 (MQ=255)
                          cGTCTGCTGGAGGTGCTGCATGGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTg                                                                                                                                                                                                                         <  1:377508/251‑1 (MQ=255)
                          cGTCTGCTGGAGGTGCTGCATGGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGGGGAGGCCCTAGATCTCTCGCCGCCAATGCTTg                                                                                                                                                                                                                         >  1:396691/1‑251 (MQ=255)
                                  ggAGGTGCTGCATGGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCCCGCGCAGGGGACGGCCTTAGATCTCTCGCCGCCCATGCTTGTTCAGGcc                                                                                                                                                                                                                 >  1:189541/1‑250 (MQ=255)
                                   gAGGTGCTGCATGGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAAAGTCACGCGCAGGGGAAGGCCCTAGATCTCTCGCCGCCAATGCCTGTTccagcac                                                                                                                                                                                                                >  1:39722/1‑245 (MQ=255)
                                   gAGGTGCTGCATGGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACAGGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCAc                                                                                                                                                                                                                <  1:86591/251‑1 (MQ=255)
                                        gccacagtgaacctgtaaataaACAAGCCATTGCAGGTGCTTTTGGTCGGGCAGCCGCGCACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCTTTTCTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTTGACGCGGGTTGTGGACCTGGCTGTATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCGCGCCGCCAAGGCGTGGTCAGGCACGCCa                                                                                                                                                                                                            <  2:439134/233‑1 (MQ=255)
                                           gcATGGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGaa                                                                                                                                                                                                         <  1:532904/250‑1 (MQ=255)
                                              gggCAACGTTAAAAAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACGGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGAt                                                                                                                                                                                                     <  2:73807/250‑1 (MQ=255)
                                                gCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGc                                                                                                                                                                                                   <  1:266315/251‑1 (MQ=255)
                                                  aaCGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGGACGTCACGCGCAGGGGGCGGGCCTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGGAGGATGCCg                                                                                                                                                                                                 >  1:128228/1‑251 (MQ=255)
                                                    gggtaaaaaaccAAGCCAGTGCAGCGGAATGTGGCCGGGCAGCCGCACACTATGGGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAAGGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTGGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGAGCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCa                                                                                                                                                                                               <  2:561543/240‑1 (MQ=255)
                                                    cGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATACAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAAAGGCACGCGCAGGTGGCGGCCTTAGATTCATCGCCGCCAAAGCTTGTTTAGGCACGCCAGCAGGATGACGcc                                                                                                                                                                                               >  1:498135/1‑250 (MQ=255)
                                                      gTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGGACGTCCCGCGCAGGGGACGGCCTTAGATCTCTCGCCGCCCAAGCTTGTTCAGGCACCCCCGAAAGGTGCCGccg                                                                                                                                                                                              >  1:247528/1‑248 (MQ=255)
                                                         aataaACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGGACGTCACGCGCCGGGGACGGCCTTAGAAATCTCCCCGCCCACGCCTGTTCCGGGACGCCAGAAGGAAGCCGCAgcac                                                                                                                                                                                           >  1:376357/1‑247 (MQ=255)
                                                         aaaaaaCAAGCCATTGCAGCGGCATTGGGTCGGGCGGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCa                                                                                                                                                                                          <  2:546035/248‑1 (MQ=255)
                                                            aaaCAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCCTAGATCTCTCGCCGCCAATGATTGTTCACGCACGCCCGAAAGATGCCGCAGACCAttt                                                                                                                                                                                       >  1:487555/1‑250 (MQ=255)
                                                              aCAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATc                                                                                                                                                                                     <  1:244331/251‑1 (MQ=255)
                                                               cAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCt                                                                                                                                                                                    <  1:345711/251‑1 (MQ=255)
                                                                aaGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCGACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTGGCCGCCGATGCGTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTg                                                                                                                                                                                   <  1:332908/251‑1 (MQ=255)
                                                                   ccATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACCCCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCg                                                                                                                                                                                <  2:270183/251‑1 (MQ=255)
                                                                      ttGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGACACACCCAGGTGAAGGACCTACATCTCTCTCCCCCCATGCTTGTTCAGGCACGCCCGAAGGATGACGCCAACCCTTATCTGGCgggg                                                                                                                                                                             >  1:652977/1‑250 (MQ=255)
                                                                                aTTTGGTCGGGCAGCAGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGGACGCCACCCGCCGGGGGCGGCCCTAGATACCTCCCCGCCCATGCTTGTTCAGGCCCGCCCGAAGGAGGCCGCAGACCACTTTCTTGCGGgaggattccaa                                                                                                                                                                    >  2:284955/1‑242 (MQ=255)
                                                                                 tttggccgggccaccgaacgctaggAGCAACAGGCAATTCTACACCGCCAGGGGTCTGACGCCTTACTGGCAAGGCTTCCACAGCTGAAAGACACCCTCGTCCGGCACGCGGGTTTTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCATACCAGGATCTGGCGGGACATATCGAAt                                                                                                                                                                   <  2:280906/226‑1 (MQ=255)
                                                                                 gtgggccggccagccccctacTATGAGCAACATGCAGATCTACAGCGCGAGAGGGCGGACGCCTTCCGGGCAATGTTTCCCCAGGGAAAACCCCCCCCCGGAGTGGAGGCGGGTGTGGGACCTGTCTGGAGGGGCCGCCTCTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGTACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGAAt                                                                                                                                                                   <  2:464603/231‑1 (MQ=255)
                                                                                    ggTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGAATccc                                                                                                                                                                <  1:198753/250‑1 (MQ=255)
                                                                                    ggTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCCCAGGGGAAGGCCCTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCACAAAGAGGCCGCAGACCATTATTTTGCGGGGGGTATCGAATTCCt                                                                                                                                                               >  1:565388/1‑251 (MQ=255)
                                                                                     gTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGAATCCCTg                                                                                                                                                              <  1:129800/251‑1 (MQ=255)
                                                                                             gCCGCACACTATGAGCAACATGCATATCTACAGCGCCAGAGTGCGGACGCCTTACTGGCAATGTTTCCACGGCGTAAATCCCCCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACTTCACGCGCAGGGGACGGCCTTAGATCTTTCGCGGCCAATGCTTGTTCAGGCACGCCAGAATGATGCCTCAGTCCATTATCGGGCGGGAGATATCGAATCCCTGCCGTTAGc                                                                                                                                                      <  1:73590/251‑1 (MQ=255)
                                                                                                gCACACTATGAGCAACATGCAGATCTACAGCGCCAGCGAGCTGCCGCCTTACTGGCAATGCTTCCACAGCGTAAATACCCCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGGAAGTCACGCGGAGGGGAAGGCCCTAGATCTTTCGACGCCCATGCTCGTTCCGGGACGCCAGAAGGATGCCGCAGGCCAATATCTGGCGGGGGAGATCGGATCCCTGccctttgagac                                                                                                                                                   >  1:251995/1‑242 (MQ=255)
                                                                                                          ggaaacttGCAGATCTACGGCGCCGGGGTGCTGCCCCCTTCGTGGAAATGTTTCCCCAGCGAAAATACCCCCACGTCCGGGCCGGGGTTTGTGGCCCGGGCTGGAGGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTc                                                                                                                                         <  2:193806/244‑1 (MQ=255)
                                                                                                           aaaacatgcatacctacgcccccAGAGGGCGGACGCCTTACTGGCAATTCTTCCACAGCTTAAATACCCCCACGTACGGGAGGCGGGTTGGGGCCGGGGCTGAATGGGCCGCCACGGGCGGGAACGTCACGCGCAGGGGACGTCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGAGGCCGCAGACCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCg                                                                                                                                        <  2:396691/230‑1 (MQ=255)
                                                                                                              acaggCAGATCGACAGCGCAAGAGTGCTGCCGCCTTACTGGCAGTGCTCCCACGGCGAAAACCCCCCCCCGTCGGGGACGCGGGTGGTGGACCGGGCTGGAGGGGCCGCCACGGGCGGGAACGTCACGCGCAGGTGACGGCCTGAGATCTCGCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCGTTATCTGGCGGGAGATATGAAATCCCTGCCGTTAGCGACTGCGACGTTCGATc                                                                                                                                     <  2:608975/247‑1 (MQ=255)
                                                                                                              aCATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTACACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGGCACGCGCCGGGGGCGGGCTTAGAACTCTCGCCCCCCATGCTTTTTCCGGCACGCCCCAAGGGAGCCGCCGACCCCTATCTGGCGGGGGGTATCGCATCCccgccggttgacactgacacgttcgatt                                                                                                                                     >  2:49090/1‑229 (MQ=255)
                                                                                                               cATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGGAGGGGAAGGCCCTAGGTCCCTCGCCGCCCATGCTTGTTCAGGCACCCCAGAAGGAAGCCGCCGACCATTATTTGTCGGGGGGTATCGAATCGCTGCCGGTAGCGAccgctaaggtcggtct                                                                                                                                    >  1:94414/1‑236 (MQ=255)
                                                                                                                   cAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCCCGCGCAGGGGACGGCCTTAGATCTCTCCCCCCCAATGCTTGTTCAGGCCCGCCCGACGGATGCCGCAGACCATTATCTTGCGGGGGGTATCGAATCCCTGCCGGTAGCGACTGCGACGTTTGATCTTGcc                                                                                                                                >  1:341668/1‑250 (MQ=255)
                                                                                                                     gATCTACAGCGACAGAGTGCTGACGCATTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGGACCTCCAGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCCATGCTTGTTCAGGCACGCCCGAAGGATGCCGCAGACCATTATCTTGCGGGAGATATAGAAACCCTGGCGTTAGCGACTGCGACGCTCGATCTTGCATg                                                                                                                              >  1:517851/1‑251 (MQ=255)
                                                                                                                                    agTGCTGACGCCTTACTGGAAATGCTTCCACAGCGTAAATACACCCCCGTCCGGGACGCGGGTTGTGGCCCGGGCTGGAGGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCgcagt                                                                                                               <  1:187150/251‑1 (MQ=255)
                                                                                                                                              ccTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACCCGCCGGTGACGGCCTTAGCTCTCTCGCCGCCCATGCTTGTTCAGGCACGCCAGAAGGATGCCGCCGACCCTTATCTGGCGGGAGATATCGAATCCCTGCCGTTTGCGGCTGCGGCGGTCGATCTTTCATGGGGCCATCTCGCCGTGCAGTGGTGc                                                                                                     >  1:448923/1‑251 (MQ=255)
                                                                                                                                                              cccACAGCGTAAATACCCCCACGTACGGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCGGTAATTTATCCACgg                                                                                     <  1:122075/250‑1 (MQ=255)
                                                                                                                                                                                ccACGTACGGGACGCGGGTTGTGGCCCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCGGTAATTTATCCACGGCACTCCGCGAGCTGTATc                                                                   <  1:340164/251‑1 (MQ=255)
                                                                                                                                                                                   cGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGGGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCCGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCCTTGAGCAATCTCGCAGTGCAGGGGTGCGGGAATTTATCCCCCGCACCCCGCGAGCTGGATCggg                                                                >  1:160872/1‑251 (MQ=255)
                                                                                                                                                                                       gtGGACGGGGGTGGTGGACCTGGTGGGATGATCCGCCACTGGCGGGAACGTCACGCGCAGGGGACGGCGTTAGATCTCGCGCCGCCAATGCTTGTTCAGGCACGCAAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCGGTAATTTATCCACGGCACTCCGCGAGCTGTATCGggtggt                                                            <  1:197782/250‑1 (MQ=255)
                                                                                                                                                                                         ggccgcgGGTTGTGGCCCGGGCTGGATGGGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCCCGCCAGAAGGATCCCGCAGACCATTTTCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCGGTAATTTATCCACGGCACTCCGCGAGCTGTATCGggtggt                                                            <  1:449117/246‑1 (MQ=255)
                                                                                                                                                                                         ggACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCCCGCGCCGGGGAAGGCCCTAGGTCTCTCGCCGCCCATGCTTTTTTCGGCACCCCAGAAGGATGCCGCAGACCATTATCTTGCGGGGGATATCTAAACACTGCCGTTTGCGACCGCGGCGTTCGGACTTGCATGGAGCAAACTCCCAGTGCCGTGGGGGGGGAATTTTTTCCCCGCACTCCCCGAGCTGTATTGGggggtgg                                                          >  1:578961/1‑250 (MQ=255)
                                                                                                                                                                                            cgcgGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGGGGCGGCCTTAGATCTCTCGCCGGCAATGCTTGTTCCGGCACGCCCGAAGGGTGTCGCCGACCCTTATCTGGGGGGGGAGATAGGATCCCTGCCCGTAGCGGCCTCGACGTTCGATCTTGCCTGGAGCAATCTCGCAGTGCGGGGGGGCGGGAATTTTTCCAAGGCACCCCGCCTGCTGTATCggggggggggcc                                                       >  2:522965/1‑242 (MQ=255)
                                                                                                                                                                                             gcgGGTTGTGGACCTGGCGGGAGGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCGCAGTGCGGTGGTGCGGTAATTTATCCACGGCACTCCGCGAGCTGTATCGGGTGGTGCGccc                                                     <  1:230788/251‑1 (MQ=255)
                                                                                                                                                                                                  gtgggggcccgggtggaagggcccgccACGGGGGGGACCGCCACGCGCGGGGGACGGCCTAAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGCCCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCGGTAATTTATCCACGGCACTCCGCGAGCTGTATCGGGTGGTGCGCCCCaaa                                                 <  2:400703/229‑1 (MQ=255)
                                                                                                                                                                                                          ccTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGGGACGGGCCTAGATCTCTCGCCGCCCATGCTTTTTCAGGCACGCCAGAAGGGGGCCGCCGACCATTTTCTGGCGGGGGAAATCGACTCCCCGCCGGTCGCGGCTGCGACGTTCGGTCTTTCGTGGAGGAATCTCCCCGTGCCGTGGGGCGGGGATTTTTTCACCGCACTCCCCGAGCGGTATCGGGTGGGGCGGCCCCAAGGGCgggg                                          >  2:69565/1‑248 (MQ=255)
                                                                                                                                                                                                               gtgaaTGAGCCGCCACTGGCGGGAACGTCACGCGCGGGTGCCGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCGGTAATTTATCCACGGCACTCCGCGAGCTGTATCGGGTGGTGCGCCCCAAAGGCGTGGTCGCGt                                    <  1:374009/247‑1 (MQ=255)
                                                                                                                                                                                                                       gccgccACGGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCGGTAATTTATCCACGGCACTCCGCGAGCTGTATCGGGTGGTGCGCCCCAAAGGCGTGGTCGCGTTTACCACg                            <  1:284955/251‑1 (MQ=255)
                                                                                                                                                                                                                          cccACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCGGTAATTTATCCACGGCACTCCGCGAGCTGTATCGGGTGGTGCGCCCCAAAGGCGTGGTCGCGTTTACCACGCTg                         <  1:168128/250‑1 (MQ=255)
                                                                                                                                                                                                                              cTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCCGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATCTCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCGGTAATTTATCCACGGCACTCCGCGAGCTGTATCGGGTGGGGCGCCCCAAAGGCGTGGGCGCGTTTACCACGCTGGTGc                     >  2:283614/1‑251 (MQ=255)
                                                                                                                                                                                                                                ggCGGGCACGACACGCGCAGGTGACGGCTTTAGATCTCTCGCCGCCAATGCTTGTGCAGGCACGCCAGAAGGATGCCGAAGACCATTATCGGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCGGTAATTTATCCACGGCACTCCGCGAGCTGTATCGGGTGGTGCGCCCCAAAGGCGTGGTCGCGTTTACCACGCTGGTGCAg                   <  1:439807/251‑1 (MQ=255)
                                                                                                                                                                                                                                ggCGGGAACGTCACGCGCAGGTGACGGCCTGAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCGGTAATTTATCCACGGCGCTCCGCGAGCTGTATCGGGTGGTGCGCCCCAAAGGCGTGGTCGCGTTTACCACGCTGGTGCAg                   <  1:107201/251‑1 (MQ=255)
                                                                                                                                                                                                                                ggCGGGAAAGTCACGCGCAGGCGACGGGCTTAGATCTATGGCCGCCAATGCTTGTCAAGGCACGACCGAAGGATGCCGCCGACCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGCAGTTCGATCTTGCATCGAGCAATCTCGCAGTGCAGTGGCGCGGTAATTTATCCACGGCACTCCGCGAGCTGTATCGGGTGGTGTGCCCCCAAGGCGTGGGCGCGttctcccccctggtgccg                   >  2:403167/1‑235 (MQ=255)
                                                                                                                                                                                                                                     gaacgcacggcgcGGGTGCCGGCATTAGACCCCCCGCCGCCAAGGCTTGTACAGGCCCGCCAGAAGGATGCCGCAGCCCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCGGTAATTTATCCACGGCACTCCGCGAGCTGTATCGGGTGGTGCGCCCCAAAGGCGTGGTCGCGTTTACCACGCTGGTGCAGGGATc              <  2:247528/242‑1 (MQ=255)
                                                                                                                                                                                                                                     gAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGCCCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCGGTAATTTATCCACGGCACTCCGCGAGCTGTATCGGGTGGTGCGCCCCAAAGGCGTGGTCGCGTTTACCACGCTGGTGCAGGGATc              <  1:522965/251‑1 (MQ=255)
                                                                                                                                                                                                                                      aaCGTCACGCGCCGGTGACGGCCTAAGATCTCTCGCCGCCAATGATTGTTCAGGCACGCCAGAAGGATGCCGCCGACCATTATCTGGCGGGCGATATCGCATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCGGTAATTTCTCCACGGAACTCCGCGAGCTGTCTCGGGTGGTGCGCCCCAAAGGCGTGGTCGGGTTTGCCACGCTGGTGCAGGGAt               >  2:39440/1‑249 (MQ=255)
                                                                                                                                                                                                                                        cGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAAGGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCGGTAATTTATCCACGGCACTCCGCGAGCTGTATCGGGTGGTGCGCCCCAAAGGCGTGGTCGCGTTTACCACGCTGGTGCAGGGATCGtt           <  1:373405/251‑1 (MQ=255)
                                                                                                                                                                                                                                                 cgggTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGAACCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCGGTAATTTATCCACGGCACTCCGCGAGCTGTATCGGGTGGTGCGCCCCAAAGGCGTGGTCGCGTTTACCACGCTGGTGCAGGGATCGTTACCCGAACt  <  1:5981/249‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                        
CTGCGCATGAAATGCAGGATATCGACCGTCTGCTGGAGGTGCTGCATGGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCGGTAATTTATCCACGGCACTCCGCGAGCTGTATCGGGTGGTGCGCCCCAAAGGCGTGGTCGCGTTTACCACGCTGGTGCAGGGATCGTTACCCGAACT  >  NC_000913/811477‑811967

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