Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 647,210 A→C 22.2% G100G (GGT→GGG citG ← 2‑(5''‑triphosphoribosyl)‑3'‑dephosphocoenzyme‑A synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913647,2100AC22.2% 136.5 / 7.2 75G100G (GGT→GGGcitG2‑(5''‑triphosphoribosyl)‑3'‑dephosphocoenzyme‑A synthase
Reads supporting (aligned to +/- strand):  ref base A (42/12);  new base C (2/18);  total (44/31)
Fisher's exact test for biased strand distribution p-value = 1.43e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCGCCGGTAAGGCCAAGCTGTTGGTACAACCGTTGACCTGCCGTCAGTTGTGAATTATTGGTACGCAGTTCGCGATCGGTCAGGCCACGGCAGAAACTTGCCGCCGTAGAACAAACGGTTGTTGGCGTTACCGGTTGGTTGAGTTGAAGCAAACGGCCAATTGCCGCACATAGCAGCCCTAAAGAAAAAATGCTGCCTTTATGCGTGTTTACGCCCGCAGTGGCGCGGAACATATCACCTTCGCAAGCCATACCAATTGGGCGTAATCCGTGGAGTACCGCTTCTGGTGCCATTTCCGCACTACAGGCACCAAATTCAATGAAACGGGGTAGCCAGCCCTGAATCGCCAGCGCGCTGCGGTGGAAATCTTCCAGCGCCATATCTTTGTGCGCACCGCAGTTAATGCGATCCACGAGGCCTGGTTTCGGTGACAGATTGACTTCAGTCAGCATGGCGCGCCAGCCCAGCAGGG  >  NC_000913/646973‑647445
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                            
tGCGCCGGTAAGGCCAAGCTGTTGGTACAACCGTTGACCTGCCGTCAGTTGTGAATTATTGGTACGCAGTTCGCGATCGGTCAGGCCACGGCAGAAACTTGCCGCCGTAGAACAAACGGTTGTTGGCGTTACCGGTTGGTTGAGTTGAAGCAAACGGCCAATTGCCGCACATAGCAGCCCTAAAGAAAAAATGCTGCCTTTATGCGTGTTTACGCCCGCAGTGGCGCGGAACATATCa                                                                                                                                                                                                                                              >  2:842482/1‑238 (MQ=255)
tGCGCCGGTAAGGCCAAGCTGTTGGTACAACCGTTGACCTGCCGTCAGTTGTGAATTATTGGTACGCAGTTCGCGATCGGTCAGGCCACGGCAGAAACTTGCCGCCGTAGAACAAACGGTTGTTGGCGTTACCGGTTGGTTGAGTTGAAGCAAACGGCAAATTGCACCCCATAGCAGCTCCAAAGAAAAAATGCTGCCCTTATACGGTTTTACCCCCCCAGGGGGCAGGAACAAATCa                                                                                                                                                                                                                                              >  1:497746/1‑238 (MQ=255)
                               cGTTGACCTGCCGTCAGTTGTGAATTATTGGTACGCAGTTCGCGATCGGTCAGGCCACGGCAGAAACTTGCCGCCGTAGAACAAACGGTTGTTGGCGTTACCGGTTGGTTGAGTTGAAGCAAACGGCCAATTGCCGCACATAGCAGCCCTAAAGAAAAAATGCTGCCTTTATGCGTGTTTACGCCCGCAGTGGCGCGGAACATATCACCTTCGCAAGCCATACCAATTGGGCGTAATcc                                                                                                                                                                                                             >  2:716680/1‑239 (MQ=255)
                               cGTTGACCTGCCGTCAGTTGTGAATTATTGGTACGCAGTTCGCGATCGGTCAGGCCACGGCAGAAACTTGCCGCCGTAGAACAAACGGTTGTTGGCGTTACCGGTTGGTTGAGTTGAAGCAAACGGCCAATTGCCGCACATAGCAGCCCTAAAGAAAAAATGCTGCCTTTATGCGGGTTTACGCCCGCAGTGCGGCGGAAAAAATCAACTTCCCAAACCATACAAAATgggggaaaacg                                                                                                                                                                                                             >  1:340413/1‑231 (MQ=255)
                                 ttGACCTGCCGTCAGTTGTGAATTATTGGTACGCAGTTCGCGATCGGTCAGGCCACGGCAGAAACTTGCCGCCGTAGCACGAACGGTTGTTGGCGTTACCGGTTGGTTGAGTTGAAGCAAACGGCCAATTGCCGCACATAGCAGCCCTAAAGAAAAAATGCTGCCTTTATGCGTGTTTACGCCCGCAGTGGCGCGGAGCATATCACCTTCGAATGCCCTACCAATTGGGCGTAATCCgg                                                                                                                                                                                                           >  2:511785/1‑238 (MQ=255)
                                 ttGACCTGCCGTCAGTTGTGAATTATTGGTACGCAGTTCGCGATCGGTCAGGCCACGGCAGAAACTTGCCGCCGTAGAACAAACGGTTGTTGGCGTTACCGGTTGGTTGAGTTGAAGCAAACGGCCAATTGCCGCACATAGCAGCCCTAAAGAAAAAATGCTGCCTTTATGCGTGTTTACGCCCGCAGTGGCGCGGAACATATCACCTTCGCAAGCCATACCAATTGGGCGTAATCCGt                                                                                                                                                                                                           >  2:762749/1‑239 (MQ=255)
                                 ttGACCTGCCGTCAGTTGTGAATTATTGGTACGCAGTTCGCGATCGGTCAGGCCACGGCAGAAACTTGCCGCCGTAGAACAAACGGTTGTTGGCGTTACCGGTTGGTTGAGTTGAAGCAAACGGCCAATTGCCGCACATAGCAGCCCTAAAGAAAAAATGCTGCCTTTATGCGGGTTTACACCCGCAGTGGCGCGGAACATATCACCTTCGCAAGCCATACCACTTGGGCGTAATCCGt                                                                                                                                                                                                           >  2:441269/1‑239 (MQ=255)
                                   gACCTGCCGTCAGTTGTGAATTATTGGTACGCAGTTCGCGATCGGTCAGGCCACGGCAGAAACTTGCCGCCGTAGAACAAACGGTTGTTGGCGTTACCGGTTGGTTGAGTTGAACCAAACGGCCAATTGCCGCACATAGCAGCCCTAAAGAAAAAATGCTGCCATTCTCCGGGTTTACGCCCGCATAGGCGCGGAACATATCACCTTCGCAAACCCTACCAATTGGGGGGAAACCgggg                                                                                                                                                                                                         >  1:663062/1‑239 (MQ=255)
                                        gCCGTCAGTTGTGAATTATTGGTACGCAGTTCGCGATCGGTCAGGCCACGGCAGAAACTTGCCGCCGTAGAACAAACGGTTGTTGGCGTTACCGGTTGGTTGAGTTGAAGCAAACGGCCAATTGCCGCCCATAGCAGCCCTAAAGAAAAAATGCTGCCTTTATGCGTGTTTACGCCCGAAGGGGCGCGGAACATATACCCTTCGCAAGCCCAACCAATGTGGCGGAACCCGTggggac                                                                                                                                                                                                     >  1:425561/1‑234 (MQ=255)
                                        gCCGTCAGTTGTGAATTATTGGTACGCAGTTCGCGATCGGTCAGGCCACGGCAGAAACTTGCCGCCGTAGAACAAACGGTTGTTGGCGTTACCGGTTGGTTGAGTTGAAGCAAACGGCCAATTGCCGCACATAGCAGCCCTAAAGAAAAAATGCTGCCTTTATGCGTGTTTACGCCCGCAGTGGCGCGGACCATATCACCTTCGCAAGCCGTACCAATTGGGGGGAATCCGGGGGGTa                                                                                                                                                                                                     >  2:760836/1‑238 (MQ=255)
                                        gCCGTCAGTTGTGAATTATTGGTACGCAGTTCGCGATCGGTCAGGCCACGGCAGAAACTTGCCGCCGTAGAACAAACGGTTGTTGGCGTTACCGGTTGGTTGAGTTGAAGCAAACGGCCAATTGCCGCACATAGCAGCCCTAAAGAAAAAATGCTGCCTTTATGCGTGTTTACGCCCGCAGTGGCGCGGAAAATATCACCTTCGCAACCCATCCCAATTTGGCGTAATCCGTGGAGTa                                                                                                                                                                                                     >  2:385248/1‑238 (MQ=255)
                                        gCCGTCAGTTGTGAATTATTGGTACGCAGTTCGCGATCGGTCAGGCCACGGCAGAAACTTGCCGCCGTAGAACAAACGGTTGTTGGCGTTACCGGTTGGTTGAGTTGAAGCAAACGGCCAATTGCCGCACATAGCAGCCCTAAAGAAAAAATGCAGCCTTTATGCGTGTTTACGCCCGCCGTGGCGCGGAACATATCACCTTCGCAAGCCATACCAATTGGGCGTAATCCGTGGAGTa                                                                                                                                                                                                     >  2:439761/1‑238 (MQ=255)
                                        gCCGTCAGTTGTGAATTATTGGTACGCAGTTCGCGATCGGTCAGGCCACGGCAGAAACTTGCCGCCGTAGAACAAACGGGTGTTGGCGTTACCGGTTGGCTGAGTTGAAGCAAACGGCCACTTGCCGCACATAGCAGCCCTAAAGCAAAAATGCTGCCTGTTTGCGTGTTTACGCCCCACGTGCCCCGAAACATACCACCTTGCCAAGccaaccaaaattggcgaaaaacaggggagac                                                                                                                                                                                                   >  1:71211/1‑216 (MQ=255)
                                         ccGTCAGTTGTGAATTATTGGTACGCAGTTCGCGATCGGTCAGGCCACGGCAGAAACTTGCCGCCGTAGAACAAACGGTTGTTGGCGTTACCGGTTGGTTGAGTTGAAGCAAACGGCCAATTGCCGCACATAGCAGCCCTAAAGAAAAAATGCTGCCTTTATGCGTTGTTACGCCCGCAGTGGCCCGGAACATATCACCTTCGCAAGCCATACCCATTGGGCTGAATCCGTGGCGTAc                                                                                                                                                                                                    >  2:732197/1‑238 (MQ=255)
                                         ccGTCAGTTGTGAATTATTGGTACGCAGTTCGCGATCGGTCAGGCCACGGCAGAAACTTGCCGCCGTAGAACAAACGGTTGTTGGCGTTACCGGTTGGTTGAGTTGAAGCAAACGGCCAATTGCCGCACATAGCAGCCCTAAAGAAAAAATGCTGCCTTTATGCGTGTTTACGCCCGCAGTGCCGCGGAACATATCACCTTCGCAAGCCATACCAATTGGGCGTAATCCGTGGAGTAc                                                                                                                                                                                                    >  2:661796/1‑238 (MQ=255)
                                                ttGTGAATTATTGGTACGCAGTTCGCGATCGGTCAGGCCACGGCAGAAACTTGCCGCCGTAGAACAAACGGTTGTTGGCGTTACCGGTTGGTTGAGTTGAAGCAAACGGCCAATTCCCGCACATAGCAGCCCTAAAGAAAAAATGCTGCCTTTATGCGTGTTTCGGCCCGCCGTGGCCCGGCACAATTAACCTTCCCAACCCATACAAATTGGGCGAAACCCGTGGacaacccgtc                                                                                                                                                                                               >  1:169446/1‑227 (MQ=255)
                                                ttGTAAATTAGTGGTCCGCGTTTCGGGATGGGCCGGGCCCCGGCAGAACCTGGCCGCCGTAGAACAAACGGTTGTGGGCGTTACCGGTTGGTTGAGTTGAAGCAAACGGCCAATTGCCGCACATAGCAGCCCTAAAGAAAAAATGCTGCCTTTATGCGTGTTTACGCCCGCAGTGGCGCGGAACATATCACCTTCGCAAGCCATACCAATTGGGCGTAATCCGTGGAGTACCGCTTCt                                                                                                                                                                                             <  2:477573/238‑1 (MQ=255)
                                                   tgAATTATTGGTACGCAGTTCGCGATCGGTCAGGCCACGGCAGAAACTTGCCGCCGTAGAACAAACGGTTGTTGGCGTTACCGGTTGGTTGAGTTGAAGCAAACGGCCAATTGCCGCACATAGCAGCCCTAAAGAAAAAATGCTGCCTTTATGCGGGTTTACGCCCGCAGTGGCGCGGAACATATCACCTTCGCAAGCCATACCAATTGGGCGTAATCCGTGGAGTACCGCTTCTggg                                                                                                                                                                                          >  2:228082/1‑237 (MQ=255)
                                                            ggTACGCAGTTCGCGATCGGTCAGGCCACGGCAGAAACTTGCCGCCGTAGAACAAACGGTTGTTGGCGTTACCGGTTGGTTGAGTTGAAGCAAACGGCCAATTGCCGCACATAGCAGCCCTAAAGAAAAAATGCTGCCTTTATGCGTGTTTACGCCCGCAGTGGCGCGGAACATATCACCTTCGCAAGCCATACCAATTGGGCGTAATCCGTGGAGTACCGCTTCTGGTGCCATTTcc                                                                                                                                                                                 >  2:692515/1‑238 (MQ=255)
                                                             gtccgattttcactgtgtggcgcgCCGCGGCAAAAGTTTGCCCCCGTAAAAAAAACGGTTTGGGGGGTGTCCGGTTGGTTGTATTGAAGGAACCGGCCAATTCCCGACCATACCACCCCTAAAGAAAAAAGGCTGCCTTTATGCGTGTTTACGCCCGCAGTGGCGCGGAACATATCCCCTCCGCAAGCCATACCAATTGGGCGTAATCCGTGGAGTACCGCTTCTGGTGCCATTTCCGc                                                                                                                                                                               <  1:434700/216‑1 (MQ=255)
                                                               aCGCAGTTCGCGATCGGTCAGGCCACGGCAGAAACTTGCCGCCGTAGAACAAACGGTTGTTGGCGTTACCGGTTGGTTGAGTTGAAGCAAACGGCCAATTGCCGCACATAGCAGCCCTAAAGAAAAAATGCTGCCTTTATGCGTGTTTACGCCCGCAGTGGCGCGGAACATATCACCTTCGCAACCCATCCAAATTGGGCGGAATCCGTGGAGTACGGCtcccggagcacttgcagcat                                                                                                                                                                             >  2:665021/1‑220 (MQ=255)
                                                               aCGCAGTTCGCGATCGGTCAGGCCACGGCAGAAACTTGCCGCCGTAGAACAAACGGTTGTTGGCGTTACCAGTTGGTTGAGTTGAAGCAAACGGCCAATTGCCGCACATAGCAGCCCTAAAGAAAAAATGCTGCCTTTATGCGTGTTTACGCCCCCAGTGGCGCGGAACATATCCCCTTATCAATCCAACCCAATTGGGGGAAATCCGGGGAGTACCGCTTCAGGGGCCCttcacccc                                                                                                                                                                              >  1:331033/1‑232 (MQ=255)
                                                                    ggtcgagtggggTCAGGCCCGGGCAGAAACTTGCCGCCGTAGAACAACCGGTTGTTGGCGTTACCGGTTGGTTGAGTTGAAGCAAACGGCCAATTGCCGCACAAAGCAGCCCTAAAGAAAAAAGGCTGCCTTTATGCGTGTTTACGCCCGCAGTGGCGCGGAACATATCACCTTCGCAAGCCATACCAATTGGGCGTAATCCGTGGAGTACCGCTTCTGGTGCCATTTCCGCACTACa                                                                                                                                                                         <  2:321176/228‑1 (MQ=255)
                                                                             gggccagcccccggaaaaattctGCCGCCTAAAACAAAAGGTGTTTGGCGGTTCCCGGTGGTGTGATGTGAACCAACGGGCAAATTCCCGCACAACCCACCCCAAAAAAAAAAATGCGCCTTTTTTGCGTGTTTCCGCCCGCAGTGGCGCGGAACATACCCCCTTCGCAAGCCATACCAATTGGGCGTAATCCGTGGAGTACCGCTTCTGGTGCCATTTCCGCACTACAGGCACCaaa                                                                                                                                                                <  1:572743/216‑1 (MQ=255)
                                                                               gTCAGGCCACGGCAGAAACTTGCCGCCGTAGAACAAACGGTTGTTGGCGTTACCGGTTGGTTGAGTTGAAGCAAACGGCCAATTGCCGCACATAGCAGCCCTAAAGAAAAAATGCTGCCTTTATGCGTGTTTCCGCCCGCAGTGGCGCGGAACATATCACCTTCGCAAGCCCTACCAATTGGGCGTAATCCGTGGAGTACCGCTTCTGGGGCCATTTCCCCCCTACAGGCACCAAAtca                                                                                                                                                             >  1:226881/1‑237 (MQ=255)
                                                                                      cggggaaaaactttGCCCCCGAAAAAAAAACGGTTGTTGGGTTTACCGGTTGGTTGAGTTGAAGCAACCGGCCAATTGCCCCACAAAGCGGCCCTAAAAAAAAAATGCTGCCTTTATGCGTTTTTACGCCCGCAGTGGCGCGGAACATATCCCCTTCGCAAGCCATACCAATTGGGCGTAATCCGTGGAGTACCGCTTCTGGTGCCATTTCCGCACTACAGGCACCAAATTCAATGaa                                                                                                                                                       <  1:762749/226‑1 (MQ=255)
                                                                                      cccggataaactttGCCGCCGGAGAACAAGCGGTTGTTGGCTTTACCGGTTGGTTGAGTTGAAGCAAACGGCCAATTGCCGCACATAACAGCCCTAAAGAAAAAATGCTGCCTTTATGCGTGTTTACGCCCGCAGTGGCGCGGAACATATCACCTTCGCAAGCCATACCAATTGGGCGTAATCCGTGGAGTACCGCTTCTGGTGCCATTTCCGCACTACAGGCACCAAATTCAATGaa                                                                                                                                                       <  2:816991/226‑1 (MQ=255)
                                                                                        cGGCAGAAACTTGCCGCCGTAGAACAAACGGTTGTTGGCGTTACCGGTTGGTTGAGTTGAAGCAAACGGCCAATTGCCGCACATAGCAGCCCTAAAGAAAAAATGCTGCCTTTATGCGTGTTTACGCCCGCAGTGGCGCGGAACATATCACCTTCGCAAGCCATACCAATTGGGCGTAATCCGTGGAGTACCGCTTCTGGTGCCATTTCTGCACTCCAGGCACCAAATTCAATGAAACg                                                                                                                                                    >  2:62586/1‑239 (MQ=255)
                                                                                              aaaCTTGCCGCCGTAGAACAAACGGTTGTTGGCGTTACCGGTTGGTTGAGTTGAAGCAAACGGCCAATTGCCGCACATAGCAGCCCTAAAGAAAAAATGCTGCCTTTATGCGTGTTTACGCCCGCAGTGGCGCGGAACATATCACCTTCGCAAGCCATACCAATTGGGCGTAATCCGTGGAGTACCGCTTCTGGTGCCATTTCCGCACTACAGGCACCAAATTCAATGAAACGGGGTa                                                                                                                                               >  2:737567/1‑238 (MQ=255)
                                                                                               accTTGCCGCAGTAAAAAAACGTTTTTTGGGCGTTACCGGTTGGTTGAGTGAAACCAAACGGCCAATTGCCGCACATAGCACCCCAAAAAAAAAAATGCTGCCTTTATGCGTGTTTACGCCCGGAGTGGCGCGGAACATATCCCCTTCGCAAGCCATACCAATTGGGCGTAATCCGTGGAGTACCGCTTCTGGTGCCATTTCCGCACTACAGGCACCAAATTCAATGAAACGGGGTagc                                                                                                                                             <  1:225295/237‑1 (MQ=255)
                                                                                                               aCAAACGTTTGTTGGCTGTACCGGTTGGTTGAGTTGAAGCAAACGGCCAATTGCCGCACATAGCAGCCCTAAAGAAAAAATGCTGCCTTTAGGCGTGTTTACGCCCGCAGTGGCGCGGAACATATCCCCTTCGCAAGCCATACCAATTGGGCGTAATCCGTGGAGTACCGCTTCTGGTGCCATTTCCGCACTACAGGCACCAAATTCAATGAAACGGGGTAGCCAGCCCTGAATCGCCa                                                                                                                             <  2:425561/239‑1 (MQ=255)
                                                                                                                   aCGGTTGTTGGCGTTACCGGTTGGTTGAGTTGAAGCAAACGGCCAATTGCCGCACATAGCAGCCCTAAAGAAAAAATGCTGCCTTTATGCGTGTTTACGCCCGCAGTGGCGCGGAACATATCACCTTCGCAAGCCATACCAATTGGGCGTAATCCGTGGAGTACCGCTTCTGGTGCCATTTCCGCCCTACAGGCACCAAATTCAATGAAACGGGGTAGCCAGCCCTGAATCGCCAacg                                                                                                                          >  2:722945/1‑238 (MQ=255)
                                                                                                                       ttgttgGCGTTCCCGGTGGGTTGAGTTGAACCAAACGGCCAATTGCCGCACATAGCGCCCTAAAAGAAAAAATGCTGCCTTTATGCGTGTTTACGCCCGCGGTGGCGCGGAACATATCCCCTTCGCAAGCCATACCAATTGGGCGTAATCCGTGGAGTACCGCTTCTGGTGCCATTTCCGCACTACAGGCACCAAATTCAATGAAACGGGGTAGCCAGCCCTGAATCGCCAGCGCGCt                                                                                                                      <  2:455329/238‑1 (MQ=255)
                                                                                                                           ttgggttttcacgggtttttggtttgaAACAAACGGCAATTGCCCGCAAATGGCAGCCCAAAAAAAAAAAGTCTCCTTTTATGCGTTTTTGCGCCGGCAGGGGCGCGGAACATATCACCTGCGCAAGCCATACCAATTGGGCGTAATCCGTGGAGTACCGCTTCTGGTGCCATTTCCGCACTACAGGCACCAAATTCAATGAAACGGGGTAGCCAGCCCTGAATCGCCAGCGCGCTGCg                                                                                                                   <  1:387273/216‑1 (MQ=255)
                                                                                                                                         ttggttgATTTGAAGCAAACGGCCAATTGCCGCAAATGGCAGCCCTAAAGAAAAAATGCTGCCTTTATGCGTGTTTACGCCCGCAGTGGCGCGGAACATATCACCTTCGCAAGCCATACCAATTGGGCGTAATCCGTGGAGTACCGCTTCTGGTGCCATTTCCGCACTACAGGCACCAAATTCAATGAAACGGGGTAGCCAGCCCTGAATCGCCAGCGCGCTGCGGTGGAAATCTTcc                                                                                                      <  2:792312/238‑1 (MQ=255)
                                                                                                                                          tggttgAGTGAAAGCAAACGGCAAATTGCCGCACATGGCAGCCTTAAAGAAAAAATGCTGCCTTTAGGCGTGTTTCCGCCCGAAGTGGCGCGGAACATATCACCTTCGCAAGCCATACCAATTGGGCGTAATCCGTGGAGTACCGCTTCTGGTGCCATTTCCGCACTACAGGCACCAAATTCAATGAAACGGGGTAGCCAGCCCTGAATCGCCAGCGCGCTGCGGTGGAAATCTTCCAg                                                                                                    <  2:308378/239‑1 (MQ=255)
                                                                                                                                            gggggtttgaAACAACCGCCCAATTGCCCCAAAAAGCGCCCTAAAAAAAAAACGCCCGCCTTTAGGCGTTTTTCCGCCGGGGGGGGCCCGGAACATATCCCCTCCGAAAGCCATACCAATTGGGCGTAATCCGTGGAGTACCGCTTCTGGTGCCATTTCCGCACTACAGGCACCAAATTCAATGAAACGGGGTAGCCAGCCCTGAATCGCCAGCGCGCTGCGGTGGAAATCTTCCAgcg                                                                                                  <  1:454308/233‑1 (MQ=255)
                                                                                                                                                      aGCAAACGGCCAATTGCCGCACATAGCAGCCCTAAAGAAAAAATGCTGCCTTTATGCGTGTTTACGCCCGCAGTGGCGCGGAACATATCACCTTCGCAAGCCATACCAATTGGGCGTAATCCGTGGAGTACCGCTTCTGGTGCCATTTCCGCCCTACAGGCACCAAATTCAATGAAACGGGGTAGCCAGCCCTGAATCGCCAGCGCGCTGCGGGGGAAATCTTCCAGCCCCATATCttt                                                                                        >  2:223043/1‑239 (MQ=255)
                                                                                                                                                           aCGGCCAATTGCCGCACATAGCAGCCCTAAAGAAAAAATGCTGCCTTTATGCGTGTTTACGCCCGCAGTGGCGCGGAACATATCACCTTCGCAAGCCATACCAATTGGGCGTAATCCGTGGAGTACCGCTTCTGGTGCCATTTCCGCACTACAGGCACCAAATTCAATGAAACGGGGTAGCCAGCCCTGAATCCCCAGCGCCCCGCGGGGGGAAACTTCCCACGCCAAATCttttggc                                                                                    >  1:8336/1‑234 (MQ=255)
                                                                                                                                                               gccaattccccacaaaggccgcCCAAAAAAAAAAAGGCTGCCTTTATGGGTGTTTACGCCCGGAGGGGCGGGGAACATTCCCCCTTCCCAAGCCATACCAATTGGGCGTAATCCGTGGAGTACCGCTTCTGGTGCCATTTCCGCACTACAGGCACCAAATTCAATGAAACGGGGTAGCCAGCCCTGAATCGCCAGCGCGCTGCGGTGGAAATCTTCCAGCGCCATATCTTTGTGCGCAc                                                                                <  1:737567/219‑1 (MQ=255)
                                                                                                                                                               ccAATTGCCGCACATAGCAGCCCTAAAGAAAGAATGCTGCCTTTATGCGTGTTTACGCCCGCAGGGGCCCGGAACATATCCCCTCCGCAAGCCATACCAATTGGGCGTAATCCGTGGAGTACCGCTTCTGGTGCCATTTCCGCACTACAGGCACCAAATTCAATGAAACGGGGTAGCCAGCCCTGAATCGCCAGCGCGCTGCGGTGGAAATCTTCCAGCGCCATATCTTTGTGCGCAc                                                                                <  2:743199/238‑1 (MQ=255)
                                                                                                                                                                     gccgccAATAGCAGCCCTAAAGAAAAAAGGCGTCCTTTATGCGGTTTTACGCCCGCAGTGGCGCGGAACATTTCCCCTTCGCAAGCCATACCAATTGGGCGTAATCCGTGGAGTACCGCTTCTGGTGCCATTTCCGCACTACAGGCACCAAATTCAATGAAACGGGGTAGCCAGCCCTGAATCGCCAGCGCGCTGCGGTGGAAATCTTCCAGCGCCATATCTTTGTGCGCACCGCAGtt                                                                         <  2:832434/239‑1 (MQ=255)
                                                                                                                                                                          cacaTAGCAGCCCTAAAGAAAAAATGCTGCCTTTATGCGTGTTTACGCCCGCAGTGGCGCGGAACATATCACCTTCGCAAGCCATACCAATTGGGCGTAATCCGTGGAGTACCGCTTCTGGTGCCATTTCCGCACTACAGGCACCAAATTCAATGAAACGGGGTAGCCAGCCCTGAATCGCCAGCGCGCTTCGGCGGAAATCTTCCAGCGCCATATCTTGTTGCAGCCCGCAGTTAAt                                                                      >  2:633753/1‑238 (MQ=255)
                                                                                                                                                                             ccacccccccTAAAAAAAAAATGCTTCTTTGTGACGTGTTCACGCCCGCGTGGGCCCGAACACTATCACCTCGGCAAGCCATACCAATTGGGCGTAATCCGTGGAGTACCGCTTCTGGTGCCATTTCCGCACTACAGGCACCAAATTCAATGAAACGGGGTAGCCAGCCCTGAATCGCCAGCGCGCTGCGGTGGAAATCTTCCAGCGCCATATCTTTGTGCGCACCGCAGTTAATGCGa                                                                  <  1:439761/232‑1 (MQ=255)
                                                                                                                                                                               agcagcCCTAAAGAAAAAATGCTGCCTTTATGCGTGTTTACGCCCGCAGTGGCGCGGAACATATCACCTTCGCAAGCCATACCAATTGGGCGTAATCCGTGGAGTACCGCTTCTGGTGCCATTTCCGCACTACAGGCACCAAATTCAATGAAACGGGGTAGCCAGCCCTGAATCCCCAGAGCGCTGCGGTGGAAATCTTCAAGCGCCATATCTTTTTGCGCACCCCGGTTAAAGCGATc                                                                >  1:156945/1‑239 (MQ=255)
                                                                                                                                                                                    cccTAAAGAAAAAATGCTGCCTTTATGCGTGTTTACGCCCGCAGTGGCGCGGAACATATCACCTTCGCAAGCCATACCAATTGGGCGTAATCCGTGGAGTACCGCTTCTGGTGCCATTTCCGCACTACAGGCACCAAATTCAATGAAACGGGGTAGCCAGCCCTGAATCGCCAGCGCGCTCCGGTGGAAATCTTCCAGCGCCCTATCTTTGTGCGCCCCGCAGTTAATGCGATCCACGa                                                           >  2:88450/1‑239 (MQ=255)
                                                                                                                                                                                    cccTAAAGAAAAAATGCTGCCTTTATGCGTGTTTACGCCCGCAGTGGCGCGGAACATATCACCTTCGCAAGCCATACCAATTGGGCGTAATCCGTGGAGTACCGCTTCTGGTGCCATTTCCGCACTACAGGCACCAAATTCAATGAAACGGGGTAGCCAGCCCTGAATCGCCAGCGCGCCGCGGTGGAAATCCTCCACCGCCATATCTTTGTGCGCACCGCAGTTAATGCGATCCACGa                                                           >  2:824664/1‑239 (MQ=255)
                                                                                                                                                                                     ccTAAAGAAAAAATGCTGCCTTTATGCGTGTTTACGCCCGCAGTGGCGCGGAACATATCACCTTCGCAAGCCATACCAATTGGGCGTAATCCGTGGAGTACCGCTTCTGGTGCCATTTCCGCACTACAGGCACCAAATTCAATGAAACGGGGTAGCCAGCCCTGAATCGCCAGCGCGCTGCGGTGGAAATCTTCCAGCGCCATATCTTTGTGCGCACCGCCGTTAATGCGATCCACGAg                                                          >  2:860830/1‑239 (MQ=255)
                                                                                                                                                                                      ccaaaaaaaaaaaGTCTGCCTTTATGCGTGTTTCCGCCGTGAGTGGCGCGGAACATTCACCCTTCGCAAGCCATACCAATTGGGCGTAATCCGTGGAGTACCGCTTCTGGTGCCATTTCCGCACTACAGGCACCAAATTCAATGAAACGGGGTAGCCAGCCCTGAATCGCCAGCGCGCTGCGGTGGAAATCTTCCAGCGCCATATCTTTGTGCGCACCGCAGTTAATGCGATCCACg                                                            <  2:162612/234‑1 (MQ=255)
                                                                                                                                                                                       tAAAGAAAAAATGCTGCCTTTATGCGTGTTTACGCCCGAAGGGGCGCGGAACATATCCCCTTCGCAAGCCATACCAATTGGGCGTAATCCGTGGAGTACCGCTTCTGGTGCCATTTCCGCACTACAGGCACCAAATTCAATGAAACGGGGTAGCCAGCCCTGAATCGCCAGCGCGCTGCGGTGGAAATCTTCCAGCGCCATATCTTTGTGCGCACCGCAGTTAATGCGATCCACGAGGc                                                        <  2:340413/239‑1 (MQ=255)
                                                                                                                                                                                          aaGAAAAAATGCTGCCTTTATGCGTGTTTACGCCCGCAGTGGCGCGGAACATATCACCTTCGCAAGCCATACCAATTGGGCGTAATCCGTGGAGTACCGCTTCTGGTGCCATTTCCGCACTACAGGCACCAAATTCAATGAAACGGGGTAGCCAGCCCTGAATCGCCAGCGCGCTGCGGTGGAAATCTTCCAGCGCCATATCTTTGTGCGCACCGCAGTTAATGCGATCCACGAgg                                                         >  2:482205/1‑236 (MQ=255)
                                                                                                                                                                                          aaGAAAAAATGCTGCCTTTATGCGTGTTTACGCCCGCAGTGGCGCGGAACATATCACCTTCGCAAGCCATACCAATTGGGCGTAATCCGTGGAGTACCGCTTCTGGTGCCATTTCCGCACTACAGGCACCAAATTCAATGAAACGGGGTAGCCAGCCCTGAATCGCCAGCGCGCTGCGGTGGAAATCTTCCAGCGCCATATCTTTGTGCGCACCGCAGTTAATGCGATCCACGAgg                                                         >  2:860283/1‑236 (MQ=255)
                                                                                                                                                                                          aaGAAAAAATGCTGCCTTTATGCGTGTTTACGCCCGCAGTGGCGCGGAACATATCACCTTCGCAAGCCATACAAATTGGGCGTAATCCGTGGAGTACCGCTTCTGGTGCCATTTCTGCACTACAGGCACCAAATTCAATGAAACGGGGTAGCCAGCCCTGAATCGCCAGCGCGCTGCGGTGGAAATCTTCCAGCGCCATATCTTTTTGCGCCCCGCAATTAATGCGTTCCCCGCgg                                                         >  2:749476/1‑236 (MQ=255)
                                                                                                                                                                                          aaGAAAAAATGCTGCCTTTATGCGTGTTTACGCCCGCAGTGGCGCGGAACATAGCACCTTCGCAAGCCATACCAATTGGGCGTAATCCGTGGAGTACCGCTTCTGGTGCCATTTCCGCACTACAGGCACCAAATTCAATGAAACGGGGTAGCCAGCCCTGAATCGCCAGCGCGCTTCGGTGGAAAACTTCCAGCGCCATATCTTTGTTAGCACCGCAGTTAAGGCGATcacctagg                                                         >  2:749464/1‑229 (MQ=255)
                                                                                                                                                                                           aGAAAAAATGCGTCTTTTATGCGTGTGTCCGCCCGCAGTGGCGCGGACCATATCCCCTGCGCAAGCCATACCAATTGGGCGTAATCCGTGGAGTACCGCTTCTGGTGCCATTTCCGCACTACAGGCACCAAATTCAATGAAACGGGGTAGCCAGCCCTGAATCGCCAGCGCGCTGCGGTGGAAATCTTCCAGCGCCATATCTTTGTGCGCACCGCAGTTAATGCGATCCACGAGGCCTg                                                     <  2:690689/239‑1 (MQ=255)
                                                                                                                                                                                              aaaaaTGTCTGCTTTTTTGGGGTTTTCCGCCCGGGGGGGCGCGAAACATTTCCCCTGCGAAAGCCATACCAATTGGGCGTAATCCGGGGAGTACCGCTTCTGGTGCCATTTCCGCACTACAGGCACCAAATTCAATGAAACGGGGTAGCCAGCCCTGAATCGCCAGCGCGCTGCGGTGGAAATCTTCCAGCGCCATATCTTTGTGCGCACCGCAGTTAATGCGATCCACGAGGCCTGGt                                                   <  1:66419/239‑1 (MQ=255)
                                                                                                                                                                                               aaaaTGCTGCCTGTAAGCGTTTTTGCGCCCGCAGTGGCGAGGAACATATCCCCTTCGAAAGCCATCCCAATTGGGCGTAATCCGGGGAGTACCGCTTCTGGTGCCATTTCCGCACTACAGGCACCAAATTCAATGAAACGGGGTAGCCAGCCCTGAATCGCCAGCGCGCTGCGGTGGAAATCTTCCAGCGCCATATCTTTGTGCGCACCGCAGTTAATGCGATCCACGAGGCCTGGttt                                                 <  2:408654/239‑1 (MQ=255)
                                                                                                                                                                                                      ctgcCTTTATGCGTGTTTACGCCCGCAGTGGCGCGGAACATATCACCTTCGCAAGCCATACCAATTGGGCGTAATCCGTGGAGTACCGCTTCTGGTGCCATTTCCGCACTACAGGCACCAAATTCAATGAAACGGGGTAGCCAGCCCTGAATCGCCAGCGCGCTGCGGTGGAAATCTTCCAGCGCCATCTCTTTTGGCCCACCGCAGTGAATGCAATCCACGAGGCCTGGTTTCGGGGa                                           >  1:378994/1‑239 (MQ=255)
                                                                                                                                                                                                         ccTTTATGCGTGTTTACGCCCGCAGTGGCGCGGAACATATCACCTTCGCAAGCCATACCAATTGGGCGTAATCCGTGGAGTACCGCTTCTGGTGCCATTTCCGCACTACAGGCACCAAATTCAATGAAACGGGGTAGCCAGCCCTGAATCGCCAGCGCGCTGCGGTGGAAATCTTCCAGCGCCATATCTTTGTGCGCACCGCAGTTAATGCGATCCCCGAGGCCTGGTTTCGGTGGCAg                                        >  2:625842/1‑239 (MQ=255)
                                                                                                                                                                                                             tATGCGTTTTTACGCCCGCAGGGGCGCGGAACATATCCCCTGCGCAAGCCATACCAATTGGGCGTAATCCGTGGAGTACCGCTTCTGGTGCCATTTCCGCACTACAGGCACCAAATTCAATGAAACGGGGTagc                                                                                                                                             <  1:334546/134‑1 (MQ=255)
                                                                                                                                                                                                                 ggttttttccgccgggtgtgggccgGAACTTTTCCCCTTCGCAATCCCTTCCAATTTGGGGTAATCCGTGGAGTACCGCTTCTGGGGCCATTTCCGCACTACAGGCACCAAATTCAATGAAACGGGGTAGCCAGCCCTGAATCGCCAGCGCGCTGCGGTGGAAATCTTCCAGCGCCATATCTTTGTGCGCACCGCAGTTAATGCGATCCACGAGGCCTGGTTTCGGTGACAGATTGACt                                 <  1:228082/216‑1 (MQ=255)
                                                                                                                                                                                                                 ggttgtTTACGCCGGCAGTGGCGCGGAACATATCCCCATCGCAAGCCATACCAATTGGGCGTAATCCGTGGAGTACCGCTTCTGGTGCCATTTCCGCACTACAGGCACCAAATTCAATGAAACGGGGTAGCCAGCCCTGAATCGCCAGCGCGCTGCGGTGGAAATCTTCCAGCGCCATATCTTTGTGCGCACCGCAGTTAATGCGATCCACGAGGCCTGGTTTCGGTGACAGATTGACt                                 <  2:497746/236‑1 (MQ=255)
                                                                                                                                                                                                                     gtTTACGCCCGCAGTGGCGCGGAACATATCACCTTCGCAAGCCATACCAATTGGGCGTAATCCGTGGAGTACCGCTTCTGGTGCCATTTCCGCACTACAGGCACCAAATTCAATGAAACGGGGTAGCCAGCCCTGAATCGCCAGCGCGCTGCGGTGGAAATCTTCCAGCGCCATATCTTTGTGCGCACCGCAGTTAATGCGATCCACGAGGCCTGGTTTCGGTGACAGATTGACTtca                              >  2:145306/1‑238 (MQ=255)
                                                                                                                                                                                                                      tttACGCCCGCAGTGGCGCGGAACATATCACCTTCGCAAGCCATACCAATTGGGCGTAATCCGTGGAGTACCGCTTCTGGTGCCATTTCCGCACTACAGGCACCAAATTCAATGAAACGGGGTAGCCAGCCCTGAATCGCCAGCGCGCTGCGGTGGAAATCTTCAAGCGCCCACTCTTTGTGCGCCCCCCAGTTAATGCGACCACCGAGGCCTGGTTTCGGTGAAAGATTGGATtcaga                            >  1:74017/1‑238 (MQ=255)
                                                                                                                                                                                                                      tttACGCCCGCAGTGGCGCGGAACATATCACCTTCGCAAGCCATACCAATTGGGCGTAATCCGTGGAGTACCGCTTCTGGTGCCATTTCCGCACTACAGGCACCAAATTCAATGAAACGGGGTAGCCAGCCCTGAATCGCCACCGCGCTGCGGTGGAAATCTTCCAGCGCCCTATCTTTTTGCGCACCGCAGTTAATGCGACCCCCGAGGCCTGGTTTCGTGGACAAATTGCCttacgc                            >  1:313344/1‑235 (MQ=255)
                                                                                                                                                                                                                        tACGCCCGCAGTGGCGCGGAACATATCACCTTCGCAAGCCATACCAATTGGGCGTAATCCGTGGAGTACCGCTTCTGGTGCCATTTCCGCACTACAGGCACCAAATTCAATGAAACGGGGTAGCCAGCCCTGaa                                                                                                                                   >  2:334546/1‑134 (MQ=255)
                                                                                                                                                                                                                        tACGCCCGCAGTGGCGCGGAACATATCACCTTCGCAAGCCATACCAATTGGGCGTAATCCGTGGAGTACCGCTTCTGGTGCCATTTCCGCACTACAGGCACCAAATTCAATGAAACGGGGTAGCCAGCCCTGAATCGCCAGCGCGCTGCGGTGGAAATCTTCCAGCGCCATATCTTTGTGCGCACCGCAGTTAATGCGATCCACGAGGCCTGGTTTCGGTGACAGATTGACTtccgtca                          >  2:323917/1‑239 (MQ=255)
                                                                                                                                                                                                                              cGCAGTGGCGCGGAACATATCACCTTCGCAAGCCATACCAATTGGGCGTAATCCGTGGAGTACCGCTTCTGGTGCCATTTCCGCACTACAGGCACCAAATTCAATGAAACGGGGTAGCCAGCCCTGAATCGCCAGCGCGCTGCGGTGGAAATCTTCCAGCGCCATATCTTTGTGCGCACCGCAGTTAATGCGATCCACGAGGCCTGGTTTCGGTGACAGATTGACTTCAGTCAGCATgg                    <  2:156945/239‑1 (MQ=255)
                                                                                                                                                                                                                                cAGTGGCGCGGAACATATCACCTTCGCAAGCCATACCAATTGGGCGTAATCCGTGGAGTACCGCTTCTGGTGCCATTTCCGCACTACAGGCACCAAATTCAATGAAACGGGGTAGCCAGCCCTGAATCGCCAGCGCGCTGCGGTGGAAATCTTCCAGCGCCATATCTTTGTgcg                                                                                   >  1:642938/1‑174 (MQ=255)
                                                                                                                                                                                                                                cAGTGGCGCGGAACATATCACCTTCGCAAGCCATACCAATTGGGCGTAATCCGTGGAGTACCGCTTCTGGTGCCATTTCCGCACTACAGGCACCAAATTCAATGAAACGGGGTAGCCAGCCCTGAATCGCCAGCGCGCTGCGGTGGAAATCTTCCAGCGCCATATCTTTGTgcg                                                                                   <  2:642938/174‑1 (MQ=255)
                                                                                                                                                                                                                                     tgagccgaaattatCACCTTAGCAAGCCATACCAATTGGGCTTAAGCCGTGGAGTTCCGCTTCTGGTGCCATTTCCGCACTACAGGCACCAAATTCAATGAAACGGGGTAGCCAGCCCTGAATCGCCAGCGCGCTGCGGTGGAAATCTTCCAGCGCCATATCTTTGTGCGCACCGCAGTTAATGCGATCCACGAGGCCTTGTTTCGGTGACAGATTGACTTCAGTCAGCATGgcgcgc               <  2:511776/227‑1 (MQ=255)
                                                                                                                                                                                                                                      cgcgGAACATATCACCTTCGCAAGCCATACCAATTGGGCGTAATCCGTGGAGTACCGCTTCTGGTGCCATTTCCGCACTACAGGCACCAAATTCAATGAAACGGGGTAGCCAGCCCTGAATCGCCAGCGCGCTGCGGTGGAAATCTTCCAGCGCCATATCTTTGTGCGCACCGCAGTTAATGCGATCCACGAGGCCTGGTTTCGGTGACAGATTGACTTCAGTCAGCATGGCGCGccc             >  2:238460/1‑237 (MQ=255)
                                                                                                                                                                                                                                            aCATATCCCCTTCGAAAGCCATACCAATTGGGCGTAATCCGTGGAGTACCGGTTCTGGTGCCATTTCCGCACTACAGGCACCAAATTCAATGAAACGGGGTAGCCAGCCCTGAATCGCCAGCGCGCTGCGGTGGAAATCTTCCAGCGCCATATCTTTGTGCGCACCGCAGTTAATGCGATCCACGAGGCCTGGTTTCGGTGACAGATTGACTTCAGTCAGCATGGCGCGccagcccag       <  2:226881/238‑1 (MQ=255)
                                                                                                                                                                                                                                             cATATCACCTTCGCAAGCCATACCAATTGGGCGTAATCCGTGGAGTACCGCTTCTGGTGCCATTTCCGCACTACAGGCACCAAATTCAATGAAACGGGGTAGCCAGCCCTGAATCGCCAGCGCGCTGCGGTGGAAATCTTCCAGCGCCATATCTTTGTGCGCACCGCAGTTACTGCGATCCACGAGGCCTGGTTTCGGTGACAGATTGACTTCAGTCAGCATGGCGCGCCAGCCcagca     >  2:32247/1‑239 (MQ=255)
                                                                                                                                                                                                                                             cATATCACCTTCGCAAGCCATACCAATTGGGCGTAATCCGTGGAGTACCGCTTCTGGTGCCATTTCCGCACTACAGGCACCAAATTCAATGAAACGGGGTAGCCAGCCCTGAATCGCCAGCGCGCTGCGGTGGAAATCTTCCAGCGCCATATCTTTGTGCGCACCGCAGTTAATGCGATCCACGAGGCCTGGTTTCGGTGACAGATTGACTtcagtcag                         <  2:719090/219‑1 (MQ=255)
                                                                                                                                                                                                                                             cATATCACCTTCGCAAGCCATACCAATTGGGCGTAATCCGTGGAGTACCGCTTCTGGTGCCATTTCCGCACTACAGGCACCAAATTCAATGAAACGGGGTAGCCAGCCCTGAATCGCCAGCGCGCTGCGGTGGAAATCTTCCAGCGCCATATCTTTGTGCGCACCGCAGTTAATGCGATCCACGAGGCCTGGTTTCGGTGACAGATTGACTtcagtcag                         >  1:719090/1‑219 (MQ=255)
                                                                                                                                                                                                                                                 atCTCCTTCGCACGCCCTTCCAATGGGGCGTAATCGGTGGAGGACCGCTTCTGGTGCCATTTCCGCACTACAGGCACCAAATTCAATGAAACGGGGTAGCCAGCCCTGAATCGCCAGCGCGCTGCGGTGGAAATCTTCCAGCGTCATATCTTTGTGCGCACCGCAGTTAATGCGATCCACGAGGCCTGGTTTCGGTGACAGATTGACTTCAGTCAGCATGGCGCGCCAGCCCAGCAggg  <  2:702327/239‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                            
TGCGCCGGTAAGGCCAAGCTGTTGGTACAACCGTTGACCTGCCGTCAGTTGTGAATTATTGGTACGCAGTTCGCGATCGGTCAGGCCACGGCAGAAACTTGCCGCCGTAGAACAAACGGTTGTTGGCGTTACCGGTTGGTTGAGTTGAAGCAAACGGCCAATTGCCGCACATAGCAGCCCTAAAGAAAAAATGCTGCCTTTATGCGTGTTTACGCCCGCAGTGGCGCGGAACATATCACCTTCGCAAGCCATACCAATTGGGCGTAATCCGTGGAGTACCGCTTCTGGTGCCATTTCCGCACTACAGGCACCAAATTCAATGAAACGGGGTAGCCAGCCCTGAATCGCCAGCGCGCTGCGGTGGAAATCTTCCAGCGCCATATCTTTGTGCGCACCGCAGTTAATGCGATCCACGAGGCCTGGTTTCGGTGACAGATTGACTTCAGTCAGCATGGCGCGCCAGCCCAGCAGGG  >  NC_000913/646973‑647445

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