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

TGCGCATGAAATGCAGGATATCGACCGTCTGCTGGAGGTGCTGCATGGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCGGTAATTTATCCACGGCACTCCGCGAGCTGTATCGGGTGGTGCGCCCCAAAGGCGTGGTCGCGTTTACCACGCTGGTGCAGGGATCGTTAC  >  NC_000913/811478‑811960
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                          
ggcgcATGAAATGCAGGATATCGACCGTCTGCTGGAGGTGCTGCAGGGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGcc                                                                                                                                                                                                                                                    <  2:218836/250‑1 (MQ=255)
              aGGATATCGACCGTCTGCTGGAGGTGCTGCATGGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCCCGCGCAGGGGACGGCCATTGATCTCTcgcc                                                                                                                                                                                                                                    >  1:354344/1‑251 (MQ=255)
               ggATATCGACCGTCTGCTGGAGGTGCTGCAGGGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTcgccg                                                                                                                                                                                                                                   <  2:444372/251‑1 (MQ=255)
                    ccacccGTCTGCGGGAGGTGCGGCGTGGCATGGTTTAAAAACAAAGCCATTGCGGCGGAATTGGGTGGGGCACCCGCCCACTATGAGCAACATCCAGATCTACAGCGCCGGTGTGCTGCCGCCTTACGGGCAATCCTTCCACGGCGTAAATCCCCCCACGTACGTGGCGCGGGTTGTGGCCCTGGTTGGATGCCCCGCCACTGGGGGGAACGTCACGCGCAGGTGACGGCCTTAGTTCTCGCGCCGCCAAt                                                                                                                                                                                                                             <  2:434357/247‑1 (MQ=255)
                       aCCGTCTGCTGGAGGTGCTGCATGGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGGTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGGACGTCACGCGCCGGTGCGGGCCTTAAATTTCTcgccccccaagct                                                                                                                                                                                                                           >  2:304376/1‑242 (MQ=255)
                         cGTCTGCTGGAGGTGCTGCATGGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGGACGTCACGCGCAGGGGACGGCCTTAGATCTCTCGCCGCCAATGcctg                                                                                                                                                                                                                        >  1:111443/1‑248 (MQ=255)
                          gTCTGCTGGAGGTGCTGCATGGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCCCGCGCCGGGGACGGCCtttgagctcctcccccccatgcctgt                                                                                                                                                                                                                        >  2:344133/1‑227 (MQ=255)
                                 ggAGGTGCTGCATGGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAgg                                                                                                                                                                                                                 <  1:426397/249‑1 (MQ=255)
                                   aGGTGCTGCATGGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACCGCCCCCCCCGGGGAAGGCCTTAGATCCCTTGCCGCCCAAGCTTGTTCCGGCAcc                                                                                                                                                                                                            >  2:244251/1‑250 (MQ=255)
                                     gTGCTGCATGGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGGGACGGCCTTAGATCTCTCGCCGCCCAAGCTTGTTCAGGCACGcc                                                                                                                                                                                                          >  1:105303/1‑251 (MQ=255)
                                      cgcttcctgcccaCGGTTCATAACCAAGCCTTAGCATCGGCATTTGGCCCGGCAGCCGCACCTTACGAGCACCACGCCTGTTTACAGCGCCCGAGTGCTGCCTCCTACCTGGCAATGTTTCCACACCGTAAATACCGCCACGTACTTGTCGCGGGTTGTGGACCTGGCGGGATGAGCCGCCACTGGCGGGTTCGTCACGCGCAGGTGACCGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCa                                                                                                                                                                                                         <  2:119131/239‑1 (MQ=255)
                                         tgcATGGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGGACGGCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACCCCAGa                                                                                                                                                                                                       >  1:412570/1‑250 (MQ=255)
                                          gcATGGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTAGGGCAGCCGCACACTATGAGCAACATGCAGATCTAAAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCAACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGGCACGCGCAGGGGACGGCCGTAAATCCCTCGCCGCCAATGCTTGTTCAGGGACGcccga                                                                                                                                                                                                       >  1:452559/1‑246 (MQ=255)
                                          gcATGGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGGCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCGCTGGCGGGAACGTCACGGGCAGGTGACTGTCTGAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGa                                                                                                                                                                                                       <  1:201497/249‑1 (MQ=255)
                                          ccatgccaccgtttAATAACCAAGCCATTGCAGCGGCATTTGGTCGGCCAGCCGCACACTATGAGCAACATGCAGATCTCCAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGTACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAg                                                                                                                                                                                                     <  1:266676/239‑1 (MQ=255)
                                             gggCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCCCACTAGGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAgg                                                                                                                                                                                                   <  2:75263/248‑1 (MQ=255)
                                                         ataaACAAGCCATTGCAGCGGCATTTGGTCGGGAAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACGCCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGCACGGAACGCGCAGGGGCCGGCCACAGATTTCCTGCCGGCCACTCTTGTTCAGGCACGCCCGAAGGaatgatcaggaga                                                                                                                                                                                      >  1:167976/1‑238 (MQ=255)
                                                          aaaaCAAGCCATTGCGGCGGCATTTGGTCGGGCAGCCGCCCACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACCCCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCAtt                                                                                                                                                                                  <  2:105303/250‑1 (MQ=255)
                                                           aaaCAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGGACGTCCCGCGCCGGGGACGGCCTTTGAACTCTCCCCGCCAAAGCTTGTTTAGGCACGCCAGAAGGATGCCGCAGACCAttt                                                                                                                                                                                 >  1:233784/1‑250 (MQ=255)
                                                            aaCAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTAt                                                                                                                                                                                <  1:376294/251‑1 (MQ=255)
                                                               aaGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGGGACGGCCTTAGATCTCTCCCCCCCCATGCTTGGTCAGGCCCGCCCGAAGGATGCCGCCGAACATTTTCTg                                                                                                                                                                             >  1:390064/1‑251 (MQ=255)
                                                                agtccATGGCAGCGGCTTTGGGTCGGGCAGCCGCACACGATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGACTTTCGGGCAATGCTTCCACAGCGTAAATACCCCCTCGTACTGGACGCGGGTTGGGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGTATGCCGCAGACCATTATCTg                                                                                                                                                                             <  2:481387/248‑1 (MQ=255)
                                                                     aTTGCAACCGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCCGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGGACGTCACGAGCCGGTGAACGACTTAGATCTCTCGCCCCCCACGCCTGTTCCAGCGCGCCAGAAGCGAGCCGCAGACCCTTTTCTGGCggg                                                                                                                                                                        >  2:390437/1‑251 (MQ=255)
                                                                          aGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCCCAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAAAGGCACGCGCCGGGGACGGGCTTAGAACTCTCCCCGCCCAAGCTTTTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGgagct                                                                                                                                                                    >  1:40976/1‑248 (MQ=255)
                                                                              gCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGCGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCCCGCCCCGGGGACCGCCTTTGATCTCCCCCCCCCCACGCCTTTTCAGGGCAGCCAGAAAGATCACGCAGACCACTATCTgggggggtgtaaccg                                                                                                                                                               >  2:207104/1‑242 (MQ=255)
                                                                                  ttGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGAATc                                                                                                                                                            <  1:192373/250‑1 (MQ=255)
                                                                                   tGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGAATcc                                                                                                                                                           <  1:463102/250‑1 (MQ=255)
                                                                                       cGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGAATCCCTGc                                                                                                                                                       <  1:24157/250‑1 (MQ=255)
                                                                                                 cacaCTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCAAGCGCCGGGGACGGGCCTTGAAATCTCGCCGCCAATGCTTTTTCAGGCACGCCAGAAGGATCCCCCAGACCCTTATTTGGCGGGGGATATCGAATCCCtccccgtagcgcct                                                                                                                                            >  2:404119/1‑238 (MQ=255)
                                                                                                    acTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATg                                                                                                                                                                                                                            <  1:165082/173‑1 (MQ=255)
                                                                                                       aTGAGCAACATGAAGATCTACAGCGCCAGAGTGCTGACGCCTTTCTGGCAATGCTGCCCCAGCGTAAATACCCCCCCGTACTGGACGGGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCTTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGAGGCCGCAGACCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACg                                                                                                                                    <  1:245718/251‑1 (MQ=255)
                                                                                                        tGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGGCACGCGCAGGGGGCGGGCCTTGAGCTCTCGCCCCCAATGCTTGTTCCAGGACGCCCGCAGGGTGCCCCAGCCCACTTTTCGGCGGGGGATATCGGAACCCTGCCCTTAGCGACAGCgccgg                                                                                                                                     >  2:364925/1‑247 (MQ=255)
                                                                                                          aGCAACATGCAGCTCAACAGCGCCAGAGGGCTGACCCCTAACTGGCATTGCTTCCACAGCGAAAAGACCCCCACGTATGGGACGCGGGTTGGGGACCTGGCGGGAGGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTc                                                                                                                                 <  2:342126/251‑1 (MQ=255)
                                                                                                                  tccctccaccacgctcttccactccaCGCCTTCCTGCCAAGGCTTCCACAGCTAAAATTCCCCCCCGTTCGGGCCGCGGGTTGTGGCCCTGGCGGGATGAGCCGCCACGGGCGGGAACGTCACGCGCAGGGGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTg                                                                                                                          <  2:336136/225‑1 (MQ=255)
                                                                                                                    aGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACCCCCACGTACTGGACGCGGGTTGTGGCCCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCAt                                                                                                                       <  1:112192/251‑1 (MQ=255)
                                                                                                                     gATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGGGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGGGATATCGAAACCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATg                                                                                                                      >  1:101393/1‑251 (MQ=255)
                                                                                                                           cAGCGCCAGAGTGCGGACGCCTTACGGGCAATGCTCCCACAGCGTAAAACCACCCACGTACGGGACGCGGGTGGTGGGCCTGGCTGGAGGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGACCTCTCGCCGCCAATGCGTGTTCAGGCTCGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGAAGCCCGGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCaa                                                                                                                <  2:239462/251‑1 (MQ=255)
                                                                                                                               ccaagagTGCTGCCGCCTTACGGGCAATGCTTCCACAGCGTAAATACCCCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGAGGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAAtctc                                                                                                            <  1:332875/248‑1 (MQ=255)
                                                                                                                                ccggagggCTGACGCTTACTGGGCAATGCTTCCACACCGAAAACACCCCCACGTACTGGAGCGGGTTGGGGGACCGGCCGGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTGGTTCAGGAACGCCAGAAGGATGCCGCAGACCATTATCTGGTGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCg                                                                                                           <  2:136397/244‑1 (MQ=255)
                                                                                                                                           aCGCCTTCCTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGGGAAGGCCCTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGGTGCCGCAGACCATTTTCTGGCGGGGGATATCGAATCCCTGCCGTTAGCGACTGCGACGGTCGGTCTTg                                                                                                                          >  1:336136/1‑225 (MQ=255)
                                                                                                                                             gCCTTTCTGGCAGTGCTTCCCCACCGTAAACACCCCCCCTTTCTGGACGGGGTTTGGGGACCTGGGTGGTGGGGCCGCCATGGGCGGGACCGTCACGCGCAGGTGACGGCCTTAGCTCCCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTGTCTGGCGGGAGATATCGAATCCCTGCCGTGAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCGCAGTGCAttggtg                                                                                              <  2:448484/251‑1 (MQ=255)
                                                                                                                                              ccTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAAAGTCACGCGCAGGTGAAGGCCTTAGAGCTCTCGCCGCCAATGCTTGTTTAGGCCCGCCAGAAGGATGCCGCAGACCAATATCTGGCGGGGGATATCGAAACACTGCCGTTAGCGGCTGCGGACGTCGATCTTTGAAGGAGCAATCTCGCCGTGCAggggtgg                                                                                             >  1:144177/1‑250 (MQ=255)
                                                                                                                                                    tGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAAAGTCACGCGCAGGGGACGGCCTTAGATCTCTCGCCCCCAATGCTTGTTCAGGCACGCCAGAAGGGTGCCGCAGACCACTATCTTGCGGGGGATATCGAATCCCTGCCGTTTACGGCTGCGACCGTTGATCTTTCATGGAGCCAtcttccagatcaggggggggggaat                                                                                       >  2:183328/1‑229 (MQ=255)
                                                                                                                                                        cAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGGACGTCACGCGCCGGCGAAGGCCTTAGGTCTCTCGCCGCCAATGCTTGTTCAGGCACGCCCGAAGGATGCCGCAGACCCTTATCTGGCGGGGGATATCCAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTTCCTGGAGCAATCTCGCAGTGCCGTGGGGCGGTCAtttt                                                                                    >  1:459075/1‑250 (MQ=255)
                                                                                                                                                          aggctcccACAGCTTAAATACCCCCACGTGCGGGAGGCGGGTTGGGGACCGGGCTGGATGAGCCGCCCCTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGATATAGCGAATCCCTTCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCGGTAATTTATc                                                                                  <  2:459075/245‑1 (MQ=255)
                                                                                                                                                                   caGCGTTAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCCGACCATTATCTGGCGGGAGATATCGAATCCCTGCCGGTAGCGACTTCGACGTTCGATCTTGCATGGAGCAATCTCGCAGTGCAGGGGGGCGGTAATTTATCCACCGCAcc                                                                         >  1:128019/1‑250 (MQ=255)
                                                                                                                                                                   caGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGGGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGGGATATCGAATCCCTGCCGTTAGCGACTGCGGCGTTCGAACTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCGGGAATTTATCCACGGCAcc                                                                         >  1:98013/1‑250 (MQ=255)
                                                                                                                                                                      cccgcctctctcccgatccttggccgcgGGTGGTGGCCCGGGCTGGAGGAGCCGCTTCTGGCGGGAACGTCCCGCTCAGGGGCCGGCCTTAGACCCCCCGCCGCCAAGGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGAAAGCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCGGTAATTTATCCACGGCACTcc                                                                       <  2:207324/227‑1 (MQ=255)
                                                                                                                                                                        aaaatccccccACGTCCTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGAGGCCGCAGACCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCGGTAATTTATCCACGGCACTCcgcg                                                                    <  1:75012/243‑1 (MQ=255)
                                                                                                                                                                           aTACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCCATGCTTGTTCAGGCCCGCCAGAAGGGTGCCGCAGGCCATTATCTGGCGGGAGATATCGAATCCCCGCCGGTAGCGACTGCGACGTTCGATCTTGCATGGAGCCATCTCGCAGGGCCGTGGGGCGGGAATTTATCCACGGCACCCcgcgggc                                                                 >  1:19171/1‑248 (MQ=255)
                                                                                                                                                                                      tACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGGGACGGCCTTAGATCTCTCGCCGCCCATGCTTGTTCAAGCACGCCAGAAGGATGCCGCCGACCCTTATCTGGCGGGAGATATCGAATCCCCGCCGTTAGCGACTGCGACGTTCGATCTTGCCTGGGGCAATCTCGCAGTGCAGTGGGGCGGGAATTTATCCACGGCACTCCCCGCGCTGTATCggggg                                                      >  1:256925/1‑249 (MQ=255)
                                                                                                                                                                                        cTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAAAGTCACGCGCAGGGGACGGCCTTCGATCTCTCGGCGCCAATGCTTGTTCCAGCACGCCAGAAGGATGCCGCAGACCATTATATGGCGGGGGATATCGAATACCTGCCGTTAGCGACTGCGACGTTGGGTCTTGCATGGAGCAATCTCGCCGTGCAGTGGGGCGGTAATTTATCCACGGGACTcc                                                                       >  1:207324/1‑232 (MQ=255)
                                                                                                                                                                                                cgGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGGCACGCGCAGGGGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCCCGCCAGAAGGATGCCCCAGACCACTATCTGGCGGGAGATATCGCATCCCCGCCGTTAGCGACTGCGACGTTCGATCTTTCATGGAGCAATCTCGCAGTGCAGGGGGGCGGTAATTTATCCCCCGCACTCCCCGAGCTGTCTCGGGGGGTGCGcccc                                             >  1:64284/1‑251 (MQ=255)
                                                                                                                                                                                                 gggTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCCCGCGCCAGGGAGGGGCTTAGATCTCTCGCCGCCACTGCTTGTTCAGGCACGCCAGAAGGATGCCGCCGACCATTATCTGACGGGAGATATCGAATCCCCGCCGTTAGCGAATGCGACGGTCGATCTTTCCTGGCGCAATCTCCCAGTGCAGTGGGGCGGTAAATTATCCACGGCACCCCGCACGCTGTATCGGGGGGGGCGccccc                                            >  1:188218/1‑250 (MQ=255)
                                                                                                                                                                                                    ttGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAAAGTCCAGCGCAGGCGGCGGCCTTAGATCCCTCGCCGCCCATGCCTGTTTAGGCACGCCCGAAAGATGCCGCCAGCCCCTATCTGGGGGGCGAAAACGGATACCTGCCGTTAACGACTCCGGCGGTCGGTCTCGCAAAGAGCAATATCGCCATGGAGGGGTGCGGTAATTTTTTCACGGCCCGCCGCGGGCTGTATCGGGTGGGGAGcccccaa                                          >  1:215791/1‑247 (MQ=255)
                                                                                                                                                                                                       tgGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCCGGTGGCGGCCGTTAGTCTCTCGCCGCCCATGCTTGTTCCGGCACGCCCGACGGGTGCCGCCGACCATTATCTGTCGGGGGAGATCGCATCCCTGCCGGTAGCGGCTGCGACGGTCGATCTTGCCTTGGGCAATCCCGCCCGGCAGGGGGGCGGTAATTTATCCCCGGGACCCCCCGAGCTGTATTGGGTGGGGCGcccccacgggg                                      >  2:215550/1‑244 (MQ=255)
                                                                                                                                                                                                                        gcccgccACTGGGGGGAACTCCACGCGCAGTGGACGGCATTAGATCTTCCGCCGCCAATGCTTGTTCAGGCACGACAGGAGGATGCCGCAGCCAATTATCTGGCGGGAGATATGGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCGGTAATTTATCCACGGCACTCCGCGAGCTGTATCGGGTGGTGCGCCCCAAAGGCGTGGTCGCGTTTACCAc                      <  2:80602/249‑1 (MQ=255)
                                                                                                                                                                                                                                 tGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGGGACGTTCGATCTTGCATGGAGCAATCTCGCAGGGCAGTGGGGCGGTAATTTATCCACGGCACTCCGCGAGCGGTATCGGGGGGGGCGCCCCAAAGGCGTGGGCGCGCTTACCACGCTGGTGcc             >  2:94972/1‑250 (MQ=255)
                                                                                                                                                                                                                                 gggCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTTTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCGGTAATTTATCCACGGCACTCCGCGAGCTGTATCGGGTGGTGCGCCCCAAAGGCGTGGTCGCGTTTACCACGCTGGTGCa             <  1:63825/250‑1 (MQ=255)
                                                                                                                                                                                                                                      ggAACGTCGCGCGCAGGGGAGGCCTGTAGCTCTCCCGCCGCCAATGCTTGGTCAGGCACGCCAGAAGGATGCCGCAGACAATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCGGTAATTTAGCCACGGCACTCCGCGAGCTGTATCGGGTGGTGCGCCCCAAAGGCGTGGTCGCGTTTACCACGCTGGTGCAGGGAt        <  2:169911/251‑1 (MQ=255)
                                                                                                                                                                                                                                       gAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCGGTAATTTATCCACGGCACTCCGCGAGCTGTATCGGGTGGTGCGCCCCAAAGGCGTGGTCGCGTTTACCACGCTGGTGCAGGGATc       <  1:404119/251‑1 (MQ=255)
                                                                                                                                                                                                                                       aaaCGTCACGCGCAGGTGACGGCCTTGAACCCCTCGCCGCCAATGTGTGTTAAGGCACGCCAGAAGGATGCCGAAGACCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCGGTAATCTATCCACGGCACTCCGCGAGCTGTATCGGGTGGTGCGCCCCAAAGGCGTGGTCGCGTTTACCACGCTGGTGCAGGGATc       <  2:128019/250‑1 (MQ=255)
                                                                                                                                                                                                                                        aacgccccgcgCAGGTGCCGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCGGTAATTTATCCACGGCACTCCGCGAGCTGTATCGGGTGGTGCGCCCCAAAGGCGTGGTCGCGTTTACCACGCTGGTGCAGGGATCg      <  2:165243/244‑1 (MQ=255)
                                                                                                                                                                                                                                            tCACGCGCAGGTGACGACCTTAGATCTCTCGCCGCCAATACTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCCTTATCTGGCGGGAGATATCGTAAACCTGCCGTTAGCGACTGCGCCGTTCGTTCTTGCATGGAGCAATCTCGCCGTGCAGTGGTGCTGTAATTTATCCACGGCACTCCGCGAGCTGTCTCGGGTGGTGCGACCCAAAGGCGAGGTCGCGTTTACCACTCTGGCGCAGGGCTCGTTAc  >  2:368611/1‑251 (MQ=255)
                                                                                                                                                                                                                                            cccggcgcAGGTGAGGGCCTAAGACCCCTCGCCGCCAATGCTGGTCCAGGCACGCCAGAAGGATGCGGCAGACCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCGGTAATTTATCCACGGCACTCCGCGAGCTGTATCGGGTGGTGCGCCCCAAAGGCGTGGTCGCGTTTACCACGCTGGTGCAGGGATCGTTAc  <  2:19171/247‑1 (MQ=255)
                                                                                                                                                                                                                                             cccgcgCAGGTGACGGCCTTAGACCCCTCGCCGCCAATGCTTGTCCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCGGTAATTTATCCACGGCACTCCGCGAGCTGTATCGGGTGGTGCGCCCCAAAGGCGTGGTCGCGTTTACCACGCTGGTGCAGGGATCGt     <  1:174346/245‑1 (MQ=255)
                                                                                                                                                                                                                                                gcgcAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCGGTAATTTATCCACGGCACTcc                                                                       <  1:68820/178‑1 (MQ=255)
                                                                                                                                                                                                                                                gcgcAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCGGTAATTTATCCACGGCACTcc                                                                       >  2:68820/1‑178 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                          
TGCGCATGAAATGCAGGATATCGACCGTCTGCTGGAGGTGCTGCATGGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCGGTAATTTATCCACGGCACTCCGCGAGCTGTATCGGGTGGTGCGCCCCAAAGGCGTGGTCGCGTTTACCACGCTGGTGCAGGGATCGTTAC  >  NC_000913/811478‑811960

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