Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 2,152,044 A→C 18.9% V360G (GTG→GGG)  yegI ← protein kinase YegI

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0077792,152,0440AC18.9% 120.5 / 2.1 76V360G (GTG→GGG) yegIprotein kinase YegI
Reads supporting (aligned to +/- strand):  ref base A (45/14);  new base C (6/11);  total (51/25)
Fisher's exact test for biased strand distribution p-value = 2.84e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TACGCTTTGCTTGTCAGACTGCCAATAATCCAGATAGCCGCCAGCACAGGAACCAAAATAATATAACGCGGTTCTGCCTGAAGGCCGCAAAGCAGCAGCGATAATGCTGACAGTGCGATCTCAAGCAGAATGATGTATTCGCGCCGTAACAGGCCTAAAGGAAGCGGCCTGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGACTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCACACCACGGGCAGTCAGTCAAATGAGCGGGGTAAACATGCATTGCTGAAACGATACATTTCTTTAACTGTTGGCGTAGAGAGTCCAGTGCCGCTACCCACGCTTTAGCCGTCGGACGCCCGGTTGCCACGCCACTTTCCG  >  NC_007779/2151810‑2152281
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                   
tACGCTTTGCTTGTCAGACTGCCAATAATCCAGATAGCCGCCAGCACAGGAACCAAAATAATATAACGCGGTTCTGCCTGAAGGCCGCAAAGCAGCAGCGATAATGCTGACAGTGCGATCTCAAGCAGAATGATGTATTCGCGCCGTAACAGGCCTAAAGGAAGCGGCCTGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGGGCTTATGAAGCCCTCACCa                                                                                                                                                                                                                                                  >  2:1627583/1‑238 (MQ=255)
     tttGCTTGTCAGACTGCCAATAATCCAGATAGCCGCCAGCACAGGAACCAAAATAATATAACGCGGTTCTGCCTGAAGGCCGCAAAGCAGCAGCGATAATGCTGACAGTGCGATCTCAAGCAGAATGATGTATTCGCGCCGTAACAGGCCTAAAGGAAGCGGCCTGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGccc                                                                                                                                                                                                                                            >  1:1195849/1‑239 (MQ=255)
     tttGCTTGTCAGACTGCCAATAATCCAGATAGCCGCCAGCACAGGAACCAAAATAATATAACGCGGTTCTGCCTGAAGGCCGCAAAGCAGCAGCGATAATGCTGACAGTGCGATCTCAAGCAGAATGATGTATTCGCGCCGTAACAGGCCTAAAGGAAGCGGCCTGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGGGGTCCTGACACCATCCCCacatg                                                                                                                                                                                                                                              >  2:670461/1‑233 (MQ=255)
     tttGCTTGTAAGACTGCCAATAATCCAGATAGCCGCCAGCACAGGAACCAAAATAATATAACGCGGTTCTGCCTGAAGGCCGCAAAGCAGCAGCGATAATGCTGACAGTGCGATCTCAAGCAGAATGATGTATTCGCGCCGTAACAGGCCTAAAGGAAGCGGCCTGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGccc                                                                                                                                                                                                                                            <  1:1293205/239‑1 (MQ=255)
               aGACTGCCAATAATCCAGATAGCCGCCAGCACAGGAACCAAAATAATATAACGCGGTTCTGCCTGAAGGCCGCAAAGCAGCAGCGATAATGCTGACAGTGCGATCTCAAGCAGAATGATGTATTCGCGCCGTAACAGGCCTAAAGGAAGCGGCCTGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGGGCCCATGACACCCACCCCATCACCCCCACTTTct                                                                                                                                                                                                                                   >  2:2129172/1‑237 (MQ=255)
                 aCTGCCAATAATCCAGATAGCCGCCAGCACAGGAACCAAAATAATATAACGCGGTTCTGCCTGAAGGCCGCAAAGCAGCAGCGATAATGCTGACAGTGCGATCTCAAGCAGAATGATGTATTCGCGCCGTAACAGGCCTAAAGGAAGCGGCCTGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGGGCTACTGACGCCATCACCAACGCCCCGACTTTCGCCa                                                                                                                                                                                                                                >  1:498591/1‑239 (MQ=255)
                  cTGCCAATAATCCAGATAGCCGCCAGCACAGGAACCAAAATAATATAACGCGGTTCTGCCTGAAGGCCGCAAAGCAGCAGCGATAATGCTGACAGTGCGATCTCAAGCAGAATGATGTATTCGCGCCGTAACAGGCCTAAAGGAAGCGGCCTGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGACTTTCGc                                                                                                                                                                                                                                  <  1:414042/236‑1 (MQ=255)
                   tGCCAATAATCCAGATAGCCGCCAGCACAGGAACCAAAATAATATAACGCGGTTCTGCCTGAAGGCCGCAAAGCAGCAGCGATAATGCTGACAGTGCGATCTCAAGCAGAATGATGTATTCGCGCCGTAACAGGCCTAAAGGAAGCGGCCTGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGAATTTCGCCAGc                                                                                                                                                                                                                              >  1:1378584/1‑239 (MQ=255)
                   tGCCAATAATCCAGATAGCCGCCAGCACAGGAACCAAAATAATATAACGCGGTTCTGCCTGAAGGCCGCAAAGCAGCAGCGATAATGCTGACAGTGCGATCTCAAGCAGAATGATGTATTCGCGCCGTAACAGGCCTAAAGGAAGCGGCCTGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGGGCCAGCGGACGCCACACCCTCGCCCAAACttttcgccg                                                                                                                                                                                                                               >  1:1050942/1‑232 (MQ=255)
                   tGCCAATAATCCAGATAGACGCCAGCACAGGAACCAAAATAATATAACGCGGTTCTGCCTGAAGGCCGCAAAGCAGCAGCGATAATGCTGACAGTGCGATCTCAAGCAGAATGATGTATTCGCGCCGTAACAGGCCTAAAGGAAGCGGCCTGCCAATCGGTTAGAAATGATAGGGTAATGGCAGGTGCAATGCTGGCGGGGCCACTGCGACCATCACCCTTTCCCAAACttttcccca                                                                                                                                                                                                                               >  2:1881859/1‑232 (MQ=255)
                           atCCAGATAGCCGCCAGCACAGGAACCAAAATAATATAACGCGGTTCTGCCTGAAGGCCGCAAAGCAGCAGCGATAATGCTGACAGTGCGATCTCAAGCAGAATGATGTATTCGCGCCGTAACAGGCCTAAAGGAAGCGGCCTGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGACTTTCGCCAGCACAAAAt                                                                                                                                                                                                                     >  1:295488/1‑238 (MQ=255)
                                     ccgccAGCACAGGAACCAAAATAATATAACGCGGTTCTGCCTGAAGGCCGCAAAGCAGCAGCGATAATGCTGACAGTGCGATCTCAAGCAGAATGATGTATTCGCGCCGTAACAGGCCTAAAGGAAGCGGCCTGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCCAACCTTCGCCACCCTCAAATTCCCGCCggtg                                                                                                                                                                                                         >  1:572570/1‑239 (MQ=255)
                                     ccgccAGCACAGGAACCAAAATAATATAACGCGGTTCTGCCTGAAGGCCGCAAAGCAGCAGCGATAATGCTGACAGTGCGATCTCAAGCAGAATGATGTATTCGCGCCGTAACAGGCCTAAAGGAAGCGGCCTGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGGGCTACTGACGCCATCACCATCGCCCACACTTTCGCCAGCACAAAATTCCCGCCgggg                                                                                                                                                                                                         >  1:1877352/1‑237 (MQ=255)
                                                  aaCCAAAATAATATAACGCGGTTCTGCCTGAAGGCCGCAAAGCAGCAGCGATAATGCTGACAGTGCGATCTCAAGCAGAATGATGTATTCGCGCCGTAACAGGCCTAAAGGAAGCGGCCTGCCAGTCGGTTGGAAATGATCGGGAAATGGCAGTTGCAATGCTGGCGGTGCCACAGAAGGCATCACCATAGCCCAGACTTTCGCCAGCACAAAAATGCCGCCCGGGGGAATGAActcc                                                                                                                                                                                             >  1:91379/1‑237 (MQ=255)
                                                         ataataTAACGCGGTTCTGCCTGAAGGCCGCAAAACAGCAGCGATAATGCTGACAGTGCGATCTCAAGCAGAATGATGTATTCGCGCCGTAACAGGCCTAAAGGAAGCGGCCTGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCATCCACGCCCAGACTTTCGCCAGCACAAAATTACCGCCGGTGGGAATGACCTCTTCGCCgag                                                                                                                                                                                    >  1:1772224/1‑239 (MQ=255)
                                                             tataACGCGGTTCTGCCTGAAGGCCGCAAAGCAGCAGCGATAATGCTGACAGTGCGATCTCAAGCAGAATGATGTATTCGCGCCGTAACAGGCCTAAAGGAAGCGGCCTGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCCCCCTCGCCCACACTTTCGCCCGCACAAAATACCCCCCCGTGGGAATGAACCCTTTGCCCCGATCa                                                                                                                                                                                >  2:495448/1‑239 (MQ=255)
                                                                  cgcgGTTCTGCCTGAAGGCCGCAAAGCAGCAGCGATAATGCTGACAGTGCGATCTCAAGCAGAATGATGTATTCGCGCCGTACCAGGCCTAAAGGAAGCGGCCTGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTTATGACGCCATCACCAACGCCCCGACTTTTGCCACCACAAACTTTCCGCCGGAGGGAATAACCCTCTCGCCGAGATCAATaaa                                                                                                                                                                           >  1:1924469/1‑239 (MQ=255)
                                                                  cgcgGTTCTGCCTGAAGGCCGCAAAGCAGCAGCGATAATGCTGACAGTGCGATCTCAAGCAGAATGATGTATTCGCGCCGTAACAGGCCTAAAGGAAGCGGCCTGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCCGACTTTCGCCAGCACAAAATTCCCCCCGGTGGTAATGACCTCTTCGCCGGGATCAATaaa                                                                                                                                                                           >  1:1723638/1‑239 (MQ=255)
                                                                  cgcgGTTCTGCCTGAAGGCCGCAAAGCAGCAGCGATAATGCTGACAGTGCGATCTCAAGCAGAATGATGTATTCGCGCCGTAACAGGCCTAAAGGAAGCGGCCTGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCCACTGACGCCATCACCATCGCCCAGACTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCCATaaa                                                                                                                                                                           >  1:342625/1‑239 (MQ=255)
                                                                     ggTTCTGCCTGAAGGCCGCAAAGCAGCAGCGATAATGCTGACAGTGCGATCTCAAGCAGAATGATGTATTCGCGCCGTAACAGGCCTAAAGGAAGCGGCCTGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCCGACTTTCGCCAGCACAAAATTCCCCCCGGGGGTAATGACCCCTTCGCCGAGATCAATAAAATa                                                                                                                                                                        >  1:301312/1‑239 (MQ=255)
                                                                       ttCTGCCTGAAGGCCGCAAAGCAGCAGCGATAATGCTGACAGTGCGATCTCAAGCAGAATGATGTATTCGCGCCGTAACAGGCCTAAAGGAAGCGGCCTGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCCTCACCATCGCCCAGACTTTCGCCAGCACAAAAATCCCGCCGGGGGTAATGACCTCTTCGCCGAGATCAATAAAATAGa                                                                                                                                                                      >  1:684861/1‑239 (MQ=255)
                                                                       ttCTGCCTGAAGGCCGCAAAGCAGCAGCGATAATGCTGACAGTGCGATCTCAAGCAGAATGATGTATTCGCGCCGTAACAGGCCTAAAGGAAGCGGCCTGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGACTTTCGCCAGCACAAAATTCCCGCCGGGGGGAATGACCTCTTCCCCGAGATCAATAAAATAGa                                                                                                                                                                      >  1:474271/1‑239 (MQ=255)
                                                                       ttCTGCCTGAAGGCCGCAAAGCAGCAGCGATAATGCTGACAGTGCGATCTCAAGCAGAATGATGTATTCGCGCCGTAACAGGCCTAAAGGAAGCGGCCTGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGGGCTACTGACGCCATCACCATCGCCCAGACTTTCGCCAGCACAAAATTCCCCCCGGGGGGAAAGACCTCTTCGCCCAGATCAAAAAAATa                                                                                                                                                                        >  1:69872/1‑237 (MQ=255)
                                                                            ccTGAAGGCCGCAAAGCAGCAGCGATAATGCTGACAGTGCGATCTCAAGCAGAATGATGTATTCGCGCCGTAACAGGCCTAAAGGAAGCGGCCTGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACCGACGCCATCACCCTCCCCCAGACTTTCGCCCACACAAAATTCCCGCCCGGGGGAATGACCTCTTCCCCGAGATCAATAAAATAGATaaa                                                                                                                                                                  >  1:1774505/1‑237 (MQ=255)
                                                                             cTGAGGGCCGAAAAGCAGCGGCGATAATGCTGACAGGGCGATCTCAAGCAGAATGTTGTATTCGCGCCGTAACAGGCCTAAAGGAAGCGGCCTGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGACTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAacac                                                                                                                                                               <  1:447491/239‑1 (MQ=255)
                                                                               gccagcagCAAAGCGGCAGCGAAATGGCGGACTGGGGGATCTCAAGAAGAAGGTTGTTTCCCCGCGGTACAGGGCCTAAAGGGGCCGGCCTGCCGGGCGTTTGGAAATGATCGGGTAATGGCAGTTCAAATGCGGGCGGTGCTCCTGCCCCCATCCCCATCCCCCATCCTTCCCCCACAAAAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAacac                                                                                                                                                               <  2:948342/230‑1 (MQ=255)
                                                                                 aGGCCGCAAAGCAGCAGCGATAATGCTGACAGTGCGATCTCAAGCAGAATGATGTATTCGCGCCGTAACAGGCCTAAAGGAAGCGGCCTGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGACTTTCGCCAGCACAAAAATCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCCATAAAATAGATAAAACCTTg                                                                                                                                                           >  1:148405/1‑239 (MQ=255)
                                                                                    ccGCAAAGCAGCAGCGATAATGCTGACAGTGCGATTTCAAGCAGAATGATGTATTCGCGCCGTAACAGGCCTAAAGGAAGCGGCCTGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGGGCTACTGACGCCCTCACCATCGCCCAGACTTTTGCCAGCACAAAAATCCCGCCGGGGGGAAAGACCCCTTCGCCGGGAACAATAAAAAAGAAAACACCTTTAtt                                                                                                                                                       >  2:447491/1‑239 (MQ=255)
                                                                                      gCAAAGCAGCAGCGATAATGCTGACAGTGCGATCTCAAGCAGAATGATGTATTCGCGCCGTAACAGGCCTAAAGGAAGCGGCCTGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCa                                                                                                                                                                                                                                                  >  1:1603446/1‑152 (MQ=255)
                                                                                      gCAAAGCAGCAGCGATAATGCTGACAGTGCGATCTCAAGCAGAATGATGTATTCGCGCCGTAACAGGCCTAAAGGAAGCGGCCTGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCa                                                                                                                                                                                                                                                  <  2:1603446/152‑1 (MQ=255)
                                                                                             agcagcGATAATGCTGACAGTGCGATCTCAAGCAGAATGATGTATTCGCGCCGTAACAGGCCTAAAGGAAGCGGCCTGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGACTTTCGCCAGCACAAAATTCCCGCCGGGGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAAAACCTTGATTATCCAgcg                                                                                                                                               >  1:1805718/1‑238 (MQ=255)
                                                                                               cagcGATAATGCTGACAGTGCGATCTCAAGCAGAATGATGTATTCGCGCCGTAACAGGCCTAAAGGAAGCGGCCTGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGACTTTCGCCAGCACAAAATTCCCGCCGGGGGGAATGACCTCTTCGCCCAGATCAATAAAATAGATAACACCTTGATTATccctttc                                                                                                                                              >  1:948342/1‑232 (MQ=255)
                                                                                                agcGATAATGCTGACAGTGCGATCTCAAGCAGAATGATGTATTCGCGCCGTAACAGGCCTAAAGGAAGCGGCCTGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGACTTTCGCCAGCAAAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATCCAgcgc                                                                                                                                              <  1:363391/236‑1 (MQ=255)
                                                                                                       aTGCTGACAGTGCGATCTCAAGCAGAATGATGTATTCGCGCCGTAACAGGCCTAAAGGAAGCGGCCTGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGGGGTACGTAACCCATCACCAATGCCCAGACTTTTGCCCGCAAAAAAATCCCGCCGGGGGGAAAAGACCCTTTCCCGGGAACAAAAAAAAAAAAAACACCCTTATTTCCCAGCGGAcacccacg                                                                                                                                     >  2:31727/1‑234 (MQ=255)
                                                                                                         gCTGACAGGGCGACCTCAAGCAGAATGAGTTTTTCGCGCCTTAACAGGCCTAAAGGAAGCGGCCTGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCCCCATCGCCCAGACTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGc                                                                                                                                                                                        <  1:490555/187‑1 (MQ=255)
                                                                                                           tGACAGTGCGATCTCAAGCAGAATGATTTTTTCGCGCCGTAACAGGCCTAAAGGAAGCGGCCTGCCAGTCGGTTGAAAAGGATCGGGTAAGGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGACTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCACACCACGGGcag                                                                                                                                 <  1:369562/238‑1 (MQ=255)
                                                                                                             aCAGTGCGATCTCAAGCAGAATGATGTATTCGCGCCGTAACAGGCCTAAAGGAAGCGGCCTGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGGTACTGACCACCTCCCCAACGCCCCAACTTTCGCCCGCACAAAAAATCCCCCGGCGGGAACAGACTCCTCCCCCAGAGCGAAAAAAATGAAAACACCCTTATTATCCAACGCACCCCCCCGGcagtca                                                                                                                              >  2:1255465/1‑239 (MQ=255)
                                                                                                             aCAGTGCGATCTCAAGCAGAATGATGTATTCGCGCCGTAACAGGCCTAAAGGAAGCGGCCTGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCCTCACCATCGCCCAGACTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATc                                                                                                                                                    >  1:581691/1‑217 (MQ=255)
                                                                                                             aCAGTGCGATCTCAAGCAGAATGATGTATTCGCGCCGTAACAGGCCTAAAGGAAGCGGCCTGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACCACCTCACCCTCCCCCCGACTTTCGCCAGCACAAAAATCCCGCCGGTGGTAAAGACCTCTTCGCCGCGATCAATAAAATAGAAAACACCTTGATTATCCAGCGCACACCACGGGcagtc                                                                                                                               >  2:1298626/1‑238 (MQ=255)
                                                                                                                 tGCGATCTCAAGCAGAATGATGTATTCGCGCCGTAACAGGCCTAAAGGAAGCGGCCTGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCCGACTTTCGCCAGCCCAAAATTCCCCCCGGGGGTAATTACCCCTTCGCCCAGGTCCATAAAAAAAATAAACCCTTGATTTTCCAACCCACACCCCGCGCAgtcagtca                                                                                                                          >  2:1034939/1‑239 (MQ=255)
                                                                                                                 tGCGATCTCAAGCAGAATGATGTATTCGCGCCGTAACAGGCCTAAAGGAAGCGGCCTGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGGGCTACTGACCCCCCCCCCATCTCCCAGACTTTTGCCCGCCCAAAAATTCCCCCCGGGGGAAAGACCCCTTTCCCCAGATAAATAAAATATATAAACCACTGATTTTCCACCGCACACCCCAGGcaaccacata                                                                                                                          >  2:9010/1‑231 (MQ=255)
                                                                                                                           aGCAGAATGATGTATTCGCGCCGTAACAGGCCTAAAGGAAGCGGCCTGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGACTTTCGCCAGCACAAAATTCCCGCCGGTGGGAAAGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCACACCACGGGCAGTCAGTCAAATGAGCggg                                                                                                                >  1:126673/1‑239 (MQ=255)
                                                                                                                                  tgatgTATTCGCGCCGTAACAGGCCTAAAGGAAGCGGCCTGCCAGTCGGTTGGAAATGATCGCGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCCGACTTTCGCCAGCACAAAATTCCCCCCGGGGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCACACCACGGGCAGGCAGGCAAATGAGCGGGGTaaaaa                                                                                                         >  1:12060/1‑237 (MQ=255)
                                                                                                                                     tgTATTCGCGCCGTAACAGGCCTAAAGGAAGCGGCCTGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGGTACTGACGCCATCACCa                                                                                                                                                                                                                                                  >  2:174196/1‑105 (MQ=255)
                                                                                                                                     tgTATTCGCGCCGTAACAGGCCTAAAGGAAGCGGCCTGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCa                                                                                                                                                                                                                                                  <  1:174196/105‑1 (MQ=255)
                                                                                                                                        aTTCGCGCCGTAACAGGCCTAAAGGAAGCGGCCTGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGGGGTAATTAAGCCCACACCCACGCCCAGACTTTCCCCAACCCAAAATTTCCCCCCGGGGGAAAGACCCCTCTCCCCGGATCAATAAAAAAAATAACAACCTTGTTTTACATCGACCACCACGGGGAAGACGTAAAAAGAGCGGGGGaaaaaagccat                                                                                                   >  2:1733740/1‑231 (MQ=255)
                                                                                                                                          tCGCGCCGTAACAGGCCTAAAGGAAGCGGCCTGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGACTTTCGCCAGCACAAAATTCCCGCCGGGGGGAATGACCTCTTCGCCGAGATCAATAAAATAGATAaca                                                                                                                                                                >  1:1645261/1‑177 (MQ=255)
                                                                                                                                          tCGCGCCGTAACAGGCCTAAAGGAAGCGGCCTGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTAACGACACCCTCACCCTCCCCCAGACTTTTCCCAGCACAAAAATCCCGCCGGGGGTCATGACCCCCTTCCCGGGAGCAATAAAAAAAATAACACCCTGATTATCCAGCCCAACCCACGGGCAGGCAGACAAAAGAGCCGGGTAAACaaacaat                                                                                                   >  2:711590/1‑231 (MQ=255)
                                                                                                                                              gcCGTAACAGGCTAAAAGGAAGGGGCCTGCCAGTCGGTTGAAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCCCCATCGCCCAGACTTTCGCCAGCAAAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGAt                                                                                                                                                                     <  1:894847/169‑1 (MQ=255)
                                                                                                                                              gcCGTAACAGGCCTAAAGGAAGCGGCCTGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCCACACCATCGCCCAGACTTTCGCCAGCACAAAATTCCCGCCGGGGGGAATGACCTCTTCGCCGCGCTCAATAAAATAGATAAAACCTTGATTATCCAGCGCACACCACGGGCAGGCAGTCAAAAGAGCGGGGGAAAAAAGCATTGCCGaa                                                                                             >  2:677933/1‑239 (MQ=255)
                                                                                                                                              gcCGTAACAGGCCTAAAGGAAGCGGCCTGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCGTCTTACGCCCCCCCCAACCCCCAGAACTTTGCCCGCACAAAATTTCCGCCCGGGGGAATTACCCTTTCGCCCGGGTCCAAAAAAAAAATAAAACCCTGAATTATCCAGGCCCACCACGGGGCCCCAGACCAAAGAGAGGGGGAAAAATGCAtaccgtaa                                                                                             >  2:1515225/1‑232 (MQ=255)
                                                                                                                                              gcCGAAACAGGCCTAAAGGAAGCGCCCTGCCAGTCGTTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGCCGCCATCACCATCGCCCAGACTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAacac                                                                                                                                                               <  1:392216/174‑1 (MQ=255)
                                                                                                                                               ccGTAACGGACAAAAGGGAAGCGCCGTGCCGGTCGGTTGAAAATGATCGGGTAATGGCAGTTGAAATGCTGGCGGTCCTACTGCCCCCATCCCCATCCCCCAGACTTCCGCCAACAAAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCACACCACGGGCAGTCAGTCAAATGAGCGGGGTAAACATGCATTGCTGaaa                                                                                            <  2:879823/239‑1 (MQ=255)
                                                                                                                                                 gTAACAGGCCTAAAGGAAGCGGCCTGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGACTTTCGCCAGCACAAAATTCCCGCCGGTGGGAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCACACCACGGGCAGTCAGGCAAATGAGCGGGGTAAACATGCATTGCTGAAACg                                                                                          >  1:205638/1‑239 (MQ=255)
                                                                                                                                                         gccgaaGGAACGCCGCCCGCCAGTCGGTTGAAAATGTTCGGGTAATGGCAGTGTCAATGCGGCCGGTGCTACTCACCCCATCACCATCGCCAAGACTTGCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCACACCACGGGCAGTCAGTCAAATGAGCGGGGTAAACATGCATTGCTGAAACGATACAt                                                                                    <  1:362585/234‑1 (MQ=255)
                                                                                                                                                            aaaGGAAGCGGCCTGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGACTTTCGCCAGCACAAAATTCCCGCCGGGGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCACACCACGGGCAGTCAGTCAAATGAGCGGGGGAAACATGCATTGCTGAAACGATACATTTCtt                                                                               >  1:232988/1‑239 (MQ=255)
                                                                                                                                                                   aggggCCTGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGACTTTCGCCAGCAAAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCACACCACGGGCAGTCAGTCAAATGAGCGGGGTAAACATGCATTGCTGAAACGATACATTTCTTTAACTg                                                                         <  1:495448/236‑1 (MQ=255)
                                                                                                                                                                    gcgggcggCCAGGCGGTTGAAAAGGATCGAGTAGTGGCATTCGCAGTGCTGGCGGTGCAAACCACGCCATCACCACCCCCCAGACTTTCGACACACAAAAATACCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCACACCACGGGCAGTCAGTCAAATGAGCGGGGTAAACATGCATTGCTGAAACGATACATTTCTTTAACTg                                                                         <  2:296855/238‑1 (MQ=255)
                                                                                                                                                                           gCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCCCCATCGCCCAGACTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCACACCACGGGCAGTCAGTCAAATGAGc                                                                                                                   <  1:696701/189‑1 (MQ=255)
                                                                                                                                                                            ccAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGACTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCACACCACGGGCAGTCAGTCAAATGAGCGGGGTAAACATGCATTGCTGAAACGATACATTTCTTTAACTGTTGGCGTag                                                                >  1:1768063/1‑239 (MQ=255)
                                                                                                                                                                             cattCGTTTGGAAATGACGGGGTAATGGCAGTTGCAATGCGGGCGGTGCTCCTGACCCCATCCCCATCGCCCAGACTTCCGCCAGCAAAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCACACCACGGGCAGTCAGTCAAATGAGCGGGGTAAACATGCATTGCTGAAACGATACATTTCTTTAACTGTTGGCGTaga                                                               <  1:31727/236‑1 (MQ=255)
                                                                                                                                                                             cAGTCGGTTGAAAATGATCGGGTAATGGCAGTTGCAATCCTGGCGGTGCTACTGCCGCCATCCCCATCGCCCAGACTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCACACCACGGGCAGTCAGTCAAATGAGCGGGGTAAACATGCATTGCTGAAACGATACATTTCTTTAACTGTTGGCGTaga                                                               <  1:677933/239‑1 (MQ=255)
                                                                                                                                                                                     tGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGACTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATaaaa                                                                                                                                                                          >  1:1477401/1‑126 (MQ=255)
                                                                                                                                                                                      ggAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGGCAATGAAGCCATCACCATCGCCCAGACTTTTGCCAGCACAAAAATCCCGCCGGGGGGAATGAACCCTTCGCCGAGGTCAATAAAAAAGATAACACCTTGATTATCCAACGCCCACCACCGGAAGTCAATCAAAAGAGCGGGGGAAAAAATCATTGCTGAAACGATAACTTTCTTTAACTGATGGCGGAGAGAATCCAgg                                                      >  2:633083/1‑238 (MQ=255)
                                                                                                                                                                                          aTGATCGGGTAATGGCAGTTGCAATGCTGGGGGTGCTACTGACGCCATCCCCATCGCCCAGACTTTCGCCAGCAAAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCACACCACGGGCAGTCAGTCAAATGAGCGGGGTAAACATGCATTGCTGAAACGATACATTTCTTTAACTGTTGGCGTAGAGAGTCCAGTGCCg                                                  <  1:1852050/239‑1 (MQ=255)
                                                                                                                                                                                          aTGATCGGGTAATGGCAGTTGCAATGCGGGCGGTGCTACTGACGCCATCACCATCGCCAAGACTTTCGCCAGCAAAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCACACCACGGGCAGTCAGTCAAATGAGCGGGGTAAACATGCATTGCTGAAACGATACATTTCTTTAACTGTTGGCGTAGAGAGTCCAGTGCCg                                                  <  1:603478/239‑1 (MQ=255)
                                                                                                                                                                                           tGATGGGGTAATGGCGTTTGAAATGCTGGCAGTGCTCCGGACCCCATCCCCATCCCCCAGTCTTTCGCCACCAAAAATTCCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCACACCACGGGCAGTCAGTCAAATGAGCGGGGTAAACATGCATTGCTGAAACGATACATTTCTTTAACTGTTGGCGTAGAGAGTCCAGTGCCGc                                                 <  1:1255465/239‑1 (MQ=255)
                                                                                                                                                                                                 ggTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCAGCGCCCAGCCTTTCGCCAGCAAAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCACACCACGGGCAGTCAGTCAAATGAGCGGGGTAAACATGCATTGCTGAAACGATACATTTCTTTAACTGTTGGCGTAGAGAGTCCAGTGCCGCTAccc                                            <  1:1502665/238‑1 (MQ=255)
                                                                                                                                                                                                                    cTGGCGGTGCTACTGACGCCATCACCATCGCCCAGACTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCACACCACGGGCAGTCAGTCAAATGAGCGGGGTAAACATGCATTGCTGAAACGATACATTTCTTTAACTGTTGGCGTAGAGAGTc                                                          >  1:184949/1‑205 (MQ=255)
                                                                                                                                                                                                                       gCGGTGCTACTGACGCCATCACCATCGCCCAGACTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCACACCACGGGCAGTCAGTCAAATGAGCGGGGTAAACATGCATTGCTGAAACGATACATTTCTTTAACTGTTGGCGTAGAGAGTCCAGTGCCGCTACCCACGCTTTAGCCGTCGGACGCCCg                     >  1:1881980/1‑239 (MQ=255)
                                                                                                                                                                                                                              ctccTGACGCCATCCCCCTCGCCCGGCTTTCGCCCACAAAAATATCGCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCACACCACGGGCAGTCAGTCAAATGAGCGGGGTAAACATGCATTGCTGAAACGATACATTTCTTTAACTGTTGGCGTAGGGAGTCCAGTGCCGCTACCCACGCTTTAGCCGTCGGACGCCCGGTTgc                <  2:12060/235‑1 (MQ=255)
                                                                                                                                                                                                                               tccggACCCCACCCCCACCGCCCAGCCTTTCCCCAGCAAAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCACACCACGGGCAGTCAGTCAAATGAGCGGGGTAAACATGCATTGCTGAAACGATACATTTCTTTAACTGTTGGCGTAGAGAGTCCAGTGCCGCTACCCACGCTTTAGCCGTCGGACGCCCGGTTgcca              <  1:9010/235‑1 (MQ=255)
                                                                                                                                                                                                                                    aCGCCATCACCATCGCCCAGACTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCACACCACGGGCAGTCAGTCAAATGAGCGGGGTAAACATGCATTGCTGAAACGATACATTTCTTTAACTGTTGGCGTaga                                                               >  2:338523/1‑185 (MQ=255)
                                                                                                                                                                                                                                    aCGCCATCACCATCGCCCAGACTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCACACCACGGGCAGTCAGTCAAATGAGCGGGGTAAACATGCATTGCTGAAACGATACATTTCTTTAACTGTTGGCGTaga                                                               <  1:338523/185‑1 (MQ=255)
                                                                                                                                                                                                                                      gCCATCACCATCGCCCAGACTTTCGCCAGCACAAAATTCCCGCCGGTGGAAATGACCTCTGCGCCGAGCTCAATAAAATAGATAACACCTTGATTATCCAGCGCACACCACGGGCAGTCAGTCAAATGAGCGGGGTAAACATGCATTGCTGAAACGATACATTTCTTTAACTGTTGGCGTAGAGAGTCCAGTGCCGCTACCCACGCTTTAGCCGTCGGACGCCCGGTTGCCACGCCACt       >  2:677716/1‑239 (MQ=255)
                                                                                                                                                                                                                                         aTCACCATCGCCCAGACTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCACACCACGGGCAGTCAGTCAAATAAGCGGGGTAAACATGCATTGCTGAAACGATACATTTCTTTAACTGTTGGCGTAGAGAGTCCAGTGCCGCTACCCACGCTTTAGCCGTCGGACGCCCGGTTGCCACGCCACttt     >  2:1090136/1‑238 (MQ=255)
                                                                                                                                                                                                                                           caccacCGCCCAGACTTTCGCCACCACAAAATCCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCACACCACGGGCAGTCAGTCAAATGAGCGGGGTAAACATGCATTGCTGAAACGATACATTTCTTTAACTGTTGGCGTAGAGAGTCCAGTGCCGCTACCCACGCTTTAGCCGTCGGACGCCCGGTTGCCACGCCACTTTCCg  <  1:1034939/239‑1 (MQ=255)
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                   
TACGCTTTGCTTGTCAGACTGCCAATAATCCAGATAGCCGCCAGCACAGGAACCAAAATAATATAACGCGGTTCTGCCTGAAGGCCGCAAAGCAGCAGCGATAATGCTGACAGTGCGATCTCAAGCAGAATGATGTATTCGCGCCGTAACAGGCCTAAAGGAAGCGGCCTGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGACTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCACACCACGGGCAGTCAGTCAAATGAGCGGGGTAAACATGCATTGCTGAAACGATACATTTCTTTAACTGTTGGCGTAGAGAGTCCAGTGCCGCTACCCACGCTTTAGCCGTCGGACGCCCGGTTGCCACGCCACTTTCCG  >  NC_007779/2151810‑2152281

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