Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,125,362 T→G 23.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,3620TG23.1% 95.6 / 6.0 57H483P (CAC→CCC) priAPrimosome factor n' (replication factor Y)
Reads supporting (aligned to +/- strand):  ref base T (16/25);  new base G (9/7);  total (25/32)
Fisher's exact test for biased strand distribution p-value = 3.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

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

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