Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913811,7210TG35.1% 88.5 / 11.1 74V67G (GTG→GGG) bioCmalonyl‑acyl carrier protein methyltransferase
Reads supporting (aligned to +/- strand):  ref base T (8/40);  new base G (24/2);  total (32/42)
Fisher's exact test for biased strand distribution p-value = 1.39e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GCGCATGAAATGCAGGATATCGACCGTCTGCTGGAGGTGCTGCATGGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCGGTAATTTATCCACGGCACTCCGCGAGCTGTATCGGGTGGTGCGCCCCAAAGGCGTGGTCGCGTTTACCACGCTGGTGCA  >  NC_000913/811479‑811949
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                      
gcgcATGAAAGGCAGGATATCGACCGTCTGCTGGAGTTGCTGCATGGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCt                                                                                                                                                                                                                                <  2:401105/251‑1 (MQ=255)
    aTGAAATGCAGGATATCGACCGTCTGCTGGAGGTGCTGCATGGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGGGACGGCCtt                                                                                                                                                                                                                               >  1:3190/1‑248 (MQ=255)
    aTGAAATGCAGGATATCGACCGTCTGCTGGAGGTGCTGCATGGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCGGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGa                                                                                                                                                                                                                            <  1:326496/251‑1 (MQ=255)
           gcgggattgcGACCGTCTGTTGGAGGTGCTGCATGGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTTTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGAtctctc                                                                                                                                                                                                                      <  2:229119/241‑1 (MQ=255)
           gCAGGATATCGACCGTCTGCTGGAGGTGCTGCATGGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTcg                                                                                                                                                                                                                     <  1:161306/251‑1 (MQ=255)
           gCAGGATACCGACCGTTAGTTGGATGTGCGTCATGGCAACGGGTAATAACAAAGCCATTGCAGCGGCATTTGGTCGTGCAGCCGCACACTAGGATAAACAGGCAGAGCTACAGCGCCAGAGGGCTGCCGCCTTACGTGCAATGCTTCCACAGCGTAAATACCCCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGTAACGTCACGCGCAGGTGACGTCCTTAGATCTCTcg                                                                                                                                                                                                                     <  2:310876/251‑1 (MQ=255)
               gATATCGACCGTCTGCTGGAGGTGCTGCATGGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGGACGGCCCGCGCCGGGGACGGACTTAGATCTCTcgccgc                                                                                                                                                                                                                 >  2:104428/1‑251 (MQ=255)
                atatCGACCGTCTGCTGGAGGTGCTGCATGGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGCACGTCCCGCGCCGGGGACGGCCTTAGATCTCTCgccgcc                                                                                                                                                                                                                >  1:451895/1‑251 (MQ=255)
                  atcgacggtctgtgtgAGGGGCTGCATGGAAACGGTAAATAAACAAGCCATTGCAGCGGCATTTGGGCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGCCGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCa                                                                                                                                                                                                               <  1:289828/235‑1 (MQ=255)
                     gACCGTCTGCTGGAGGTGCTGCATGGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGGACGGCACGGGGAGGGGACGGCCTTTGATCTCTCcccgcccaagc                                                                                                                                                                                                           >  2:305013/1‑246 (MQ=255)
                          tCTGCTGGAGGTGCTGCATGGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCCATGCCTGtt                                                                                                                                                                                                      >  1:359593/1‑251 (MQ=255)
                           ctgctgGAGGTGCTGCATGGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCCGGGGGCGGCCTTAGATCTCTCGCCGCCAAAGCTTGTTc                                                                                                                                                                                                     >  2:425050/1‑251 (MQ=255)
                             gctgGAGGTGCTGCATGGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATGCACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGATCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCa                                                                                                                                                                                                    <  1:314724/250‑1 (MQ=255)
                                ggAGGTGCTGCATGGCACCGGTTAATAAACAAGCCATGGCAGCGGCAGTTGGTCGGGCAGCCGCACACGATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCGCGCCGCCAATGCTTGTTCAgg                                                                                                                                                                                                  <  1:167189/249‑1 (MQ=255)
                                 gAGGTGCTGCATGGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACCCCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCAc                                                                                                                                                                                               <  1:180590/251‑1 (MQ=255)
                                 gAGGTGCTGCATGGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGGCCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCAc                                                                                                                                                                                               <  1:95650/251‑1 (MQ=255)
                                 gAGGGGCTGCATGGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCAc                                                                                                                                                                                               <  1:352276/251‑1 (MQ=255)
                                  gggTGCTGCATGGCAACTGTTAAAAAACAAGCCATTGCATCGTCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGCCGCCTTACTGGCAATGCTTCCTCAGCTTAAATACACCCACGTACTGGACGCGGGTGGTGGACTTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCATGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACg                                                                                                                                                                                              <  1:408849/250‑1 (MQ=255)
                                   ggggcggcATGGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGGGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACCCCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGc                                                                                                                                                                                             <  2:36823/245‑1 (MQ=255)
                                     tgctgcATGGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGGTCCGGGCCGGGGACGGGCTTATATATCTCGCCGCCAATGCTTGTTCAGGCACGccc                                                                                                                                                                                           >  1:336184/1‑250 (MQ=255)
                                          gatggcactgttaaataaACAAGCCATTGCAGCGGCATTTGGCCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCGGACGCCTTACGGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAgg                                                                                                                                                                                      <  2:409410/238‑1 (MQ=255)
                                           aTGGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGGCACGCGCAGGGGACGGCCTTAGGTCTCTCGCCGCCAATGCTTGTTTCGGCACGCCAGAAgg                                                                                                                                                                                      >  1:411395/1‑250 (MQ=255)
                                              gcaccggtaaataacaaaGCCCTTGCAGCGGCAGTTGGTGGGGCAGCCGCACACTAGGACCAACAGGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTCCCCCAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGc                                                                                                                                                                                  <  2:94477/235‑1 (MQ=255)
                                                    gTTAATAACCACGCCATTGCAGCGGCAGTTGGTCGGGCAGCCGCCCACTATGAGCAACAGGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGa                                                                                                                                                                            <  2:119218/251‑1 (MQ=255)
                                                     ttAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGGCAAGCGCCGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCCGCAGGATGGCGCAGAc                                                                                                                                                                           >  1:308098/1‑251 (MQ=255)
                                                        ataaACAAGACATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGGCACGCGCAGGGGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTTAGGCACGCCCGCAGGATGCCGccaacca                                                                                                                                                                         >  1:100864/1‑244 (MQ=255)
                                                             aaaGCCATTGCGGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGCCGCCTTACGGGCAATGCTTCCACAGCGTAAATACACCCCCGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCt                                                                                                                                                                   <  2:251905/250‑1 (MQ=255)
                                                               aGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCCCGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAAAGCTTGGTCAGGCACGCCAGAAGGATGCCGCCGACCATTATCTgg                                                                                                                                                                 >  1:379219/1‑251 (MQ=255)
                                                                gCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCCGGGGACGGCCTTAGATCTCTCGCCGCCCATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGc                                                                                                                                                                >  1:119049/1‑251 (MQ=255)
                                                                  cATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTCCTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGCGGCGGGTTGTGGCCCTGTCTGGATGAGCCGCCACTGGCGGGAACGTCACGAGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCgg                                                                                                                                                              <  1:196804/251‑1 (MQ=255)
                                                                   aTTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCTTTACGGGCAATGCTTCCACAGCGTAAATACCCCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCggg                                                                                                                                                             <  1:413222/251‑1 (MQ=255)
                                                                   aTTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCAAGCGCAGGGGACGGCCTTAGATCTCTCGCCGCCAATGATTGTTCAGGCACACCCGAAGGGAGCCGCAGACCAATATTCGGCggg                                                                                                                                                             >  2:32989/1‑251 (MQ=255)
                                                                      gcatcggtaTTTGGCCGGGCAGCCGCTCACTATGAGCAACATGCAGATCTACAGTGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGCAAATACACCCACGTACTGGACGCGGGTTGTGGACCGGGCGGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGgaga                                                                                                                                                          <  1:500809/243‑1 (MQ=255)
                                                                         gCGGGATTTGGTCGGGCAGCCTCACACTATGAGCAACATGCAGATCTACAGCGCCATAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCGCGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGGCGGTCTTAGATCTCTCGCCGCCAATGCTTTTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGatat                                                                                                                                                       <  1:331285/251‑1 (MQ=255)
                                                                         ccGGCGTTTGGTCGGGCTGCCCCACCCGACGAGAAAAATGCAGATCTCCGGCGCGGGAGGTCTGCCGCCTTACTGGCAATGCGTCCACAGGTGAAAGGCACCCACGTGGTGGGGGCGTTTTGGGGCCCTGGCTGGAGGTGCCGCCACGGGCGGGGACGTAACGCGCAGGTGACGGACTTAGTACTGTCGCTTCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGTCCAGTTTCTTGTGGGAGatat                                                                                                                                                       <  2:293906/250‑1 (MQ=255)
                                                                            gCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCCCGCGCAGGGGACGGCCTTAGATCCTTCGCCGCCCATGCTTGTTCCGGCCCGCCAGAAGGATGCCGCCGACCATTATCTGGCggggggtattg                                                                                                                                                     >  1:403565/1‑242 (MQ=255)
                                                                                   gTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGAGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCCCGCGCAGGGGACGGCCTTAGATCTCTCGCCGCCCATGCTTGTTCAGGCCCGCCCAACGGATGCCCCAGACCATTATAGGGCGGGGGAAATCGAATCCCTg                                                                                                                                            >  1:267683/1‑251 (MQ=255)
                                                                                    ccGGGCAGCCGCCCATTAGGAGCAACATGCAGATCTACAGCGCCAGAGTGCGTACGCCTTACTGGCAATGCTTCCACAGCGTAAATACCCCCACGTACTGGACGCGGGTTGGGGACCGGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTGAGATCTCGCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGAATCCCTGc                                                                                                                                           <  2:451895/250‑1 (MQ=255)
                                                                                      gggCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCCCGCGCCGGGGACGGCCTTAGATCTCTCGCCGCCCATGCTTGTTCAGGCACGCCCCAAGGATGCCGCCGACCATTATCTGGCGGGGGATATCGAATCCccgccg                                                                                                                                         >  1:309018/1‑251 (MQ=255)
                                                                                        gCAGCCGCCCACTAGGAGCAACATGCAGATCAACGCGGCCGGGGTCCGGACGCTTTCCTGGAATGGCTTCACCAGCGAAAATACCCCCACGTACTGGGCGCGTGGTTGGGCCCGGGTTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGGGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGAGGCCGCAGCCCATTATCTGGCGGGAGATATCGAATCCCTGCCGtt                                                                                                                                       <  2:195267/251‑1 (MQ=255)
                                                                                             cGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGGACGTCACGCGCAGGGGACGGCCCTAGATCTCCCGCCGCCCATGCTTGTTCAGGGACGCCCGCAGGATGGCGCAGACCCTTATCTGGCggggggtatcaactcccctccgttggcg                                                                                                                                    >  2:50331/1‑225 (MQ=255)
                                                                                               cacccgaagAGCAACATGCAGACCTACAGCGCCAGATGGCTGCCGCCTTACTGGCAATGCTTCCCCAGCGAAAATACCCCCCCGTACTGGACGCGGGTTGTGGCCCGGGCTGGATGGGCCGCCACTGGCGGGAACGTCACTCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACt                                                                                                                                <  1:270200/243‑1 (MQ=255)
                                                                                               cacaCTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGGACGTCACGCGCCGGGGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGGTGCCGCCGACCATTATCTGGCGGGAGATAACGAATCCCTGcccgttgcgccc                                                                                                                                >  2:123878/1‑241 (MQ=255)
                                                                                                 cacTATGAGCAAAATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATCCCCCCCCGTCCTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGc                                                                                                                              <  1:415641/251‑1 (MQ=255)
                                                                                                  acTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCCATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGGACGTCACGCGCAGGTGACCGCCCTAGGTCTCTCGCCGCCAATGCTTGTTTAGGCACGCCAGAAGGATGCCGCCGACCATTATCTGGCGGGAGAAATCGAATCCCTGCCGGTAGCGACTGCg                                                                                                                             >  1:147095/1‑251 (MQ=255)
                                                                                                   cTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAAAGTCACGCGCAGGGGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACCCCCGAAGAGTGCCGCAGACCAATATCTGGCGGGAGGTATCGCATCCCTGCCGGTAGCGGCTGCgg                                                                                                                            >  1:85950/1‑250 (MQ=255)
                                                                                                     gtaagcaacttgaaaatttacggcgcCTGAGTGCTTCCGCGTAACTGGCATTGCTTCCAAAGGGAAAAAACCCCCACGTACGGAAGGGGGGTTGGGGCCTTGGTTGGATGGGCCGCCTCGTGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCCCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGAGGCCGCAGACCATTACCGGGCGGGAAATAGCGAACCCCTGCCGTTAGCGCCTGCGACg                                                                                                                          <  2:1739/229‑1 (MQ=255)
                                                                                                     aTGAGCAACATGCAGATCTACACCGCCAGAGTGCTGACGCATTACTGGCAATGCTGCCACAGCGTAAATCCACCCACGTACTGGACGCGGGTTGGGGACCTGGCTGGATGGGCCGCCACTGGCGGGAACGTCACGCGCAGTTGACGGCCTTAGACCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACg                                                                                                                          <  1:133610/251‑1 (MQ=255)
                                                                                                            aCATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCCCCCCAGAAGGATGCCGCAGACCCTTATCTGGCGGGGGATATTGACTCCCTGCCGTTAGCGACTGCGCCGTTCGATc                                                                                                                  >  1:310804/1‑251 (MQ=255)
                                                                                                                 cAGATCTACAGCGCAAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACCCCCACGTACGGGACGCGGTTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCa                                                                                                             <  1:414749/251‑1 (MQ=255)
                                                                                                                   gcgctccggcgcCAGTGTGCTGACGCCTTACTGGGAATGCTCCCACGGCTTAAATCCCCCCACGTCCTGGACGCGGGTGGTGGACCGGGCGGGATGGTCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGAGCTCTCGCTGCCAAGGCTGGTTCAGGCACGCCAGAAGGATGCCGCAGACCTTTATCTGGCGGTAGATATCGAATCCCTGCCGGTAGCGACTGCGACGTTCGATCTTGCATg                                                                                                           <  2:463120/243‑1 (MQ=255)
                                                                                                                            cgcCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAAAGTCACGCGGAGGGGACGGCCTTAGATCTCTCGCCGCCCAACATTGTTCAGGCACGCCAGAAGGATGCCCCCGACCATTATCTGGCGGGAGATATCGAATCCCCGCCGTTAGCGGCTGCGCCCTTCGGTCTTGccaggggccatcc                                                                                                  >  1:77923/1‑239 (MQ=255)
                                                                                                                                 agagTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGGGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCgc                                                                                               >  1:138916/1‑250 (MQ=255)
                                                                                                                                          gacgcttttcgtgataggcgtCCACAGCGTAAATCCCCCCACGGACGGGACGCGGTTGGGGCCCGTGGCTGGATGGGCCGCCACGGGCGGGAACGTCACGCTCAGGTGACGTCCTTAGATCACTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCGGGCGGGAGATATCGAATCCCTTCCGTTAGCGACTTCGACGTGCGATCTTGCATGTAGCAATCTCGcagtgcagtg                                                                                      <  2:300717/231‑1 (MQ=255)
                                                                                                                                            gCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGGACGTCACGCGCAGGGGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGAATCCCCGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCGCAGTGCAggggtg                                                                                   >  1:313734/1‑251 (MQ=255)
                                                                                                                                                 aCTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAAAGCCCCGCGCCGGTGACGGCCCTAGAGCTCTCGCCGCCCATGCTTGGTCCGGCACCCCAGAAAGATGCCCCCGACCCATATCTTGCGGGGGATATCGAATCCCTGCCGTTTACGTCTGCGCCGTTCGATCTTGCGTTGAGCAAtcttgccgtggaggggaggggtt                                                                              >  1:263358/1‑231 (MQ=255)
                                                                                                                                                    gtcAATGCTTCCACCGCGTAAATACACCCACGTCCTGGTCGCGGTTTGTGGCCCGGGCTGGATGAGCCGCCACTGGCGGGAACGGCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCGGTAAtt                                                                           <  1:237847/249‑1 (MQ=255)
                                                                                                                                                      cAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCAAGCGCAGGGGACGGCCTTAGATTTCTCGCCGCCCATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATTTTGCGGGGAAATTCGAAACCCCGCCCTTAGCGACTGCGGAGTTCGAGCTTTCATTGAGCAATATCGCAGGGCAGTGGGGCGGGAAtttt                                                                         >  1:446422/1‑250 (MQ=255)
                                                                                                                                                                      agaaaatacccccctgtacgggccgcgGTTTGTGGCCCGGGCTGGATGCGCCGCCTCTGGGGGGAACGTCACGCGCGGGGGCCGGCCTTAGATCTCTCGCCGCCAATGCTGGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCGGTAATTTATCCACGGCACTCcg                                                           <  2:263358/228‑1 (MQ=255)
                                                                                                                                                                          aaTACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGGGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCCGAAGGATGCCGCAGACCATTATCTGGCGGGGGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCGGGAATTTATCCACGGGACTCcgcgcg                                                       >  1:417992/1‑249 (MQ=255)
                                                                                                                                                                            tACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGGGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCGCAGTGCAGTGGGGCGGTAATTTTTCCACCGCACTCCGCGAGCt                                                     >  1:173966/1‑251 (MQ=255)
                                                                                                                                                                            aacccccACTTACTGGACGCGGGTTGTGGCCCTGGCGGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCGGTAATTTATCCACGGCACTCCGCGAGCt                                                     <  1:14520/247‑1 (MQ=255)
                                                                                                                                                                              cacCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGGACGTCACGCGCAGGGGAAGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCCTTATCTGGCGGGAGATATCCAATCCCTGCCGTTAGCGACTGCGGCCTTCGATCTTGCCTGGAGCAATTTCGCAGGGCAGTTGGGCGGTAATTTATCCACGGCACTCCGCGAGCTGt                                                   >  1:134585/1‑251 (MQ=255)
                                                                                                                                                                                cctcgttccTGGTCGCGTGTTGTGGATCGGGCTGGATTGCCCGCCGCTGGCGGAAACGCCACGCGCAGGTGACGGCCTTAGATCCCGCGCCGCCAATGCTTGTTAAGGCACGCAAGAAGGATACCGCAGACCATTTTCGGTCGGGAGATATCGAATCCCTGCCGGTAGCGACTGCGACGTTCGAGCTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCGGTAATTTATCCGCGGCACTCCGCGAGCTGTAt                                                 <  2:22363/243‑1 (MQ=255)
                                                                                                                                                                                    cGTACTGGACGCGGGTTGTGGACCTGGCTGGAGGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCGGTAATTTATCCACGGCACTCCGCGAGCTGTATCggg                                             <  1:346360/251‑1 (MQ=255)
                                                                                                                                                                                        cTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGGTCTCTCGCCGCCCATGCTTGTTCAAGGACGCCCGAAGGATGCCGCAGACCATTATTTGGCGGGGGATATCGAATCCCTGCCGTTAGCGACTGCGGCGTTTCGTCTTTCATTGAGCAATCTCTCTGTGCAGGGGGGGGGGAAATTTTCCCCGGGCCTCCGCGAGCTGTATCggggggg                                         >  2:93936/1‑247 (MQ=255)
                                                                                                                                                                                          ggccgggggTTGGGGACCGGGGTGGTGGAGCCGCCACTGGCGGGAACTTCACGCGCAGGTGACGGCCTTAGATCCCCCGCCGCCAATGCTGGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCGGGCGGGAGATATCGATTCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCGGTAATTTATCCACGGCACTCCGCGAGCTGTATCGggtggt                                         <  1:84715/243‑1 (MQ=255)
                                                                                                                                                                                            aCGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAAAGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGGAATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGGTCTTGCATGGAGCAATCTCGCCGTGCCGTGGTGCGGTAATTTATCCCCGGCACCCCGCGGGCTGTATCGggtggggggc                                     >  1:381909/1‑246 (MQ=255)
                                                                                                                                                                                             cgggggTTGTGGCCCGGGGTGGAGGAGCCGCCCCTGGCGGGACCGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCGGTAATTTATCCACGGCACTCCGCGAGCTGTATCGGGTGGTGCGcc                                    <  2:379219/248‑1 (MQ=255)
                                                                                                                                                                                             cgggggTTGTGGACCTGGTTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCGGTAATTTATCCACGGCACTCCGCGAGCTGTATCGGGTGGTGCGcc                                    <  1:139526/248‑1 (MQ=255)
                                                                                                                                                                                             cgcgGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCGGTAATTTATCCACGGCACTCCGCGAGCTGTATCGGGTGGTGCGcc                                    <  1:104428/251‑1 (MQ=255)
                                                                                                                                                                                                gggTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTGAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCGGTAATTTATCCACGGCACTCCGCGAGCTGTATCGGGTGGTGCGCCCCa                                 <  1:123878/251‑1 (MQ=255)
                                                                                                                                                                                                       tgGACCTGGCTGGATGAGCCGCCACTGGCGGCAAAGTCACGCGCAGGGGACGGCCTTAGATCTCTCGCCGCCCATGCTTGTTCAGGCACGCCCGAAGGATGCCGGAGGCCATTATCTGGCGGGAGATATCGAATCCCCGCCGTTAGGGGCTGCGGCGTTCGATCTTGCATGGAGACATCTCGCAGTGCAGTGGGGCGGTAATTTATTCACGGCACGCCGCGAGCGGGATAGGGGGGGGCGCCCCAAAGGCg                           >  1:168452/1‑251 (MQ=255)
                                                                                                                                                                                                                                  ggCGGGACCGTCACGCGGGGGTGAGGGCCTTAGATCTCCCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCTCAGACCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCGGTAATTTATCCACGGCACTCCGCGAGCTGTATCGGGTGGTGCGCCCCAAAGGCGTGGTCGCGTTTACCACGCTGGTGCa  <  2:364562/250‑1 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                      
GCGCATGAAATGCAGGATATCGACCGTCTGCTGGAGGTGCTGCATGGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCGGTAATTTATCCACGGCACTCCGCGAGCTGTATCGGGTGGTGCGCCCCAAAGGCGTGGTCGCGTTTACCACGCTGGTGCA  >  NC_000913/811479‑811949

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