Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 149,436 A→G 9.2% S56P (TCC→CCC)  panB ← 3‑methyl‑2‑oxobutanoate hydroxymethyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913149,4360AG 9.2% 199.1 / 2.9 81S56P (TCC→CCC) panB3‑methyl‑2‑oxobutanoate hydroxymethyltransferase
Reads supporting (aligned to +/- strand):  ref base A (25/43);  new base G (7/4);  total (34/47)
Fisher's exact test for biased strand distribution p-value = 1.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCACGTTCGGTCAGCATTTGTACGGTTTCTACCAGCCACTCACCGCCTTCAATTTTGACCATGTTAGCACCGGCACGCATAACCGTTGCGGCGTTTTCGAAGGCTTGTTCCGGCGTGGCATACGCCATAAACGGCAGGTCAGCCAGCAGCAGGCAGTTTGGTGCGCCGCGACGTACGGCGGCAGTGTGGTAGGCGATATCGGCAACGGTAACTGGCAGGGTGGAGTCGTGCCCCTGAACCGTCATGCCCAGCGAATCGCCCACCAGCATGACGTTAAGCCCTTCATCAGCAAAGAGTTTGGCGAAGCTATAGTCATAAGCGGTGATGGTCGCGAAACGTTTTTTTTCCTGTTTGTACTTCTGCAGTAAGGAGATGGTGGTCGGTTTCATAACGTATCCTGATAAATTGATGTTGTGCTGTCTGGCATTTTATCAGTCACATTGGTGGGG  >  NC_000913/149212‑149661
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                  
ggCACGTTCGGTCAGCATTTGTACGGTTTCTACCAGCCACTCACCGCCTTCAATTTTGACCATGTTAGCACCGGCACGCATAACCGTTGCGGCGTTTTCGAAGGCTTGTTCCGGCGTGGCATACGCCATAAACGGCAGGTCAGCCAGCAGCCGGCGGTTTGGTGCGCCGCCACGTACGGCGGCAGGCTGGGAGGCGGAATCGGCAACGGTGACTGGGAGGGGGGGGt                                                                                                                                                                                                                                  >  1:185074/1‑227 (MQ=255)
ggCACGTTCGGTCAGCATTTGTACGGTTTCTACCAGCCACTCACCGCCTTCAATTTTGACCATGTTAGCACCGGCACGCATAACCGTTGCGGCGTTTTCGAAGGCTTGTTCCGGCGTGGCATACGCCATAAACGGCAGGTCAGCCAGCAGCAGGCCGTTTGGTGCGCCGCGACGTACGGCGGCAGTGTGGGAGGCGATAGCGGCAACGGGAACGTGCAGGGGGGAGt                                                                                                                                                                                                                                  >  2:612564/1‑227 (MQ=255)
 gCACGTTCGGTCAGCATTTGTACGGCTTCTACCAGCCACTCACCGCCTTCCATTTTGACCATGTTAGCACCGGCACGCATAACCGTTGCGGCGTTTTCGAAGGCTTGTTCCGGCGTGGCATACGCCATAAACGGCAGGTCAGCCAGCAGCAGGCGGTTTGGTGCGCCGCGACGAACGGCGGCCGCTGGGTACGCGCTACCTGCAACGGCGACCGGCCGGGCGGGGTCg                                                                                                                                                                                                                                >  1:201679/1‑228 (MQ=255)
 gCACGTTCGGTAAGCATTTGTACGGTTTCTCCCAGCCACTCACCGCCTTCAATTTTGCCCATGTTAGCACCGGCACGCATAACCGTGGCGGCGTTTTCGAAGGCTTGGTCCGGCGTGGCATACGCCATAAACGGCAGGTCAGCCAGCAGCAGGCAGTTTGGTGCGCCGCGACGTACGGCGGCAGTGTGGTAGGCGATATCGGCAACGGTAACTGGCAGGGTGGAGTCg                                                                                                                                                                                                                                <  2:558221/228‑1 (MQ=255)
           tcacaaTTTGTACGTGTTCTACCGGCCACCCACCTCCTTCAATTTTGACCATGTTAGCACCGGCACGCATAACCTTTGCGGCGTTTTCGAAGGCTTGGTCCGGCGTGGCATACGCCATAAACGGCAGGTCAGCCAGCAGCAGGCAGTTTGGTGCGCCGCGACGTACGGCGGCAGTGTGGTAGGCGATATCGGCAACGGTAACTGGCAGGGTGGAGTCGTGCCCCTGaa                                                                                                                                                                                                                      <  2:405064/223‑1 (MQ=255)
           tCAGCTTTTTTCCGGTTTCGCCCCGCCACTCTCCGCCTTCAATTTTGACCCTGTTAGCCCCGGCACGCATAACCGTTGCGGCGTTTTCGAAGGCTTGTTCCGGCGTGGCATACGCCATAAACGGCAGGTCAGCCAGCAGCAGGCAGTTTGGTGCGCCGCGACGTACGGCGGCAGTGTGGTAGGCGATATCGGCAACGGTAACTGGCAGGGTGGAGTCGTGCCCCTGaa                                                                                                                                                                                                                      <  1:177705/228‑1 (MQ=255)
                           ttCTACCAGCCACTCACCGCCTTCAATTGTGACCAGGTTAGCACCGGCACGCATAACCGTTGCGGCGTTTTCGAAGGCTTGTTCCGGCGTGGCATACGCCATAAACGGCAGGTCAGCCAGCAGCAGGCAGTTTGGTGCGCCGCGACGTACGGCGGCAGTGTGGTAGGCGATATCGGCAACGGTAACTGGCAGGGTGGAGTCGTGCCCCTGAACCGTCATGCCCAGCg                                                                                                                                                                                                       <  2:66168/227‑1 (MQ=255)
                            tCTACCAGCCACTCACCGCCTTCAATTTTAACCATGTTAGCACCGGCACGCATAACCGTTGCGGCGTTTTCGAAGGCTTGTTCCGGCGTGGCATACGCCATAAACGGCAGGTCAGCCAGCAGCAGGCAGTTTGGTGCGCCGCGACGTACGGCGGCAGTGTGGTAGGCGATATCGGCACCGGTAACGGGCAGGGTGGAGTGGTGCCCCGGAACCGTCATGCCCAGCGa                                                                                                                                                                                                      >  2:684843/1‑227 (MQ=255)
                                 cAGCCACTCCCCGCCTTCAATTTTGACCATGTTAGCACCGGCACGCATAACCGTTGCGGCGTTTTCGAAGGCTTGTTCCGGCGTGGCATACGCCATAAACGGCAGGTCAGCCAGCAGCAGGCAGTTTGGTGCGCCGCGACGTACGGCGGCAGTGTGGTAGGCGATATCGGCAACGGTAACTGGCAGGGTGGAGTCGTGCCCCTGAACCGTCATGCCCAGCGAATCGcc                                                                                                                                                                                                <  2:43703/228‑1 (MQ=255)
                                    ccACTCACCGCCTTCAATTTTGACCATGTTAGCACCGGCACGCATAACCGTTGCGGCGTTTTCGAAGGCTTGTTCCGGCGTGGCATACGCCATAAACGGCAGGTCAGCCCGCCGCAGGCAGTTTGGTGCGCCGCGACGTCCGGCGGCAGGGTGGTAGGCGAAATCGGCAACGGCAACTGGCAGGGCGGAGTCGTGCCCCTGAACCGCCCGGCCCAGCGAATCGccccc                                                                                                                                                                                             >  1:506397/1‑226 (MQ=255)
                                    ccACTCACCGCCTTCAATTTTGACCATGTTAGCACCGGCACGCATAACCGTTGCGGCGTTTTCGAAGGCTTGTTCCGGCGTGGCATACGCCATAAACGGCAGGTCAGCCAGCAGCAGGCAGGTTGGTGCGCCGCGACGTACGGCGGCAGTGGGGTAGGCGATGTTGGCGACGGGAACTGGCGGGGGGGGAGCGGGCCCCTGAACCCTCAGGCCCCGCAAAACccccca                                                                                                                                                                                             >  1:597036/1‑226 (MQ=255)
                                    ccACTCACCGCCTTCAATTTTGACCATGTTAGCACCGGCACGCATAACCGTTCCGGCGTTTTCGAAGGCTTGTTCCGGCGTGGCATACGCCATAAACGGCAGGTCAGCCAGCAGCAGGCAGGTTGGTGCGCCCCGACGCACGGCGGCAGGGTGGTAGGCGAAATCTGACACGGGAACGGGCAGGGGGGGGGCGGGCCACCGAACCGtgcgccccacgcaatcccacac                                                                                                                                                                                             >  1:224319/1‑207 (MQ=255)
                                         ccccGGCTTTAATTTTGACCCTGTTTGCCCCGTGACGCATACCCGTGGCGGCGTTTTCGAAGGCTTGTTCCGGCGTGGCATGCGCCATAAACGGCAGGGCAGCCAGCAGCAGGAAGTTTGGTGCGCCGCGACGTACGGCGGCAGTGTGGTAGGCGATATCGGCAACGGTAACTGGCAGGGTGGAGTCGTGCCCCTGAACCGTCATGCCCAGCGAATCGCCCACCAGCa                                                                                                                                                                                        <  1:140987/226‑1 (MQ=255)
                                         cACCGCCTTCAATTTTGGCCATGTTAGCACCGGCACGCATAACCGTTGCGGCGTTGTCGAAGGCTTGTTCCGGCGTGGCATACGCCATAAACGGCAGGTCAGCCAGCAGCAGGCAGTTTGGTGCGCCGCGACGTACGGCGGCAGTGTGGTAGGCGATATCGGCAACGGTAACTGGCAGGGTGGAGTCGTGCCCCTGAACCGTCATGCCCAGCGAATCGCCCACCAGCa                                                                                                                                                                                        <  2:113801/228‑1 (MQ=255)
                                           ccGCCTTCAATTTTGACCATGTTTGCCCCGGCACGCATAACCGTTGCGGCGTTTTCGAAGGCTTGGTCCGGCGTGGCATACGCCATAAACGGCAGGTCAGCCAGCAGCAGGCAGTTTGGTGCGCCGCGACGTACGGCGGCAGTGTGGTAGGCGATATCGGCAACGGTAACTGGCAGGGTGGAGTCGTGCCCCTGAACCGTCATGCCCAGCGAATCGCCCACCAGCATg                                                                                                                                                                                      <  2:597036/228‑1 (MQ=255)
                                           ccGCCTTCAATGTTGACCATGTTGGCACCGGCACGCATAACCGTTGCGGCCTTTTCGAAGGCTTGTTCCGGCGTGGCATACGCCGTAAACGGCAGGTCAGCCAGCAGCAGGCAGTTTGGTGCGCCGCGACGTACGGCGGCAGTGTGGTAGGCGATATCGGCAACGGTAACTGGCAGGGTGGAGTCGTGCCCCTGAACCGTCATGCCCAGCGAATCGCCCACCAGCATg                                                                                                                                                                                      <  2:289445/228‑1 (MQ=255)
                                              ccTTCAATTTTGACCATGTTAGCACCGGCACGCATAACCGTTGCGGCGTTTTCGAAGGCTTGTTCCGGCGTGGCATACGCCATAAACGGCAGGTCAGCCAGCAGCCGGCAGTTTGGTGCGCCGCGACGTACGGCGGCAGTGTGGTAGGCGATATCGGCAACGGGACCTGGCAGGGTGGAGTCGTGCCCCTGAACCGTCATGTCCAGCGAACCGCCCACCAGCATGACg                                                                                                                                                                                   >  2:641715/1‑228 (MQ=255)
                                                 tCAATTTTGACCATGTTAGCCCCGGCACGCATAACCGTTGCGGCGTTTTCGAAGGCTTGGTCCGGCGTGGCATACGCCATAAACGGCAGGTCAGCCAGCAGCAGGCAGTTTGGTGCGCCGCGACGTACGGCGGCAGTGTGGTAGGCGATATCGGCAACGGTAACTGGCAGGGTGGAGTCGTGCCCCTGAACCGTCATGCCCAGCGAATCGCCCACCAGCATGACGTTa                                                                                                                                                                                <  2:54154/228‑1 (MQ=255)
                                                                ttAGCCCCGGCTCGCATAACCGGTGCGGCGTTTTCGAAGGCTTGTTCCGGCGTGGCATACGCCATAAACGGCAGGTCAGCCAGCAGCAGGCAGTTTGGTGCGCCGCGACGTACGGCGGCAGTGTGGTAGGCGATATCGGCAACGGTAACTGGCAGGGTGGAGTCGTGCCCCTGAACCGTCATGCCCAGCGAATCGCCCACCAGCATGACGTTAAGCCCTTCATCAGCa                                                                                                                                                                 <  2:229456/228‑1 (MQ=255)
                                                                         gCACGCATAACCGTTGCGGCGTTTTCGAAGGCTTGTTCCGGCGTGGCATACGCCATAAACGGCAGGTCAGCCAGCAGCAGGCAGTTTGGTGCGCCGCGACGTACGGCGGCAGTGTGGTAGGCGATATCGGCAACGGTAACTGGCAGGGGGGAGGCGTGCCCCTGGAACGTCCTGCCCCGGCAATCGCCCACCAGCATGGCGTTAAGCCCCTCCTCAACAAAAAGGTTg                                                                                                                                                        >  1:482885/1‑228 (MQ=255)
                                                                         gCACGCATAACCGTTGCGGCGTTTTCGAAGGCTTGTTCCGGCGTGGCATACGCCATAAACGGCAGGTCAGCCAGCAGCAGGCAGTTTGGTGCGCCGCGACGTACGGCGGCAGCGTGGTAGGCGATATCGGCAACGGTAACTGGCAGGGTGGAGTCGTGCCCCTGAACCGTCATGCCCAGCGAATCGCCCACCAGCATGACGTTAAGCCCTTCATCAGCAAAAAGCTTg                                                                                                                                                        >  2:278855/1‑228 (MQ=255)
                                                                         gCACGCATAACCGTTGCGGCGTTTTCGAAGGCTTGTTCCGGCGTGGCATACGCCATAAACGGCAGGTCAGCCAGCAGCAGGCAGTTTGGTGCGCCGCGACGAACGGCGGCAGTGTGGTAGGCGATATCGGCAACGGTAACTGGCAGGGGGGAGTCGTGCCCCTGAACCGCCATGCCCAGCGAATCGCCCACCAGCATGACGGTAAGCCCTTCATCAGCAAAACGTTTg                                                                                                                                                        >  2:595109/1‑228 (MQ=255)
                                                                         gCACGCATAACCGTTGCGGCGTTTTCGAAGGCTTGTTCCGGCGTGGACTCCGCCATAAACGGCAGGTCAGCCACCAGCAGCGCGTTTGGGGCGCCGCGACGTACGGCGGCAGTGTGGTGGGCGATATCGGCAACGGTGCCTGGACGGGGGGGGCCGTGCCCCTGGACCGTCGTGCCCCGCGACGCGCCCCCCCGCATGACCGTAAGCCCCtcatcatccaacgacctt                                                                                                                                                        >  1:576780/1‑216 (MQ=255)
                                                                         gCACGCATAACCGCTGCGGCGTTTTCGAAGGCTTGTTCCGGCGTGGCATACGCCATAAACGGCAGGTCAGCCAGCAGCAGGCAGTTTGGTGCGCCGCGACGTACGGCGGCAGTGCGGTAGGCGATATCGGCAACGGAAACTGGCAGGGTGGAGTCGGGCCCCGGAACCGGCATGCCCCGGGGAACGCACCCCAGCAGGCCGTTAAGCCCTCTATAAGCCAAGAGTTTg                                                                                                                                                        >  1:105364/1‑228 (MQ=255)
                                                                          cACGCATAACCGTTGCGGCGTTTTCGAAGGCTTGTTCCGGCGTGGCATACGCCATAAACGGCAGGTCAGCCAGCAGCAGGCAGTTTGGTGCGCCGCGACGTACGGCGGCAGTGTGGTAGGCGATATCGGCAACGGTAACTGGCAGGGTGGAGTCGTGCCCCTGAACCGTCATGCCCAGCGAATCGCCCACCAGCATGACGTTAAGCCCTTCATCAGCAAAGAGTTTgg                                                                                                                                                       <  2:651875/228‑1 (MQ=255)
                                                                          cACGCATAACCGTTGCGGCGTTTTCGAAGGCTTGTTCCGGCGTGGCATACGCCATAAACGGCAGGTCAGCCAGCAGCAGGCAGTTTGGTGCGCCGCGACGTACGGCGGCAGTGTGGGAGGCGATATCGGCAACGGTAACTGGCAGGGGGGCGTCGTGCCCCTGAACCGTCATGCCCCGCGAATCGGCCACCAGCATGACGTTAAGCCCTTCATCAGCAAAGAGTTTg                                                                                                                                                        >  2:238707/1‑227 (MQ=255)
                                                                           ccGCATAACCGTTGGGGCGTTTTCCGAGGCTTGTTCCGGCGTGGCATACGCCATAAACGGCAGGTCAGCCAGCAGCAGGTAGTTTGGTGCGCCGCGACGTACGGCGGCAGTGTGGTAGGCGATATCGGCAACGGTAACTGGCAGGGTGGAGTCGTGCCCCTGAACCGTCATGCCCAGCGAATCGCCCACCAGCATGACGTTAAGCCCTTCATCAGCAAAGAGTTTGGc                                                                                                                                                      <  2:185074/227‑1 (MQ=255)
                                                                           acgcagacgctgggcggcgTTTGCGAGGGTGTGTCCCGGCGGGGCATACGCGATAAACGGCAGGGCAGGCAGCAGCAGGCAGTTTGGTGCGCCGCGACGTACGGCGGCAGTGGGGTAGGCGATATCGGCAACGGTAACTGGCAGGGTGGAGGCGTGCCCCTGAACCGTCATGCCCAGCGAATCGCCCACCAGCATGACGTTAAGCCCTTCATCAGCAAAGAGTTTGGc                                                                                                                                                      <  1:612564/215‑1 (MQ=255)
                                                                             gCATAACCGTTGCGGCGTGTGCGAGGGCTTGGTCCGGCGTGGCAGACGCCATAAGCGGGAGGTCAGCCAGAAGCAGGCAGTTTGGTGCGCCGCGACGTACGGCGGCAGTGTGGTAGGCGATATCGGCAACGGTAACTGGCAGGGTGGAGTCGTGCCCCTGAACCGTCATGCCCAGCGAATCGCCCACCAGCATGACGTTAAGCCCTTCATCAGCAAAGAGTTTGGCGa                                                                                                                                                    <  2:329267/228‑1 (MQ=255)
                                                                              gacccccGTTGCGGCGTTTTCGAAGGCTTGTTCCGGCGGGGCATACGCCAGAAACGGCAGGTCAGCCAGCAGCAGGCAGTTTGGTGCGCCGCGACGTACGGCGGCAGTGTGGTAGGCGATATCGGCAACGGTAACTGGCAGGGTGGAGTCGTGCCCCTGAACCGTCATGCCCAGCGAATCGCCCACCAGCATGACGTTAAGCCCTTCATCAGCAAAGAGTTTGGCGa                                                                                                                                                    <  2:1712/222‑1 (MQ=255)
                                                                                    cttggcggggTTTTCGAAGGGTTGTTCCGGGGTGGCATACGGCCTAAGCGGCAGGTCAGCCGGCAGCAGGAAGTTTGGTGCGCCGCGACGTACGGCGGCAGTGTGGTAGGCGATATCGGCAACGGTAACTGGCAGGGTGGAGTCGTGCCCCTGAACCGTCATGCCAAGCGAATCGCCCACCAGCATGACGTTAAGCCCTTCATCAGCAAAGAGTTTGGCGAAGCtata                                                                                                                                             <  1:581014/224‑1 (MQ=255)
                                                                                     gTTGCCGCGTTTTCGAAGGCTTGTTCCGGCGTGGCCTGCGCCATAAACGGCAGGTCAGGCAGCAGCAGGCAGTTTGGTGCGCCGCGACGTACGGCGGCAGTGTGGTAGGCGATATCGGCAACGGTAACTGGCAGGGTGGAGTCGTGCCCCTGAACCGTCATGCCCAGCGAATCGCCCACCAGCATGACGTTAAGCCCTTCATCAGCAAAGAGTTTGGCGAAGCtata                                                                                                                                             <  2:589278/227‑1 (MQ=255)
                                                                                         cggcgTTTTAGAAGTGTAGTTGCGGGTGTGGATCCCCCAGAAACGGCAGGTCAGCCGGCAGCCGGCAGTTGGGGTCGCCGCGACGTACGGGGGCAGTGTGGTAGGCGATATCGGCAACGGTAACTGGCAGGGTGGAGTCGTGCTCCTGAACCGTCATGCCCAGCGAATCGCCCACCAGCATGACGTTAAGCCCTTCATCAGCAAAGAGTTTGGCGAAGCTATAGTCAt                                                                                                                                        <  1:165121/228‑1 (MQ=255)
                                                                                           gcgTTTTCGAAGGCTTGTTCCGGCGTGGCATACGCCATAAACGGCAGGTCAGCCAGCAGCAGGCAGTTTGGTGCGCCGCGACGTACGGCGGCAGTGTGGTAGGCGATATCGGCAACGGTAACTGGCAGGGTGGAGTCGTGCCCCTGAACCGTCATGCCCAGCGAATCGCCCACAAGCATGACGGTAAGCCCTTCATCAGCAAAGAGTTTGGCGAAGCCATAGTCAGaa                                                                                                                                      >  2:77611/1‑228 (MQ=255)
                                                                                                ttCGAAGGCTTGTTCCGGCGTGGCATACGCCATAAACGGCAGGTCAGCCAGCAGCAGGCAGTTTGGTGCGCCGCGACGTACGGCGGCAGGGGGGTAGGCGATATCGGCAACGGAAACTGGCCGGGCGGCGGCGTGCCCCTGAACCTTCATGCCCAGCGCAACGCCCACCAGCCTGACGTTAAGCCCTTCATCCGCAAACAGATTTGCGCAGCTATCGTCATAAACggg                                                                                                                                 >  1:506665/1‑227 (MQ=255)
                                                                                                  cGAAGGCTTGTTCCGGCGTGGCATACGCCATAAACGGCAGGTCAGCCAGCAGCAGGCAGTTTGGTGCGCCGCGACGTACGGCGGCAGTGTGGTAGGCGATATCGGCAACGGTAACTGGCAGGGTGGAGTCGTGCCCCTGAACCGTCATGCCCAGCGAATCGCCCACCAGCATGACGTTAAGCCCTTCATCAGCAAAGAGTTTGGCGAAGCTATAGTCATAAGCGGTGa                                                                                                                               <  2:542885/228‑1 (MQ=255)
                                                                                                           gTTCCGGCGTGGCATACGCCATAAACGGCAGGTCAGCCAGCAGCAGGCAGTTTGGTGCGCCGCGACGTACGGCGGCAGTGTGGTAGGCGATATCGGCAACGGTAACTGGCAGGGTGGAGTCGTGCCCCTGAACCGTCATGCCCAGCGAATCGCCCACCAGCATGACGTTAAGCCCTTCATCAGCAAAGAGTTTGTGGAAGCTATAGTCATAACCGGGGATGGTcccgt                                                                                                                      >  2:671688/1‑227 (MQ=255)
                                                                                                           gTTCCGGCGTGGCATACGCCATAAACGGCAGGTCAGCCAGCAGCAGGCAGTTTGGTGCGCCGCGACGTACGGCGGCAGTGTGGTAGGCGATATCGGCAACGGTAACTGGCAGGGTGGAGTCGTGCCCCTGAACCGTCATGCCCAGCGAATCGCCCACCAGCATGACGTTAAGCCCTTCATCAGCAAAGAGTTTGGCGAAGCTATCGTCATAAGCGGGGATGGTCCCGa                                                                                                                      >  2:376127/1‑228 (MQ=255)
                                                                                                           gTTCCGGCGTGGCATACGCCATAAACGGCAGGTCAGCCAGCAGCAGGCAGTTTGGTGCGCCGCGACGTACGGCGGCAGTGGGGTAGGCGATATCGGCAACGGTAACTGGCAGGGTGGGGTCGTGCCCCTGAACCGTCATGCCCAGCGAATCGCCCACCAGCATGGCGTTAGACCCTTCCTCAGCAAAGAGTTTGGCGAAACCATAATCaaaacggggtaaggaggcga                                                                                                                      >  1:112694/1‑212 (MQ=255)
                                                                                                                 gCGTGGCATACGCCATAAACGGCAGGTCAGCCAGCAGCAGGCAGTTTGGTGCGCCGCGACGTACGGCGGCAGTGTGGTAGGCGATATCGGCAACGGTAACTGGCAGGGTGGAGTCGTGCCCCTGAACCGTCATGCCCAGCGAATCGCCCACCAGCATGACGTTAAGCCCTTCATCAGCAAAGAGTTTGGCGAAGCTATAGTCATAAGCGGTGATGGTCGCGAAACGtt                                                                                                                >  2:448172/1‑228 (MQ=255)
                                                                                                                                  aaaCGGCAGGTCAGCCAGCAGCAGGCAGTTTGGTGCGCCGCGACGTACGGCGGCAGTGTGGTAGGCGATATCGGCAACGGTAACTGGCAGGGGGGAGTCGGGCCCCTGAACCGCCATGCCCAGCGAATCGCCCACCAGCATGACGTTAAGCCCTTCCTACGCAAAGAGTTGGGCGAAGCTATAGTCATAAGGGGTGATGGGCGCGAAACGGTTTTTTTCCTGGTTGTa                                                                                                >  1:654711/1‑228 (MQ=255)
                                                                                                                                   aaCGGCAGGTCAGCCAGCAGCAGGCAGTTTGGTGCGCCGCGACGTACGGCGGCAGTGTGGTAGGCGATATCGGCAACGGTAACTGGCAGGGTGGAGTCGTGCCCCTGAACCGTCATGCCCAGCGAATCGCCCACCAGCATGACGTTAAGCCCTTCATCAGCAAAGAGTTTGGCGAAGCTATAGTCATAAGCGGTGATGGTCGCGAAACGTTTTTTTTCCTGTTTGTa                                                                                                <  2:224319/227‑1 (MQ=255)
                                                                                                                                        cAGGTCAGGCAGGATGAAGCAGGTTGGTGCGCCGCGACGTACGGCGGCAGTGTGGTAGGCGATATCGGCAACGGTAACTGGCAGGGTGGAGTCGTGCCCCTGAACCGTCATGCCCAGCGAATCGCCCACCAGCATGACGTTAAGCCCTTCATCAGCAAAGAGTTTGGCGAAGCTATAGTCATAAGCGGTGATGGTCGCGAAACGTTTTTTTTCCTGTTTGTACTTCTg                                                                                          <  2:139717/228‑1 (MQ=255)
                                                                                                                                        cAGGTCAGCCCGCAGCAGGCAGTTTGGGTCGCCGCGACGTACGGCGGCAGTGGGGTAGGCGATATCGGCAACGGTAACTGGCAGGGTGGAGTCGTGCCCCTGAACCGTCATGCCCAGCGAATCGCCCACCAGCATGACTTTAAGCCCTTCATCAGCAAAGAGTTTGGCGAAGCTATAGTCATAAGCGGTGATGGTCGCGAAACGTTTTTTTTCCTGTTTGTACTTCTg                                                                                          <  2:691213/228‑1 (MQ=255)
                                                                                                                                          ggTCAGCCAGCAGCAGGCAGTTTGGTGCGCCGCGACGTACGGCGGCAGTGTGGTAGGCGATATCGGCAACGGTAACTGGCAGGGTGGAGTCGTGCCCCTGAACCGGCATGCCCAGCGAATCGCCCACCAGCATGACGTTAAGCCCTTCATCAGCAAAGAGTTTGGCGAAGCTATAGTCATAAGCGGTGATGGTCGCGAACCGTTTTTTTTCCTGTTTGTACTTCTGcg                                                                                        >  2:592260/1‑227 (MQ=255)
                                                                                                                                           gtgagccagcAGCAGGGAAGTGGGTGCGCCGCGACGTACGGCGGGAGTTTGGTAGGCGATATCGGCCACGGTAACTGGCAGGGTGGAGTCGTGCCCCTGACCCGTCATGCCCAGCGAATCGCCCACCAGCATGACGTTAAGCCCTTCATCAGCAAAGAGTTTGGCGAAGCTATAGTCATAAGCGGTGATGGTCGCGAAACGTTTTTTTTCCTGTTTGTACTTCTGCAg                                                                                       <  2:676851/225‑1 (MQ=255)
                                                                                                                                           ggtagccgggAGCGGGCAGTTGGGGCCGCCGCGACGTACGGCGGCATTGCGGTAGGCGACATCGGCAACGGTAACTGGCAGGGTTGAGGCGTGCCCCTGAACCGTCATGCCCTGCGAATCGCCCACCAGCATGACGTTAAGTCCTTCATCAGCAAAGAGTTTGGCGAAGCTATAGTCATAAGCGGTGATGGTCGCGAAACGTTTTTTTTCCTGTTTGTACTTCTGCAg                                                                                       <  1:609448/225‑1 (MQ=255)
                                                                                                                                           gccggccaggAGCCGGCAGGTTGGTGCGCCGCGACGTACGGCGGCAGTGTGGTAGGCGATATCGGCAACGGTAACTGGCAGGGTGGAGTCGTGCCCCTGAACCGTCATGCCCAGCGAATCGCCCACCAGCATGACGTTAAGCCCTTCATCAGCAAAGAGTTTGGCGAAGCTATAGTCATAAGCGGTGATGGTCGCGAAACGTTTTTTTTCCTGTTTGTACTTCTGCAg                                                                                       <  2:654711/224‑1 (MQ=255)
                                                                                                                                            tcaggcgggggcgggaggTTTGGTGGGCCGGGGCGGAGGGGGGCCGGTGGGGAGGCAAGATCGGCAACGGTACCTGGCAGGGTGGGTTCGTGCCCCTGAACCGTCATGCCCAGCGAATCGCCCACCAGCATGACGTTAAGCCCTTCATCAGCAAAGAGTTTGGCGAAGCTATAGTCATAAGCGGTGATGGTCGCGAAACGTTTTTTTTCCTGTTTGTACTTCTGCAGt                                                                                      <  1:684843/211‑1 (MQ=255)
                                                                                                                                                ccagcAGCAGGCAGTTTGGTGCCCCGCGACGTGCGGCGGCAGTGTGGTAGGCGATATCGGCAACGGTAACTGGCAGGGTGGAGTCGTGCCCCTGAACCGTCATGCCCAGCGAATCGCCCACCAGCATGACGTTAAGCCCTTCATCAGCAAAGAGTTTGGCGAAGCTATAGTCATAAGCGGTGATGGTCGCGAAACGTTTTTTTTCCTGTTTGTACTTCTGCAGTAAgg                                                                                 <  2:258151/228‑1 (MQ=255)
                                                                                                                                                          gcagTTTGGTGCGCCGCGACGTACGGCGGCAGTGTGGTAGGCGATAACGGCAACCGTAACTGGCAGGGTGGAGTCGTGCCCCTGAACCGTCATGCCCAGCTAATCGCACACCAGCATGACGTTAAGCCCTTCATCAGCAAAGAGTTTGGCGAAGCTATAGTCATAAGCGGTGATGGTCGCGAAACGTTTTTTTTCCTGTTTGTACTTCTGCAGTAAGGAGATggtggt                                                                       <  2:576780/228‑1 (MQ=255)
                                                                                                                                                           aggTTGGGTGCGCGGGAGGTACGGCGGCAGTGTGGTAGGCGATATCGGCAACGGTAACTGGCAGGGTGGAGTCGTGCCCCTGAACCGTCATGCCCAGCGAATCGCCCACCAGCATGACGTTAAGCCCTTCATCAGCAAAGAGTTTGGCGAAGCTATAGTCATAAGCGGTGATGGTCGCGAAACGTTTTTTTTCCTGTTTGTACTTCTGCAGTAAGGAGATggtggt                                                                       <  2:201679/224‑1 (MQ=255)
                                                                                                                                                            agTTTGGTGCGCCGCGACGTACGGCGGCAGTGTGGTAGGCGATATCGGCAACGGTAACTGGCAGGGTGGAGTCGTGCCCCTGAACCGTCATGCCCAGCGAATCGCCCACCAGCATGACGTTAAGCCCTTCATCAGCAAAGAGTTTGGCGAAGCTATAGTCATAAGCGGTGATGGTCGCGAAACGTTTTTTTTCCTGTTTGTACTTCTGCAGTAAGGAGATggtggt                                                                       <  2:112694/226‑1 (MQ=255)
                                                                                                                                                            agTTTGGGGCGCCGCGACGTACGGCGGCAGTGTGGTAGGCGATATCGGCAACGGTAACTGGCAGGGTGGAGTCGTGCCCCTGAACCGTCATGCCCAGCGAATCGCCCACCAGCATGACGTTAAGCCCTTCATCAGCAAAGAGTTTGGCGAAGCTATAGTCATAAGCGGTGATGGTCGCGAAACGTTTTTTTTCCTGTTTGTACTTCTGCAGTAAGGAGATggtggt                                                                       <  2:386589/226‑1 (MQ=255)
                                                                                                                                                            agTTGGGTGCGCCGGCACGTACGGCGGCAGTTGGGTAGGCGATATCGGCAACGGTAACTGGCAGGGTGGAGTCGTGCCCCTGAACCGTCATGCCCAGCGAATCGCCCACCAGCATGACGTTAAGCCCTTCATCAGCAAAGAGTTTGGCGAAGCTATAGTCATAAGCGGTGATGGTCGCGAAACGTTTTTTTTCCTGTTTGTACTTCTGCAGTAAGGAGATggtggt                                                                       <  2:482885/226‑1 (MQ=255)
                                                                                                                                                              tttGGTGCGCCGCGACGTACGGCGGCAGTGTGGTAGGCGATATCGGCAACGGTAACTGGCAGGGTGGAGTCGTGCCCCTGAACCGTCATGCCCAGCGAATCGCCCACCAGCATGACGTTAAGCCCTtc                                                                                                                                                                        >  1:20605/1‑128 (MQ=255)
                                                                                                                                                               ttGGTGCGCCGCGACGTACGGCGGCAGTGTGGTAGGCGATATCGGCAACGGTAACTGGCAGGGTGGAGTCGTGCCCCTGAACCGTCATGCCCAGCGAATCGCCCACCAGCATGACGTTAAGCCCTtca                                                                                                                                                                       <  2:20605/128‑1 (MQ=255)
                                                                                                                                                                 ggTGGGGCGCGAGGTTGGGGGGCAGTGTGGTAGGCGATATCGGCAACGGTAACTGGCAGGGTGGAGTCGTGCCCCTGAACCGTCATGCCCAGCGAATCGCCCACCAGCATGACGTTAAGCCCTTCATCAGCAAAGAGTTTGGCGAAGCTATAGTCATAAGCGGTGATGGTCGCGAAACGTTTTTTTTCCTGTTTGTACTTCTGCAGTAAGGAGATGGTGGTCGGTTTc                                                                <  1:361383/228‑1 (MQ=255)
                                                                                                                                                                  gTGCGCCGCGACGGACGGCGGCAGTGTGGGAGGCGATATCGGCAACGGTAACTGGCAGGGTGGAGTCGTGCCCCTGAACCGTCATGCCCAGCGAATCGCCCACCAGCATGACGTTAAGCCCTTCATCAGCAAAGAGTTTGGCGAAGCTATAGTCATAAGCGGTGATGGTCGCGAAACGTTTTTTTTCCTGTTTGTACTTCTGCAGTAAGGAGATGGTGGTCGGTTTc                                                                <  2:690298/227‑1 (MQ=255)
                                                                                                                                                                    gcgcCGCGACGTACGGCGGCAGTGTGGTAGGCGATATCGGCAACGGTACCTGGCAGGGTGGAGTCGTGCCCCTGAACCGTCATGCCCAGCGAATCGCCCACCAGCATGACGTTAAGCCCCTCATCAGCCAACGGTTTGGACAAGCAATAGTCATACGCGGGGGTGGTCGCCAAACGTTTTTTTTCCTGATTTTACTTCTGCAGTAAGGAGACGGTGGCAGGTTTCTTa                                                             >  1:285278/1‑228 (MQ=255)
                                                                                                                                                                    gcgcCGCGACGTACGGCGGCAGTGTGGTAGGCGATATCGGCAACGGTAACTGGCAGGGTGGAGTCGTGCCCCGGAACCGTCATGCCCAGCGAATCGCCCACCAGCATGACGTTAAGCCCTTCATCAGCAAAGAGTTTGGCGAAGCTATAGTCATAAGCGGTGATGGTCGCGAAACGTTTTTTTTCCTGTTTGTACTTCTGCAGTAAGGAGATGGTGGTCGGTTTCATa                                                             <  1:448172/228‑1 (MQ=255)
                                                                                                                                                                      cccggGACGTCAGGCGGCAGTGGGGTAGGCGATAGCGGAAACGGTAACTGGCAGGGTGGAGTCGTGCCCCTGAACCGTCATGCCCAGCGAATCGCCCACCAGCATGACGTTAAGCCCTTCATCAGCAAAGAGTTTGGCGAAGCTATAGTCATAAGCGGTGATGGTCGCGAAACGTTTTTTTTCCTGTTTGTACTTCTGCAGTAAGGAGATGGTGGTCGGTTTCATaa                                                            <  2:674804/226‑1 (MQ=255)
                                                                                                                                                                        cgcgACGTACGGCGGCAGTGTGGTAGGCGATATCGGCAACGGTAACTGGCAGGGTGGGGTCGTGCCCCTGAACCGTCATGCCCAGCGAATCGCCCACCAGCATGACGTTAAGCCCTTCATCAGCAAAGAGTTTGGCGAAGCt                                                                                                                                                >  2:466266/1‑142 (MQ=255)
                                                                                                                                                                        cgcgACGTACGGCGGCAGTGTGGTAGGCGATATCGGCAACGGTAACTGGCAGGGTGGAGTCGTGCCCCTGAACCGTCATGCCCAGCGAATCGCCCACCAGCATGACGTTAAGCCCTTCATCAGCAAAGAGTTTGGCGAAGCTATAGTCATAAGCGGTGATGGTCGCGAAACGTTTTTTTTACTGTTTGTAATTCTGCAGTAAGGGGATGGTGGTAGGTTTCATAACGt                                                         >  2:222922/1‑228 (MQ=255)
                                                                                                                                                                         gcgACGTACGGCGGCAGTGTGGTAGGCGATATCGGCAACGGTAACTGGCAGGGTGGGGTCGTGCCCCTGAACCGTCATGCCCAGCGAATCGCCCACCAGCATGACGTTAAGCCCTTCATCAGCAAAGAGTTTGGCGAAGCta                                                                                                                                               <  1:466266/142‑1 (MQ=255)
                                                                                                                                                                           gacggtaggcggcAGTGTGGTAGGCGATACTGGAAACGGTAACTGGCAGGGTGGAGTAGTTCCCCTGAACCGTCATGCCCAGAGAATCGCCCACCAGCATGACGTTAAGCCCTTCATCAGCAAAGAGTTTGGCGAAGCTATAGTCATAAGCGGTGATGGTCGCGAAACGTTTTTTTTCCTGTTTGTACTTCTGCAGTAAGGAGATGGTGGTCGGTTTCATAACGTATc                                                      <  2:683381/221‑1 (MQ=255)
                                                                                                                                                                                    gggcAGTGGGGTAGGCGATATCGGCAACGGTAACTGGCAGGGTGGAGTCGTGCCCCTGAACCGTCATGCCCAGCGAATCGCCCACCAGCATGACGTTAAGCCCTTCATCAGCAAAGAGTTTGGCGAAGCTATAGTCATAAGCGGTGATGGTCGCGAAACGTTTTTTTTCCTGTTTGTACTTCTGCAGTAAGGAGATGGTGGTCGGTTTCATAACGTATCCTGATAAAt                                             <  2:85231/227‑1 (MQ=255)
                                                                                                                                                                                         gtgtGGTAGGCGATATCGGCAACGGTAACTGGCAGGGTGGAGTCGTGCCCCTGAACCGTCATGCCCAGCGAATCGCCCACCAGCATGACGTTAAGCCCTTCATCAGCAAAGAGTTTGGCGAAGCTCTAGTCATAAGCGGTGATGGTCGCGAAACGTTTTTTTTCCTGTTTGTACTTCTGCAGTAAGGAGATGGTGGTCGGTTTCATAACGTATCCTGTTAAATTg                                           >  2:695802/1‑225 (MQ=255)
                                                                                                                                                                                               tAGGCGATATCGGCAACGGTAACTGGCAGGGTGGAGTCGTGCCCCTGAACCGTCATGCCCAGCGAATCGCCCACCAGCATGACGTTAAGCCCTTCATCAGCAAAGAGTTTGGCGAAGCTATAGTCATAAGCGGTGATGGTCGCGAAACGTTTTTTTTCCTGTTTGTACTTCTGCAGTAAGGAGATGGTGGTCGGTTTCATAACGTATCCTGATAAATTGATGTtgtg                                   <  2:506397/227‑1 (MQ=255)
                                                                                                                                                                                               tAGGCGAGATCGGCAACGGTAACTGGCAGGGTGGAGTCGTGCTCCTGAACCGTCATGCCCAGCGAATCGCCCACCAGCATGACGTTAAGCCCTTCATCAGCAAAGAGTTTGGCGAAGCTATATTCATAAGCGGTGATGGTCGCGAAACGTTTTTTTTCCTGTTTGTACTTCTGCAGTAAGGAGATGGTGGTCGGTTTCATAACGTATCCTGATAAATTGATGTtgtg                                   <  2:178565/227‑1 (MQ=255)
                                                                                                                                                                                                 ggCGATATCGGCAACGGTAACTGGCAGGGTGGAGTGGTGCCCCTGAACCGTCATGCCCAGAGAATCGCCCACCAGCATGACGTTAAGCCCTTCATCCGCAAAGAGTTTGGCGAAGCTATAGTCATAAGCGGTGATGTTCGCGCAACGTTTTTTTTCCGGTTTGTACTTCTTAGTCAAGGGAGTGGTGGGCGGTTTCATAACGTATCCGTATaaaattgaatgtgagga                                >  1:336328/1‑214 (MQ=255)
                                                                                                                                                                                                 ggCGATATCGGCAACGGTAACTGGCAGGGTGGAGTCGTGCCCCTGAACCGTCATGCCCAGCGAATCGCCCACCAGCATGACGTTAAGCCCTTCATCAGCAAAGAGTTTGGCGAAGCTATAGTCATAAGCGGTGATGGTCGCGAAACGTTTTTTTTCCTGTTTGTACTTCTGCAGTAAGGAGATGGTGGTCGGTTTCATAACGTATCCTGATAAATTGAtgttgt                                    >  2:106410/1‑224 (MQ=255)
                                                                                                                                                                                                  gCGATATCGGCAACGGTAACTGGCAGGGTGGAGTCGTGCCCCTGAACCGTCATGCCCAGCGAATCGCCCACCAGCATGACGTTAAGCCCTTCATCAGCAAAGAGTTTGGCGAAGCTATAGTCATAAGCGGTGATGGTCGCGAAACGTTTTTTTTCCTGTTTGTACTTCTGCAGTAAGGAGATGGTGGTCGGTTTCATAACGTATCCTGATAAATTGATGTtgtg                                   <  1:106410/224‑1 (MQ=255)
                                                                                                                                                                                                      tatCGACAACGGTAACTGGCAGGGTGGAGTCGTGCCCCTGAACCGTAATGCCCAGCGAATAGCCCACCAGCATGACGTTAAGCCCTTCATCAGCAAAGAGTTTGGCGAAGCTATAGTCATAAGCGGTGATGGCcgc                                                                                                                        >  2:268082/1‑136 (MQ=255)
                                                                                                                                                                                                      gagcGGCAACGGTAACTGGCAGGGGGGAGTGGTGCCCATGAACCGTCATGCCCAGCGAATCGCCCACCAGCATGACGTTAAGCCATTCATCAGCAAAGAGTTTGGCGAAGCTATAGTCATAAGCGGTGATGGTCGCGAAACGTTTTTTTTCCTGTTTGTACTTCTGCAGTAAGGAGATGGTGGTCGGTTTCATAACGTATCCTGATAAATTGATGTTGTGCTGTc                              <  1:376127/222‑1 (MQ=255)
                                                                                                                                                                                                               cGGTTACTGGCAGGGGGGGGTCGTGCCCATGAACCGTCAGGCCCAGCGAATCGCCCACCAGCATGACGTTAAGCCCTTCAGCAGCAAAGAGTTTGGCGAAGCTATAGTCATAAGCGGTGATGGTCGCGAAACGTTTTTTTTCCTGTTTGTACTTCTGCAGTAAGGAGATGGTGGTCGGTTTCATAACGTATCCTGATAAATTGATGTTGTGCTGTCTGGCATTTTAt                   <  2:105364/227‑1 (MQ=255)
                                                                                                                                                                                                                ggTAACTGGGCGGGTGGGGTCGGGCCCCTGAACCGTCATGGCCAGCGAAGTCCCCACCAGCATGACGTTAAGCCCTTCATCAGCAAAGAGTTTGGCGAAGCTATAGTCATAAGCGGTGATGGTCGCGAAACGTTTTTTTTCCTGTTTGTACTTCTGCAGTAAGGAGATGGTGGTCGGTTTCATAACGTATCCTGATAAATTGATGTTGTGCTGTCTGGCATTTTATCa                 <  1:42112/228‑1 (MQ=255)
                                                                                                                                                                                                                          aGGGTGGAGTCGTGCCCCTGAACCGTCATGCCCAGCGAATCGCCCACCAGCATGACGTTAAGCCCTTCATCAGCAAAGAGTTTGGCGAAGCTATAGTCATAAGCGGTGATGGTCGCGAAACGTTTTTTTTCCTGTTTGTACTTCTGCAGTAAGGAGATGGTGGTCGGTTTCATAACGTATCCTGATAAATTGATGTTGTGCTGTCTGGCATTTTATc                  >  2:42112/1‑217 (MQ=255)
                                                                                                                                                                                                                            ggTGGAGTCGTGCCCCTGAACCGTCATGCCCAGCGAATCGCCCACCAGCATGACGTTAAGCCCTTCATCAGCAAACAGTTTGGCGAAGCTATAGTCATAAGCGGTGATGGTCGCGAAAAGTTTTTTTTCCTGTTTGTACTTCTGCAGTAAGGAGATGGTGGTCGGTTTCATAACGTATCCTGATAAATTGATGTTGTGCGGTCTGGCATTTTATCAGTCACATtggtg     >  2:255463/1‑228 (MQ=255)
                                                                                                                                                                                                                              aggAGTCGTGCCCCTGAACCGTCATGCCCAGCGAATCGCCCACCAGCATGACGTTAAGCCCTTCCTCCGCAAAGAGTTTGGCACACCTATAGTCATAAGCGGTGATGGTCGCGCAACGTTTTTTTTCCTGt                                                                                                     >  1:60297/2‑131 (MQ=255)
                                                                                                                                                                                                                                gAGTCGTGCCCCTGAACCGTCATGCCCAGCGAATCGCCCACCAGCATGACGTTAAGCCCTTCATCAGCAAAGAGTTTGGCGAAGCTATAGTCATAAGCGGTGATGGTCGCGAAACGTTTTTTTTCCTGTTTGTACTTCTGCAGTAAGGAGATGGTGGTCGGTTTCATAACGTATCCTGATAAATTGATGTTGTGCTGTCTGGCATTTTATCAGTCACATTGGTgggg  >  2:359634/1‑227 (MQ=255)
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                  
GGCACGTTCGGTCAGCATTTGTACGGTTTCTACCAGCCACTCACCGCCTTCAATTTTGACCATGTTAGCACCGGCACGCATAACCGTTGCGGCGTTTTCGAAGGCTTGTTCCGGCGTGGCATACGCCATAAACGGCAGGTCAGCCAGCAGCAGGCAGTTTGGTGCGCCGCGACGTACGGCGGCAGTGTGGTAGGCGATATCGGCAACGGTAACTGGCAGGGTGGAGTCGTGCCCCTGAACCGTCATGCCCAGCGAATCGCCCACCAGCATGACGTTAAGCCCTTCATCAGCAAAGAGTTTGGCGAAGCTATAGTCATAAGCGGTGATGGTCGCGAAACGTTTTTTTTCCTGTTTGTACTTCTGCAGTAAGGAGATGGTGGTCGGTTTCATAACGTATCCTGATAAATTGATGTTGTGCTGTCTGGCATTTTATCAGTCACATTGGTGGGG  >  NC_000913/149212‑149661

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: