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

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