Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,320,007 G→A H219H (CAC→CAT phnI ← carbon‑phosphorus lyase core complex subunit PhnI

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,320,0070GA95.2% 159.3 / ‑4.0 66H219H (CAC→CATphnIcarbon‑phosphorus lyase core complex subunit PhnI
Reads supporting (aligned to +/- strand):  ref base G (1/2);  new base A (36/24);  total (38/28)
Fisher's exact test for biased strand distribution p-value = 5.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCCGGGGCTTGCAGAGCGCGGTCGACCAGCGCCATCGCCATCGCTTTGCGCTCGCTCATGCCGAACACCAGCCCGTAGCCGCGCGTGAAGTGCGGCGGCTCACCCGGCGGGTCGATAAAACCGTTAACCATTTCACACTCAGTCATCAGCAGTTCGCCGACGTTCACCGCAAATCCCAGCTCTTCCGGCACAATCGACACGTCGATATAGCCGCTGCGGATCTCGCCTGCGAACGGGTGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGTCGCCGCGCATCAACTGCTGCAAACGGGAGGAGCGTGAGCAGGGGTAAACCGGCGGCGTGCGGGTGATGTCATCCGGCTGTGCGCCGCTATCCTCTTCAAACTTCGCCAGCCCCTGACGCGCCAGCAGGCTGAAAACGTGCGGCGACGGCTGCTGTTCGCTGTCGGCGGTGGTCAGCGTCG  >  NC_000913/4319771‑4320245
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                               
ctccgggCTTGCAGAGAGCAGTCGCCCAGCGCCATAGCCATCGCTTTGCGCTCGCTCATGCCGCACACCAGCCCGTAGCCGCGCGTGAAGTGCGGCGGCTCACCCGGCGGGTCGATAAAACCGTTAACAATTTAACTCACAGTCATAATCAGCTCGCCGACGTTCACCGCAACACCCAGCCCTTCCAGCCCAATCGACACGCACTTATAGCCTCGGCGGAGCTCGCCTGCGCACGGGTGATTGCGCCCgct                                                                                                                                                                                                                                   >  2:239486/5‑249 (MQ=255)
               gcgcGGTCGACCAGCGCCATCGCCATCGCTTTGCGCTCGCTCATGCCGAACACCAGCCCGTAGCCGCGCGTGAAGTGCGGCGGCTCACCCGGCGGGTCGATAAAACCGTTAACCATTTCACACTCAGTCATCAGCAGTTCGCCGACGTTCACCGCAAATCCCAGCTCTTCCGGCACAATCGACACGTCGATATAGCCGCTGCGGATCTCGCCTGCGAACGGATGATTGCGCCCGTAACCGCGTTGGGTGGa                                                                                                                                                                                                                    <  1:46237/251‑1 (MQ=255)
                  cggccGACCAGCGCCATCGCCATCGCTTTGCGCTCGCTCATGCCGAACCCCAGCCCGTAGCCGCGCGTGAAGTGCGGCGGCTCACCCGGCGGGTCGATAAAACCGTTAACCATTTCACACCCAGTCATCAGCAGTTCGCCGACGTTCACCGCAAATCCCAGCTCTTCCGGCACAATCGACACGTCGATATAGCCGCTGCGGATCTCGCCTGCGAACGGATGATTGCGCCCGTAACCGCGTTGGGTGGAGTa                                                                                                                                                                                                                 <  1:71241/247‑1 (MQ=255)
                   ggTCGACCAGCGCCATCGCCATCGCTTTGCGCTCGCTCATGCCGAACACCAGCCCGTAGCCGCGCGTGAAGTGCGGCGGCTCACCCGGCGGGTCGATAAAACCGTTAACCATTTCACACTCAGTCATCAGCAGTTCGCCGACGTTCACCGCAAATCCCAGCTCTTCCGGCACAATCGACACGTCGATATAGCCGCTGCGGATCTCGCCTGCGAACGGATGATTGCGCCCGTAACCGCGTTGGGGGGAGt                                                                                                                                                                                                                  >  1:70886/1‑249 (MQ=255)
                          cAGCGCCATCGCCATCGCTTTGCGCTCGCTCATGCCGAACACCAGCCCGTAGCCGCGCGTGAAGTGCGGCGGCTCACCCGGCGGGTCGATAAAACCGTTAACCATTTCACACTCAGTCATCAGCAGTTCGCCGACGTTCACCGCAAATCCCAGCTCTTCCGGCACAATCGACACGTCGATATAGCCGCTGCGGATCTCGCCTGCGAACGGATGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGgccagcg                                                                                                                                                                                                         <  1:75957/251‑1 (MQ=255)
                          cAGCGCCATCGCCATCGCTTTGCGCTCGCTCATGCCGAACACCAGCCCGTAGCCGCGCGTGAAGTGCGGCGGCTCACCCGGCGGGTCGATAAAACCGTTAACCATTTCACACTCAGTCATCAGCAGTTCGCCGACGTTCACCGCAAATACCAGCTCTTCCGGCCCAATCGACACGTCGATATAGCCGCTGCGGATCTCGCCTGCGAACGGATGGTTGCGCCCGTAACCGCGTTGGGTGGAGTAGgccagcg                                                                                                                                                                                                         >  1:220776/1‑251 (MQ=255)
                           aGCGCCATCGCCATCGCTTTGCGCTCGCTCATGCCGAACACCAGCCCGTAGCCGCGCGTGAAGTGCGGCGGCTCACCCGGCGGGTCGATAAAACCGTTAACCATTTCACACTCAGTCATCAGCAGTTCGCCGACGTTCACCGCAAATCCCAGCTCTTCCGGCACAATCGCCACGTCGATATAGCCGCTGCGGATCTCGCCTGCGAACGGATGATTGCGCCCGTAAACGCGTTGGGTGGAGTAGgccagcgc                                                                                                                                                                                                        >  1:43237/1‑251 (MQ=255)
                           aGCGCCATCGCCATCGCTTTGCGCTCGCTCATGCCGAACACCAGCCCGTAGCCGCGCGTGAAGTGCGGCGGCTCACCCGGCGGGTCGATAAAACCGTTAACCATTTCACACTCAGTCATCAGCAGTTCGCCGACGTTCACCGCAAATCCCAGCTCTTCCGGCACAATCGACACGTCGATATAGCCGCTGCGGATCTCGCCTGCGAACGGATGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGgccagcgc                                                                                                                                                                                                        >  1:30223/1‑251 (MQ=255)
                             cgccatcgccatCGCTTTGCGCTCGCTCATGCCGAACACCAGCCCGTAGCCGCGCGTGAAGTGCGGCGGCTCACCCGGCGGGTCGATAAAACCGTTAACCATTTCTCACTCAGTCATCAGCAGTTCGCCGACGTTCACCGCAAATCCCAGCTCTTCCGGCACAATCGACACGTCGATATAGCCGCTGCGGATCTCGCCTGCGAACGGTTGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGgccagcgcca                                                                                                                                                                                                      <  1:150910/251‑1 (MQ=255)
                               ccatcgccatcgCTTTGCGCTCGCTCATGCCGAACACCAGCCCGTAGCCGCGCGTGAAGTGCGGCGGCTCACCCGGCGGGTCGATAAAACCGTTAACCATTTCACACTCAGTCATCAGCAGTTCGCCGACGTTCACCGCAAATCCCAGCTCTTCCGGCACAATCGACACGTCGATATAGCCGCTGCGGATCTCGCCTGCGAACGGATGATTGCGCCCGTAACCGCGTTGGGGGGAGTAGgccagcgccagc                                                                                                                                                                                                    >  1:160676/1‑251 (MQ=255)
                                 atcgccaccgtTTTTCGCTCGCTCATGCCGAACACCAGCCCGTAGCCGCGCGTGAAGTGCGGCGGCTCACCCGGCGGGTCGATAAAACCGTTAACCATTTCACACTCAGTCATCAGCAGTTCGCCGACGTTCACCGCAAATCCCAGCTCTGCCGGCACAATCGACTCGTCGATATAGCCGCTGCGGATCTCGCCTGCGAACGGATGATTGCTCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCaa                                                                                                                                                                                                  <  2:323045/251‑1 (MQ=255)
                                         cggcttGCGCTCTCTTCTGCCGAACCCCATCCCGTAGCCGCGCGTGAAGTGTGGCGGCTCACCCGGCGGGTCGATAAAACCGTTAACCATGTCACACCCAGTCAGCAGCAGTTCGCCGTCGTTCACCGCAAATCCCAGCTCTTCCGGCACAATCGACACGTCTATATAGCCGCTGCGGATCTCGCCTGCGAACGGATGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCt                                                                                                                                                                                           <  2:164987/246‑1 (MQ=255)
                                          gcTTTGCGCTCGCTCATGCCGAACACCAGCCCGTAGCCGCGCGTGAAGTGCGGCGGCTCACCCGGCGGGTCGATAAAACCGTTAACCATTTCACACTCAGTCATCAGCAGTTCGCCGACGTTCACCGCAAATCCCAGCTCTTCCGGCACAATCGACACGTCGATATAGCCGCTGCGGATCTCGCCTGCGAACGGATGATTACGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTtc                                                                                                                                                                                         >  1:35666/1‑251 (MQ=255)
                                                          tGCCGAACACCAGCCCGTAGCCGCGCGAGAAGTGCGGCGGCTCACCCGGCGGGTCGATAAAACCGTTAACCATTTCACACTCAGTCATCAGCAGTTCGCCGACGTTCACCGCAAATCCCAGCTCTTCCGGCACAATCGACACGACGATATAGCCGCTGCGGATCTCGCCTGCGAACGGATGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGTCGCCGCGCATCAAc                                                                                                                                                                         >  1:301215/1‑251 (MQ=255)
                                                           gCCGAACACCAGCCCGTAGCCGCGCGTGAAGTGCGGCGGCTCACCCGGCGGGTCGATAAAACCGTTAACCATTTCACACTCAGTCATCAGCAGTTCGCCGACGTTCACCGCAAATCCCAGCTCTTCCGGCACAATCGACACGTCGATATAGCCGCTGCGGATCTCGCCTGCGAACGGATGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGTCGCCGCGCATCAAc                                                                                                                                                                         <  1:54041/250‑1 (MQ=255)
                                                             ccaactccgtcccGCAGCCGTGCGTGAAGTGCGGCGGTTCCCCGGGCGTGTCGAGAAAACCGTTAACCATGTCACACTTAGTCATCAGCAGTTCGCCGACGTTCACCGCAAATCCCAGCGCGTCCGGCACAAGCGTCCCGTCGATAGAGCCGCTGCGGATCGCGCCGGCGAACGGATGATTGCGCCCGGAACCGCGTTGGGTGTAGTAGGCCAGCGCCAGCAAATAGCCTTCTTTGCCGCGCTTCAActgc                                                                                                                                                                      <  2:232472/241‑1 (MQ=255)
                                                              gAACACCAGCCCGTAGCCGCGCGTGAAGTGCGGCGGCTCACCCGGCGGGTCGATAAAACCGTTAACCATTTCACACTCAGTCATCAGCAGTTCGCCGACGTTCACCGCAAATCCCAGCTCTTCCGGCACAATCGACACGTCGATATAGCCGCTGCGGATCTCGCCTGCGAACGGATGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGTCGCCGCGCATCCActgct                                                                                                                                                                     >  1:106906/1‑251 (MQ=255)
                                                              gAACACCAGCCCGTAGCCGCGCGTGAAGTGCGGCGGCTCACCCGGCGGGTCGATAAAACCGTTAACCATTTCACACTCAGTCATCAGCAGTTCGCCGACGTTCACCGCAAATCCCAGCTCTTCCGGCACAATCGACACGTCGATATAGCCGCTGCGGATCTCGCCTGCGAACGGATGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGTCGCCGCGCATCAActgct                                                                                                                                                                     >  1:131588/1‑251 (MQ=255)
                                                               acagccctcccgtatccgctcattaaGTGCGGCTGCTCACCCGGCTGGTCGATACAACCTTTAACCATCTCACCCTCACTCAGCAGCAGTTCGCCGACGTTCACCGCAAATCCCAGCTCATCCGGCACAATCGTCACGTCTATATAGCCGCTGCGGATCTCGCCTGCGAACGGATGATTGCGCCCGTAACCGCGTTGGGTGGAGGAGGCCAGCGCCCGCAAATAGCCTTCGTCGCCGCGCATCAActgc                                                                                                                                                                      <  2:382515/225‑1 (MQ=255)
                                                                accccAGCCCGTTGCCGCGCGTGAAGTGCGGCGGCTCACCCGGCGGGTCGATAAAACCGTTAACCATTTCACACTCAGTCATCAGCAGTTCGCCGACGTTCACCGCAAATCCCAGCTCTTCCGGCACAATCGACACGTCGATATAGCCGCTGCGGATCTCGCCTGCGAACGGATGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGTCGCCGCGCATCAACtgctgc                                                                                                                                                                   <  1:360300/248‑1 (MQ=255)
                                                                      gCCCGTAGCCGCGCGTGAAGTGCGGCGGCTCACCCGGCGGGTCGATAAAACCGTTAACCATTTCACACTCAGTCATCAGCAGTTCGCCGACGTTCACCGCAAATCCCAGCTCTTCCGGCACAATCGACACGTCGATATAGCCGCTGCGGATCTCGCCTGCGAACGGATGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGTCGCCGCGCATCAACTGCTGCAAACgg                                                                                                                                                             >  1:143917/1‑251 (MQ=255)
                                                                       cccGTAGCCGCGCGTGAAGTGCGGCGGCTCACCCGGCGGGTCGATAAAACCGTTAACCATTTCACACTCAGTCATCAGCAGTTCGCCGACGTTCACCGCAAATCCCAGCTCTTCCGGCACAATCGACACGTCGATATAGCCGCTGAGGATCTCGCCTGCGAACGGATGATTGCGCCCGTAACCGCGTTGGGGGGGGTAGGCCCGCGCCAGCAAATAGCCTTCGTCGCCGCGCATCACCTGCTGCAAAAggg                                                                                                                                                            >  1:87616/1‑251 (MQ=255)
                                                                            aGCCGCGCGTGAAGTGCGGCGGCTCACTCGGCGGGTCGATAAAACCGTTAACCATTTCACACTCAGTCATCAGCAGTTCGCCGACGTTCACCGCAAATCCCAGCTCTTCCGGCACAATCGACACGTCGATATAGCCGCTGCGGATCTCGCCTGCGAACGGATGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGTCGCCGCGCATCAACTGCTGCAAACGGgaggag                                                                                                                                                       <  1:209576/251‑1 (MQ=255)
                                                                                gcgcgTGAAGTGCGGCGGCTCACCCGGCGGGTCGATAAAACCGATAACCATTTCACACTCAGTCATCAGCAGTTCGCCGACGTTCACCGCAAATCCCAGCTCTTCCGGCACAATCGACACGTCGATATAGCCGCTGCGGATCTCGCCTGCGAACGGATGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGTCGCCGCGCATCAACTGCTGCAAACGGGAGGAGCGTg                                                                                                                                                   >  1:132160/1‑251 (MQ=255)
                                                                                           gcggcgGCTCACCCGGCGGGTCGATAAAACCGTTAACCATTTCACACTCAGTCATCAGCAGTTCGCCGACGTTCACCGCAAATCCCAGCTCTTCCGGCACAATCGACACGTCGATATAGCCGCTGCGGATCTCGCCTGCGAACGGATGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGTCGCCGCGCATCAACTGCTGCAAACGGGAGGAGCGTGAGCAGGGGTa                                                                                                                                         <  1:363349/250‑1 (MQ=255)
                                                                                               cggcTCACCCGGCGGGTCGATAAAACCGTTAACCATTTCACACTCAGTCATCAGCAGTTCGCCGACGTTCACCGCAAATCCCAGCTCTTCCGGCACAATCGACACGTCGATATAGCCGCTGCGGATCTCGCCTGCGAACGGATGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATGGCCTTCGTCGCCGCGCATCAACTGCTGCAAACGGGAGGAGCGTGAGCAGGGGTAAACcg                                                                                                                                    >  1:172277/1‑251 (MQ=255)
                                                                                                 gcTCACCCGGCGGGTCGATAAAACCGTTAACCATTTCACACTCAGTCATCAGCAGTTCGCCGACGTTCACCGCAAATCCCAGCTCTTCCGGCACAATCGACACGTCGATATAGCCGCTGCGGATCTCGCCTGCGAACGGATGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGTCGCCGCGCATCAACTGCTGCAAACGGGAGGAGCGTGAGCAGGGGTAAACcggc                                                                                                                                 >  1:216847/1‑251 (MQ=255)
                                                                                                     aCCCGGCGAGTCGATAAAACCGTTAACCATTTCACACTCAGTCATCAGCAGTTCGCGGACGATCACCGCAAATCCCAGCTCTTCCGGCACACTCGACACGTCGATATAGCCGCTGCGGATATCGCCTGCGCAAGGCTGATTGCGAACGTAACCGCGTTGGGTGGAGTAGGCCCACGCAAGCAAATAGCCTTCATCGCCGCGCACCAACTTGTGCAAACGGGAGGAGCGTGAGCAGGGGTAAACCGgcggcg                                                                                                                             >  2:307309/1‑251 (MQ=255)
                                                                                                       ccGGCGGGTCGATAAAACCGTTAACCATTTCACCCTCAGTCATCAGCAGTTCGCCGACGTTCACCGCAAATCCCAGCTCTTCCGGCACAATCGACACGTCGATATAGCCGCTGCGGATCTCGCCTGCGAACGGATGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGTCGCCGCGCATCAACTGCTGCAAACGGGAGGAGCGTGAGCAGGGGTAAACCGGCGGCGTg                                                                                                                           >  1:107296/1‑251 (MQ=255)
                                                                                                         ggcggGTCGATAAAACCGTTAACCATTTCACACTCAGTCATCAGCAGTTCGCCGACGTTCACCGCAAATCCCAGCTCTTCCGGCACAATCGACACGTCGATATAGCCGCTGCGGATCTCGCCTGCGAACGGATGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGTCGCCGCGCATCAACTGCTGCAAACGGGAGGAGCGTGAGCAGGGGTAAACCGGCGGCGTGCg                                                                                                                         <  1:176911/251‑1 (MQ=255)
                                                                                                            gggTCGATAAAACCGTTAACCATTTCACACTCAGTCATCAGCAGTTCGCCGACGTTAACCGCAAATCCCAGCTCTTCCGGCACAATCGACACGTCGATATAGCCGCTGCGGATCTCGCCTGCGAACGGATGATTGCGCCCGTAACCGCGGTGGGTGGAGTAGGCCAGCGCCAGCACATAGCCTTCGTCGCCGCGCATCAACTGCTGCAAACGGGAGGAGCGTGAGCAGGGGTAACCCGGCGGCGTGCGGGt                                                                                                                      >  1:366063/1‑251 (MQ=255)
                                                                                                              gtctataaagacgtcaaCCATCTCACGCTCAGTCACCAGCAGTTCGCCGGCGTTCACCGTAAAACCCAGCTCTTCCGGCACAATCGACGCGGCGATTTAGCCGCTGCGGAGCTCGCCTGCGAACGGATGATTGCGCCCGTAACCGCGTTGGGTGGAGGAGGCCAGCGTCAGCAAATAGCCTTCGTCGCCGCGCATCAACTGCTGCAAACGGGAGGAGCGTGAGCAGGGGTAAATCGGCGGCGTGCGGGTGa                                                                                                                    <  2:369529/236‑1 (MQ=255)
                                                                                                                       aCCGTTAACCATTTCACACTCAGTCAGCAGCAGTTCGCCGACGTTCACCGCAAATCCCAGCTCTCCCGGCACAATCGACACGTCGATATAGCCGCTGCGGATCTCGCCTGCGAACGGATGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGTCGCCGCGCATCAACTGCTGCAAACGGGAGGAGCGTGAGCAGGGGTAAACCGGCGGCGTGCGGGTGATGTCATCCg                                                                                                           <  2:84654/251‑1 (MQ=255)
                                                                                                                        ccGTTAACCATTTCACACTCAGTCATCAGCAGTTCGCCGACGTTCACCGCAAATCCCAGCTCTTCCGGCACAATCGACACGTCGATATAGCCGCTGCGGATCTCGCCTGCGAACGGATGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGTCGCCGCGCATCAACTGCTGCAAACGGGAGGAGCGTGAGCAGGGGTAAACCGGCGGCGTGCGGGTGATGTCATCCgg                                                                                                          >  1:252235/1‑251 (MQ=255)
                                                                                                                         cGTTAACCATTTCACACTCAGTCATCAGCAGTTCGCCGACGTTCACCGCAAATCCCAGCTCTTCCGGCACAATCGACACGTCGATATAGCCGCTGCGGATCTCGCCTGCGAACGGATGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGTCGCCGCGCATCAACTGCTGCAAACGGGAGGAGCGTGAGCAGGGGAAAACCGGCGGCGTGCGGGTGAGGTCATCCGGc                                                                                                         >  2:120644/1‑251 (MQ=255)
                                                                                                                          gtcaactaTTTCACACACAGTCAGCAGCAGCTCGCCGCCTTTCACCGCAAAGTCCAGCTCTTCCGGCGCAATTGACACGTCGATATAGCCGCTGCGGTGCTCGCCCGTGAAAGGCTGATTGCGGCCGTAACCGCGTGGGGTGGAGTAGGCCAGCGCCAGCAAATAGCTGTCGTCGCCGCGCATCAACTGCTGCTAACGGGAGGAGGGTGAGCAGTGGTAAACCGGCGGCGTGCGGGTGATGTCATCCGGca                                                                                                        <  2:82272/244‑2 (MQ=255)
                                                                                                                             aaCCATTTCACACTCAGTCATCAGCAGTTCGCCGACGTTCACCGCAAATCCCCGCTCTTCCGGCACAATCGACACGTCGATATAGCCGCTGCGGATCTCGCCTGCGAACGGATGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGACAGCTACTAGCCTTCGTCGCCGCTCATCAATCGCTGCAATCGGGAGGAACGATAGCGGGGGTAAACCGTGCGCGTGCgggacgtgtacacctgctgtg                                                                                                     >  2:199623/1‑233 (MQ=255)
                                                                                                                               ccATTTCACACTCAGTCATCAGCAGTTCGCCGACGTTCACCGCAAATCCCAGCTCTTCCGGCACAATCGACACGTCGATATAGCCGCTGCGGATCTCGCCTGCGAACGGATGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGTCGCCGCGCATCAACTGCTGCAAACGGGAGGAGCGTGAGCAGGGGTAAACCGGCGGCGTGCGGGTGATGTCATCCGGCTGTgcg                                                                                                   >  1:222664/1‑251 (MQ=255)
                                                                                                                                 aTTTCACACTCAGTCATCAGCAGTTCGCCGACGTTCACCGCAAATCCCAGCTCTTCCGGCACAATCGACACGTCGATATAGCCGCTGCGGATCTCGCCTGCGAACGGATGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGTCGCCGCGCATCAACTGCTGCAAACGGGAGGAGCGTGAGCAGGGGTAAACCGGCGGCGTGCGGGTGATGTCATCCGGCTGTGcgcc                                                                                                 >  1:64092/1‑251 (MQ=255)
                                                                                                                                    tCACACTCAGTCATCAGCAGTTCGCCGACGTTCACCGCAAATCCCAGCTCTTCCGGCACAATCGACACGTCGATATAGCCGCTGCGGATCTCGCCTGCGAACGGATGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGTCGCCGCGCATCAACTGCTGCAAACGGGAGGAGCGTGAGCAGGGGTAAACCGGCGGCGTGCGGGTGATGTCATCCGGCTGTGCGCCGCt                                                                                              >  1:3206/1‑251 (MQ=255)
                                                                                                                                       caccactgcatctgctgTTCGCCGCCGTTCACCGCAAAGCCCATATAGTCCGGCTCAAACGACTCGTCGAGATAGCCGCGGCGGTTCACGCCTGCGAACGGATGATGTCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCTTCGCCGCGCATCAACTTCTGCAAACGTGAGGAGCGGGAGCAGTGGTAAACCGGCGGCGGGCGGGTGATGTCAGCCGGCTGTGCACCGCAATc                                                                                           <  2:220776/235‑1 (MQ=255)
                                                                                                                                               catcaGCAGTTCGCCGACGTTCACCGCAAATCCCAGCTCTTCCGGCACAATCGACACGTCGATATAGCCGCTGCGGATCTCGCCTGCGAACGGATGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGTCGCCGCGCATCAACTGCTGCAAACGGGAGGAGCGTGAGCAGGGGTAAACCGGCGGCGTGCGGGTGATGTCATCCGGCTGTGCGCCGCTATCCTCTTCaa                                                                                   <  1:120707/251‑1 (MQ=255)
                                                                                                                                                atcaGCAGTTCGCCGAAGTTCACCGCAAATCCCAGCTCTTCCTGCACCATCGACACGTCGCTATAGCCGCTGCGGATCTCGCCTGCGAACGGATGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGTCGCCGCGCATCACCTGCTGCAAACGGGAGGAGCGTGAGCAGGGGTAAAACGGCGCAGGGGGAGTGATGTAATCCGGCTGTGCGCCGCTATCCTATTCaac                                                                                  >  2:18010/1‑250 (MQ=255)
                                                                                                                                                      tgTTCGCCGACGTTCACCGCAAATCCCAGCTCTTCCGACACAATCGACACGTCGATATAGCCGCTGCGGATCTCGCCTGCGAACGGATGATTGAGCCAGTAACCGCGTTGGGGGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGTCGCCGCGCATCAGCTGCTGCAAACGGGAGGAGCGTGAGCAGGGGTAAACCAGCGGCGCGCGGGTGATGTCATCCGGCTGTGCGCCGCTATCCTCttacaacttggc                                                                            >  2:4460/2‑241 (MQ=255)
                                                                                                                                                             cGACGTTCACCGCAAATCCCAGCTCTTCCGGCACAATCGACACGTCGATATAGCCGCTGCGGATCTCGCCTGCGAACGGATGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGTCGCCGCGCATCAACTGCTGCAAACGGGAGGAGCGTGAGCAGGGGTAAACCGGCGGCGTGCGGGTGATGTCATCCGGCTGTGCGCCGCTATCCTCTTCAAACTTCGCCAGcccc                                                                     >  1:147217/1‑251 (MQ=255)
                                                                                                                                                               actttCACCGCAAATCCCAGCTCTTCCGGCCCAATCGACACGTCGATATAGCCGCTGCGGATCTCGCCTGCGAACGGATGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGTCGCCGCGCATCAACTGCTGCAAACGGGAGGAGCGTGAGCAGGGGTAAACCGGCGGCGTGCGGGTGATGTCTGCCGGCTGTGCGCCGCTATCCTCGTCAAACTTCGCCAGCCCCTg                                                                   <  2:235139/248‑1 (MQ=255)
                                                                                                                                                               aCGTTCACCGCAAATCCCAGCTCTTCCGGCACAATCGACACGTCGATATAGCCGCTGCGGATCTCGCCTGCGAACGGATGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGTCGCCGCGCATCAACTGCTGCAAACGGGAGGAGCGTGAGCAGGGGTAAACCGGCGGCGTGCGGGTGATGTCATCCGGCTGTGCGCCGCTATCCTCTTCAAACTTCGCCAGCCCCTg                                                                   >  1:130758/1‑251 (MQ=255)
                                                                                                                                                                cGTTCACCGCAAATCCCAGCTCTTCCGGCACAATCGACACGTCGATATAGCCGCTGCGGATCTCGCCTGCGAACGGATGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGTCGCCGCGCATCAACTGCTGCAAACGGGAGGAGCGTGAGCAGGGGTAAACCGGCGGCGTGCGGGTGATGTCATCCGGCTGTGCGCCGCTATCCTCTTCAAACTTCGCCAGCCCCTg                                                                   >  1:173611/1‑250 (MQ=255)
                                                                                                                                                                    cACCGCAAATCCCAGCTCTTCCGGCACAATCGACACGTCGATATAGCCGCTGCGGATCTCGCCTGCGAACGGATGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGTCGCCGCGCATCAACTGCTGCAAACGGGAGGAGCGTGAGCAGGGGTAAACCGGCGGCGTGCGGGTGATGTCATCCGGCTGTGCGCCGCTATCCTCTTCAAACTTCGCCAGCCCCTGAcgcg                                                              >  1:234701/1‑251 (MQ=255)
                                                                                                                                                                    cACCGCAAATCCCAGCTCTTCCGGCACAATCGACACGTCGATATAGCCGCTGCGGATCTCGCCTGCGAACGGATGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGTCGCCGCGCATCAACTGCTGCAAACGGGAGGAGCGTGAGCAGGGGTAAACCGGCGGCGTGCGGGTGATGTCATCCGGCTGTGCGCCGCTATCCTCTTCAAACTTCGCCAGCCCCTGAcgc                                                               <  1:283464/250‑1 (MQ=255)
                                                                                                                                                                               ccAGCTATTCCGGCACAATCGCCTCGTGGATATAGCCGCTGCGGATCACGCCTGCGACCGGAGGATGGCGCCCGTAGCGGCGTTGGGTGGAGTAGGCCAGCGCCAGAAAATATCCTTCGTCGCCGCGCATCAACTGCTGCATACGGGAGGAGCGTGAGCAGGGGTAAACCTGCGGCGTGCGGGTGATGTCATCCGGCTGGGCGCCGCTATCCTCTGCAAACTTCGCCAGCCCCTGACGCGCCAGCAGGCtt                                                   <  2:316578/251‑2 (MQ=255)
                                                                                                                                                                                            cacaATCGACGCGACGATATAGCCGCTGCGGATCTCGCCGGCGAACGGATGAGTGCGCCCGTAACCGCGTTGGGTGGCGTAGGCAAGCGCCAGCAAAGAGCATACGTCGCCGCGCACCAACTGCTGCAAACGGGAGGAGCTTGATCAGGGGTAAACCGGCGGCGGGCGGGTGATGTCATCCGGCTGTGCGCCGCTATCCTCTTCAAACTTCGCCAGCCCCTGACGCGCCAGCAGGCTGAAAACGTgcggcg                                      <  2:159716/251‑1 (MQ=255)
                                                                                                                                                                                                aTCGACACGTCGATATAGCCGCTGCGGATCTCGCCTGCGAACGGATGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGTCGCCGCGCATCAACTGCTGCAAACGGGAGGAGCGTGAGCAGGGGTAAACCGGCGGCGTGCGGGTGATGTCATCCGGCTGTGCGCCGCTATCCTCTTCAAACTTCGCCAGCCCCTGACGCGCCAGCAGGCTGAAAACGTGCGGCGACgg                                  >  1:97368/1‑251 (MQ=255)
                                                                                                                                                                                                aTCGACACGTCGATATAGCCGCTGCGCATATCGCCTGCGAACGGATGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGTCGACGCGCATCAACTGCTGCAAACGGGAGGAGCGTGAGCAGGGGTAAACCGGCGGCGTGCGGGTGATGTCATCCGGCTGTGCGCAGCTATCCTCTTCAAACTTCGCCAGCCCCTGACGCGCCAGCAGGCTGAAAACATGCGGCGACgg                                  >  1:317348/1‑251 (MQ=255)
                                                                                                                                                                                                   gacacgcagaaataGCCGCTGCGGATCTCGCCTGCGCACGGATGATTGCGCCCGTACCCGCGATGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGCCGACGCGCATCAGCTGCTGCAAACGGAAGGAGCGTGCGCCGGGGTAAACCGGCGGCGACCGGGGGATGTCATCCGGCTGAGCGCCGCCATACTCTTCAAACTTCGCCCGCCCCTGTCGCTCCAGCAGGCGGCAAACGTACGACGACGgctg                               >  2:39385/12‑251 (MQ=255)
                                                                                                                                                                                                    acacGTCGATATAGCCGCTGCGGATCTCGCCTGCGAACGGATGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGTCGCCGCGCATCAACTGCTGCAAACGGGAGGAGCGTGAGCAGGGGTAAACCGGCGGCGTGCGGGTGATGTCATCCGGCTGTGCGCCGCTATCCTCTTCAAACTTCGCCAGCCCCTGACGCGCCAGCAGGCTGAAAACGTGCGGCGACGgctgc                              >  1:314785/1‑251 (MQ=255)
                                                                                                                                                                                                     cacGTCGATATAGCCGCTGCGGATCTCGCCTGCGAACGGATGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGTCGCCGCGCATCAACTGCTGCAAACGGGAGGAGCGTGAGCAGGGGTAAACCGGCGGCGTGCGGGTGATGTCATCCGGCTGTGCGCCGCTATCCTCTTCAAACTTCGCCAGCCCCTGACGCGCCAGCAGGCTGAAAACGTGCGGCGACGgctgct                             <  1:248834/251‑1 (MQ=255)
                                                                                                                                                                                                      acGTCGATATAGCCGCTGCGGATCTCGCCTGCGAACGGATGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGTCGCCGCGCATCAACTGCTGCAAACGGGAGGAGCGTGAGCAGGGGTAAACCGGCGGCGTGCGGGTGATGTCATCCGGCTGTGCGCCGCTATCCTCTTCAAACTTCGCCATCCCCTGACGCCCCAGCAGGCTGAAAACGTGCGGCGACGGctgctg                            >  2:277472/1‑251 (MQ=255)
                                                                                                                                                                                                           gATATAGCCGCTGCGGATCTCGCCTGCGAACGGATGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGTCGCCGCGCATCAACTGCTGCAAACGGGAGGAGCGTGAGCAGGGGTAAACCGGCGGCGTGCGGGTGATGTCATCCGGCTGTGCGCCGCTATCCTCTTCAAACTTCGCCAGCCCCTGACGCGCCAGCAGGCTGAAAACGTGCGGCGACGGCTGCTGTTCGc                       >  1:256070/1‑251 (MQ=255)
                                                                                                                                                                                                               agctccctgcgtactcctccggCGAACGGAGGATTGCCCCCGTAACCGCTTGGGGGGGAGTAGGCCAGTGCCACCAAATAGCCTTCGTGCCCGCGCATCAACTGCTGCAACCGGGAGGAGCGTGATCAGTGGTAAACCGGCGGCGTGCGGGTGATTTCATCCGGCTGTGCTCCGCTATCCTCTTCAAACTTCGCCAGCCCCTGACTCGCCAGCAGGCTGAAAACGTGCGGCGACGGCTGCTGTTCGCTGTc                   <  2:216847/230‑1 (MQ=255)
                                                                                                                                                                                                                aGCCGCTGCGGAGCTCGCCTGCGAACGGATGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGTCGCCGCGCATCAACTGCTGCAAACGGGAGGAGCGTGAGCAGGGGTAAACCGGCGGCGTGCGGGTGATGTCATCCGGCTGTGCGCCGCTATCCTCTTCAAACTTCGCCAGCCCCTGACGCGCCAGCAGGCTGAAAACGTGCGGCGACGGCTGCTGTTCGCTGTc                   <  1:329908/250‑1 (MQ=255)
                                                                                                                                                                                                                   cGCTGCGGATCTCGCCTGCGAACGGATGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGTCGCCGCGCATCAACTGCTGCAAACGGGAGGAGCGTGAGCAGGGGTAAACCGGCGGCGTGCGGGTGATGTCATCCGGCTGTGCGCCGCTATCCTCTTCAAACTTCGCCAGCCCCTGACGCGCCAGCAGGCTGAAAACGTGCGGCGACGGCTGCTGTTCGCTGTcggcg               >  1:106072/1‑251 (MQ=255)
                                                                                                                                                                                                                    gcggcggctctggtcTGCGAACGGGTGATTGCCCCCTTAGCCGCGGTGGGTGGAGCAGGCGAGCGCCAGCAAATAGCCTGCGCCGCCGCGCATCAACTGCTGCAAACGGGAGGAGCGTTAGCATGGTTAAACCGGCGGCGTGCGGGGGAGGTGAGCCGGCTGTGCGCCGCTATCTTCTTCAACCTTCGCCAGCCCCTTACGCGCCAGCAGGCTGAAAACGTGCGTCGACGGCTGCTGTTCGCTGTcggcgg              <  2:347653/237‑1 (MQ=255)
                                                                                                                                                                                                                         ggATCTCGCCTGCGAACGGATGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGTCGCCGCGCATCAACTGCTGCAAACGGGAGGAGCGTGAGCAGGGGTAAACCGGCGGCGTGCGGGTGATGTCATCCGGCTGTGCGCCGCTATCCTCTTCAAACTTCGCCAGCCCCTGACGCGCCAGCAGGCTGAAAACGTGCGGCGACGGCTGCTGTTCGCTGTCGGCGGTGGTc         >  1:68225/1‑251 (MQ=255)
                                                                                                                                                                                                                               cGCCTGCGAACGGATGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGTCGCCGCGCATCAACTGCTGCAAACGGGAGGAGCGTGAGCAGGGGTAAACCGGCGGCGTGCGGGTGATGTCATCCGGCTGTGCGCCGCTATCCTCTTCAAACTTCGCCAGCCCCTGACGCGCCAGCAGGCTGAAAACGTGCGGCGACGGCTGCTGTTCGCTGTCGGCGGTTGTCa        <  1:317683/246‑1 (MQ=255)
                                                                                                                                                                                                                                gCCTGCGAACGGGTGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGTCGCCGCGCATCAACTGCTGCAAACGGGAGGAGCGTGAGCAGGGGTAAACCGGCGGCGTGCGGGTGATGTCATCCGGCTGTGCGCCGCTATCCTCTTCAAACTTCGCCAGCCCCTGACGCGCCAGCAGGCTGAAAACGTGCGGCGACGGCTGCTGTTCGCTGTCGGCGGTGGTCAGCGTCg  <  1:239486/251‑1 (MQ=255)
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                               
TCCGGGGCTTGCAGAGCGCGGTCGACCAGCGCCATCGCCATCGCTTTGCGCTCGCTCATGCCGAACACCAGCCCGTAGCCGCGCGTGAAGTGCGGCGGCTCACCCGGCGGGTCGATAAAACCGTTAACCATTTCACACTCAGTCATCAGCAGTTCGCCGACGTTCACCGCAAATCCCAGCTCTTCCGGCACAATCGACACGTCGATATAGCCGCTGCGGATCTCGCCTGCGAACGGGTGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGTCGCCGCGCATCAACTGCTGCAAACGGGAGGAGCGTGAGCAGGGGTAAACCGGCGGCGTGCGGGTGATGTCATCCGGCTGTGCGCCGCTATCCTCTTCAAACTTCGCCAGCCCCTGACGCGCCAGCAGGCTGAAAACGTGCGGCGACGGCTGCTGTTCGCTGTCGGCGGTGGTCAGCGTCG  >  NC_000913/4319771‑4320245

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