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

GTAATGCAACCAGCGTCACATCCGGGAAATGGTGACCTTTCGCCAGCATTTGTGTACCAATCAAAATCCGCGCGCCGCCGCGATGTACTTCTGCCAGTTGCTGTTCCAGCGCCCCTTTGCGGCTGGTGGTATCGCGGTCGATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGGACCAGGTGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCACCGGACGCTGACTGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCGTGGCACAGCAGTGCAGGCGCAAAGCCACGGCGGTTAAGAAAGAGAATGACCTGGTTATCAGCCTGTAAATGCTGGCGCATACGAGTGAT  >  NC_000913/4125141‑4125588
                                                                                                                                                                                                                                |                                                                                                                                                                                                                                    
gttaggCACCCAGCGTCACCTCCGGGAAAGGGTGCCCTTCCGGCACCATTTGGTTACCAATCAAAATCCGCGCGCCGCCGCGATGTTCTTCTGCCATTTGCGGTCCCAGCGCCCCTTTGCGGCTGGTGGTATCGCGGTCGATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGGACCAGGTGCGTgg                                                                                                                                                                                                                                <  1:360359/223‑1 (MQ=255)
      cacccAGCGCCACATCCGGGAAATGGTGACCTTTCGCCAGCATTTGGGTACCAATCAAAATCCGCGCGCCGCCGGGATGTACTTCTGCCAGTTGCTGTTCCAGCGCCCCTTTGCGGCTGGTGGTATCGCGGTCGATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGGACCAGGTGCGTGGAACCGc                                                                                                                                                                                                                          <  2:564440/225‑1 (MQ=255)
       aaCCAGCGTCACATCCGGGAAATGGTGACCTTTCGCCAGCATTTGTGTACCAATCAAAATCCGCGCGCCGCAGCGATGTACTTCTGCCAGTTGCTGTTCCAGCGCCCCTTTGCGGCTGGCGGTATCGCGGTCGATACGAGAAATGGGCACGCCGGGGACCAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCAAGCCCCACGGGGACCAGGGGCGTGGAccc                                                                                                                                                                                                                            >  2:459307/1‑225 (MQ=255)
         ccAGCGTCACATCCGGGAAATGGTGACCTTTCGCCAGCATTTGTGTACCAATCAAAATCCGCGCGCCGCCGCGATGTACTTCTGCCAGTTGCTGTTCCAGCGCCCCTTTGCGGCTGGTGGTATCGCGGTCGATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGGGTCTGTTCAAGGTGTTCGGTGCCCAGGCCCCCGGGGGACAGGGGCGGGGACACGCAgg                                                                                                                                                                                                                       >  1:398903/1‑228 (MQ=255)
         ccAGCGTCACATCCGGGAAATGGTGACCTTTCGCCAGCATTTGTGTACCAATCAAAATCCGCGCGCCGCCGCGATGTACTTCTGCCAGTTGCTGTTCCAGCGCCCCTTTGCGGCTGGTGGTATCGCGGTCGATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGGACCAGGTGCGTGGAACCGCAg                                                                                                                                                                                                                        >  2:9991/1‑227 (MQ=255)
         ccAGCGTCACATCCGGGAAATGGTGACCTTTCGCCAGCATTTGTGTACCAATCAAAATCCGCGCGCCGCCGCGATGTACTTCTGCCAGTTGCTGTTCCAGCGCCCCTTTGCGGCTGGTGGTATCGCGGTCGATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGGCTGTTCAAGCTGTTCGGGGGCCCGCCCCCAGGGGACCAGGGGCGCGGAAACCCAgg                                                                                                                                                                                                                       >  1:563286/1‑228 (MQ=255)
             cGTCACATCCGGGAAATGGTGACCTTTCGCCAGCATTTGTGTACCAATCAAAATCCGCGCGCCGCCGCGATGTACTTCTGCCAGTTGCTGTTCCAGCGCCCCTTTGCGGCTGGTGGTATCGCGGTCGATACGAGAAATGGGCACGCCGGGGAACAACGGCGGGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGGACCAGGTGCGTGGAACCGCAGGAAgg                                                                                                                                                                                                                   >  2:180653/1‑228 (MQ=255)
                                     tttCGCCCGCTTTTGTGTACCAATCAAAATCCGCGCGCCGCCGCGATGTACTTCTGCCAGTTGCTGTTCCAGCGCCCCTTTGCGGCTGGTGGTATCGCGGTCGATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGGACCAGGTGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCACCGGACGCt                                                                                                                                                                                           <  2:570159/227‑1 (MQ=255)
                                     atccGCCAGCATTTGTGCACCCATCAAAATCCGCGCGCCGCCGCGATGTGCTTCTGCCAGTTGCTGTTCCAGCGCCCCTTTGCGGCTGGTGGTATCGGGGTCGATTCGAGAAAGGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGGACCAGGTGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCACCGGACGCt                                                                                                                                                                                           <  2:571037/224‑1 (MQ=255)
                                      ttCGCCAGCATTGGTGTCCCAATCAAAATCCGCGCGCCGGCGCGATGTACTTCGGCCAGTTGCTGTTCCAGCGCCCCTTTGCGGCTGGTGGTATCGCGGTCGATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGGACCAGGTGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCACCGGACGCTGa                                                                                                                                                                                         <  2:269163/228‑1 (MQ=255)
                                            aGCATTTGGGTACCAATCAAAATCCGCGCGCCGGCGCGATGTACTTCGGCCAGTTGCTGTTCCAGCGCCCCTTTGCGGCTGGTGGTATCGCGGTCGATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGGTCGGTGCCCAGCCCCACGGGGACCAGGTGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCACCGGACGCTGACTGTc                                                                                                                                                                                    <  2:398903/227‑1 (MQ=255)
                                              cATTTGTGTACCAATCAAAATCCGCGCGCCGCCGCGATGTACTTCTGCCAGTTGCTGTTCCAGCGCCCCTTTGCGGCTGGTGGTATCGCGGTCGATACGAGAAATGGGCACGCCGGGGCACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGGGCCCAGCCCCACGGGGACCAGGGGCGTGGAAACGCAGGAAGGGCCCGGGCGGGGCACCGGACGCGGAGTGTcaaa                                                                                                                                                                                 >  1:269163/1‑226 (MQ=255)
                                              cATTTGTGTACCAATCAAAATCCGCGCGCCGCCGCGATGTACTTCTGCCAGTTGCTGTTCCAGCGCCCCCTTGCGGCTGGTGGTATCGCGGTCGATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTCCGGTGCCCCGCCCCCCGGGGACCAGGTGCGTGGAACCCCCGGAAGGGCCCTGGCCCGGCCCCGGAccgtgaacgtacca                                                                                                                                                                                 >  1:253499/1‑215 (MQ=255)
                                                      tACCAATCAAAACCCCGGCGCCGGCGGGGGGTGCTTTCGCAAGTTGCTGTTCCAGCGCCCTTTGGCGGCGGGGGGTATCGCGGTCGATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGGACCAGGTGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCACCGGACGCTGACTGTCACAGTGGTGGc                                                                                                                                                                        <  1:536837/228‑1 (MQ=255)
                                                        ccAATCAAAATCCGCGCGCCGCCGCGATGTACTTCTGCCAGTTGCTGTTCCAGCGCCCCTTTGCGGCTGGTGGTATCGCGGTCGATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCCGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGGACCAGGGGCGAGGAACCGCCGGAAGGGCACTGGGGCGGGCCCGGAAGCGGCCTGGCACAGGGGGGGGAg                                                                                                                                                                      >  1:575139/1‑228 (MQ=255)
                                                        ccAATCAAAATCCGCGCGCCGCCGCGATGTACTTCTGCCAGTTGCTGTTCCAGCGCCCCTTTGCGGCTGGTGGTATCGCGGTCGATACGAGAAATGGGCACGCCGGGGAACAACGGCCCGAGCGTCTGTTCAAGCTGTTCGGGGGCCACCCCAACGGGGCCCAGGGGCGTGGAACCCCAGGAAGGGCCATGCGCGGGGCCCCGGCCCTGACTGGACAAGCGGTGGCAg                                                                                                                                                                      >  1:247182/1‑228 (MQ=255)
                                                                      cgcgcCGCCGCGATGTACTTCTGCCAGTTGCTGTTCCAGCGCCCCTTTGCGGCTGGTGGTATCGCGGTCGATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCCCGGGGACCAGGGGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCCCCGGACGCTAGCTGGCACAGGGGGGGCCACGCAgcgtgggcag                                                                                                                                                        >  1:140999/1‑219 (MQ=255)
                                                                      cgcgcCGCCGCGATGTACTTCTGCCAGTTGCTGTTCCAGCGCCCCTTTGCGGCTGGTGGTATCGCGGTCGATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCCCGGGGACCAGGGGCGTGGAACCCCAGGAAGGGCCCCGGCCCGGCCCCGGGCGCTTACGGCACAAGTGGGGGCAgcgccaggggggggg                                                                                                                                                        >  1:593581/1‑217 (MQ=255)
                                                                      cgcgcCGCCGCGATGTACTTCTGCCAGTTGCTGTTCCAGCGCCCCTTTGCGGCTGGTGGTATCGCGGTCGATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGGACCAGGTGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCACCGGACGCTGACTGTAACAGCGGTGGCAGCGCATGTGGtgctg                                                                                                                                                        >  2:135446/1‑228 (MQ=255)
                                                                      cccgccgcccgCGTTGTTCTTTGGCCGGTTGCTTGTCCCGGGCCCCTTTGGGGCTGGTGGTATCGCGGGCGATTCGAGAAAGGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGGACCAGGTGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCACCGGACGCTGACTGTCACAGTGGTGGCAGCGCAGATGGtgct                                                                                                                                                         <  1:420566/226‑1 (MQ=255)
                                                                                cggATGTTCTGCTGCCAGTTTCTGTACCGGCTCCCTTTTGCGGCTGGTGTACTCGCGGTAGAGCTGGGAAATGGGTCCGCCCGGGACCAACGGCGCGGGCTCCTGTTCGAGCGGTTCGGTGCCCAGCCCCACGGGGACCAGGTGCTTGGAACCGCAGGAAGGGCACTGGCGCGGCACCGGACGCTGACTGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATg                                                                                                                                               <  1:473338/226‑1 (MQ=255)
                                                                                  tatgtctttCTGCCCGTTGCTGTTCCGTCGCCCCTTTGCGGCTGGTGGTCTCGCGGTGGATCCGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGGACCAGGTGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCACCGGACGCTGACTGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCa                                                                                                                                             <  1:870/221‑1 (MQ=255)
                                                                                   gtttttctttttcccgTTGCTGTTCCCGCGCCCCTTTGGGGTTGGTGGTGTCGCGGTCGGTACGAGAAATGGGCACGCCGGGCAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGGACCAGGTGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCACCGGACGCTGACTGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCa                                                                                                                                             <  1:557521/213‑1 (MQ=255)
                                                                                   aTGTACTTCTGCCTGTTGCGGTTCCAGCGCCCCGTTGCGGCTGGTGGTATCGCGGTCGATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGGACCAGGTGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCACCGGACGCTGACTGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCa                                                                                                                                             <  2:561460/227‑1 (MQ=255)
                                                                                    tGTATTTCTGCCATTTGTTGTTCCCGGCCCCCTTTGCGGCTGGTGGTATCGCGGGCGATACGAGAAATGGGCACCGCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTGCGGTGCCCAGCCCCACGGGGACCAGGTGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCACCGGACGCTGACTGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGc                                                                                                                                           <  1:9991/228‑1 (MQ=255)
                                                                                     gTACTTCGCCAGGTTGCTGTGCCAGCGCCCCTTTGCGGGTGGGGGGATCGGGGGCGATACGAGAAAGGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTGCAAGCTGTTCGGTGCCCAGCCCCACGGGGACCAGGTGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCACCGGACGCTGACTGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGc                                                                                                                                           <  1:6865/227‑1 (MQ=255)
                                                                                                aGTTGCTGTTCCAGCGCCCCTTTGCGGCTGGTGGTATCGCGGTCGATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGGACCAGGTGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCACCGGACGCTGACTGTCACAGTGGTGGCAGCGCCGATGGTGCTGCGCCTGTTGAAGCGTGTAGGAGTGa                                                                                                                               >  2:387554/1‑228 (MQ=255)
                                                                                                    tGCTGTTCCAGCGCCCCTTTGCGGCTGGTGGTATCGCGGTCGATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGGACCAGGTGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCACCGGACCCTGACTGTCACAGTGGGGGCAGCGCAGATGGGGCTGCGCCCGATGCACCGCGTAGAAGGGATCg                                                                                                                            >  2:650289/1‑228 (MQ=255)
                                                                                                                gcCCCTTTGGGGCTGGGTGTATCGGGGTCGATACGAGAAATGGGCACGCCGGGGAAACACGGCGCGAGCGTCTGTTCAAGCTGGTCGGTGCCCAGCCCCACGGGGACCAGGTGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCACCGGACGCTGACTGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCa                                                                                                                 <  2:408368/227‑1 (MQ=255)
                                                                                                                   ccTTTGCGGCTGGTGGTATCGCGGTCGATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCCCGGGGACCAGGTGCGGGGAACCGCAGGGAGGGCACTTGCCCGGCACCGGCGGCTGCCTGTCACGGTGGCGGCAGCGCGGATGGTGCTGCCCCTCAGTCAGCGGGGAGGAGTGATCGCAACGGGGGGAtggg                                                                                                             >  1:648683/1‑225 (MQ=255)
                                                                                                                      ttGCGGCTGGGGGTATCGCGGTCGATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGGACCAGGTGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCACCGGACGCTGACTGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGc                                                                                                           <  2:454046/227‑1 (MQ=255)
                                                                                                                               ggggtATCGCGGTCGATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGGACCAGGTGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCACCGGACGCTGACTGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCca                                                                                                 <  2:593581/226‑1 (MQ=255)
                                                                                                                                         ggTCGATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGGACCAGGTGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCACCGGACGCTGACTGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGCCGTGGc                                                                                       >  2:392444/1‑228 (MQ=255)
                                                                                                                                          gTCGATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGGACCAGGTGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCACCGGACGCTGACTGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCCCAGTCGTGGca                                                                                      >  2:92536/1‑228 (MQ=255)
                                                                                                                                          gTCGATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGGACCAGGTGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCACCGGACGCTGACTGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGACGTGGca                                                                                      >  2:439282/1‑228 (MQ=255)
                                                                                                                                             gATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGGACCAGGTGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCACCGGACGCTGACTGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGACGGGGcaca                                                                                    >  2:602303/1‑227 (MQ=255)
                                                                                                                                             gATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGGACCAGGTGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCACCGGACGCTGACTGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGGGTAGTAGGGATCGCAACGTGGGCATTAGGCAATCCAGCCCCAGGCGGGGcaaga                                                                                   >  1:97985/1‑225 (MQ=255)
                                                                                                                                                         gggCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGGACCAGGGGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCACCGGACGCTGACTGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCGTGGCACAGCAGTGCAGgcgc                                                                       >  2:575856/1‑228 (MQ=255)
                                                                                                                                                                      aacaacGGGGCGAGCGGCTGTTCAAGCTGTTCGGTGCCCAGCACCACGGGGACCAGGTGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCACCGGACGCTGACTGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCGTGGCACAGCAGTGCAGGCGCAAAGCCAcggc                                                            <  2:563286/226‑1 (MQ=255)
                                                                                                                                                                         aacGGCGCGGGCGGCTGTTCAAGCGGGTCGGGGCCCAGCCCCAGGGGGACCAGGTGGGTGGAACCGCGGGAGGGGCACGGGCGCGGCACCGGACGCTGACTGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCGTGGCACAGCAGTGCAGGCGCAAAGCCACGGCGGtt                                                        <  1:156138/227‑1 (MQ=255)
                                                                                                                                                                                     ggtgtttCATGCTGTTCGGGGCCCAGCCCCACGGGGACCAGGGGGGTGGAACCGCAGGAAGGGCACTGGCGCGGCACCGGACGCTGACCGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCGTGGCACAGCAGTGCAGGCGCAAAGCCACGGCGGTTAAGAAAGAGAATg                                           <  1:459307/223‑1 (MQ=255)
                                                                                                                                                                                       ttGTTCAAGCGTTGCGGGGCCCAGCCCCACGGGGACCAGGGGCGGGGAAACGCAGGAAGGGCACTGGCGCGGCACCGGACGCTGACGGGCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCGTGGCACAGCAGTGCAGGCGCAAAGCCACGGCGGTTAAGAAAGAGAATGAc                                         <  1:602303/227‑1 (MQ=255)
                                                                                                                                                                                       tggtccacgCTGTTCGGGGCCCGGCCCCACGGGGACCAGGTGCGTGGAACCGCAGGAAGGGCAGTGGCGCGGCACCGGACGCTGACTGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCGTGGCACAGCAGTGCAGGCGCAAAGCCACGGCGGTTAAGAAAGAGAATGAc                                         <  1:553702/220‑1 (MQ=255)
                                                                                                                                                                                               ggtggtCGGTGCCCACCCCCCCGGGGAACGGGTGGGTGGAACCGCAGGAGGGGCACTGGGGCGGCCCCGGTCGCTGACTGTCACAGGGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCGTGGCACAGCAGTGCAGGCGCAAAGCCACGGCGGTTAAGAAAGAGAATGACCTGGTTAt                                 <  1:24167/223‑1 (MQ=255)
                                                                                                                                                                                               agggTTCGGGGCTCGGCCCCACGGGGACGAGGGTCGTGGAACCGGAGGAAGGGCGCTGGCGCGGCACCGGAGCGTGACTGTCACAGTGGTGGCAGCGCAGATGGGGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCGTGGCACAGCAGTGCAGGCGCAAAGCCACGGCGGTTAAGAAAGAGAATGACCTGGTTAt                                 <  1:135446/225‑1 (MQ=255)
                                                                                                                                                                                                cTGTTCGGTGCCCAGCCCCCAGGGGACCAGGTGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCACCGGACGCTGACTGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCGTGGCACAGCAGTGCAGGCGCAAAGCCACGGCGGTTAAGAAAGAGAATGACCTGGTTAt                                 <  2:97985/227‑1 (MQ=255)
                                                                                                                                                                                                 tggtCGGGGGCCAGCCCCACGGAGACCAGCGTGGGGGAACCGAAGGAAGGGCAGTGGCGCGGCACCGGACGCTGACTGGCACAGTGGGGGCGGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGAGCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCGTGGCACAGCAGTGCCGGCGCAAAGCCACGGCGGTTAAGAAAGAGAATGACCTGGTTATc                                <  1:33763/224‑1 (MQ=255)
                                                                                                                                                                                                  gTTCGGTGCCCAGCCCGCCGGGGCCCAGGTGGGGGGAACCGGAGGGAGGGCACTGGGGCGGCACCGGACGCTGAATGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATCCGGCAATCCAGCCACAGTCGTGGCACAGCAGTGCAGGCGCAAAGCCACGGCGGTTAAGAAAGAGAATGACCTGGTTATCAg                              <  1:636599/228‑1 (MQ=255)
                                                                                                                                                                                                    tCGGTGCCCAGCCCCACGGGGACCAGGTGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCACCGGACGCTGACTGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGTGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCGTGGCACAGCAGTGCAGGCGCAAAGCCACGGCGGTTACGAAAGGGCATGCACTCGTTTTCAGcc                            >  2:564409/1‑228 (MQ=255)
                                                                                                                                                                                                    tCGGTGCCCAGCCCCACGGGGACCAGGTGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCACCGGACGCTGACTGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCACTCAAGCCACAGTCGTGGCACAGCAGTGCAGGCGCCAAGCCACGGCGGTTcagtaaaggaagagcccggtagtcagcc                            >  2:640012/1‑206 (MQ=255)
                                                                                                                                                                                                    tCGGTGCCCAGCCCCACGGGGACCAGGTGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCACCGGACGCTGACTGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCGTGGCACAGCAGTGCAGGCGCAAAGCCACGGCGGTCAAGAAagaggcaaaacagggtattcacc                            >  1:539842/1‑209 (MQ=255)
                                                                                                                                                                                                       gggCCCAGCCCCACGGGGCCCAGGGGGGTGGGACCGGAGGAAGGGCACTGGGGGGGGACCGGAAGCTGACTGGCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGGAGGGATCGCAAAGTGGGCATTCGGCAATCCAGCCACAGTCGTGGCACAGCAGTGCAGGCGCAAAGCCACGGCGGTTAAGAAAGAGAATGACCTGGTTATCAGCCTGt                         <  1:387554/226‑1 (MQ=255)
                                                                                                                                                                                                             aGCCCCACGGGCACCAGGTTCGGGGAACCGCAGGAAGGGCCCGTGCGCGGCACCGGACGCTGACTGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCGTGGCACAGCAGTGCAGGCGCAAAGCCACGGCGGTTAAGAAAGAGAATGACCTGGTTATCAGCCTGTAAATGc                   <  2:209114/228‑1 (MQ=255)
                                                                                                                                                                                                                cccACGGGGACCAGGGGCGTGGAACCGGAGGAAGGGCACTGGCGCGGCACCGGACGCTGACTGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCGTGGCACAGCAGTGCAGGCGCAAAGCCACGGCGGTTAAGAAAGAGAATGACCTGGTTATCAGCCTGTAAATGCTg                 <  2:66217/227‑1 (MQ=255)
                                                                                                                                                                                                                  cccGGGGACCAGGGGCGTGGAACCGCAGGACGGGCACTGGCGCGGCCCCGGACGCTGACTGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCGTGGCACAGCAGTGCAGGCGCAAAGCCACGGCGGTTAAGAAAGAGAATGACCTGGTTATCAGCCTGTAAATGCTGgcg              <  1:466530/226‑1 (MQ=255)
                                                                                                                                                                                                                         aCCAGGTGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCACCGGACGCTGACTGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCGTGGCACAGCAGTGCAGGCGCAAAGCCACGGCGGTTTAGAAAAGGAATAACCCGGTTATCCGCCTTTaaaacgcgggccctacg        >  2:330108/1‑213 (MQ=255)
                                                                                                                                                                                                                           cAGGGGCGTGGAACGGCAGGAAGGGCACTGGCGCGGCACCGGACGCTGACTGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCGTGGCACAGCAGTGCAGGCGCAAAGCCACGGCGGTTAAGAAAGAGAATGACCTGGTTATCAGCCTGTAAATGCTGGCGCATACGa       <  2:13106/226‑1 (MQ=255)
                                                                                                                                                                                                                               gtgggtggAACCGCGGGAGGGGCACTGGCGGGGCACCGGACGCTGACTGTCACAGGGGTGGCAGCGCAGAGGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCGTGGCACAGCAGTGCAGGCGCAAAGCCACGGCGGTTAAGAAAGAGAATGACCTGGTTATCAGCCTGTAAATGCTGGCGCATACGAGTGAt  <  1:469469/228‑1 (MQ=255)
                                                                                                                                                                                                                                |                                                                                                                                                                                                                                    
GTAATGCAACCAGCGTCACATCCGGGAAATGGTGACCTTTCGCCAGCATTTGTGTACCAATCAAAATCCGCGCGCCGCCGCGATGTACTTCTGCCAGTTGCTGTTCCAGCGCCCCTTTGCGGCTGGTGGTATCGCGGTCGATACGAGAAATGGGCACGCCGGGGAACAACGGCGCGAGCGTCTGTTCAAGCTGTTCGGTGCCCAGCCCCACGGGGACCAGGTGCGTGGAACCGCAGGAAGGGCACTGGCGCGGCACCGGACGCTGACTGTCACAGTGGTGGCAGCGCAGATGGTGCTGCGCCTGATGCAGCGTGTAGTAGTGATCGCAACGTGGGCATTCGGCAATCCAGCCACAGTCGTGGCACAGCAGTGCAGGCGCAAAGCCACGGCGGTTAAGAAAGAGAATGACCTGGTTATCAGCCTGTAAATGCTGGCGCATACGAGTGAT  >  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: