Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,125,339 T→G 16.0% T491P (ACC→CCC)  priA ← Primosome factor n' (replication factor Y)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,125,3390TG16.0% 141.1 / 2.6 75T491P (ACC→CCC) priAPrimosome factor n' (replication factor Y)
Reads supporting (aligned to +/- strand):  ref base T (24/34);  new base G (7/9);  total (32/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATGGTGACCTTTCGCCAGCATTTGTGTACCAATCAAAATCCGCGCGCCGCCGCGATGTACTTCTGCCAGTTGCTGTTCCAGCGCCCCTTTGCGGCTGGTGGTATCGCGGTCGATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGGACCAGGTGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCACCGGACGCTGACTGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCGTGGCACAGCAGTGCAGGCGCAAAGCCACGGCGGTTAAGAAAGAGAATGACCTGGTTATCAGCCT  >  NC_000913/4125115‑4125563
                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                 
aGAAAACAGCCCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATGGTGACCTTTCGCCAGCATTTGTGTTCCAATCAAAATCCGCGCGCCGCCGCGATGTACTTCTGCCAGTTGCTGTTCCAGCGCCCCTTTGCGGCTGGTGGTATCGCGGTCGATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGcc                                                                                                                                                                                                                                <  1:595133/228‑1 (MQ=255)
  aaaaaggagccgggtccctgtccgGTAATGCAACCAGGGTCCCGTCGGGAAAATGGTGACCTTTCGCCAGCATTTGGCTACCAATCAAAACCCGCGCGCCGCCGCGAGGTTCTTCTGCCAGTTGCTGTTCCAGCGCCCCTTTGCGGCTGGTGGTATCGCGGTCGATACGTGAAATGGGCGCGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCa                                                                                                                                                                                                                              <  2:553012/209‑1 (MQ=255)
          gcgcCGTCCACGTCCAGTAATGCACCCAGCGTCACATCCGGGAAATGGTGCCCTTTTGCCAGCATTTGTGTACCAATCAAAATCCGCGCGCCGCCGCGATGTACTTCTGCCAGTTGCTGTTCCAGCGCCCCTTTGCGGCTGGTGGTATCGCGGTCGATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACg                                                                                                                                                                                                                      <  1:603475/228‑1 (MQ=255)
              cgtccacgtccaGTAATGCAACCAGCGTCACATCCGGGAAATGGTGACCTTTCGCCAGCATTTGTGTACCAATAAAAATCCGCGAGCGGCCGCGATGTACTTCTGCCAGTTGCTGTTCCAGCGCCCCTTTGCGGCTGGTGGAATGCCGGTCGTTACCGGAAATGGGCACTCCCGGGCAAAACGGCCGGGGCGCTCGGTCAAGCTGGTGGGCGACCACCCCCACCGGGAc                                                                                                                                                                                                                >  2:230901/1‑229 (MQ=255)
                tccacgtccaGTAATGCAACCAGCGTCACATCCGGGAAATGGTGACCTTTCGCCAGCATTTGTGTACCAATCAAAATCCGCGCGCCGCCGCGATGTACTTCTGCCAGTTGCTGTTCCAGCGCCCCTTTGCGGCTGGTGGGATCGCGGTCGATACGAGAAATGGGCACGCCGGGGAACCACGCCGCGGGGGTCTGTTCAAACGGTTCGGGGCCCAGCCCCCCggggcgac                                                                                                                                                                                                               >  2:549629/1‑225 (MQ=255)
                 ccccgcccgttatggacaCCAGCGCCCCATCCGGGAAATGGGGACCTTTCGCCAGCATTTGTGTACCAATCAAAATCCGCGCGCCGCGGGGATGTACTTCTGCCAGTTGCTGTTCCAGCGCCCCTTTGCGGCTGGTGGTATCGCGGTCGATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGGACCAg                                                                                                                                                                                                             <  2:696355/212‑1 (MQ=255)
                 cccacgcccgTAATTCAACCACGGCCACAGCCGGGAAAGGGTGCCCTTGTGCCAGCATTTGTGTACCAAGCAAAATCCGCGCGCCGCCGCGATGTACTTCTGCCAGTTGCTGTTCCAGCGCCCCTTTGCGGCTGGTGGTATCGCGGTCGATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGGACCAg                                                                                                                                                                                                             <  2:188604/220‑1 (MQ=255)
                              tGCAACCAGCGTCACATCCGGGAAATGGTGACCTTTCGCCAGCATTTGTGTACCAATCAAAATCCGCGCGCCGCCGCGATGTACTTCTGCCAGTTGCTGTTCCAGCGCCCCTTTGCGGCTGGTGGTATCGCGGTCGATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGGGCGTCTGTTCAAGCTGTTCGGGGCCCGGCCCCCCGGGGGCACGGTGGGGGGAAccc                                                                                                                                                                                                >  1:314478/1‑228 (MQ=255)
                               gCAACCAGCGTCACATCCGGGAAATGGTGACCTTTCGCCAGCATTTGTGTACCAATCAAAATCCGCGCGCCGCCGCGATGTACTTCTGCCAGTTGCTGTTCCAGCGCCCCTTTGCGGCTGGTGGTATCGCGGTCGATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCCCGGGGACCAGGTGCGTGGAACCg                                                                                                                                                                                                >  1:117250/1‑228 (MQ=255)
                                cAACCAGCGTCACATCCGGGAAATGGTGACCTTTCGCCAGCATTTGTGTACCAATCAAAATCCGCGCGCCGCCGCGATGTACTTCTGCCAGTTGCTGTTCCAGCGCCCCTTTGCGGCTGGTGGTATCGCGGTCGATACGAGAAATGGGCACGCCGGGGGAGAACGGCGCGAGCGTCCGTTCAAGCTGGTCGGTGGCCAACCCCCCGGGGGCCAGGGGCGGGGCAccacc                                                                                                                                                                                              >  1:320056/1‑226 (MQ=255)
                                   ccAGCGTCACATCCGGGAAATGGTGACCTTTCGCCAGCATTTGTGTACCAATCAAAATCCGCGCGCCGCCGCGATGTACTTCTGCCAGTTGCTGTTCCAGCGCCCCTTTGCGGCTGGTGGTATCGCGGTCGATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGGCCAGCCCCCCGGGGACAAGGGACGTGGAACCGCAggg                                                                                                                                                                                           >  1:753200/1‑228 (MQ=255)
                                     aGCGTCACATCCGGGAAATGGTGACCTTTCGCCAGCATTTGTGTACCAATCAAAATCCGCGCGCCGCCGCGATGTACTTCTGCCAGTTGCTGTTCCAGCGCCCCTTTGCGGCTGGTGGAATCGCGGCCGATACGCGAAATGGGCACGCCGGGGAAAAACGGGGCGCGCGTCTTGTTAAGCGTTACGGGGCCCACCCCCACGGGGCACCCGGGCGGGGGACCGCacagaa                                                                                                                                                                                         >  2:448306/1‑224 (MQ=255)
                                      ccgtccaaTCCGGGAAAGTGTGCCCTTTCTCCAGCATTTGTGGCCCAACCAAACCCCGCCGGCGGGCGGGATGTACTGCTGCCAGTTGCTGTGCCAGCGCCCCTTTGCGGCGGGGGGTATCGCGGTCGATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGGACCAGGTGCGTGGAACCGCAGGa                                                                                                                                                                                           <  2:320056/219‑1 (MQ=255)
                                       cGTCACATCCGGGAAATGGTTACCTTTCGCCAGCATTTGTGTACCAATCAAAATCCGCGCGCCGCCGCGATGTACTTCTGCCAGTTGCTGTTCCAGCGCCCCTTTGCGGCTGGTGGTATCGCGGTCGATACGAGAAATGGGCACGCCGGGGAACAACGGCCGGAGGGCCGGTTCCAGCTGTTCGGTGCCCACCCCCCCGCGGCCCGGGGGCGGGGAAcacccgggaggg                                                                                                                                                                                       >  2:405624/1‑218 (MQ=255)
                                       cGTCACATCCGGGAAATGGTGACCTTTCGCCAGCATTTGTGTACCAATCAAAATCCGCGCGCCGCCGCGATGTACTTCTGCCAGTTGCTGTTCCAGCGCCCCTTTGCGGCTGGTGGTATCGCGGTCGATACGAGAAATGGGAACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGGGCCCAGCCCCACGGGGACCAGGGGCGGGGAACCGCAGGAAagg                                                                                                                                                                                       >  1:291864/1‑229 (MQ=255)
                                         tCACATCCGGGAAATGGTGACCTTTCGCCAGCATTTGTTTACCAATCAAAATCCGCGCGCCGCCGCGATGTACTTCTGCCAGTTGCTGTTCCAGCGCCCCTTTGCGGCTGGTGGTATCGCGGGCGATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGGGCGGCTGTTCAAACTGTTCGGTGCCCAGCCCCCCGGGGGCCAggggggggggacacgggaaggggc                                                                                                                                                                                      >  2:545400/1‑206 (MQ=255)
                                         tCACATCCGGGAAATGGTGACCTTTCGCCAGCATTTGTGTACCAATCAAAATCCGCGCGCCGCCGCGATGTACTTCTGCCAGTTGCTGTTCCAGCGCCCCTTTGCGGCTGGTGGTATCGCGGTCGATACGAGAAATGGGCCCGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGGGCCCAGCCCCCCGGGGACCAGGGGCGCGGAACCGCAGGAAGGGc                                                                                                                                                                                      >  1:181454/1‑228 (MQ=255)
                                                     aaTGGTGGCCTTTCCCCAGCATTTGGGTACCAATCAAAAGCCGCGCGCCGCCGCGATGGACTTCTGCCAGTTGCTGTTCCAGCGCCCCTTTGCGGCTGGTGGTATCGCGGTCGATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGGACCAGGTGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCAc                                                                                                                                                                         <  1:549629/229‑1 (MQ=255)
                                                         tgcgctcttCGCCAGCATTTTGTGACAAACCGAAACCCCCGTCCCCCCGGGATGTAGTCTGCCCAGTGGTTGTACCAGGGCCCCTGTGCGGCTGGGGGTATCGCGGTCGATGCGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGGTCAAGCTGTTCGGTGACCAGCCCCACGGGGACCAGGTGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCACCTGa                                                                                                                                                                     <  2:456046/222‑1 (MQ=255)
                                                         gTGACCTTTCGCCAGCATTTGTGTACCAATCAAAATCCGCGCGCCGCCGCGATGTACTTCTGCCAGTTGCTGTTCCAGCGCCCCTTTGCGGCTGGTGGTATCGCGGTCTATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGGACCCGGTGCGTGGAACCGCAGGAAGGGCACTGACGCGCGACCgg                                                                                                                                                                      >  1:665345/1‑228 (MQ=255)
                                                            aCCTTTCGCCAGCATTTGTGTACCAATCAAAATCCGCGCGCCGCCGCGATGTACTTCTGCCAGTTGCTGTTCCAGCGCCCCTTTGCGGCTGGTGGTATCGCGGTCGATACGGCAAATGGGCCCGCCGGGGCACAACGGCCCGAGCGCCTGTTCAAGCTCTTCGGGGCCCACCCCCCCGGGGGCCAGGTGCGCTGCACCGCCGGAAGGGCCCTGCGCGGGAcccggaccc                                                                                                                                                                 >  2:113782/1‑222 (MQ=255)
                                                                   ggcatctttttTGGACCAACCAAAAGCCGCGCGCCGCCGCGATGTTCTTCTGCCAGTTGCTGTTCCCGGGCCCCTTTGCGGCTGGGGGTATCGCGGTCGATACGAGAAAGGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGGACCAGGTGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCACCGGACGCTGACTGt                                                                                                                                                           <  2:659650/222‑1 (MQ=255)
                                                                                   ccAATCAAAATCCGCGCGCCGCCGCGATGTACTTCTGCCAGTTGCTGTTCCAGCGCCCCTTTGCGGCTGGTGGTATCGCGGTCGATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGcc                                                                                                                                                                                                                                >  2:158943/1‑146 (MQ=255)
                                                                                   ccAATCAAAATCCGCGCGCCGCCGCGATGTACTTCTGCCAGTTGCTGTTCCAGCGCCCCTTTGCGGCTGGTGGTATCGCGGTCGATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGcc                                                                                                                                                                                                                                <  1:158943/146‑1 (MQ=255)
                                                                                    cAATCAAAATCCGCGCGCCGCCGCGATGTACTTCTGCCAGTTGCTGTTCCAGCGCCCCTTTGCGGCTGGTGGTATCGCGGTCGATACGAGAAATGGGCACGCCGGGGACCAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGGACCAGGCGCGTGGAACCGCAGGAAGGGCACTGGCCCGGCACCGGCCCCTGCCGGGCCCAGTGGGGGcgcgc                                                                                                                                           >  1:692692/1‑225 (MQ=255)
                                                                                      aTCAAAATCCGCGCGCCGCCGCGATGTACTTCTGCCAGTTGCTGTTCCAGCGCCCCTTTGCGGCTGGTGGTATCGCGGTCGATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGGCCCACCGGGGCCCGGGGCGTGGAACCGCAGGCAGGGCAACGGCGCGGCCCCGGCCGCCGAcgtccaccgtcgggggagcgga                                                                                                                                         >  2:705918/1‑208 (MQ=255)
                                                                                      aTCAAAATCCGCGCGCCGCCGCGATGTACTTCTGCCAGTTGCTGTTCCAGCGCCCCTTTGCGGCTGGTGGTATCGCGGTCGATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACg                                                                                                                                                                                                                      >  1:261685/1‑153 (MQ=255)
                                                                                      aTCAAAATCCGCGCGCCGCCGCGATGTACTTCTGCCAGTTGCTGTTCCAGCGCCCCTTTGCGGCTGGTGGTATCGCGGTCGATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGGGCCAGGTGCGTGGAACCGCAGGAAGGGCACCGGCCCGGCACCGGCCGCTGCCCGGCACAGGGGGGGGAccgcc                                                                                                                                         >  2:731227/1‑228 (MQ=255)
                                                                                      aTCAAAATCCCCGCGCCGGCGGGAGGGACTTCGTGCAGTTGCTGTTCCCGCGCTCCTTTGCGGCTGGTGGTATCGCGGGCGATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGGACCAGGTGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCACCGGACGCTGACTGTCACAGTGGTGGCAGCGCa                                                                                                                                         <  2:564170/229‑1 (MQ=255)
                                                                                         aaaaTCCGCGCGCCGCCGCGATGTACTTCTGCCAGTTGCTGTTCCAGCGCCCCTTTGCGGCTGGTGGTATCGCGGTCGATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGGACCAGGTGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCACCGGACGCTGACTGTCACAGTGGTGGCAGCGCAGa                                                                                                                                       >  1:23059/1‑228 (MQ=255)
                                                                                                 cgcgcCGCCGCGATGTACTTCTGCCAGTTGCTGTTCCAGCGCCCCTTTGCGGCTGGTGGTATCGCGGTCGATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGGCCCAGGTGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCACCGGACGCCGCCTGTCACAGTGGTGGCAGCGCAGCTGGggctgc                                                                                                                              >  1:522759/1‑229 (MQ=255)
                                                                                                          gggagttACTTCTGCCCGTTGCTGTTCCAGGGGCCCTTGGCGGCTGGTTGTATCGCGGTCGATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGGACCAGGTGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCACCGGACGCTGACTGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGc                                                                                                                     <  1:414161/223‑1 (MQ=255)
                                                                                                           ggtgtgtcTTGTGCCGTTTGCGGTGCCAGCGCCCCTTGGGGGGTGGTGGTATCGGGGTCGATACGAGAAATGGGCCCGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGGACCAGGTGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCACCGGACGCTGACTGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCa                                                                                                                    <  2:117250/222‑1 (MQ=255)
                                                                                                           gggggtctgtgggCCGGTTGGTGTTCCGTCGCCCCTTGGCGGCTGGGGGTGTCGCGGTGGTTCCGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGGACCAGGTGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCACCGGACGCTGACTGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCa                                                                                                                    <  2:402258/217‑1 (MQ=255)
                                                                                                                                        gcgcCCCTTTGCGGCTGGTGGTATCGCGGTCGATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGGACCAGGTGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCACCGGACGCTGACTGTCACAGGGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGGGTAGTAGGGATGGCAACGTGGGcg                                                                                        >  1:259742/1‑228 (MQ=255)
                                                                                                                                           ccccTTTGCGGCTGGTGGTATCGCGGTCGATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGGCCCAGGTGCGTGGAACCGCAGGAAGGGCACTCGCCCGGCCCCGGACGCCGACTCTCACCGGGGGGGCCGCCCGGGTGCGGCTGCGGCGTATGCCGCGTTGTGGAGTGACCTCCACGTGGGCAtt                                                                                      >  2:22113/1‑228 (MQ=255)
                                                                                                                                                 gtgggCTGGGGGTTCCGCGGGCGAGGAGTGAAATGGGCCCGCCGGGGAACAACGGCGGGGGCGTCTGGTCAAGCTGTGCGGGGCCCAGCCCCACGGGGACCAGGTGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCACCGGACGCTGACTGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAAt                                                                               <  2:663032/226‑1 (MQ=255)
                                                                                                                                                  gCGGGTGGTGGGTTCGGGGTCGATACGAGAAAGGGCCACGCCGGGGAACATCGGCGCGGGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCGCGGGGACCAGGTGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCCCCGGACGCTGACTGTCACAGTGGTGGCCGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATc                                                                              <  2:314478/229‑1 (MQ=255)
                                                                                                                                                       tggtggTATCGCGGTCGATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGGGCCCAGCCCCACGGGGACCAGGGGCGTGGAACCGCAGGAAGGGCACCGGCCCGGGCCCGGACGCTGACTGTCCCAGGGGGGGCAGCCCAGATGGGGCGGCCCGTCATGCCGCGTGGAGTAGGGATTGCAACGAGGGCATTTgggcacacagga                                                                         >  2:14409/1‑219 (MQ=255)
                                                                                                                                                        ggtggtATCGGGGGCGATACGAGAAATGGGCCCGCCGGGGAACAGCGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGGACCAGGTGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCACCGGACGCTGACTGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCca                                                                        <  1:104844/229‑1 (MQ=255)
                                                                                                                                                        ggtggtATCGCGTGCGAATCGGAAAACGGTCACGCCGGGGAAAAACGGGGCGTGGGTCTGTTCAAGGTGGTCGGTGCCCAGCCCCCGCGGGACCAGGTGCGTGGAACCGCAGGAAGGGCCCTGGCGCGGCCCCGGACGCTGACTGTCCCAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCca                                                                        <  2:205601/229‑1 (MQ=255)
                                                                                                                                                         gtggtATCGCGGTCGATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGGACCAGGTGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCACCGGACGCTGACTGTCACAGTGGTGGCAGCGCAGAGGGTGCTGGGCCTGATGCAGCGGGTGGTGGTGATCGCAACGTGGGCATTCGGCAATCCAGCcac                                                                       >  1:262276/1‑229 (MQ=255)
                                                                                                                                                         gtggtAGCGGGGTCGATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGGACCAGGTGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCACCGGACGCTGACTGTCACAGTGGGGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCcac                                                                       <  1:607104/229‑1 (MQ=255)
                                                                                                                                                         ggggtggcggggtggaggcGAGAAATGGGCACGCCGGGGAACAACGGCGCGGGCGTCTGTTCAAGCTGGTCGGTGCCCAGCCCCACGGGGAGCAGGTGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCACCGGACGCTGACTGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCcac                                                                       <  2:243969/211‑1 (MQ=255)
                                                                                                                                                             tggcggggtggtgaCGGGAAAGGGGCACGGCGGGGAACGAGGGGCGGGGCGGCTGTCCAAGCGGGTCGGTGCCCGGCCCCACGGGGACCAGTTGCGTGGAACCGGAGGAGGGGCATGGGCGGGGCACCGGACGCTGATTGTCACAGTGTTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTTGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGt                                                                    <  2:254435/215‑1 (MQ=255)
                                                                                                                                                              ttCGCGGTCGATTCGAGAAAGGGGCACGCCGGGGACCAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGGACCAGGTGCTTGGAACCGCAGGAAGGGCACTGGCGCGGCACCGGACGCTGACTGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCg                                                                  <  2:559222/228‑1 (MQ=255)
                                                                                                                                                              aTCGCGGTCGATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGGACCAGGTGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCACCGGACGCTGACTGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTCGTAGCGGTCGCAACGGGGACATTCGGCAATCCAGCCACAGTCg                                                                  >  1:685374/1‑229 (MQ=255)
                                                                                                                                                               tgtcggtggttagagtgaaggggCACGGCGGGGAACACCGGCGGGAGGGTTGGTTCAAGCGGGTCGGTGCCCAGCCCCGCGGGGAACAGGTGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCTCCGGACGCTGAATGTCACAGTGGTGGCAGCGCAGATGGGGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCg                                                                  <  2:344628/208‑1 (MQ=255)
                                                                                                                                                               tgcgggatgATAGGGGAAAGGGGCCCGTCGGGGAACAACGGGGCGAGGGTCGGTTCAAGGTGTTCGGGGCCCAGCCCCCGGGGGTCCAGGTGCATGGAACCGCAGGAAGGGCTCGGGCGCGGCACCGGACGCTGACTGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCg                                                                  <  2:80928/220‑1 (MQ=255)
                                                                                                                                                                 gggttcggtgggagaCAGGGGCCCGCGGGGGGACACGGGGGCGGGCCTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGGACCGGGTGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCACCGGACGCTGACTGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCGTgg                                                               <  2:291864/218‑1 (MQ=255)
                                                                                                                                                                    gTCGATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGGACCAGGTGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCACCGGCCGCTGACTGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCCGCGTGTCGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGcccc                                                                       >  1:243969/1‑216 (MQ=255)
                                                                                                                                                                    gTCGATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGGACCAGGTGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCACCGGACGCTGACTGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGCGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCGGGGcac                                                            >  1:503660/1‑229 (MQ=255)
                                                                                                                                                                                 aTGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGGACCAGGTGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCACCGGACGCTGACTGGCACAGTGGGGGCGGCGCAGATGGTGCGGCGCCTGAGGCAGCGTGTAGTAGGGATCGCAACGAGGGCATTCGGCAATCCAGCCACAGTGGTGGCCcagcagagagggcg                                               >  2:723258/1‑221 (MQ=255)
                                                                                                                                                                                   gtgCACGCCGGGGAACAACGGCGCGAGCGTCTGGTCAAGCTGTTCGGTGCCCAGCCCCACGGGGACCAGGTGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCACCGGACGCTGACTGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTAATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCGTGGCACAGCAGTGCAGGCGCa                                             <  1:450920/227‑1 (MQ=255)
                                                                                                                                                                                        cGCCGGGGAACAGAGGCGTGAGTGTCGGTTAAAGCTGTTCGGTGCCCAGCCCCACGGGGACCAGGTGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCACCGGACGCTGACTGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCGTGGCACAGCAGTGCAGGCGCAAAGc                                         <  1:456108/228‑1 (MQ=255)
                                                                                                                                                                                           cGGGGAACAGCGGGGCGAGCTTCTGTTGAGGTTGTGCGGGGCCCAGCCCCACGGGGACCAGGTGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCACCGGACGCTGACTGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCGTGGCACAGCAGTGCAGGCGCAAAGCCAcg                                     <  1:571936/229‑1 (MQ=255)
                                                                                                                                                                                              ggAACCAGGGCGCGGGCGTCTGTTCAGGCGGTTCGGGGCCCAGCCCCACGGGGACCAGGTGCGTGGAACCGGAGGAAGGGCACTGGCGCGGCACCGGACGCTGACTGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCGTGGCACAGCAGTGCAGGCGCAAAGCCAcggc                                   <  2:262276/228‑1 (MQ=255)
                                                                                                                                                                                                 acgagggcggcgagggttGTTCAAGCTGTTCGGTGCCCAGTCCTAGGGGGCCCAGGGGGGGGGAACCGCAGGAAGTGCAGTGGCGCGGGACCGGACGGTGACTGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCGTGGCACAGCAGTGCAGGCGCAAAGCCACGGCGGtt                               <  1:189502/212‑1 (MQ=255)
                                                                                                                                                                                                 acgacGGCGCGGGGGTGTGTTCAAGCTGGTCGGGGCCCAGCCCCACGGGGCACAGGGGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCACCGGACGCTGACTGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCGTGGCACAGCAGTGCAGGCGCAAAGCCACGGCGGtt                               <  2:23059/226‑1 (MQ=255)
                                                                                                                                                                                                     cGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGGACCAGGTGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCACCGGACGCTGACTGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCGTGGCACAGCAGTGCAGGCGCAAAGCCACGGCGGTTtag                            >  1:115100/1‑228 (MQ=255)
                                                                                                                                                                                                     cGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGGACCAGGTGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCACCGGACGCTGACTGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGGGATCGCAACGTGGGCATTCGGCAATCCAGCCCCAGTCGGTGCACAGCGGTGCAGGGCCAAACCAAAGGCGGGTcaga                           >  2:27380/1‑229 (MQ=255)
                                                                                                                                                                                                     cGGCGAGGGCGTGTGGTCAGGTAGTGCGGGCCCCAGCCCCCGGGGGACCAGGTGCGTGGAACCGCAGTAAGGGCACTGGGGCGGCACCGGACGCTGACTGTCCCAGTGGTGGGAGCGCAGGTGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCGTGGCACAGCAGTGCAGGCGCAAAGCCACGGCGGTTaaga                           <  2:181129/229‑1 (MQ=255)
                                                                                                                                                                                                        cgcgAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGGACCAGGTGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCACCGGACGCTGACTGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACCGCCGTGGCACAGCAGTGCAGGCGCAAAGCAAATGCGGTTtagg                           >  1:539431/1‑225 (MQ=255)
                                                                                                                                                                                                           gAGCGTCTGTTCAAGGTGTTCGGTGCCCAGCCCCAGGGGTACCAGGTGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCACCGGACGCTGACGGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCGTGGCACAGCAGTGCAGGCGCAAAGCCACGGCGGTTAAGAAAgaga                     <  1:507436/229‑1 (MQ=255)
                                                                                                                                                                                                            agtggttGTTCAAAGTGTTGGGGGCCCGGCCCCACGGGGACCAGGTGCGTGGAACCGCAGTAAGGGCACGGGCGCGGCCCCGGACGCTGACTGTCACAGGGGTGGCGGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCGTGGCACAGCAGTGCAGGCGCAAAGCCACGGCGGTTAAGAAAGAGaa                    <  2:504908/223‑1 (MQ=255)
                                                                                                                                                                                                             gCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGGACCAGGTGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCACCGGACGCTGACTGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCACCGTGGGCATTCGGCAATCCAGCCACAGTCGTGGCACAGCAGTGCAGGCGCAAAGCCACGGCGGGTAAGAAAAAGaa                    >  1:183811/1‑228 (MQ=255)
                                                                                                                                                                                                                tgtggtcaatgtGTTCGGTGCCCCGCCCCCCGGGGCCCAGGTGCGTGGAACCGCGGGAAGGGCCCTGGCGCGGCACCGGACGCTGACTGTCCCAGTGGTGGCAGCGCAGATGGTGCGGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCGTGGCACAGCAGTGCAGGCGCAAAGCCACGGCGGTTAAGAAAGAGAATGAc                <  2:82450/218‑1 (MQ=255)
                                                                                                                                                                                                                tCTGTTCAATCTGTTCGGTGCCCAGCCCCACGGGGACCAGGTGCGTGGAACCGAAGGAAGGGCACTGGCGCGGCACCGGACGCTGACTGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCGTGGCACAGCAGTGCAGGCGCAAAGCCACGGCGGTTAAGAAAGAGAATGAc                <  1:547545/229‑1 (MQ=255)
                                                                                                                                                                                                                tCTGTTCAAGCTGGTCGGTGCCCAGCCCCACGGGGACCAGGTGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCACCGGACGCTGACTGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCGTGGCACAGCAGTGCAGGCGCAAAGCCACGGCGGTTAAGAAAGAGAATGAc                <  1:545400/229‑1 (MQ=255)
                                                                                                                                                                                                                   gTTCAAGGTGTTCGGGCCCCAGCCCCCCGGGGACCAGGTGCGTGGAACCGCAGGAAGGGCAATGGCGCGGCACCGGACGGTGACTGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCGTGGCACAGCAGTGCAGGCGCAAAGCCACGGCGGTTAAGAAAGAGAATGACCt              <  1:705434/228‑1 (MQ=255)
                                                                                                                                                                                                                         gcgtgtCGGTGCCCAGCCCCACGGGGACCAGGTGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCACCGGACGCTGACTGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCGTGGCACAGCAGTGCAGGCGCAAAGCCACGGCGGTTAAGAAAGAGAATGACCTGGTTAt        <  1:127534/223‑1 (MQ=255)
                                                                                                                                                                                                                            ggtCGGGGCCCAGCCCCACGGGGACCAGGTGCGTGGAACCGCAGGAAGGGCACTGGGTCGGCACCGGACGCTGACTGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCTGCCTCGGTCGTGGCACAGCAGTGCAGGCGCAAAGCCACGGCGGTTAAGAAAGAGAATGACCTGGTTATCAGc    <  1:113342/227‑1 (MQ=255)
                                                                                                                                                                                                                            gTTCGGTGCCCGGCCCCCGGGGGACCCGTGGCGTGGAACGGCAGGAAGGGCACTGGGGGGGCACCGGACGCTGGATGTCCCAGGGGTGGCAGCGCAGATGGTGCTGCCCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCGTGGCACAGCAGTGCAGGCGCAAAGCCACGGCGGTTAAGAAAGAGAATGACCTGGTTATCAGc    <  2:170159/229‑1 (MQ=255)
                                                                                                                                                                                                                            gTTCGGTGCCCAGCCCCACGGGGACCAGGTGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCACCGGACGCTGACTGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCGTGGCACAGCAGTGCAGGCGCAAAGCCACGGCGGTTAAGAAAGAGAATGACCTGGTTATCAGc    <  1:405624/229‑1 (MQ=255)
                                                                                                                                                                                                                               cGGTGCCCAGGCCCAAGGGGCACAGGGGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCACCGGACGCTGACGGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGACTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCGTGGCACAGCAGTGCAGGCGCAAAGCCACGGCGGTTAAGAAAGAGAATGACCTGGTTATCAGCCt  <  1:731228/228‑1 (MQ=255)
                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                 
AGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATGGTGACCTTTCGCCAGCATTTGTGTACCAATCAAAATCCGCGCGCCGCCGCGATGTACTTCTGCCAGTTGCTGTTCCAGCGCCCCTTTGCGGCTGGTGGTATCGCGGTCGATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGGACCAGGTGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCACCGGACGCTGACTGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCGTGGCACAGCAGTGCAGGCGCAAAGCCACGGCGGTTAAGAAAGAGAATGACCTGGTTATCAGCCT  >  NC_000913/4125115‑4125563

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