Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,125,363 T→G 26.1% H483P (CAC→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,3630TG26.1% 80.7 / 6.8 52H483P (CAC→CCC) priAPrimosome factor n' (replication factor Y)
Reads supporting (aligned to +/- strand):  ref base T (12/24);  new base G (9/7);  total (21/31)
Fisher's exact test for biased strand distribution p-value = 1.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTAATGCAACCAGCGTCACATCCGGGAAATGGTGACCTTTCGCCAGCATTTGTGTACCAATCAAAATCCGCGCGCCG‑CCGCGATGTACTTCTGCCA‑GTTGCTGTTCCAGCGCCCCTTTGCGGCTGGTGGTATCGCGGTCGATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGGACCAG‑GTGCGTGGAACCGCAGGAAGGGCACTG‑GCGCGGCACCGGACGCTGACTGTCAC‑AGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCGTGGCACAGCAGTGCAGGCGCAAAGCCACGGCGGTTAAGAAAGAGAATGACCTGGTTATCAGCCTGTAAATGCTGGCGCATACGAGTGAT  >  NC_000913/4125142‑4125589
                                                                                                                                                                                                                                |                                                                                                                                                                                                                                    
gttaggCACCCAGCGTCACCTCCGGGAAAGGGTGCCCTTCCGGCACCATTTGGTTACCAATCAAAATCCGCGCGCCG‑CCGCGATGTTCTTCTGCCA‑TTTGCGGTCCCAGCGCCCCTTTGCGGCTGGTGGTATCGCGGTCGATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGGACCAG‑GTGCGTgg                                                                                                                                                                                                                                <  2:360359/223‑1 (MQ=255)
      cacccAGCGCCACATCCGGGAAATGGTGACCTTTCGCCAGCATTTGGGTACCAATCAAAATCCGCGCGCCG‑CCGGGATGTACTTCTGCCA‑GTTGCTGTTCCAGCGCCCCTTTGCGGCTGGTGGTATCGCGGTCGATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGGACCAG‑GTGCGTGGAACCGc                                                                                                                                                                                                                          <  1:564440/225‑1 (MQ=255)
       aaCCAGCGTCACATCCGGGAAATGGTGACCTTTCGCCAGCATTTGTGTACCAATCAAAATCCGCGCGCCG‑CAGCGATGTACTTCTGCCA‑GTTGCTGTTCCAGCGCCCCTTTGCGGCTGGCGGTATCGCGGTCGATACGAGAAATGGGCACGCCGGGGACCAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCAAGCCCCACGGGGACCAG‑GGGCGTGGAccc                                                                                                                                                                                                                            >  1:459307/1‑225 (MQ=255)
         ccAGCGTCACATCCGGGAAATGGTGACCTTTCGCCAGCATTTGTGTACCAATCAAAATCCGCGCGCCG‑CCGCGATGTACTTCTGCCA‑GTTGCTGTTCCAGCGCCCCTTTGCGGCTGGTGGTATCGCGGTCGATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGGGTCTGTTCAAGGTGTTCGGTGCCCAGGCCCCCGGGGGACAG‑GGGCGGGGACACGCAgg                                                                                                                                                                                                                       >  2:398903/1‑228 (MQ=255)
         ccAGCGTCACATCCGGGAAATGGTGACCTTTCGCCAGCATTTGTGTACCAATCAAAATCCGCGCGCCG‑CCGCGATGTACTTCTGCCA‑GTTGCTGTTCCAGCGCCCCTTTGCGGCTGGTGGTATCGCGGTCGATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGGACCAG‑GTGCGTGGAACCGCAg                                                                                                                                                                                                                        >  1:9991/1‑227 (MQ=255)
         ccAGCGTCACATCCGGGAAATGGTGACCTTTCGCCAGCATTTGTGTACCAATCAAAATCCGCGCGCCG‑CCGCGATGTACTTCTGCCA‑GTTGCTGTTCCAGCGCCCCTTTGCGGCTGGTGGTATCGCGGTCGATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGGCTGTTCAAGCTGTTCGGGGGCCCGCCCCCAGGGGACCAG‑GGGCGCGGAAACCCAgg                                                                                                                                                                                                                       >  2:563286/1‑228 (MQ=255)
             cGTCACATCCGGGAAATGGTGACCTTTCGCCAGCATTTGTGTACCAATCAAAATCCGCGCGCCG‑CCGCGATGTACTTCTGCCA‑GTTGCTGTTCCAGCGCCCCTTTGCGGCTGGTGGTATCGCGGTCGATACGAGAAATGGGCACGCCGGGGAACAACGGCGGGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGGACCAG‑GTGCGTGGAACCGCAGGAAgg                                                                                                                                                                                                                   >  1:180653/1‑228 (MQ=255)
                                     tttCGCCCGCTTTTGTGTACCAATCAAAATCCGCGCGCCG‑CCGCGATGTACTTCTGCCA‑GTTGCTGTTCCAGCGCCCCTTTGCGGCTGGTGGTATCGCGGTCGATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGGACCAG‑GTGCGTGGAACCGCAGGAAGGGCACTG‑GCGCGGCACCGGACGCt                                                                                                                                                                                           <  1:570159/227‑1 (MQ=255)
                                     atccGCCAGCATTTGTGCACCCATCAAAATCCGCGCGCCG‑CCGCGATGTGCTTCTGCCA‑GTTGCTGTTCCAGCGCCCCTTTGCGGCTGGTGGTATCGGGGTCGATTCGAGAAAGGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGGACCAG‑GTGCGTGGAACCGCAGGAAGGGCACTG‑GCGCGGCACCGGACGCt                                                                                                                                                                                           <  1:571037/224‑1 (MQ=255)
                                      ttCGCCAGCATTGGTGTCCCAATCAAAATCCGCGCGCCG‑GCGCGATGTACTTCGGCCA‑GTTGCTGTTCCAGCGCCCCTTTGCGGCTGGTGGTATCGCGGTCGATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGGACCAG‑GTGCGTGGAACCGCAGGAAGGGCACTG‑GCGCGGCACCGGACGCTGa                                                                                                                                                                                         <  1:269163/228‑1 (MQ=255)
                                            aGCATTTGGGTACCAATCAAAATCCGCGCGCCG‑GCGCGATGTACTTCGGCCA‑GTTGCTGTTCCAGCGCCCCTTTGCGGCTGGTGGTATCGCGGTCGATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGGTCGGTGCCCAGCCCCACGGGGACCAG‑GTGCGTGGAACCGCAGGAAGGGCACTG‑GCGCGGCACCGGACGCTGACTGTc                                                                                                                                                                                    <  1:398903/227‑1 (MQ=255)
                                              cATTTGTGTACCAATCAAAATCCGCGCGCCG‑CCGCGATGTACTTCTGCCA‑GTTGCTGTTCCAGCGCCCCTTTGCGGCTGGTGGTATCGCGGTCGATACGAGAAATGGGCACGCCGGGGCACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGGGCCCAGCCCCACGGGGACCAG‑GGGCGTGGAAACGCAGGAAGGGCCCGG‑GCGGGGCACCGGACGCGGAGTGTcaaa                                                                                                                                                                                 >  2:269163/1‑226 (MQ=255)
                                                      tACCAATCAAAACCCCGGCGCCG‑GCGGGGGGTGCTTTCGCAA‑GTTGCTGTTCCAGCGCCCTTTGGCGGCGGGGGGTATCGCGGTCGATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGGACCAG‑GTGCGTGGAACCGCAGGAAGGGCACTG‑GCGCGGCACCGGACGCTGACTGTCAC‑AGTGGTGGc                                                                                                                                                                        <  2:536837/228‑1 (MQ=255)
                                                        ccAATCAAAATCCGCGCGCCG‑CCGCGATGTACTTCTGCCA‑GTTGCTGTTCCAGCGCCCCTTTGCGGCTGGTGGTATCGCGGTCGATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCCGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGGACCAG‑GGGCGAGGAACCGCCGGAAGGGCACTG‑GGGCGGGCCCGGAAGCGGCCTGGCAC‑AGGGGGGGGAg                                                                                                                                                                      >  2:575139/1‑228 (MQ=255)
                                                        ccAATCAAAATCCGCGCGCCG‑CCGCGATGTACTTCTGCCA‑GTTGCTGTTCCAGCGCCCCTTTGCGGCTGGTGGTATCGCGGTCGATACGAGAAATGGGCACGCCGGGGAACAACGGCCCGAGCGTCTGTTCAAGCTGTTCGGGGGCCACCCCAACGGGGCCCAG‑GGGCGTGGAACCCCAGGAAGGGCCATGCGCGGGGCCCCGG‑CCCTGACTGGACA‑AGCGGTGGCAg                                                                                                                                                                      >  2:247182/1‑228 (MQ=255)
                                                                      cgcgcCG‑CCGCGATGTACTTCTGCCA‑GTTGCTGTTCCAGCGCCCCTTTGCGGCTGGTGGTATCGCGGTCGATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCCCGGGGACCAG‑GGGCGTGGAACCGCAGGAAGGGCACTG‑GCGCGGCCCCGGACGCTAGCTGGCAC‑AGGGGGGGCCACGCAgcgtgggcag                                                                                                                                                        >  2:140999/1‑219 (MQ=255)
                                                                      cgcgcCG‑CCGCGATGTACTTCTGCCA‑GTTGCTGTTCCAGCGCCCCTTTGCGGCTGGTGGTATCGCGGTCGATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCCCGGGGACCAG‑GGGCGTGGAACCCCAGGAAGGGCCCCG‑GCCCGGCCCCGGGCGCTTAC‑GGCACAAGTGGGGGCAgcgccaggggggggg                                                                                                                                                        >  2:593581/1‑217 (MQ=255)
                                                                      cgcgcCG‑CCGCGATGTACTTCTGCCA‑GTTGCTGTTCCAGCGCCCCTTTGCGGCTGGTGGTATCGCGGTCGATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGGACCAG‑GTGCGTGGAACCGCAGGAAGGGCACTG‑GCGCGGCACCGGACGCTGACTGTAAC‑AGCGGTGGCAGCGCATGTGGtgctg                                                                                                                                                        >  1:135446/1‑228 (MQ=255)
                                                                      cccgccgcccgCGTTGTTCTTTGGCCG‑GTTGCTTGTCCCGGGCCCCTTTGGGGCTGGTGGTATCGCGGGCGATTCGAGAAAGGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGGACCAG‑GTGCGTGGAACCGCAGGAAGGGCACTG‑GCGCGGCACCGGACGCTGACTGTCAC‑AGTGGTGGCAGCGCAGATGGtgct                                                                                                                                                         <  2:420566/226‑1 (MQ=255)
                                                                                cggATGTTCTGCTGCCA‑GTTTCTGTACCGGCTCCCTTTTGCGGCTGGTGTACTCGCGGTAGAGCTGGGAAATGGGTCCGCCCGGGACCAACGGCGCGGGCTCCTGTTCGAGCGGTTCGGTGCCCAGCCCCACGGGGACCAG‑GTGCTTGGAACCGCAGGAAGGGCACTG‑GCGCGGCACCGGACGCTGACTGTCAC‑AGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATg                                                                                                                                               <  2:473338/226‑1 (MQ=255)
                                                                                  tatgtctttCTGCCC‑GTTGCTGTTCCGTCGCCCCTTTGCGGCTGGTGGTCTCGCGGTGGATCCGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGGACCAG‑GTGCGTGGAACCGCAGGAAGGGCACTG‑GCGCGGCACCGGACGCTGACTGTCAC‑AGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCa                                                                                                                                             <  2:870/221‑1 (MQ=255)
                                                                                   aTGTACTTCTGCCT‑GTTGCGGTTCCAGCGCCCCGTTGCGGCTGGTGGTATCGCGGTCGATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGGACCAG‑GTGCGTGGAACCGCAGGAAGGGCACTG‑GCGCGGCACCGGACGCTGACTGTCAC‑AGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCa                                                                                                                                             <  1:561460/227‑1 (MQ=255)
                                                                                    tGTATTTCTGCCA‑TTTGTTGTTCCCGGCCCCCTTTGCGGCTGGTGGTATCGCGGGCGATACGAGAAATGGGCACCGCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTGCGGTGCCCAGCCCCACGGGGACCAG‑GTGCGTGGAACCGCAGGAAGGGCACTG‑GCGCGGCACCGGACGCTGACTGTCAC‑AGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGc                                                                                                                                           <  2:9991/228‑1 (MQ=255)
                                                                                     gTACTTC‑GCCAGGTTGCTGTGCCAGCGCCCCTTTGCGGGTGGGGGGATCGGGGGCGATACGAGAAAGGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTGCAAGCTGTTCGGTGCCCAGCCCCACGGGGACCAG‑GTGCGTGGAACCGCAGGAAGGGCACTG‑GCGCGGCACCGGACGCTGACTGTCAC‑AGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGc                                                                                                                                           <  2:6865/227‑1 (MQ=255)
                                                                                                a‑GTTGCTGTTCCAGCGCCCCTTTGCGGCTGGTGGTATCGCGGTCGATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGGACCAG‑GTGCGTGGAACCGCAGGAAGGGCACTG‑GCGCGGCACCGGACGCTGACTGTCAC‑AGTGGTGGCAGCGCCGATGGTGCTGCGCCTGTTGAAGCGTGTAGGAGTGa                                                                                                                               >  1:387554/1‑228 (MQ=255)
                                                                                                    tGCTGTTCCAGCGCCCCTTTGCGGCTGGTGGTATCGCGGTCGATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGGACCAG‑GTGCGTGGAACCGCAGGAAGGGCACTG‑GCGCGGCACCGGACCCTGACTGTCAC‑AGTGGGGGCAGCGCAGATGGGGCTGCGCCCGATGCACCGCGTAGAAGGGATCg                                                                                                                            >  1:650289/1‑228 (MQ=255)
                                                                                                                gcCCCTTTGGGGCTGGGTGTATCGGGGTCGATACGAGAAATGGGCACGCCGGGGAAACACGGCGCGAGCGTCTGTTCAAGCTGGTCGGTGCCCAGCCCCACGGGGACCAG‑GTGCGTGGAACCGCAGGAAGGGCACTG‑GCGCGGCACCGGACGCTGACTGTCAC‑AGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCa                                                                                                                 <  1:408368/227‑1 (MQ=255)
                                                                                                                   ccTTTGCGGCTGGTGGTATCGCGGTCGATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCCCGGGGACCAG‑GTGCGGGGAACCGCAGGGAGGGCACTT‑GCCCGGCACCGGCGGCTGCCTGTCAC‑GGTGGCGGCAGCGCGGATGGTGCTGCCCCTCAGTCAGCGGGGAGGAGTGATCGCAACGGGGGGAtggg                                                                                                             >  2:648683/1‑225 (MQ=255)
                                                                                                                      ttGCGGCTGGGGGTATCGCGGTCGATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGGACCAG‑GTGCGTGGAACCGCAGGAAGGGCACTG‑GCGCGGCACCGGACGCTGACTGTCAC‑AGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGc                                                                                                           <  1:454046/227‑1 (MQ=255)
                                                                                                                               ggggtATCGCGGTCGATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGGACCAG‑GTGCGTGGAACCGCAGGAAGGGCACTG‑GCGCGGCACCGGACGCTGACTGTCAC‑AGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCca                                                                                                 <  1:593581/226‑1 (MQ=255)
                                                                                                                                         ggTCGATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGGACCAG‑GTGCGTGGAACCGCAGGAAGGGCACTG‑GCGCGGCACCGGACGCTGACTGTCAC‑AGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGCCGTGGc                                                                                       >  1:392444/1‑228 (MQ=255)
                                                                                                                                          gTCGATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGGACCAG‑GTGCGTGGAACCGCAGGAAGGGCACTG‑GCGCGGCACCGGACGCTGACTGTCAC‑AGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCCCAGTCGTGGca                                                                                      >  1:92536/1‑228 (MQ=255)
                                                                                                                                          gTCGATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGGACCAG‑GTGCGTGGAACCGCAGGAAGGGCACTG‑GCGCGGCACCGGACGCTGACTGTCAC‑AGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGACGTGGca                                                                                      >  1:439282/1‑228 (MQ=255)
                                                                                                                                             gATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGGACCAG‑GTGCGTGGAACCGCAGGAAGGGCACTG‑GCGCGGCACCGGACGCTGACTGTCAC‑AGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGACGGGGcaca                                                                                    >  1:602303/1‑227 (MQ=255)
                                                                                                                                             gATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGGACCAG‑GTGCGTGGAACCGCAGGAAGGGCACTG‑GCGCGGCACCGGACGCTGACTGTCAC‑AGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGGGTAGTAGGGATCGCAACGTGGGCATTAGGCAATCCAGCCCCAGGCGGGGcaaga                                                                                   >  2:97985/1‑225 (MQ=255)
                                                                                                                                                         gggCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGGACCAG‑GGGCGTGGAACCGCAGGAAGGGCACTG‑GCGCGGCACCGGACGCTGACTGTCAC‑AGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCGTGGCACAGCAGTGCAGgcgc                                                                       >  1:575856/1‑228 (MQ=255)
                                                                                                                                                                      aacaacGGGGCGAGCGGCTGTTCAAGCTGTTCGGTGCCCAGCACCACGGGGACCAG‑GTGCGTGGAACCGCAGGAAGGGCACTG‑GCGCGGCACCGGACGCTGACTGTCAC‑AGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCGTGGCACAGCAGTGCAGGCGCAAAGCCAcggc                                                            <  1:563286/226‑1 (MQ=255)
                                                                                                                                                                         aacGGCGCGGGCGGCTGTTCAAGCGGGTCGGGGCCCAGCCCCAGGGGGACCAG‑GTGGGTGGAACCGCGGGAGGGGCACGG‑GCGCGGCACCGGACGCTGACTGTCAC‑AGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCGTGGCACAGCAGTGCAGGCGCAAAGCCACGGCGGtt                                                        <  2:156138/227‑1 (MQ=255)
                                                                                                                                                                                     ggtgtttCATGCTGTTCGGGGCCCAGCCCCACGGGGACCAG‑GGGGGTGGAACCGCAGGAAGGGCACTG‑GCGCGGCACCGGACGCTGACCGTCAC‑AGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCGTGGCACAGCAGTGCAGGCGCAAAGCCACGGCGGTTAAGAAAGAGAATg                                           <  2:459307/223‑1 (MQ=255)
                                                                                                                                                                                       ttGTTCAAGCGTTGCGGGGCCCAGCCCCACGGGGACCAG‑GGGCGGGGAAACGCAGGAAGGGCACTG‑GCGCGGCACCGGACGCTGACGGGCAC‑AGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCGTGGCACAGCAGTGCAGGCGCAAAGCCACGGCGGTTAAGAAAGAGAATGAc                                         <  2:602303/227‑1 (MQ=255)
                                                                                                                                                                                       tggtccacgCTGTTCGGGGCCCGGCCCCACGGGGACCAG‑GTGCGTGGAACCGCAGGAAGGGCAGTG‑GCGCGGCACCGGACGCTGACTGTCAC‑AGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCGTGGCACAGCAGTGCAGGCGCAAAGCCACGGCGGTTAAGAAAGAGAATGAc                                         <  2:553702/220‑1 (MQ=255)
                                                                                                                                                                                               ggtggtCGGTGCCCACCCCCCCGGGGAACGG‑GTGGGTGGAACCGCAGGAGGGGCACTG‑GGGCGGCCCCGGTCGCTGACTGTCAC‑AGGGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCGTGGCACAGCAGTGCAGGCGCAAAGCCACGGCGGTTAAGAAAGAGAATGACCTGGTTAt                                 <  2:24167/223‑1 (MQ=255)
                                                                                                                                                                                               agggTTCGGGGCTCGGCCCCACGGGGACGAG‑GGTCGTGGAACCGGAGGAAGGGCGCTG‑GCGCGGCACCGGAGCGTGACTGTCAC‑AGTGGTGGCAGCGCAGATGGGGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCGTGGCACAGCAGTGCAGGCGCAAAGCCACGGCGGTTAAGAAAGAGAATGACCTGGTTAt                                 <  2:135446/225‑1 (MQ=255)
                                                                                                                                                                                                cTGTTCGGTGCCCAGCCCCCAGGGGACCAG‑GTGCGTGGAACCGCAGGAAGGGCACTG‑GCGCGGCACCGGACGCTGACTGTCAC‑AGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCGTGGCACAGCAGTGCAGGCGCAAAGCCACGGCGGTTAAGAAAGAGAATGACCTGGTTAt                                 <  1:97985/227‑1 (MQ=255)
                                                                                                                                                                                                 tggtCGGGGGCCAGCCCCACGGAGACCAGCGTG‑GGGGAACCGAAGGAAGGGCAGTG‑GCGCGGCACCGGACGCTGACTGGCAC‑AGTGGGGGCGGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGAGCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCGTGGCACAGCAGTGCCGGCGCAAAGCCACGGCGGTTAAGAAAGAGAATGACCTGGTTATc                                <  2:33763/224‑1 (MQ=255)
                                                                                                                                                                                                  gTTCGGTGCCCAGCCCGCCGGGGCCCAG‑GTGGGGGGAACCGGAGGGAGGGCACTG‑GGGCGGCACCGGACGCTGAATGTCAC‑AGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATCCGGCAATCCAGCCACAGTCGTGGCACAGCAGTGCAGGCGCAAAGCCACGGCGGTTAAGAAAGAGAATGACCTGGTTATCAg                              <  2:636599/228‑1 (MQ=255)
                                                                                                                                                                                                    tCGGTGCCCAGCCCCACGGGGACCAG‑GTGCGTGGAACCGCAGGAAGGGCACTG‑GCGCGGCACCGGACGCTGACTGTCAC‑AGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGTGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCGTGGCACAGCAGTGCAGGCGCAAAGCCACGGCGGTTACGAAAGGGCATGCACTCGTTTTCAGcc                            >  1:564409/1‑228 (MQ=255)
                                                                                                                                                                                                       gggCCCAGCCCCACGGGGCCCAG‑GGGGGTGGGACCGGAGGAAGGGCACTG‑GGGGGGGACCGGAAGCTGACTGGCAC‑AGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGGAGGGATCGCAAAGTGGGCATTCGGCAATCCAGCCACAGTCGTGGCACAGCAGTGCAGGCGCAAAGCCACGGCGGTTAAGAAAGAGAATGACCTGGTTATCAGCCTGt                         <  2:387554/226‑1 (MQ=255)
                                                                                                                                                                                                             aGCCCCACGGGCACCAG‑GTTCGGGGAACCGCAGGAAGGGCCCGT‑GCGCGGCACCGGACGCTGACTGTCAC‑AGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCGTGGCACAGCAGTGCAGGCGCAAAGCCACGGCGGTTAAGAAAGAGAATGACCTGGTTATCAGCCTGTAAATGc                   <  1:209114/228‑1 (MQ=255)
                                                                                                                                                                                                                cccACGGGGACCAG‑GGGCGTGGAACCGGAGGAAGGGCACTG‑GCGCGGCACCGGACGCTGACTGTCAC‑AGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCGTGGCACAGCAGTGCAGGCGCAAAGCCACGGCGGTTAAGAAAGAGAATGACCTGGTTATCAGCCTGTAAATGCTg                 <  1:66217/227‑1 (MQ=255)
                                                                                                                                                                                                                  cccGGGGACCAG‑GGGCGTGGAACCGCAGGACGGGCACTG‑GCGCGGCCCCGGACGCTGACTGTCAC‑AGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCGTGGCACAGCAGTGCAGGCGCAAAGCCACGGCGGTTAAGAAAGAGAATGACCTGGTTATCAGCCTGTAAATGCTGgcg              <  2:466530/226‑1 (MQ=255)
                                                                                                                                                                                                                           cAG‑GGGCGTGGAACGGCAGGAAGGGCACTG‑GCGCGGCACCGGACGCTGACTGTCAC‑AGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCGTGGCACAGCAGTGCAGGCGCAAAGCCACGGCGGTTAAGAAAGAGAATGACCTGGTTATCAGCCTGTAAATGCTGGCGCATACGa       <  1:13106/226‑1 (MQ=255)
                                                                                                                                                                                                                               gtgggtggAACCGCGGGAGGGGCACTG‑GCGGGGCACCGGACGCTGACTGTCAC‑AGGGGTGGCAGCGCAGAGGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCGTGGCACAGCAGTGCAGGCGCAAAGCCACGGCGGTTAAGAAAGAGAATGACCTGGTTATCAGCCTGTAAATGCTGGCGCATACGAGTGAt  <  2:469469/228‑1 (MQ=255)
                                                                                                                                                                                                                                |                                                                                                                                                                                                                                    
GTAATGCAACCAGCGTCACATCCGGGAAATGGTGACCTTTCGCCAGCATTTGTGTACCAATCAAAATCCGCGCGCCG‑CCGCGATGTACTTCTGCCA‑GTTGCTGTTCCAGCGCCCCTTTGCGGCTGGTGGTATCGCGGTCGATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGGACCAG‑GTGCGTGGAACCGCAGGAAGGGCACTG‑GCGCGGCACCGGACGCTGACTGTCAC‑AGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCGTGGCACAGCAGTGCAGGCGCAAAGCCACGGCGGTTAAGAAAGAGAATGACCTGGTTATCAGCCTGTAAATGCTGGCGCATACGAGTGAT  >  NC_000913/4125142‑4125589

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