Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,439,013 T→G 6.7% I84L (ATT→CTT)  ytfJ ← putative transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,439,0130TG 6.7% 149.2 / 2.4 50I84L (ATT→CTT) ytfJputative transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (22/24);  new base G (2/2);  total (24/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01

AGCAGGTCCATCACTTGCTGCACCTCTTCCGGTGTGAGCGCCCCGTCTTTGGCCCATTGCACGCGACCATCCTTATCCAGCACCACTACAGCGGAACTTTCTTCATCCAGCTGCCAGGCACCGAGCGCGACGCCATTGCTATCGACAATAAACTGCGACCAGGGATAAAGCTTTTTATTACTCTCCAGACTGCTGCGCACAAACATCCCTGAACCCGGAATCGCGTCGTCGGTGTTAACAATGGTGGTGGTCTGGTAACGATCGTGCGGTAACTTCGCTGATTTAATCGCTTCAATCAGCGTCGCGTTTTTCTCTTTTGCAGAGGTGCGACCAGCAATATGTTGCAGTACTCGCACTTTTCCCACTAACTGCGCGCTGTTCCAGGTTTTGTAGCTAAACTGATCTTTATCAAGCACCAACTCGCCCCGATCGGTAATGCCAATCGGCGGCACTCGCTGACCGGTTTCGAACTGATGTGCGGAAGCCATCATCGGCAACAGCAGGCAGGTGAGTGCCAGAATCTTGCGTAGGGTCATGGTGTTTCCTTCTTATGATATGCAGGTGATCCGACCACTTGGGTCTGAG  >  NC_000913/4438727‑4439311
                                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                                                          
aagaggtgcagcaagtgatggaCCTCTTCCGGTGTGAGCGCCCCGTCTTTGGCCCATTGCACGCGACCATCCTTATCCAGCACCACTACAGCGGAACTTTCTTCATCCAGCTGCCAGGCACCGAGCGCGACGCCATTGCTATCGACAATAACCTGCGACCAGGGATAAAGCTTTTTATTACTCTCCAGACTGCTGCGCACAAACATCCCTGAACCCGGAATCGAGTCGTCGGTGTTAACAATGGGGGTGGTCAGGGAACGATCGTGCGGTAACTTCCCCGATTTAATCCCTCCAATAAGcc                                                                                                                                                                                                                                                                                              >  1:18047/22‑300 (MQ=255)
   ctctgccggtctgtcggcACCTCTTCCGGTGTTAGCGCCCCGGCTTTGGCCCATTGCACGCTATCATCCTTATCCAGCACCACTACAGCGGAACTTGCTTCATCCAGCTGCCAGGCACCGAGCGCGACGCCATTGCTATCGACAATAAACTGCGACCAGGGATAAAGCTTGTTATTACTCTCCAGACTGCTGCGCACAAACAGCCCTGAACCCGGAATCGCGTCGTCGGTGTTAACAATGGTGGTGGTCTGGTAACGATCGTGCGGTAACTTCGCTGATTGAAGCGATTAAATCAGCgaag                                                                                                                                                                                                                                                                                           <  1:554656/285‑4 (MQ=255)
             cTTGCTGCACCTCTTCCGGTGTGAGCGCCCCGTCTTTGGCCCATTGCACGCGACCATCCTTATCCAGCACCACTACAGCGGAACTTTCTTCATCCAGCTGCCAGGCACCGAGCGCGACGCCATTGCTATCGACAATAAACTGCGACCAGGGATAAAGCTTTTTATTACTCTCCAGACTGCTGCGCACAAACATCCCTGAACCCGGAATCGCGTCGTCGGTGTTAACAATGGTGGTGGTCTGGTAACGATCGTGCGGTAACTTCGCTGATTTAATCGCTTCAATCAGCGTCGCGTTTTtctc                                                                                                                                                                                                                                                                                 >  2:503604/1‑301 (MQ=255)
               tgctgcACCTCTTCCGGTGTGAGCGCCCCGTCTTTGGCCCATTGCACGCGACCATCCTTATCCAGCACCACTACAGCGGAACTTTCTTCATCCAGCTGCCAGGCACCGAGCGCGACGCCATTGCTATCGACAATAAACTGCGACCAGGGATAAAGCTTTTTATTACTCTCCAGACTGCTGCGCACAAACATCCCTGAACCCGGAATCGCGTCGTCGGTGTTAACAATGGTGGTGGTCTGGTAACGATCGTGCGGTAACTTCGCTGATTGAAGCGATTAAATCAGCgaagagatcggaagaa                                                                                                                                                                                                                                                                               >  2:555370/1‑286 (MQ=255)
                ctttccctactcgtcgctcttcctttctttCTTTGGCCCATTGCACGCGACCATCCTTATCCAGCACCACTACAGCGGAACTTTCTTCATCCAGCGGCCAGGCACCGAGCGCGACGCCATTGCTATCGACAATAAACTGCGACCAGGGATAAAGCTTTTTATTACTCTCCAGACTGCTGCGCACAAACATCCCTGAACCCGGAATCGCGTCGTCGGTGTTAACAATGGTGGTGGTCTGGTAACGATCGTGCGGTAACTTCGCTGATTTAATCGCTTCAATCAGCGTCGCGTTTTTCTCttt                                                                                                                                                                                                                                                                              <  1:386112/272‑1 (MQ=255)
                   gtcccTCTTCCGGTGTGAGCGCCCCGTCTTTGGCCCATTGCACGCGACCATCCTTATCCAGCACCACTACAGCGGAACTTTCTTCATCCAGCTGCCAGGCACCGAGCGCGACGCCATTGCTATCGACAATAAACTGCGACCAGGGATAAAGCTTTTTATTACTCTCCAGACTGCTGCGCACAAACATCCCTGAACCCGGAATCGCGTCGTCGGTGTTAACAATGGTGGTGGTCTGGTAACGATCGTGCGGTAACTTCGCTGATTTAATCGCTTCAATCAGCGTCGCGTTTTTCTCTTTTGc                                                                                                                                                                                                                                                                           <  1:196700/298‑1 (MQ=255)
                   gcACCTCTTCCGGTGTTAGCGTCCCGTCTTTGGCCCATTGCACGCGACCTGCCTTATCCAGCACCACTACAGCGGAACTTTCTTCATCCAGCTGCCAGGCACCGAGCGCGACGCCATTGCTATCGACAATAAACTGCGACCAGGGATAAAGCTTTTTATTACTCTCCAGACTGCTGCGCACAAACATCCCTGAACCCGGAATCGCGTCGTCGGTGTTAACAATGGTGGTGGTCTGGTAACGATCGTGCGGTAACTTCGCTGATTTAATCGCTTCAATCAGCGTCGCGTTTTTCTCTTTTGc                                                                                                                                                                                                                                                                           <  1:582695/301‑1 (MQ=255)
                      ccTCTTCCGGTGTGAGCGCCCCGTCTTTGGCCCATTGCACGCGACCATCCTTATCCAGCACCACTACAGCGGAACTTTCTTCATCCAGCTGCCAGGCACCGAGCGCGACGCCATTGCTATCGACAATAAACTGCGACCAGGGATAAAGCTTTTTATTACTCTCCAGACTGCTGCGCACAAACATCCCTGAACCCGGAATCGCGTCGTCGGTGTTAACAATGGTGGTGGTCTGGTAACGATCGTGCGGTAACTTCGCTGATTTAATCGCTTCAATCAGCGTCGCGTTTTTCTCTTTTGCaga                                                                                                                                                                                                                                                                        >  2:125197/1‑301 (MQ=255)
                          ttCCGGTGTGAGCGCCCCGTCTTTGGCCCATTGCACGCAACCATCCTTATCCAGCACCACTACAGCGGAACTTTCTTCATCCAGCTGCCAGGCACCGAGCGCGACGCCATTGCTATCGACAATAAACTGCGACCAGGGATAAAGCTTTTTATTACTCTCCAGACTGCTGCGCACAAACATCCCTGAACCCGGAATCGCGTCGTCGGTGTTAACAATGGTGGTGGTCTGGTAACGATCGTGCGGTAACTTCGCTGATTTAATCGCTTCAATCAGCGTCGCGTTTTTCTCTTTTGCAAAGGTg                                                                                                                                                                                                                                                                    >  2:144960/1‑301 (MQ=255)
                                            gTCTTTGGCCCATTGCACGCGACCATCCTTATCCAGCACCACTACAGCGGAACTTTCTTCATCCAGCTGCCAGGCACCGAGCGCGACGCCATTGCTATCGACAATAAACTGCGACCAGGGATAAAGCTTTTTATTACTCTCCAGACTGCTGCGCACAAACATCCCTGCACCCGGAATCGCGTCGTCGGTGTTAACAATGGTGGTGGTCTGGTAACGATCGTGCGGTAACTTCGCTGATTTAATCGCTTCAATCAGCGTCGCGTTTTTCTCtttagatcggaagagcacacgtcagaactcc                                                                                                                                                                                                                                                  >  2:386626/1‑273 (MQ=255)
                                              ttttGGCCCATTGCACGCGACCATCCTTATCCAGCACCACTACAGCGGAACTTTCTTCATCCAGCTGCCAGGCACCGAGCGCGACGCCATTGCTATCGACAATAAACTGCGACCAGGGATAAAGCTTTTTATTACTCTCCAGACTGCTGCGCACAAACATCCCTGAACCCGGAATCGCGTCGTCGGTGTTAACAATGGTGGTGGTCTGGTAACGATCGTGCGGTAACTTCGCTGATTTAATCGCTTCAATCAGCGTCGCGTTTTTCTCTTTTGCAGAGGTGCGACCAGCAATATGTTGCAg                                                                                                                                                                                                                                                <  2:301689/300‑1 (MQ=255)
                                                            ccgcgACCCTCCTTATCCAGCACCACTACAGCGGAACTTTCTTCATCCAGCTGCCAGGCACCGAGCGCGACGCCATTGCTATCGACAATAAACTGCGACCAGGGATAAAGCTTTTTATTACTCTCCAGACTGCTGCGCACAAACATCCCTGAACCCGGAATCGCGTCGTCGGTGTTAACAATGGTGGTGGTCTGGTAACGATCGTGCGGTAACTTCGCTGATTTAATCGCTTCAATCAGCGTCGCGTTTTTCTCTTTTGCAGAGGTGCGACCAGCAATATGTTGCAGTACTCGCACTTTTc                                                                                                                                                                                                                                  <  1:125073/300‑1 (MQ=255)
                                                                      ccTTATCCAGGACCACTACAGCGGAACTTTCTTCATCCAGCTGCCAGGCACCGAGCGCGACGCCATTGCTATCGACAATAAACTGCGACCAGGGATAAAGCTTTTTATTACTCTCCAGACTGCTGCGCACAAACATCCCTGAACCCGGAATCGCGTCGTCGGTGTTAACAATGGTGGTGGTCTGGTAACGATCGTGCGGTAACTTCGCTGATTTAATCGCTTCAATCAGCGTCGCGTTTTTCTCTTTTGCAGAGGTGCGACCAGCAATATGTTGCAGTACTCGCACTTTTCCCACTAACTg                                                                                                                                                                                                                        <  1:502932/301‑1 (MQ=255)
                                                                       cgtttCCATCACTCCTCCAGCGGGACTTTCTTCTTCCGGCTGTCAGGCAACGAGCGCGAGGCCATTGCTATCGACAATAAACTGCGACCAGGGATAAAGCTTTTTATTACTCTCCAGACTGCTGCGCACAAACATCCCTGAACCCGGAATCGCGTCGTCGGTGTTAACAATGGTGGTGGTCTGGTAACGATCGTGCGGTAACTTCGCTGATTTAATCGCTTCAATCAGCGTCGCGTTTTTCTCTTTTGCAGAGGTGCGACCAGCAATATGTTGCAGTACTCGCACTTTTCCCACTAACTgc                                                                                                                                                                                                                       <  1:629476/297‑1 (MQ=255)
                                                                       cTTATCCAGCACCACTACAGCGGAACTTTCTTCATCCAGCTGCCAGGCACCGAGCGCGACGCCATTGCTATCGACAATAAACTGCGACCAGGGATAAAGCTTTTTATTACTCTCCAGACTGCTGCGCACAAACATCCCTGAACCCGGAATCGCGTCGTCGGTGTTAACAATGGTGGTGGTCTGGTAACGATCGTGCGGTAACTTCGCTGATTTAATCGCTTCAATCAGCGTCGCGTTTTTCTCTTTTGCAGAGGTGCGACCAGCAATATGTTGCAGTACTCGCACTTTTCCCACTAACTgc                                                                                                                                                                                                                       >  2:420769/1‑301 (MQ=255)
                                                                                 accacTACGGCGGAACTTTTTTCATCCAGTTGCCGGGCACCGGGCGCGACGCCATTGCTATCGACAATAAAATGCGACCAGGGATAAAGCTTTTTATTACTCTCCAGACTGCTGCGCACAAACATCCCTGAACCCGGAATCGCGTCGTCGGTGTTAACAATGGTGGTGGTCTGGTAACGATCGTGCGGTAACTTCGCTGATTTAATCGCTTCAATCAGCGTCGCGTTTTTCTCTTTTGCAGAGGTGCGACCAGCAATATGTTGCAGTACTCGCACTTTTCCCACTAACTGCGCGCTGTTcc                                                                                                                                                                                                             <  1:120930/301‑1 (MQ=255)
                                                                                   cacTACAGCGGAACTTTCTTCATCCAGCTGCCAGGCACCGAGCGCGACGCCATTGCTATCGACAATAAACTGCGACCAGGGATAAAGCTTTTTATTACTCTCCAGACTGCTGCGCACAAACATCCCTGAACCCGGAATCGCGTCGTCGGTGTTAACAATGGTGGTGGTCTGGTAACGATCGTGCGGTAACTTCGCTGATTTAATCGCTTCAATCAGCGTCGCGTTTTTCTCTTTTGCAGAGGTGCGACCAGCAATATGTTGCAGTACTCGCACTTTTCCCACTAACTGCGCGCTGTTCCAg                                                                                                                                                                                                           <  2:577551/301‑1 (MQ=255)
                                                                                       aCAGCGGAACTTTCTTCATCCAGCTGCCAGGCACCGAGCGCGACGCCATTGCTATCGACAATAAACTGCGACCAGGGATAAAGCTTTTTATTACTCTCCAGACTGCTGCGCACAACCATCCCTGAACACGGAATCGCGTCGTCGGTGTTAACAATGGTGGTGGTCTGGTAACGATCGTGCGGTAACTTCGCTGATTTAATCGCTTCAATCAGCGTCGCGTTTTTCTCTTTTGCAGAGGTGCGACCAGCAATATGTTGCAGTACTAGCACTTTTCCAACTCAATGCGCGCAGTTCCAGGttt                                                                                                                                                                                                       >  1:464414/1‑301 (MQ=255)
                                                                                          gCGGAACTTTCTTCATCCAGCTGCCAGGCACCGAGCGCGACGCCATTGCTATCGACAATAAACTGCGACCAGGGATAAAGCTTTTTATTACTCTCCAGACTGCTGCGCACAAACATCCCTGAACCCGGAATCGCGTCGTCGGTGTTAACAATGGTGGTGGTCTGGTAACGATCGTGCGGTAACTTCGCTGATTTAATCGCTTCAATCAGCGTCGCGTTTTTCTCTTTTGCAGAGGTGCGACCAGCAATATGTTGCAGTACTCGCACTTTTCCCACTAACTGCGCGCTGTTCCAGGTTTTga                                                                                                                                                                                                    >  2:206898/1‑300 (MQ=255)
                                                                                                                       cccGAGCGCGATGCCATTGCTATCGACAATAAACTGCGACCAGGGATAAAGCTTTTTATTACTCTCCAGACTGCTGCGCACAAACATCCCTGAACCCGGAATCGCGTCGTCGGTGTTAACAATGGTGGTGGTCTGGTAACGATCGTGCGGTAACTTCGCTGATTTAATCGCTTCAATCAGCGTCGCGTTTTTCTCTTTTGCAGAGGTGCGACCAGCAATATGTTGCAGTACTCGCACTTTTCCCACTAACTGCGCGCTGTTCCAGGTTTTGTAGCTAAACTGATCTTTATCAAGCACCAAc                                                                                                                                                                       <  2:464978/300‑1 (MQ=255)
                                                                                                                       aCCGATCGCGACGCCATTGCTATCGACAATAAACTGCGACCATGGATAAAGCTTTTTATTACTCTCCAGACTGCTGCGCACAAACATCCCTGAACCCGGAATCGCGTCGTCGGTGTTAACAATGGTGGTGGTCTGGTAACGATCGTGCGGTAACTTCGCTGATTTAATCGCTTCAATCAGCGTCGCGTTTTTCTCTTTTGCAGAGGTGCGACCAGCAATATGTTGCAGTACTCGCACTTTTCCCACTAACTGCGCGCTGTTCCAGGTTTTGTAGCTAAACTGATCTTTATCAAGCACCAAc                                                                                                                                                                       <  2:18098/301‑1 (MQ=255)
                                                                                                                                                 cAATAAACTGCGACCAGGGATAAAGCTTTTTATTACTCTCCAGACTGCTGTGCACAAACATCCCTGAACCCGGAATCGCGTCGTCGGTGTTAACAATGGTGGTGGTCTGGTAACGATCGTGCGGTAACTTCGCTGATTCAATCGCTTCAATCAGCGTCGCGTTTTTCTCTTTTGCAGAGGTGCGACCAGCAATATGTTGCAGTACTCGCACTTTTCCCACTAACTGCGCGCTGTTCCAGGTTTTGTAGCCAAACTGATCTTTATCAAGCACCAACTCACCCCGATCGGTAATGCCaaccag                                                                                                                                             >  1:276772/1‑297 (MQ=255)
                                                                                                                                                       aCTGCTACCAGGGTTAAAGCTTTTTTTTACTCTCCAGACTGTTTCGCACTATCATCCCTGATCCCGGTATCGCGTCGTCGGTGTTAACAATGGTGGTGGTCTGGTAACGATCGTGCGGTAACTTCGCTGATTTAATCGCTTCAATCAGCGTCGCGTTTTTCTCTTTTGCAGAGGTGCGACCAGCAATATGTTGCAGTACTCGCACTTTTCCCACTAACTGCGCGCTGTTCCAGGTTTTGTAGCTAAACTGATCTTTATCAAGCACCAACTCGCCCCGATCGGTAATGCCAATCGGCGGCAc                                                                                                                                       <  1:206655/301‑1 (MQ=255)
                                                                                                                                                                        cacgatgctcttccgttctgACTGCTGCGCACAAACATCCCTGAACCCGGAATCGCGTCGTCGGTGTTAACAATGGTGGTGGTCTGGTAACGATCGTGCGGTAACTTCGCTGATTTAATCGCTTCAATCAGCGTCGCGTTTTTCTCTTTTGCAGAGGTGCGACCAGCAATATGTTGCAGTACTCGCACTTTTCCCACTAACTGCGCGCTGTTCCAGGTTTTGTAGCTAAACTGATCTTTATCAAGCACCAACTCGCCCCGATCGGTAATGCCAATCGGCGGCACTCGCTGACCGGTTTCGa                                                                                                                      <  1:265172/282‑1 (MQ=255)
                                                                                                                                                                                   aCTCTCCAGACTGCTGCGCACAAACATCCCTGAACCCGGAATCGCGTCGTCGGTGTTAACAATGGTGGTGGTCTGGTAACGATCGTGCGGTAACTTCGCTGATTTAATCGCTTCAATCAGCGTCGCGTTTTTCTCTTTTGCAGAGGTGCGACCAGCAATATGTTGCAGTACTCGCACTTTTCCCACTAACTGCGCGCTGTTCCAGGTTTTGTAGCTAAACTGATCTTTATCAAGCACCAACTCGCCCCGATCGGTAATGCCAATCGGCGGCACTCGCTGACCGGTTTCGAACTGATGTGCg                                                                                                           <  2:63076/301‑1 (MQ=255)
                                                                                                                                                                                           gACTGCTGCGCACAAACATCCCTGAACCCGGAATCGCGTCGTCGGTGTTAACAATGGTGGTGGTCTGGTAACGATCGTGCGGTAACTTCGCTGATTTAATCGCTTCAATCAGCGTCGCGTTTTTCTCTTTTGCAGAGGTGCGACCAGCAATATGTTGCAGTACTCGCACTTTTCCCACTAACTGCGCGCTGTTCCAGGTTTTGTAGCTAAACTGATCTTTATCAAGCACCAACTCGCCCCGATCGGTAATGCCAATCGGCGGCACTCGCTGACCGGTTTCGaagatcggaagagcacacgt                                                                                                   >  2:265440/1‑283 (MQ=255)
                                                                                                                                                                                             ctgctgCGCACAAACATCCCTGAACCCGGAATCGCGTCGTCGGTGTTAACAATGGTGGTGGTCTGGTAACGATCGTGCGGTAACTTCGCTGATTTAATCGCTTCAATCAGCGTCGCGTTTTTCTCTTTTGCAGAGGTGCGACCAGCAATATGTTGCAGTACTCGCACTTTTCCCACTAACTGCGCGCTGTTCCAGGTTTTGTAGCTAAACTGATCTTTATCAAGCACCAACTCGCCCCGATCGGTAATGCCAATCGGCGGCCATCGCTGACCGGTTTCGAACTGATGTGCGGAAGcaaacc                                                                                                 >  1:175781/1‑296 (MQ=255)
                                                                                                                                                                                                     cacaAACATCCCTGAACCCGGAATCGCGTCGTCGGTGTTAACAATGGTGGTGGTCTGGTAACGATCGTGCGGTAACTTCGCTGATTTTATCGCTTCAATCAGCGTCGCGTTTTTCTCTTTTGCAGAGGTGCGACCAGCAATATGTTGCAGTACTCGCACTTTTCCCACTAACTGCGCGCTGTTCCAGGTTTTGTAGCTAAACTGATCTTTATCAAGCACCAACTCGCCCCGATCGGTACTGCCAAACGGCGGCACACACTGACAGGTTTAGAACTAATGTACGGAAGAcatcccacgcaaa                                                                                         >  1:269080/1‑292 (MQ=255)
                                                                                                                                                                                                            aTCCCTGAACCCGGAATCGCGTCGTCGGTGTTAACAATGGTGGTGGTCTGGTAACGATCGTGCGGTAACTTCGCTGATTTAATCGCTTCAATCAGCGTCGCGTTTTTCTCTTTTGCAGAGGTGCGACCAGCAATATGTTGCAGTACTCGCACTTTTCCCACTAACTGCGCGCTGTTCCAGGTTTTGTAGCTAAACTGATCTTTATCAAGCACCAACTCGCCCCGATCGGTAATGCCAATCGGCGGCACTCGCTGACCGGTTTCGAACTGATGTGCGGAAGCCATCATCGGCAACAgcaggc                                                                                  >  1:101328/1‑301 (MQ=255)
                                                                                                                                                                                                              cccTGAACCCGGAATCGCGTCGTCGGTGTTAACAATGGTGGTGGTCTGGTAACGATCGTGCGGTAACTTCGCTGATTTAATCGCTTCAATCAGCGTCGCGTTTTTCTCTTTTGCAGAGGTGCGACCAGCAATATGTTGCAGTACTCGCACTTTTCCCACTAACTGCGCGCTGTCCCAGGTTTTGTAGCTAAACTGATCTTTATCAAGCACCAACTCGCCCCGATCGGTAATGCCAATCGGCAGCACTCGCTGACCGGTTTCGAACTGATGTGAGGAAGCAATCATCGGCAACAgcaggccg                                                                                >  1:160899/1‑299 (MQ=255)
                                                                                                                                                                                                               ctttccctacagtgcgctctttcgttcttAAAAATGGTGGTGGTCTGGTAACGATCGTGCGGTTACTTCGCTGATTTAATCGCTTCAATCAGCGTCGCGTTTTTCTCTTTTGCAGAGGTGCGACCAGCAATATGTTGCAGTACTCGCACTTTTCCCACTAACTGCGCGCTGTTCCAGGTTTTGTAGCTAAACTGATCTTTATCAAGCACCAACTCGCCCCGATCGGTAATGCCAATCGGCGGCACTCGCTGACCGGTTTCGAACTGATGTGCGGAAGCCATCATCGGCAACAGcaggcagg                                                                               <  1:546259/274‑1 (MQ=255)
                                                                                                                                                                                                               ccTGAACCCGGAATCGCGTCGTCGGTGTTAACAATGGTGGTGGTCTGGTAACGATCGTGCGGTAACTTCGCTGATTTAATCGCTTCAATCAGCGTCGCGTTTTTCTCTTTTGCAGAGGTGCGACCAGCAATATGTTGCAGTACTCGCACTTTTCCCACTAACTGCGCGCTGTTCCAGGTTTTGTAGCTAAACTGATCTTTATCAAGCACCAACTCGCCCCGATCGGTAATGCCAATCGGCGGCACTCGCTGACCGGTTTCGAACTGATGTGCGGAAGCCAACATCGGCAACAGcaggcagg                                                                               >  2:84981/1‑301 (MQ=255)
                                                                                                                                                                                                                 tGAACCCGGAATCGCGTCGTCGGTGTTAACAATGGTGGTGGTCTGGTAACGATCGTGCGGTAACTTCGCTGATTTAATCGCTTCAATCAGCGTCGCGTTTTTCTCTTTTGCAGAGGTGCGACCAGCAATATGTTGCAGTACTCGCACTTTTCCCACTAACTGCGCGCTGTTCCAGGTTTTGTAGCTAAACTGATCTTTATCAAGCACCAACTCGCCCCGATCGGTAATGCCAATCGGCGGCACTCGCTGACCGGTTTCGAACTGATGTGCGGAAGCCATCATCGGCAACAGCAGGCAGGTg                                                                             <  2:101452/301‑1 (MQ=255)
                                                                                                                                                                                                                          aaTCGCGTGTTCGGTGTTAATATTGTTGGTGGTCTGTTACCGATCTTGCGGTAACTTTGCTGATTTATTCGCTTCAATCAGCGTCGCGTTTTTCTCTTTTGCAGAGGTGCGACCAGCATTATGTTGCAGTACTCGCACTTTTCCCACTAACTGCGCGCTGTTCCAGGTTTTGTAGCTAAACTGATCTTTATCAAGCACCAACTCGCCCCGATCGGTAATGCCAATCGGCGGCACTCGCTGACCGGTTTCGAACTGATGTGCGGAAGCCATCATCGGCAACAGCAGGCAGGTGAGTGCCAGa                                                                    <  1:144834/301‑1 (MQ=255)
                                                                                                                                                                                                                              gcgTCGTCGGTGTTAACAATGGTGGTGGTCTGGTAACGATCGTGCGGTAACTTCGCTGATTTAATCGCTTCAATCAGCGTCGCGTTTTTCTCTTTTGCAGAGGTGCGACCAGCAATATGTTGCAGTACTCGCACTTTTCCCACTAACTGCGCGCTGTTCCAGGTTTTGTAGCTACACTGATCTTTATCAAGCACCAACTCGCCCCGATCGGTAATGCCAATCGGCGGCACTCGCTGACCGGTTTCGAACTGATGTGCGGAAGCCATCATCGGCAACAGCAGGCAGGTACGTGCCAGAATCt                                                                >  2:538436/1‑301 (MQ=255)
                                                                                                                                                                                                                                           tAACAATGGTGGTGGTCTGGTAACGATCGTGCGGTAACTTCGCTGCTTTAATCGCTTCAATCAGCGTCGCGTTTTTCTCTTTTGCAGATGTGCGACCAGCAATATGTTGCAGTACTCGCACTTATCCCACTAACTGCGCGCTGTTCCAGGTTTTGTAGCTAAACTGATCTTCATCAAGCACCAACTCGCCCCGATTGGTAATGCCCATCGGCGGCACTCGCGGACCGGGTTCGAACGGATGTGAGGAAGCAATCATAGGAAACAGCAGGACGGTAACTGCCAGAATCTGGAGTGAGGTCag                                                   >  1:513078/1‑300 (MQ=255)
                                                                                                                                                                                                                                           tAACAATGGTGGTGGTCTGGTAACGATCGTGCGGTAACTTCGCTGATTTAATCGCTTCAATCAGCGTCGCGTTTTTCTCTTTTGCAGAGGTGCGACCAGCAATATGTTGCAGTACTCGCACTTTTCCCACTAACTGCGCGCTGTTCCAGGTTTTGTAGCTAAACTGATCTTTATCAAGCACCAACTCGCCCCGATCGGTAATGCCAATCGGCGGCACTCGCTGACCGGTTTCGAACTGATGTGCGGAAGCCATCATCGGCAACAGCAGGCAGgagatcggaagagcaaacgtcggaactcc                                                   >  2:546932/1‑276 (MQ=255)
                                                                                                                                                                                                                                                    aggtggtCTGGTAACGATCGTGCGGTAACTTCGCTGATTTAATCGCTTCAATCAGCGTCGCGTTTTTCTCTTTTGCAGAGGTGCGACCAGCAATATGTTGCAGTACTCGCACTTTTCCCACTAACTGCGCGCTGTTCCAGGTTTTGTAGCTAAACTGATCTTTATCAAGCACCAACTCGCCCCGATCGGTAATGCCAATCGGCGGCACTCGCTGACCGGTTTCGAACTGATGTGCGGAAGCCATCATCGGCAACAGCAGGCAGGTGAGTGCCAGAACCCTTCATAAGGTCCAGGTGTTTcc                                          >  1:382278/2‑301 (MQ=255)
                                                                                                                                                                                                                                                     tgttctCTGGTAACGATCGTGCGGTAACTTCGCTGATTTAATCGCTTCAATCAGCGTCGCGTTTTTCTCTTTTGCAGAGGTGCGACCAGCAATATGTTGCAGTACTCGCACTTTTCCCACTAACTGCGCGCTGTTCCAGGTTTTGTAGCTAAACTGATCTTTATCAAGCACCAACTCGCTCCGATCGGTAATGCCAATCGGCGGCACTCGCTGACCGGTTTCGAACTGATGTGCGGAAGCCATCATCGGCAACAGCAGGCAGGTGAGTGCCAGAATCTTGCGTAGGGTCATGGTGTTTCct                                         <  1:93779/296‑1 (MQ=255)
                                                                                                                                                                                                                                                           cTGGTAACGATCGTGCGGTAACTTCGCTGATTTAATCGCTTCAATCAGCGTCGCGTTTTTCTCTTTTGCAGAGGTGCGACCAGCAATATGTTGCAGTACTCGCACTTTTCCCACTAACTGCGCGCTGTTCCAGGTTTTGTAGCTAAACTGATCTTTATCAAGCACCAACTCGCTCCGATCGGTAATGCCAATCGGCGGCACTCGCTGACCGGTTTCGAACTGATGTGCGGAAGCCATCATCGGCAACAGCAGGCAGGTGAGTGCCAGAATCTTGCGTAGGGTCATGGTGTTTCctagatcg                                   >  2:93903/1‑295 (MQ=255)
                                                                                                                                                                                                                                                            tggtactgttggTGCGGTAACTTCTCTGATTTAATCGCTTCAATCAGCGTCGCGTTTTTCTCTTTTGCAGAGGTGCGACCAGCAATATGTTGCAGTACTCGCACTTTTCCCACTAACTGCGCGCTGTTCCAGGTTTTGTAGCTAAACTGATCTTTATCAAGCACCAACTCGCCCCGATCGGTAATGCCAATCGGCGGCACTCGCTGACCGGTTTCGAACTGATGTGCGGAAGCCATCATCGGCAACAGCAGGCAGGTGAGTGCCAGAATCTTGCGTAGGGTCATGGTGTTTCCTTCTTATg                                  <  1:84857/290‑1 (MQ=255)
                                                                                                                                                                                                                                                              gTAACGATCGTGCGGTAACTTCGCTGATTTAATCGCTTCAATCAGCGTCGCGTTTTTCTCTTTTGCAGAGGTGCGACCAGCAATATGTTGCAGTACTCGCACTTTTCCCACTAACTGCGCGCTGTTCCAGGTTTTGTAGCTAAACTGATCTTTATCAAGCACCAACTCGCCCCGATCGGTAATGCCAATCGGCGGCACTCGCTGACCGGTTTCGAACTGATGTGCGGAAGCCATCATCGGCAACAGCAGGCAGGTGAGGGCCAGAATCTTGCGTAGGGGCATGGTGTTTCCTTCTTATGat                                >  2:623927/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                     tCGTGCGGTAACTTCGCTGATTTTATCGCTTCAATCAGCGTCGCGTTTTTCTCTTTTGCAGAGGTGCGACCAGCAATATGTTGCAGTACTCGCACTTTTCCCACTAACTGCGCGCTGTTCCAGGTTTTGTAGCTAAACTGATCTTTATCAAGCACCAACTCGCCCCGATCGGTAATGCCAATCGGCGGCACTCGCTGACCGGTTTCGAACTGATGTGCGGAAGCCATCATCGGCAACAGCAGGCAGGTGAGTGCCAGAATCTTGCGTAGGGTCATGGTGTTTCCTTCTTATGATATGCAgg                         <  2:269378/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                           ggTAACTTCGCTGGTTTATTCTCTTCAATCTGCTTCTCGTTTTTCTCTTTTGCAGAGGTGCGACCAGCAATATGTTGCAGTACTCGCACTTTTCCCACTAACTGCGCGCTGTTCCAGGTGTTGTAGCTAAACTGATCTTTATCAAGCACCAACTCGCCCCGATCGGTAATGCCAATCGGCGGCACTCGCTGACCGGTTTCGAACTGATGTGCGGAAGCCATCATCGGCAACAGCAGGCAGGTGAGTGCCAGAATCTTGCGTAGGGTCATGGTGTTTCCTTCTTATGATATGCAGGTGATcc                   <  1:420206/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                cgacGCTGATTGAAGCGATTAAATCAGCGTCGCGTTTTTCTCTTTTGCAGAGGTGCGACCAGCAATATGTTGCAGTACTCGCACTTTTCCCACTAACTGCGCGCTGTTCCAGGTTTTGTAGCTAAACTGATCTTTATCAAGCACCAACTCGCCCCGATCGGTAATGCCAATCGGCGGCACTCGCTGACCGGTTTCGAACTGATGTGCGGAAGCCATCATCGGCAACAGCAGGCAGGTGAGTGCCAGAATCTTGCGTAGGGTCATGGTGTTTCCTTCTTATGATATGCAGGTGATCCGACCa              >  2:259049/4‑301 (MQ=255)
                                                                                                                                                                                                                                                                                      ttaTTTACTCGCTTCCATCAGCGTCGCGTTTTTCTCTTTTGCAGAGGTGCGTCCAGCAATATGTTGCAGTACTCGCACTTTTCCCACTAACTGCGCGCTGTTCCAGGTTTTGTAGCTAAACTGATCTTTATCAAGCACCAACTCGCCCCGATCGGTAATGCCAATCGGCGGCACTCGCTGACCGGTTTCGAACTGATGTGCGGAAGCCATCATCGGCAACAGCAGGCAGGTGAGTGCCAGAATCTTGCGTAGGGTCATGGTGTTTCCTTCTTATGATATGCAGGTGATCCGACCACTTggg        <  1:623068/299‑1 (MQ=255)
                                                                                                                                                                                                                                                                                      ggATTGAAGCGATTAAATCAGCGTCGCGTTTTTCTCTTTTGCAGAGGTGCGACCAGCAATATGTTGCAGTACTCGCACTTTTCCCACTAACTGCGCGCTGTTCCAGGTTTTGTAGCTAAACTGATCTTTATCAAGCACCAACTCGCCCCGATCGGTAATGCCAATCGGCGGCACTCGCTGACCGGTTTCGAACTGATGTGCGGAAGCCATCATCGGCAACAGCAGGCAGGTGAGTGCCAGAATCTTGCGTAGGGTCATGGTGTTTCCTTCTTATGATATGCAGGTGATCCGACCACTTggg        <  1:258781/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                                        aTTTAATCGCTTCAATCAGCGTCGCGTTTTTCTCTTTTGCAGAGGTGCGACCAGCAATATGTTGCAGTACTCGCACTTTTCCCACTAACTGCGCGCTGTTCCAGGTTTTGTAGCTAAACTGATCTTTATCAAGCACCAACTCGCCCCGATCGGTAATGCCAATCGGCGGCACTCGCTGACCGGTTTCGAACTGATGTGCGGAAGCCATCATCGGCAACAGCAGGCAGGTGAGTGCCAGAATCTTGCGTAGGGTCATGGTGTTTCCTTCTTATGATATGCAGGTGATCCGACCACTTGGGTc      >  1:146040/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                          ttAATCGCTTCAATCAGCGTCGCGTTTTTCTCTTTTGCAGAGGTGCGACCAGCAATATGTTGCAGTACTCGCACTTTTCCCACTAACTGCGCGCTGTTCCAGGTTTTGTAGCTAAACTGATCTTTATCAAGCACCAACTCGCCCCGATCGGTAATGCCAATCGGCGGCACTCGCTGACCGGTTTCGAACTGATGTGCGGAAGCCATCATCGGCAACAGCAGGCAGGTGAGTGCCAGAATCTTGCGTAGGGTCATGGTGTTTCCTTCTTATGATATGCAGGTGATCCGACCACTTGGGTCTg    <  2:140456/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                            taTCGCTTCAATCAGCGTCGCGTTTTTCTCTTTTGCAGAGGGGCGACCAGCAATATGTTGCAGTACTCGCACTTTTCCCACTAACTGCGCGCTGTTCCAGGTTTTGTAGCTAAACTGATCTTTATCAAGCACCAACTCGCCCCGATCGGTAATGCCAATCGGCGGCACTCGCTGACCGGTTTCGAACTGATGTGCGGAAGCCATCATCGGCAACAGCAGGCAGGTGAGTGCCAGAATCTTGCGTAGGGTCATGGTGTTTCCTTCTTATGATATGCAGGTGATCCGACCACTTGGGTCTGAg  <  2:382774/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                                                          
AGCAGGTCCATCACTTGCTGCACCTCTTCCGGTGTGAGCGCCCCGTCTTTGGCCCATTGCACGCGACCATCCTTATCCAGCACCACTACAGCGGAACTTTCTTCATCCAGCTGCCAGGCACCGAGCGCGACGCCATTGCTATCGACAATAAACTGCGACCAGGGATAAAGCTTTTTATTACTCTCCAGACTGCTGCGCACAAACATCCCTGAACCCGGAATCGCGTCGTCGGTGTTAACAATGGTGGTGGTCTGGTAACGATCGTGCGGTAACTTCGCTGATTTAATCGCTTCAATCAGCGTCGCGTTTTTCTCTTTTGCAGAGGTGCGACCAGCAATATGTTGCAGTACTCGCACTTTTCCCACTAACTGCGCGCTGTTCCAGGTTTTGTAGCTAAACTGATCTTTATCAAGCACCAACTCGCCCCGATCGGTAATGCCAATCGGCGGCACTCGCTGACCGGTTTCGAACTGATGTGCGGAAGCCATCATCGGCAACAGCAGGCAGGTGAGTGCCAGAATCTTGCGTAGGGTCATGGTGTTTCCTTCTTATGATATGCAGGTGATCCGACCACTTGGGTCTGAG  >  NC_000913/4438727‑4439311

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