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

AGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATGGTGACCTTTCGCCAGCATTTGTGT‑ACCAATCAAAATCCGCGCGCCGCCGCGATGTACTTCTG‑CCAGTTGCTGTTCCAGCGCCCCTTTGCGGCTGGTGGTATCGCGGTCGATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGG‑ACCAGGTGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCACCGGACGCTGACTGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCGTGGCACAGCAGTGCAGGCGCAAAGCCACGGCGGTTAAGAAAGAGAATGACCTGGTTATCAGCCT  >  NC_000913/4125115‑4125563
                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                 
aGAAAACAGCCCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATGGTGACCTTTCGCCAGCATTTGTGT‑TCCAATCAAAATCCGCGCGCCGCCGCGATGTACTTCTG‑CCAGTTGCTGTTCCAGCGCCCCTTTGCGGCTGGTGGTATCGCGGTCGATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGcc                                                                                                                                                                                                                                <  1:595133/228‑1 (MQ=255)
  aaaaaggagccgggtccctgtccgGTAATGCAACCAGGGTCCCGTCGGGAAAATGGTGACCTTTCGCCAGCATTTGGCT‑ACCAATCAAAACCCGCGCGCCGCCGCGAGGTTCTTCTG‑CCAGTTGCTGTTCCAGCGCCCCTTTGCGGCTGGTGGTATCGCGGTCGATACGTGAAATGGGCGCGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCa                                                                                                                                                                                                                              <  2:553012/209‑1 (MQ=255)
          gcgcCGTCCACGTCCAGTAATGCACCCAGCGTCACATCCGGGAAATGGTGCCCTTTTGCCAGCATTTGTGT‑ACCAATCAAAATCCGCGCGCCGCCGCGATGTACTTCTG‑CCAGTTGCTGTTCCAGCGCCCCTTTGCGGCTGGTGGTATCGCGGTCGATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACg                                                                                                                                                                                                                      <  1:603475/228‑1 (MQ=255)
              cgtccacgtccaGTAATGCAACCAGCGTCACATCCGGGAAATGGTGACCTTTCGCCAGCATTTGTGT‑ACCAATAAAAATCCGCGAGCGGCCGCGATGTACTTCTG‑CCAGTTGCTGTTCCAGCGCCCCTTTGCGGCTGGTGGAATGCCGGTCGTTACCGGAAATGGGCACTCCCGGGCAAAACGGCCGGGGCGCTCGGTCAAGCTGGTGGGCGACCACCCCCACCGGG‑Ac                                                                                                                                                                                                                >  2:230901/1‑229 (MQ=255)
                tccacgtccaGTAATGCAACCAGCGTCACATCCGGGAAATGGTGACCTTTCGCCAGCATTTGTGT‑ACCAATCAAAATCCGCGCGCCGCCGCGATGTACTTCTG‑CCAGTTGCTGTTCCAGCGCCCCTTTGCGGCTGGTGGGATCGCGGTCGATACGAGAAATGGGCACGCCGGGGAACCACGCCGCGGGGGTCTGTTCAAACGGTTCGGGGCCCAGCCCCCCggggcgac                                                                                                                                                                                                               >  2:549629/1‑225 (MQ=255)
                 ccccgcccgttatggacaCCAGCGCCCCATCCGGGAAATGGGGACCTTTCGCCAGCATTTGTGT‑ACCAATCAAAATCCGCGCGCCGCGGGGATGTACTTCTG‑CCAGTTGCTGTTCCAGCGCCCCTTTGCGGCTGGTGGTATCGCGGTCGATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGG‑ACCAg                                                                                                                                                                                                             <  2:696355/212‑1 (MQ=255)
                 cccacgcccgTAATTCAACCACGGCCACAGCCGGGAAAGGGTGCCCTTGTGCCAGCATTTGTGT‑ACCAAGCAAAATCCGCGCGCCGCCGCGATGTACTTCTG‑CCAGTTGCTGTTCCAGCGCCCCTTTGCGGCTGGTGGTATCGCGGTCGATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGG‑ACCAg                                                                                                                                                                                                             <  2:188604/220‑1 (MQ=255)
                              tGCAACCAGCGTCACATCCGGGAAATGGTGACCTTTCGCCAGCATTTGTGT‑ACCAATCAAAATCCGCGCGCCGCCGCGATGTACTTCTG‑CCAGTTGCTGTTCCAGCGCCCCTTTGCGGCTGGTGGTATCGCGGTCGATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGGGCGTCTGTTCAAGCTGTTCGGGGCCCGGCCCCCCGGGG‑GCACGGTGGGGGGAAccc                                                                                                                                                                                                >  1:314478/1‑228 (MQ=255)
                               gCAACCAGCGTCACATCCGGGAAATGGTGACCTTTCGCCAGCATTTGTGT‑ACCAATCAAAATCCGCGCGCCGCCGCGATGTACTTCTG‑CCAGTTGCTGTTCCAGCGCCCCTTTGCGGCTGGTGGTATCGCGGTCGATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCCCGGGG‑ACCAGGTGCGTGGAACCg                                                                                                                                                                                                >  1:117250/1‑228 (MQ=255)
                                cAACCAGCGTCACATCCGGGAAATGGTGACCTTTCGCCAGCATTTGTGT‑ACCAATCAAAATCCGCGCGCCGCCGCGATGTACTTCTG‑CCAGTTGCTGTTCCAGCGCCCCTTTGCGGCTGGTGGTATCGCGGTCGATACGAGAAATGGGCACGCCGGGGGAGAACGGCGCGAGCGTCCGTTCAAGCTGGTCGGTGGCCAACCCCCCGGGG‑GCCAGGGGCGGGGCAccacc                                                                                                                                                                                              >  1:320056/1‑226 (MQ=255)
                                   ccAGCGTCACATCCGGGAAATGGTGACCTTTCGCCAGCATTTGTGT‑ACCAATCAAAATCCGCGCGCCGCCGCGATGTACTTCTG‑CCAGTTGCTGTTCCAGCGCCCCTTTGCGGCTGGTGGTATCGCGGTCGATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGGCCAGCCCCCCGGGG‑ACAAGGGACGTGGAACCGCAggg                                                                                                                                                                                           >  1:753200/1‑228 (MQ=255)
                                     aGCGTCACATCCGGGAAATGGTGACCTTTCGCCAGCATTTGTGT‑ACCAATCAAAATCCGCGCGCCGCCGCGATGTACTTCTG‑CCAGTTGCTGTTCCAGCGCCCCTTTGCGGCTGGTGGAATCGCGGCCGATACGCGAAATGGGCACGCCGGGGAAAAACGGGGCGCGCGTCTTGTTAAGCGTTACGGGGCCCACCCCCACGGGGCACCCGG‑GCGGGGGACCGCacagaa                                                                                                                                                                                         >  2:448306/1‑224 (MQ=255)
                                      ccgtccaaTCCGGGAAAGTGTGCCCTTTCTCCAGCATTTGTGG‑CCCAACCAAACCCCGCCGGCGGGCGGGATGTACTGCTG‑CCAGTTGCTGTGCCAGCGCCCCTTTGCGGCGGGGGGTATCGCGGTCGATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGG‑ACCAGGTGCGTGGAACCGCAGGa                                                                                                                                                                                           <  2:320056/219‑1 (MQ=255)
                                       cGTCACATCCGGGAAATGGTTACCTTTCGCCAGCATTTGTGT‑ACCAATCAAAATCCGCGCGCCGCCGCGATGTACTTCTG‑CCAGTTGCTGTTCCAGCGCCCCTTTGCGGCTGGTGGTATCGCGGTCGATACGAGAAATGGGCACGCCGGGGAACAACGGCCGGAGGGCCGGTTCCAGCTGTTCGGTGCCCACCCCCCCGCGG‑CCCGGGGGCGGGGAAcacccgggaggg                                                                                                                                                                                       >  2:405624/1‑218 (MQ=255)
                                       cGTCACATCCGGGAAATGGTGACCTTTCGCCAGCATTTGTGT‑ACCAATCAAAATCCGCGCGCCGCCGCGATGTACTTCTG‑CCAGTTGCTGTTCCAGCGCCCCTTTGCGGCTGGTGGTATCGCGGTCGATACGAGAAATGGGAACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGGGCCCAGCCCCACGGGG‑ACCAGGGGCGGGGAACCGCAGGAAagg                                                                                                                                                                                       >  1:291864/1‑229 (MQ=255)
                                         tCACATCCGGGAAATGGTGACCTTTCGCCAGCATTTGTTT‑ACCAATCAAAATCCGCGCGCCGCCGCGATGTACTTCTG‑CCAGTTGCTGTTCCAGCGCCCCTTTGCGGCTGGTGGTATCGCGGGCGATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGGGCGGCTGTTCAAACTGTTCGGTGCCCAGCCCCCCGGGG‑GCCAggggggggggacacgggaaggggc                                                                                                                                                                                      >  2:545400/1‑206 (MQ=255)
                                         tCACATCCGGGAAATGGTGACCTTTCGCCAGCATTTGTGT‑ACCAATCAAAATCCGCGCGCCGCCGCGATGTACTTCTG‑CCAGTTGCTGTTCCAGCGCCCCTTTGCGGCTGGTGGTATCGCGGTCGATACGAGAAATGGGCCCGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGGGCCCAGCCCCCCGGGG‑ACCAGGGGCGCGGAACCGCAGGAAGGGc                                                                                                                                                                                      >  1:181454/1‑228 (MQ=255)
                                                     aaTGGTGGCCTTTCCCCAGCATTTGGGT‑ACCAATCAAAAGCCGCGCGCCGCCGCGATGGACTTCTG‑CCAGTTGCTGTTCCAGCGCCCCTTTGCGGCTGGTGGTATCGCGGTCGATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGG‑ACCAGGTGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCAc                                                                                                                                                                         <  1:549629/229‑1 (MQ=255)
                                                         tgcgctcttCGCCAGCATTT‑TGTGACAAACCGAAACCCCCGTCCCCCCGGGATGTA‑GTCTGCCCAGTGGTTGTACCAGGGCCCCTGTGCGGCTGGGGGTATCGCGGTCGATGCGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGGTCAAGCTGTTCGGTGACCAGCCCCACGGGG‑ACCAGGTGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCACCTGa                                                                                                                                                                     <  2:456046/222‑1 (MQ=255)
                                                         gTGACCTTTCGCCAGCATTTGTGT‑ACCAATCAAAATCCGCGCGCCGCCGCGATGTACTTCTG‑CCAGTTGCTGTTCCAGCGCCCCTTTGCGGCTGGTGGTATCGCGGTCTATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGG‑ACCCGGTGCGTGGAACCGCAGGAAGGGCACTGACGCGCGACCgg                                                                                                                                                                      >  1:665345/1‑228 (MQ=255)
                                                            aCCTTTCGCCAGCATTTGTGT‑ACCAATCAAAATCCGCGCGCCGCCGCGATGTACTTCTG‑CCAGTTGCTGTTCCAGCGCCCCTTTGCGGCTGGTGGTATCGCGGTCGATACGGCAAATGGGCCCGCCGGGGCACAACGGCCCGAGCGCCTGTTCAAGCTCTTCGGGGCCCACCCCCCCGGGG‑GCCAGGTGCGCTGCACCGCCGGAAGGGCCCT‑GCGCGGGAcccggaccc                                                                                                                                                                 >  2:113782/1‑222 (MQ=255)
                                                                   ggcatctttttTGG‑ACCAACCAAAAGCCGCGCGCCGCCGCGATGTTCTTCTG‑CCAGTTGCTGTTCCCGGGCCCCTTTGCGGCTGGGGGTATCGCGGTCGATACGAGAAAGGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGG‑ACCAGGTGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCACCGGACGCTGACTGt                                                                                                                                                           <  2:659650/222‑1 (MQ=255)
                                                                                   ccAATCAAAATCCGCGCGCCGCCGCGATGTACTTCTG‑CCAGTTGCTGTTCCAGCGCCCCTTTGCGGCTGGTGGTATCGCGGTCGATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGcc                                                                                                                                                                                                                                >  2:158943/1‑146 (MQ=255)
                                                                                   ccAATCAAAATCCGCGCGCCGCCGCGATGTACTTCTG‑CCAGTTGCTGTTCCAGCGCCCCTTTGCGGCTGGTGGTATCGCGGTCGATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGcc                                                                                                                                                                                                                                <  1:158943/146‑1 (MQ=255)
                                                                                    cAATCAAAATCCGCGCGCCGCCGCGATGTACTTCTG‑CCAGTTGCTGTTCCAGCGCCCCTTTGCGGCTGGTGGTATCGCGGTCGATACGAGAAATGGGCACGCCGGGGACCAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGG‑ACCAGGCGCGTGGAACCGCAGGAAGGGCACTGGCCCGGCACCGGCCCCTGCCGGGCCCAGTGGGGGcgcgc                                                                                                                                           >  1:692692/1‑225 (MQ=255)
                                                                                      aTCAAAATCCGCGCGCCGCCGCGATGTACTTCTG‑CCAGTTGCTGTTCCAGCGCCCCTTTGCGGCTGGTGGTATCGCGGTCGATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGGCCCACCGGG‑GCCCGGGGCGTGGAACCGCAGGCAGGGCAACGGCGCGGCCCCGGCCGCCGAcgtccaccgtcgggggagcgga                                                                                                                                         >  2:705918/1‑208 (MQ=255)
                                                                                      aTCAAAATCCGCGCGCCGCCGCGATGTACTTCTG‑CCAGTTGCTGTTCCAGCGCCCCTTTGCGGCTGGTGGTATCGCGGTCGATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACg                                                                                                                                                                                                                      >  1:261685/1‑153 (MQ=255)
                                                                                      aTCAAAATCCGCGCGCCGCCGCGATGTACTTCTG‑CCAGTTGCTGTTCCAGCGCCCCTTTGCGGCTGGTGGTATCGCGGTCGATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGG‑GCCAGGTGCGTGGAACCGCAGGAAGGGCACCGGCCCGGCACCGGCCGCTGCCCGGCACAGGGGGGGGAccgcc                                                                                                                                         >  2:731227/1‑228 (MQ=255)
                                                                                      aTCAAAATCCCCGCGCCGGCGGGAGGGACTTCGT‑GCAGTTGCTGTTCCCGCGCTCCTTTGCGGCTGGTGGTATCGCGGGCGATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGG‑ACCAGGTGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCACCGGACGCTGACTGTCACAGTGGTGGCAGCGCa                                                                                                                                         <  2:564170/229‑1 (MQ=255)
                                                                                         aaaaTCCGCGCGCCGCCGCGATGTACTTCTG‑CCAGTTGCTGTTCCAGCGCCCCTTTGCGGCTGGTGGTATCGCGGTCGATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGG‑ACCAGGTGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCACCGGACGCTGACTGTCACAGTGGTGGCAGCGCAGa                                                                                                                                       >  1:23059/1‑228 (MQ=255)
                                                                                                 cgcgcCGCCGCGATGTACTTCTG‑CCAGTTGCTGTTCCAGCGCCCCTTTGCGGCTGGTGGTATCGCGGTCGATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGG‑CCCAGGTGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCACCGGACGCCGCCTGTCACAGTGGTGGCAGCGCAGCTGGggctgc                                                                                                                              >  1:522759/1‑229 (MQ=255)
                                                                                                          gggagttACTTCTG‑CCCGTTGCTGTTCCAGGGGCCCTTGGCGGCTGGTTGTATCGCGGTCGATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGG‑ACCAGGTGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCACCGGACGCTGACTGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGc                                                                                                                     <  1:414161/223‑1 (MQ=255)
                                                                                                           ggtgtgtcTTGTG‑CCGTTTGCGGTGCCAGCGCCCCTTGGGGGGTGGTGGTATCGGGGTCGATACGAGAAATGGGCCCGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGG‑ACCAGGTGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCACCGGACGCTGACTGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCa                                                                                                                    <  2:117250/222‑1 (MQ=255)
                                                                                                           gggggtctgtggg‑CCGGTTGGTGTTCCGTCGCCCCTTGGCGGCTGGGGGTGTCGCGGTGGTTCCGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGG‑ACCAGGTGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCACCGGACGCTGACTGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCa                                                                                                                    <  2:402258/217‑1 (MQ=255)
                                                                                                                                        gcgcCCCTTTGCGGCTGGTGGTATCGCGGTCGATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGG‑ACCAGGTGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCACCGGACGCTGACTGTCACAGGGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGGGTAGTAGGGATGGCAACGTGGGcg                                                                                        >  1:259742/1‑228 (MQ=255)
                                                                                                                                           ccccTTTGCGGCTGGTGGTATCGCGGTCGATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGG‑CCCAGGTGCGTGGAACCGCAGGAAGGGCACTCGCCCGGCCCCGGACGCCGACTCTCACCGGGGGGGCCGCCCGGGTGCGGCTGCGGCGTATGCCGCGTTGTGGAGTGACCTCCACGTGGGCAtt                                                                                      >  2:22113/1‑228 (MQ=255)
                                                                                                                                                 gtgggCTGGGGGTTCCGCGGGCGAGGAGTGAAATGGGCCCGCCGGGGAACAACGGCGGGGGCGTCTGGTCAAGCTGTGCGGGGCCCAGCCCCACGGGG‑ACCAGGTGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCACCGGACGCTGACTGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAAt                                                                               <  2:663032/226‑1 (MQ=255)
                                                                                                                                                  gCGGGTGGTGGGTTCGGGGTCGATACGAGAAAGGGCCACGCCGGGGAACATCGGCGCGGGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCGCGGGG‑ACCAGGTGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCCCCGGACGCTGACTGTCACAGTGGTGGCCGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATc                                                                              <  2:314478/229‑1 (MQ=255)
                                                                                                                                                       tggtggTATCGCGGTCGATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGGGCCCAGCCCCACGGGG‑ACCAGGGGCGTGGAACCGCAGGAAGGGCACCGGCCCGGGCCCGGACGCTGACTGTCCCAGGGGGGGCAGCCCAGATGGGGCGGCCCGTCATGCCGCGTGGAGTAGGGATTGCAACGAGGGCATTTgggcacacagga                                                                         >  2:14409/1‑219 (MQ=255)
                                                                                                                                                        ggtggtATCGGGGGCGATACGAGAAATGGGCCCGCCGGGGAACAGCGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGG‑ACCAGGTGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCACCGGACGCTGACTGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCca                                                                        <  1:104844/229‑1 (MQ=255)
                                                                                                                                                        ggtggtATCGCGTGCGAATCGGAAAACGGTCACGCCGGGGAAAAACGGGGCGTGGGTCTGTTCAAGGTGGTCGGTGCCCAGCCCCCGCGGG‑ACCAGGTGCGTGGAACCGCAGGAAGGGCCCTGGCGCGGCCCCGGACGCTGACTGTCCCAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCca                                                                        <  2:205601/229‑1 (MQ=255)
                                                                                                                                                         gtggtATCGCGGTCGATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGG‑ACCAGGTGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCACCGGACGCTGACTGTCACAGTGGTGGCAGCGCAGAGGGTGCTGGGCCTGATGCAGCGGGTGGTGGTGATCGCAACGTGGGCATTCGGCAATCCAGCcac                                                                       >  1:262276/1‑229 (MQ=255)
                                                                                                                                                         gtggtAGCGGGGTCGATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGG‑ACCAGGTGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCACCGGACGCTGACTGTCACAGTGGGGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCcac                                                                       <  1:607104/229‑1 (MQ=255)
                                                                                                                                                         ggggtggcggggtggaggcGAGAAATGGGCACGCCGGGGAACAACGGCGCGGGCGTCTGTTCAAGCTGGTCGGTGCCCAGCCCCACGGGG‑AGCAGGTGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCACCGGACGCTGACTGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCcac                                                                       <  2:243969/211‑1 (MQ=255)
                                                                                                                                                             tggcggggtggtgaCGGGAAAGGGGCACGGCGGGGAACGAGGGGCGGGGCGGCTGTCCAAGCGGGTCGGTGCCCGGCCCCACGGGG‑ACCAGTTGCGTGGAACCGGAGGAGGGGCATGGGCGGGGCACCGGACGCTGATTGTCACAGTGTTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTTGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGt                                                                    <  2:254435/215‑1 (MQ=255)
                                                                                                                                                              ttCGCGGTCGATTCGAGAAAGGGGCACGCCGGGGACCAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGG‑ACCAGGTGCTTGGAACCGCAGGAAGGGCACTGGCGCGGCACCGGACGCTGACTGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCg                                                                  <  2:559222/228‑1 (MQ=255)
                                                                                                                                                              aTCGCGGTCGATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGG‑ACCAGGTGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCACCGGACGCTGACTGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTCGTAGCGGTCGCAACGGGGACATTCGGCAATCCAGCCACAGTCg                                                                  >  1:685374/1‑229 (MQ=255)
                                                                                                                                                               tgtcggtggttagagtgaaggggCACGGCGGGGAACACCGGCGGGAGGGTTGGTTCAAGCGGGTCGGTGCCCAGCCCCGCGGGG‑AACAGGTGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCTCCGGACGCTGAATGTCACAGTGGTGGCAGCGCAGATGGGGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCg                                                                  <  2:344628/208‑1 (MQ=255)
                                                                                                                                                               tgcgggatgATAGGGGAAAGGGGCCCGTCGGGGAACAACGGGGCGAGGGTCGGTTCAAGGTGTTCGGGGCCCAGCCCCCGGGGG‑TCCAGGTGCATGGAACCGCAGGAAGGGCTCGGGCGCGGCACCGGACGCTGACTGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCg                                                                  <  2:80928/220‑1 (MQ=255)
                                                                                                                                                                 gggttcggtgggagaCAGGGGCCCGCGGGGGGACACGGGGGCGGGCCTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGG‑ACCGGGTGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCACCGGACGCTGACTGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCGTgg                                                               <  2:291864/218‑1 (MQ=255)
                                                                                                                                                                    gTCGATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGG‑ACCAGGTGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCACCGGCCGCTGACTGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCCGCGTGTCGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGcccc                                                                       >  1:243969/1‑216 (MQ=255)
                                                                                                                                                                    gTCGATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGG‑ACCAGGTGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCACCGGACGCTGACTGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGCGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCGGGGcac                                                            >  1:503660/1‑229 (MQ=255)
                                                                                                                                                                                 aTGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGG‑ACCAGGTGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCACCGGACGCTGACTGGCACAGTGGGGGCGGCGCAGATGGTGCGGCGCCTGAGGCAGCGTGTAGTAGGGATCGCAACGAGGGCATTCGGCAATCCAGCCACAGTGGTGGCCcagcagagagggcg                                               >  2:723258/1‑221 (MQ=255)
                                                                                                                                                                                   gtgCACGCCGGGGAACAACGGCGCGAGCGTCTGGTCAAGCTGTTCGGTGCCCAGCCCCACGGGG‑ACCAGGTGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCACCGGACGCTGACTGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTAATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCGTGGCACAGCAGTGCAGGCGCa                                             <  1:450920/227‑1 (MQ=255)
                                                                                                                                                                                        cGCCGGGGAACAGAGGCGTGAGTGTCGGTTAAAGCTGTTCGGTGCCCAGCCCCACGGGG‑ACCAGGTGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCACCGGACGCTGACTGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCGTGGCACAGCAGTGCAGGCGCAAAGc                                         <  1:456108/228‑1 (MQ=255)
                                                                                                                                                                                           cGGGGAACAGCGGGGCGAGCTTCTGTTGAGGTTGTGCGGGGCCCAGCCCCACGGGG‑ACCAGGTGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCACCGGACGCTGACTGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCGTGGCACAGCAGTGCAGGCGCAAAGCCAcg                                     <  1:571936/229‑1 (MQ=255)
                                                                                                                                                                                              ggAACCAGGGCGCGGGCGTCTGTTCAGGCGGTTCGGGGCCCAGCCCCACGGGG‑ACCAGGTGCGTGGAACCGGAGGAAGGGCACTGGCGCGGCACCGGACGCTGACTGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCGTGGCACAGCAGTGCAGGCGCAAAGCCAcggc                                   <  2:262276/228‑1 (MQ=255)
                                                                                                                                                                                                 acgagggcggcgagggttGTTCAAGCTGTTCGGTGCCCAGTCCTAGGGGG‑CCCAGGGGGGGGGAACCGCAGGAAGTGCAGTGGCGCGGGACCGGACGGTGACTGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCGTGGCACAGCAGTGCAGGCGCAAAGCCACGGCGGtt                               <  1:189502/212‑1 (MQ=255)
                                                                                                                                                                                                 acgacGGCGCGGGGGTGTGTTCAAGCTGGTCGGGGCCCAGCCCCACGGGG‑CACAGGGGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCACCGGACGCTGACTGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCGTGGCACAGCAGTGCAGGCGCAAAGCCACGGCGGtt                               <  2:23059/226‑1 (MQ=255)
                                                                                                                                                                                                     cGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGG‑ACCAGGTGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCACCGGACGCTGACTGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCGTGGCACAGCAGTGCAGGCGCAAAGCCACGGCGGTTtag                            >  1:115100/1‑228 (MQ=255)
                                                                                                                                                                                                     cGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGG‑ACCAGGTGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCACCGGACGCTGACTGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGGGATCGCAACGTGGGCATTCGGCAATCCAGCCCCAGTCGGTGCACAGCGGTGCAGGGCCAAACCAAAGGCGGGTcaga                           >  2:27380/1‑229 (MQ=255)
                                                                                                                                                                                                     cGGCGAGGGCGTGTGGTCAGGTAGTGCGGGCCCCAGCCCCCGGGGG‑ACCAGGTGCGTGGAACCGCAGTAAGGGCACTGGGGCGGCACCGGACGCTGACTGTCCCAGTGGTGGGAGCGCAGGTGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCGTGGCACAGCAGTGCAGGCGCAAAGCCACGGCGGTTaaga                           <  2:181129/229‑1 (MQ=255)
                                                                                                                                                                                                        cgcgAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGG‑ACCAGGTGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCACCGGACGCTGACTGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACCGCCGTGGCACAGCAGTGCAGGCGCAAAGCAAATGCGGTTtagg                           >  1:539431/1‑225 (MQ=255)
                                                                                                                                                                                                           gAGCGTCTGTTCAAGGTGTTCGGTGCCCAGCCCCAGGGGT‑ACCAGGTGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCACCGGACGCTGACGGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCGTGGCACAGCAGTGCAGGCGCAAAGCCACGGCGGTTAAGAAAgaga                     <  1:507436/229‑1 (MQ=255)
                                                                                                                                                                                                            agtggttGTTCAAAGTGTTGGGGGCCCGGCCCCACGGGG‑ACCAGGTGCGTGGAACCGCAGTAAGGGCACGGGCGCGGCCCCGGACGCTGACTGTCACAGGGGTGGCGGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCGTGGCACAGCAGTGCAGGCGCAAAGCCACGGCGGTTAAGAAAGAGaa                    <  2:504908/223‑1 (MQ=255)
                                                                                                                                                                                                             gCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGG‑ACCAGGTGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCACCGGACGCTGACTGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCACCGTGGGCATTCGGCAATCCAGCCACAGTCGTGGCACAGCAGTGCAGGCGCAAAGCCACGGCGGGTAAGAAAAAGaa                    >  1:183811/1‑228 (MQ=255)
                                                                                                                                                                                                                tgtggtcaatgtGTTCGGTGCCCCGCCCCCCGGGG‑CCCAGGTGCGTGGAACCGCGGGAAGGGCCCTGGCGCGGCACCGGACGCTGACTGTCCCAGTGGTGGCAGCGCAGATGGTGCGGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCGTGGCACAGCAGTGCAGGCGCAAAGCCACGGCGGTTAAGAAAGAGAATGAc                <  2:82450/218‑1 (MQ=255)
                                                                                                                                                                                                                tCTGTTCAATCTGTTCGGTGCCCAGCCCCACGGGG‑ACCAGGTGCGTGGAACCGAAGGAAGGGCACTGGCGCGGCACCGGACGCTGACTGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCGTGGCACAGCAGTGCAGGCGCAAAGCCACGGCGGTTAAGAAAGAGAATGAc                <  1:547545/229‑1 (MQ=255)
                                                                                                                                                                                                                tCTGTTCAAGCTGGTCGGTGCCCAGCCCCACGGGG‑ACCAGGTGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCACCGGACGCTGACTGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCGTGGCACAGCAGTGCAGGCGCAAAGCCACGGCGGTTAAGAAAGAGAATGAc                <  1:545400/229‑1 (MQ=255)
                                                                                                                                                                                                                   gTTCAAGGTGTTCGGGCCCCAGCCCCCCGGGG‑ACCAGGTGCGTGGAACCGCAGGAAGGGCAATGGCGCGGCACCGGACGGTGACTGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCGTGGCACAGCAGTGCAGGCGCAAAGCCACGGCGGTTAAGAAAGAGAATGACCt              <  1:705434/228‑1 (MQ=255)
                                                                                                                                                                                                                         gcgtgtCGGTGCCCAGCCCCACGGGG‑ACCAGGTGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCACCGGACGCTGACTGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCGTGGCACAGCAGTGCAGGCGCAAAGCCACGGCGGTTAAGAAAGAGAATGACCTGGTTAt        <  1:127534/223‑1 (MQ=255)
                                                                                                                                                                                                                            ggtCGGGGCCCAGCCCCACGGGG‑ACCAGGTGCGTGGAACCGCAGGAAGGGCACTGGGTCGGCACCGGACGCTGACTGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCTGCCTCGGTCGTGGCACAGCAGTGCAGGCGCAAAGCCACGGCGGTTAAGAAAGAGAATGACCTGGTTATCAGc    <  1:113342/227‑1 (MQ=255)
                                                                                                                                                                                                                            gTTCGGTGCCCGGCCCCCGGGGG‑ACCCGTGGCGTGGAACGGCAGGAAGGGCACTGGGGGGGCACCGGACGCTGGATGTCCCAGGGGTGGCAGCGCAGATGGTGCTGCCCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCGTGGCACAGCAGTGCAGGCGCAAAGCCACGGCGGTTAAGAAAGAGAATGACCTGGTTATCAGc    <  2:170159/229‑1 (MQ=255)
                                                                                                                                                                                                                            gTTCGGTGCCCAGCCCCACGGGG‑ACCAGGTGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCACCGGACGCTGACTGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCGTGGCACAGCAGTGCAGGCGCAAAGCCACGGCGGTTAAGAAAGAGAATGACCTGGTTATCAGc    <  1:405624/229‑1 (MQ=255)
                                                                                                                                                                                                                               cGGTGCCCAGGCCCAAGGGG‑CACAGGGGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCACCGGACGCTGACGGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGACTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCGTGGCACAGCAGTGCAGGCGCAAAGCCACGGCGGTTAAGAAAGAGAATGACCTGGTTATCAGCCt  <  1:731228/228‑1 (MQ=255)
                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                 
AGAAAACAGCGCGCCGTCCACGTCCAGTAATGCAACCAGCGTCACATCCGGGAAATGGTGACCTTTCGCCAGCATTTGTGT‑ACCAATCAAAATCCGCGCGCCGCCGCGATGTACTTCTG‑CCAGTTGCTGTTCCAGCGCCCCTTTGCGGCTGGTGGTATCGCGGTCGATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGG‑ACCAGGTGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCACCGGACGCTGACTGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCGTGGCACAGCAGTGCAGGCGCAAAGCCACGGCGGTTAAGAAAGAGAATGACCTGGTTATCAGCCT  >  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: ‑