Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,551,915 A→C 23.5% V186G (GTG→GGG)  malP ← maltodextrin phosphorylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,551,9150AC23.5% 76.0 / 4.1 52V186G (GTG→GGG) malPmaltodextrin phosphorylase
Reads supporting (aligned to +/- strand):  ref base A (24/14);  new base C (5/9);  total (29/23)
Fisher's exact test for biased strand distribution p-value = 1.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCAGTATGGTTGTCGTTTGGATAGAGAACTTTGGTCAGTTTTTCCGCATTGATGCCCTGCTGTTCGGCACGCAAGAAATCACCGTCGTTAAATTTAGTCAGATCAAACGGATGCGCGTGCGTCGCCTGCCACAGACGCAGCGGCTGCGCCACGCCATTACGATAGCCGACAACGGGGAG‑ATCCCACGCTTGACCGGTAATGGTAAACTCCGGCTCCCAGCGTCCGTCTTTCGTCACTTTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGATGCCAGTCATCCGGCGC‑TTCAACCTGTTTGCCATCG‑ACAAAAGACTGGCGGAACAAACCATATTGATAGTTCAGACCGTAACCCGTCGCAGACTGACCGACAGTTGCCATTGAGTCGAGGAAGCACGCCGCCAGACGTCCCAGACCACCGTTACCCAGCG  >  NC_000913/3551681‑3552143
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                      
gCAGTATGGTTGTCGTTTGGATAGAGAACTTTGGTCAGTTTTTCCGCATTGATGCCCTGCTGTTCGGCACGCAAGAAATCACCGTCGTTAAATTTAGTCAGATCAAACGGATGCGCGTGCGTCGCCTGCCACAGACGCAGCGGCTGCGCCACGCCATTACGATAGCCGCCAACGGGGAG‑ATCCCACGCTTGACCGGTAAATGTAAACCCCGGCTCCCAGCGTCCGTCTTTCGGCCCttt                                                                                                                                                                                                                                    >  1:120346/1‑239 (MQ=255)
gCAGTATGGTTGTCGTTTGGATAGAGAACTTTGGTCAGTTTTTCCGCATTGATGCCCTGCTGTTCGGCACGCAAGAAATCACCGTCGTTAAATTTAGTCAGATCAAACGGATGCGCGTGCGTCGCCTGCCACAGACGCAGCGGCTGCGCCACGCCATTACGATAGCCGACAACGGGGGG‑ATCCCACGCTTGACCGGTAATGGTAAACTCCGGCTCCCAGCGTCCGCCTTTCGTCACttt                                                                                                                                                                                                                                    >  2:42828/1‑239 (MQ=255)
gCAGTATGGTTGTCGTTTGGATAGAGAACTTTGGTCAGTTTTTCCGCATTGATGCCCTGCTGTTCGGCACGCAAGAAATCACCGTCGTTAAATTTAGTCAGATCAAACGGATGCGCGTGCGTCGCCTGCCACAGACGCAGCGGCTGCGCCACGCCATTACGATAGCCGACAACGGGGAG‑ATCCCACGCTTGACCGGTAATGGTAAACTCCGGCTCCCAGCGTCCGTCTTTCGTCCCttt                                                                                                                                                                                                                                    >  2:737597/1‑239 (MQ=255)
                 tGGATAGAGAACTTTGGTCAGTTTTTCCGCATTGATGCCCTGCTGTTCGGCACGCAAGAAATCACCGTCGTTAAATTTAGTCAGATCAAACGGATGCGCGTGCGTCGCCTGCCACAGACGCAGCGGCTGCGCCACGCCATTACGATAGCCGACAACGGGGAG‑ATCCCACGCTTGACCGGGAATGGTAAACTCCGGCCCCCAGCGCCCGTCTTTCGTCACTTTACCGCCAAACCCAAccc                                                                                                                                                                                                                   >  1:394958/1‑238 (MQ=255)
                        agaACTTTGGTCAGTTTTTCCGCATTGATGCCCTGCTGTTCGGCACGCAAGAAATCACCGTCGTTAAATTTAGTCAGATCAAACGGATGCGCGTGCGTCGCCTGCCACAGACGCAGCGGCTGCGCCACGCCATTACGATAGCCGACAACGGGGAG‑ATCCCACGCTTGACCGGTAATGGTAAACTCCGGCTCCCAGCGTCCGTCTTTCGTCACTTTACCGCCAATCCCTACCTGCACAt                                                                                                                                                                                                             >  2:811019/1‑238 (MQ=255)
                        agaACTTTGGTCAGTTTTTCCGCATTGATGCCCTGCTGTTCGGCACGCAAGAAATCACCGTCGTTAAATTTAGTCAGATCAAACGAATGCGCGTGCGTCGCCTGCCACAGACGCAGCGGCTGCGCCACGCCATTACGATAGCCGACAACGGGGAG‑ATCCCACGCTTGACCGGTAATGGTAAACTCCGGCTCCCAGCTTCCGTCTTTCGTCACTTTACCCCCAATCCCTACCTGCACAt                                                                                                                                                                                                             >  2:607030/1‑238 (MQ=255)
                                  tCAGTTTTTCCGCATTGATGCCCTGCTGTTCGGCACGCAAGAAATCACCGTCGTTAAATTTAGTCAGATCAAACGGATGCGCGTGCGTCGCCTGCCACAGACGCAGCGGCTGCGCCACGCCATTACGATAGCCGACAACGGGGAG‑ATCCCACGCTTGACCGGTAATGGTAAACTCCGGCTCCCAGCGTCCGTCTTTCGTCACTTTACCGCCAATCCCTACCTGCACATCCaagccttcg                                                                                                                                                                                                  >  2:270965/1‑231 (MQ=255)
                                    aGTTTTTCCGCATTGATGCCCTGCTGTTCGGCACGCAAGAAATCACCGTCGTTAAATTTAGTCAGATCAAACGGATGCGCGTGCGTCGCCTGCCACAGACGCAGCGGCTGCGCCACGCCATTACGATAGCCGACAACGGGGAG‑ATCCCACGCTTGACCGGTAATGGTAAACTCCGGCTCCCAGCGTCCGTCTTTCGTCACTTTACCGCCAATCCCTACCTGCCCATCCAGTGCTTCGtt                                                                                                                                                                                                >  2:657840/1‑239 (MQ=255)
                                                         tgctgTTCGGCACGCAAGAAATCACCGTCGTTAAATTTAGTCAGATCAAACGGATGCGCGTGCGTCGCCTGCCACAGACGCAGCGGCTGCGCCACGCCATTACGATAGCCGACAACGGGGAG‑ATCCCACGCTTGACCGGTAATGGTAAACTCCGGCTCCCAGCGTCCGTCTTTCGTCACTTTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCCCGGGTAGtt                                                                                                                                                                           >  2:497727/1‑239 (MQ=255)
                                                                   cACGCAAGAAATCACCGTCGTTAAATTTAGTCAGATCAAACGGATGCGCGTGCGTCGCCTGCCACAGACGCAGCGGCTGCGCCACGCCATTACGATAGCCGACAACGGGGAG‑ATCCCACGCTTGACCGGTAATGGTAAACTCCGGCTCCCAGCGTCCGTCTTTCGTCCCTTTACCGCCAATCCCTACCTGCACATCCAGTGCTGCGTTGTTGCGGaaaccagggggattaaacgggat                                                                                                                                                                  >  1:779716/1‑217 (MQ=255)
                                                                   cACGCAAGAAATCACCGTCGTTAAATTTAGTCAGATCAAACGGATGCGCGTGCGTCGCCTGCCACAGACGCAGCGGCTGCGCCACGCCATTACGATAGCCGACAACGGGGAG‑ATCCCACGCTTGACCGGTAATGGTAAACTCCGGCTCCCAGCGTCCGCCTTTCGTCACTTTACCCCCAATCCCTACCTGCACATCCATGGCTTCGTTGTGGCGGAACCACGGGCAGCTatctcgctg                                                                                                                                                                  >  2:222505/1‑230 (MQ=255)
                                                                        aaGAAATCACCGTCGTTAAATTTAGTCAGATCAAACGGATGCGCGTGCGTCGCCTGCCACAGACGCAGCGGCTGCGCCACGCCATTACGATAGCCGACAACGGGGAG‑ATCCCACGCTTGACCGGTAATGGTAAACTCCGGCTCCCAGCGTCCGTCTTTCGTCACTTTACCGCCAATCCCTACCTGCCCATCCCGTGCTTCGTTGTGGCGGAACAACGGGAGGTTCATGCGATGCCatac                                                                                                                                                            >  2:156689/1‑236 (MQ=255)
                                                                                      gtTAAATTTAGTCAGATCAAACGGATGCGCGTGCGTCGCCTGCCACAGACGCAGCGGCTGCGCCACGCCATTACGATAGCCGACAACGGGGAG‑ATCCCACGCTTGACCGGTAATGGTAAACTCCGGCTCCCAGCGTCCGTCTTTCGTCACTTTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACCGGTAGTTACTCGGATGCCACACATACGGCCG‑TTCa                                                                                                                                              >  2:479395/1‑238 (MQ=255)
                                                                                               aGTCAGATCAAACGGATGCGCGTGCGTCGCCTGCCACAGACGCAGCGGCTGCGCCACGCCATTACGATAGCCGACAACGGGGAG‑ATCCCACGCTTGACCGGTAATGGTAAACTCCGGCTCCCAGCGTCCGTCTTTCGTCACTTTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGCCGGAACCACGGGTAGTTACTGCGATGCCAGTCATCCGGCGC‑TTCACCCTGTTTg                                                                                                                                     >  2:170689/1‑238 (MQ=255)
                                                                                                gccAGATAAAACGGATGGCGGGGGGTCGCCGGCCAAAGGCGACGCGGCTGCGCCCCGCCATTACGATAGCCGACAACGGGGAG‑ATCCCACGCTTGACCGGTAATGGTAAACTCCGGCTCCCAGCGTCCGTCTTTCGTCACTTTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGATGCCAGTCATCCGGCGC‑TTCAACCTGTTTGcc                                                                                                                                   <  2:394958/237‑1 (MQ=255)
                                                                                                               tGCGCGTGCGTCGCCTGCCACAGACGCAGCGGCTGCGCCACGCCATTACGATAGCCGACAACGGGGAG‑ATCCCACGCTTGACCGGTAATGGTAAACTCCGGCTCCCAGCGTCCGTCTTTCGTCCCTGTACCGCCAATACCACCCTGCCCATCCAGTGCTTCGTTGTGCCGGAACCACGGTGATTTACAGCGATGCCCG‑CAACCGGCGCGTTCAACCCGTTTGcccaccccaaaaaaacg                                                                                                                    >  1:621118/1‑224 (MQ=255)
                                                                                                               tGCGCGTGCGTCGCCTGCCACAGACGCAGCGGCTGCGCCACGCCATTACGATAGCCGACAACGGGGAG‑ATCCCACGCTTGACCGGTAATGGTAAACTCCGGCTCCCAGCGTCCGTCTTTCGTCACTTTACCGCCAATCCCTACCTGCCCATCCAGTGCTTCGTTGTGGGGGAACCCCGGGTAGTTACTGCGATCCCCGTCATCCCGCGC‑TTCAACCCGTTTGCCATCC‑AAAAAAAAccg                                                                                                                   >  1:43007/1‑237 (MQ=255)
                                                                                                               tGCGCGTGCGTCGCCTGCCACAGACGCAGCGGCTGCGCCACGCCATTACGATAGCCGACAACGGGGAG‑ATCCCACGCTTGACCGGTAATGGTAAACTCCGGCTCCCAGCGTCCGTCTTTCGTCACTTTACCGCCAATCCCTACCTGCACTTCCAGGGCTTCGTTGTGGCGGAACCCAGGGTCGTTATTTCGGTTCCCTTCCTCCCGCGC‑TCTAACCCGTTGTCCGTCG‑ACAAAAAAcgc                                                                                                                   >  1:340779/1‑237 (MQ=255)
                                                                                                               tGCGCGTGCGTCGCCTGCCACAGACGCAGCGGCTGCGCCACGCCATTACGATAGCCGACAACGGGGAG‑ATCCCACGCTTGACCGGTAATGGTAAACTCCGGCTCCCAGCGTCCGTCTTTCGTCACTTTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGATGCCAGCCATCCGGCCC‑TTCAAACCGTTTGACCTCGCACaaaaacgg                                                                                                                    >  1:716271/1‑235 (MQ=255)
                                                                                                               tGCGCGTGCGTCGCCTGCCACAGACGCAGCGGCTGCGCCACGCCATTACGATAGCCGACAACGGGGAG‑ATCCCACGCTTGACCGGTAATGGTAAACTCCGGCTCCCAGCGTCCGCCTTTCGTCCCTTTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCCAGTCCAGTCATCCGGCAC‑TTCAACCCTTTTGcccatcaaacaaacagt                                                                                                                    >  1:218685/1‑224 (MQ=255)
                                                                                                                 gcgttccgcccccTGCAACAGACGCGGGGCGTGCGCCCCGCCTTTCGGAGAGCCGACAGCGGGGGG‑ATCCCACCCTTGCCCGGAAATGGAAACCTGCGGCCCCACGCCGCCGTCGTGCTTCACTTTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGATGCCAGTCATCCGGCGC‑TTCAACCTGTTTGCCATCG‑ACAAAAGACTgg                                                                                                                  <  1:156689/227‑1 (MQ=255)
                                                                                                                         tccccTGCCACAGCCCGCGCGGCTGCGCACCGCCATTCCGAAAGCCGACAACGGGGAG‑ATCCCACGTTTGACCGGTAATGGTAAACTCCGGCTCCCAGCGTCCGTCTTTCGTCCCTTTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGATGCCAGTCATCCGGCGC‑TTCAACCTGTTTGCCATCG‑ACAAAAGACTGGCGGAACaa                                                                                                          <  2:301263/235‑1 (MQ=255)
                                                                                                                          cgccccccccagcccggggcgCTGCCCCCCGCCTTTCCGAGCGCCGACAGCGGGGAG‑CCCCCCCGGTTCCCCGATAAGGGAAACCTCCGGCTCCCACCCTCCGTCTTTCGTCACTTGACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGATGCCAGTCATCCGGCGC‑TTCAACCTGTTTGCCATCG‑ACAAAAGACTGGCGGAACAAAc                                                                                                        <  1:42828/219‑1 (MQ=255)
                                                                                                                          ccctctccccagcggccggggCTGCGCCACGCCATCACTATGGCCGGCGTGGGGCAG‑CCCCCTCGTTTAGTCGGTTATGACACTCCCCGGCTCCCGACGCCCGGTTTTCTTCCGTTTGCCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGATGCCAGTCATCCGGCGC‑TTCAACCTGTTTGCCATCG‑ACAAAAGACTGGCGGAACAAAc                                                                                                        <  1:782699/220‑1 (MQ=255)
                                                                                                                               taaaggggacaggcgccgcgcccacgccaTTACGATAGCCGACAACGGGGAA‑TTCCCCCGCTTGCCCGTTAATGGTAAACTCCGGCTCCCAGCTTCCGTCTTTCGTCACTTTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGATGCCAGTCATCCGGCGC‑TTCAACCTGTTTGCCATCG‑ACAAAAGACTGGCGGAACAAACCatat                                                                                                   <  2:838477/218‑1 (MQ=255)
                                                                                                                                ccACAGACGCAGCGGCTACGCCACGCCATTACGATAGCCGACAACGGGGCG‑ATCCCACGCTTGACCGGTAATGGTAAACTCCCGCTCCCAGCGTCCGGCTTTCGTCACTTTACCGCCAATACCTACCTGCACATCCAGCGCTTCGCTGTGGCGGAACCACAGGCGGATACCGCACGTCCCGTCATTCCGCGC‑TTCACACTGGTTGCCCACG‑CCAAAAGAAAGGGGGGACAAACCacaaa                                                                                                  >  1:693631/1‑235 (MQ=255)
                                                                                                                                          gcgggtggcccccccctatttggATACCCGACACCGGGGAG‑ATCCCCCGCTGGACCGGTAAGTGTAACCCCCCGCTCCAAGGCTCCGTTTTTGCTCACTTTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGATGCCAGTCATCCGGCGC‑TTCAACCTGTTTGCCATCG‑ACAAAAGACTGGCGGAACAAACCATATTGATAGTTCa                                                                                         <  1:737597/216‑1 (MQ=255)
                                                                                                                                          aGCGGGCGCACCACGCCATTACGATAGCCGACAACGGGGAG‑ATCCCCCGCTTGACCGGTAATGGTAAACTCCGGCTCCCAGCGCCCGTCTTTCGTCACTTTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGATGCCAGTCATCCGGCGC‑TTCAACCTGTTTGCCATCG‑ACAAAAGACTGGCGGAACAAACCATATTGATAGTTCa                                                                                         <  2:576652/238‑1 (MQ=255)
                                                                                                                                          aGCGGCTGCGCCACGCCATTACGATAGCCGACAACGGGGAG‑ATCCCACGCTTGACCGGTAATGGTAAACTCCGGCTCCCAGCGTCCGTCTTTCGTCACTTTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCCCGGGGAGTTACTGCGATGCCAGTCATCCGCGGC‑TTCAACCTGTTTGCCATCG‑CAAAAAAACGGGCGGGAAAAAcacaaatgtaagattaa                                                                                        >  1:753503/1‑223 (MQ=255)
                                                                                                                                                                    gCCGACAACGGGGAG‑ATCCCACGCTTGACCGGTAATGGTAAACTCCGGCTCCCAGCGTCCGTCTTTCGTCACTTTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGATGCCAGTCATCCGGCGC‑TTCAACCTGTTTTCCATCG‑CCAAAAGACTGGCGGAACAAACCATATTTTTAGTTCCAACCCTAACCCCTCGCCGACTGaacca                                                              >  2:567557/1‑235 (MQ=255)
                                                                                                                                                                    gCCGACAACGGGGAG‑ATCCCACGCTTGACCGGTAATGGTAAACTCCGGCTCCCAGCGTCCGTCTTTCGTCACTTTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGATGCCAGTCATCCGGCGC‑TTCAACCTGTTTGCCATCG‑ACAAAAGACTGGCGGAACAAACCCTATTGATAGTTCAGACCGTAAcccccgccagcatggccgc                                                              >  2:53409/1‑223 (MQ=255)
                                                                                                                                                                    gCCGACAACGGGGAG‑ATCCCACGCTTGACCGGTAATGGTAAACTCCGGCTCCCAGCGTCCGTCTTTCGTCACTTTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGATGCCAGTCATCCGGCGC‑TTCAACCTGTTTGCCATCG‑ACAAAAGACTGGCGGAACAAACCATATTGCTAATTCAGACCGTAACCCGTCGCAGACTggacgc                                                              >  2:628929/1‑234 (MQ=255)
                                                                                                                                                                     ccgcaccgggggAA‑CTCCCCCGCTTACCCTGTTATGGAAACCCCCGGCCCCCAGGTCCCGGTTTTCCTCCCTTGACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGATGCCAGTCATCCGGCGC‑TTCAACCTGTTTGCCATCG‑ACAAAAGACTGGCGGAACAAACCATATTGATAGTTCAGACCGTAACCCGTCGCAGACTGACCGAc                                                             <  1:657840/231‑1 (MQ=255)
                                                                                                                                                                         acaaggggAGAATCCCACGCTTGCCCGGTAAGGGTAAACCCCGGCCCCCAGCTTCCGTCTTTCGTCCCTTTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGATGCCAGTCATCCGGCGC‑TTCAACCTGTTTGCCATCG‑ACAAAAGACTGGCGGAACAAACCATATTGATAGTTCAGACCGTAACCCGTCGCAGACTGACCGACAGt                                                          <  2:753503/235‑1 (MQ=255)
                                                                                                                                                                          aaCGGGGAG‑ATCCCACGCTTGACCGGTAATGGTAAACTCCGGCTCCCAGCGTCCGTCTTTCGTCACTTTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGGAGTTACTGCGATGCCAGTCATCCGGCGC‑TTCAACCTGTTTGCCATCG‑ACAAAAGACTGGCGGAACAAACCATATTGATTGTTCAAAACGTAACCCCTCCCAGAGTAACCGaaaatt                                                         >  1:577118/1‑233 (MQ=255)
                                                                                                                                                                               ggAC‑ATCCCACGCTTGGCCGGTAATGGTAACCCCCGGCTCCAAGCGTCCGTTTTTGCTCACTTTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGATGCCAGTCATCCGGCGC‑TTCAACCTGTTTGCCATCG‑ACAAAAGACTGGCGGAACAAACCATATTGATAGTTCAGACCGTAACCCGTCGCAGACTGACCGACAGTTGCCAt                                                    <  2:751553/238‑1 (MQ=255)
                                                                                                                                                                                 ag‑aTCCCACGTTTGACCGGTAATGGTAAACCCCGGCTCCCAGCGTCCGTCTTTCGTCCTTTTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGATGCCAGTCATCCGGCGC‑TTCAACCTGTTTGCCATCG‑ACAAAAGACTGGCGGAACAAACCATATTGATAGTTCAGACCGTAACCCGTCGCAGACTGACCGACAGTTGCCATTg                                                  <  2:120346/238‑1 (MQ=255)
                                                                                                                                                                                              tgccccgaaaTGGTAACCCCCGGCTCCCGGCGTCCTTCTTTCGTCACTTTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGATGCCAGTCATCCGGCGC‑TTCAACCTGTTTGCCATCG‑ACAAAAGACTGGCGGAACAAACCATATTGATAGTTCAGACCGTAACCCGTCGCAGACTGACCGACAGTTGCCATTGAGTCGAGGAAGCa                                     <  2:253684/231‑1 (MQ=255)
                                                                                                                                                                                              tgcccGGTAAGGGTAACCTCCGGCCCCAAGCGTCCTTTTTTCGTCACTTTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGATGCCAGTCATCCGGCGC‑TTCAACCTGTTTGCCATCG‑ACAAAAGACTGGCGGAACAAACCATATTGATAGTTCAGACCGTAACCCGTCGCAGACTGACCGACAGTTGCCATTGAGTCGAGGAAGCa                                     <  2:43007/236‑1 (MQ=255)
                                                                                                                                                                                              ggACCGTTAATGGAAAACTCGGGCCCCCGGGTGCCGTTTTTCTCCACTTCACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGATGCCAGTCATCCGGCGC‑TTCAACCTGTTTGCCATCG‑ACAAAAGACTGGCGGAACAAACCATATTGATAGTTCAGACCGTAACCCGTCGCAGACTGACCGACAGTTGCCATTGAGTCGAGGAAGCa                                     <  2:175304/238‑1 (MQ=255)
                                                                                                                                                                                                aCCGGTAATGGTAAACTCCGGCTCCCAGCGTCCGTCTTTCGTCACTTTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGATTCCAGTCCTCCGGCGC‑TTCAACCTGTTTGCCATCC‑ACAAAAGACTGGCCGCACAAACCATATTGTTAGTTCCTCCCGTCACCCCCCCCACACTACCCCACAAGTGCCAtttttcccgggaaacac                                    >  1:504633/1‑223 (MQ=255)
                                                                                                                                                                                                    gTAATGGTAAACTCCGGCTCCCAGCGTCCGTCTTTCGTCACTTTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGATGCCAGTCATCCGGCGC‑TTCAACCTGTTTGCCATCG‑ACCAAAGACTGGCGGAACCAACCCTATTGATAGTTCAGACCGCAACCCGTCGCCGACTGACCGACAGTTGCCCTTGAGTCGCGGAAGGAcgccg                                >  2:37927/1‑238 (MQ=255)
                                                                                                                                                                                                        tgtaaacccccGGCCCCCAGCCTCCGTCTTTCGTCCCTTTACCGCCAATCCCTACCTGCACATCCATTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGATGCCAGTCATCCGGCGC‑TTCAACCTGTTTGCCATCG‑ACAAAAGACTGGCGGAACAAACCATATTGATAGTTCAGACCGTAACCCGTCGCAGACTGACCGACAGTTGCCATTGAGTCGAGGAAGCACGCCGCCAGa                           <  2:556877/230‑1 (MQ=255)
                                                                                                                                                                                                        tGGTAAACTCCGGCCCCCAGCGTCCGTCTTTCGTCACTTTCCCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGATGCCAGTCATCCGGCGC‑TTCAACCTGTTTGCCATCG‑ACAAAAGACTGGCGGAACAAACCATATTGATAGTTCAGACCGTAACCCGTCGCAGACTGACCGACAGTTGCCATTGAGTCGAGGAAGCACGCCGCCAGa                           <  2:101090/239‑1 (MQ=255)
                                                                                                                                                                                                         ggTAAACTCCGGCTCCCAGCGTCCGTCTTTCGTCACTTTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGATGCCCGTCATCCGGCGC‑TTCAACCTGTTTGCCATCG‑ACAAAAGACTGGCGGAACAAACCat                                                                                                     >  1:251949/1‑164 (MQ=255)
                                                                                                                                                                                                         ggTAAACTCCGGCTCCCAGCGTCCGTCTTTCGTCACTTTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGATGCCAGTCATCCGGCGC‑TTCAACCTGTTTGCCATCG‑ACAAAAGACTGGCGGAACAAACCat                                                                                                     <  2:251949/164‑1 (MQ=255)
                                                                                                                                                                                                          gaaaaCTCCGCCCCCCGGTGTCCTTCTTTCGTCACTTTCCCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGATGCCAGTCATCCGGCGC‑TTCAACCTGTTTGCCATCG‑ACAAAAGACTGGCGGAACAAACCATATTGATAGTTCAGACCGTAACCCGTCGCAGACTGACCGACAGTTGCCATTGAGTCGAGGAAGCACGCCGCCAGa                           <  2:801855/235‑1 (MQ=255)
                                                                                                                                                                                                                  cGGCTCCCAGCGTCCGTCTTTCGTCACTTTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGATGCCAGTCATCCGGCGC‑TTCAACCGGTTTGCCATCG‑ACAAAAGACTGGCGGAACAAACCATATTGATAGTTCAGACCGCAAACCGACGCACACTGACCGACAGTTGCCATTTAGTCGAGGAAGAAAGCCCCCAGACGGCCCAGa                  >  1:178592/1‑238 (MQ=255)
                                                                                                                                                                                                                     cccccggggtccatcTTGCTACACTTGACCGCCAACCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGATGCCAGTCATCCGGCGC‑TTCAACCTGTTTGCCATCG‑ACAAAAGACTGGCGGAACAAACCATATTGATAGTTCAGACCGTAACCCGTCGCAGACTGACCGACAGTTGCCATTGAGTCGAGGAAGCACGCCGCCAGACGTCCCAGacca               <  2:456438/230‑1 (MQ=255)
                                                                                                                                                                                                                       cccAGCTTCCTGTTTTCGTCCCTTTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGATGCCAGTCATCCGGCGC‑TTCAACCTGTTTGCCATCG‑ACAAAAGACTGGCGGAACAAACCATATTGATAGTTCAGACCGTAACCCGTCGCAGACTGACCGACAGTTGCCATTGAGTCGAGGAAGCACGCCGCCAGACGTCCCAGACCACCg            <  2:577118/239‑1 (MQ=255)
                                                                                                                                                                                                                             ctccgtcTTTCGTCCCTTTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGATGCCAGTCATCCGGCGC‑TTCAACCTGTTTGCCATCG‑ACAAAAGACTGGCGGAACAAACCATATTGATAGTTCAGACCGTAACCCGTCGCAGACTGACCGACAGTTGCCATTGAGTCGAGGAAGCACGCCGCCAGACGTCCCAGACCACCGTTAcc       <  2:716271/237‑1 (MQ=255)
                                                                                                                                                                                                                                 gtcTTTCGTCCCTTACCCGCCAATCCCTACCTGCTCATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGATGCCAGTCATCCGGCGC‑TTCAACCTGTTTGCCATCG‑ACAAAAGACTGGCGGAACAAACCATATTGATAGTTCAGACCGTAACCCGTCGCAGACTGACCGACAGTTGCCATTGAGTCGAGGAAGCACGCCGCCAGACGTCCCAGACCACCGTTACCCAgcg  <  2:58510/239‑1 (MQ=255)
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                      
GCAGTATGGTTGTCGTTTGGATAGAGAACTTTGGTCAGTTTTTCCGCATTGATGCCCTGCTGTTCGGCACGCAAGAAATCACCGTCGTTAAATTTAGTCAGATCAAACGGATGCGCGTGCGTCGCCTGCCACAGACGCAGCGGCTGCGCCACGCCATTACGATAGCCGACAACGGGGAG‑ATCCCACGCTTGACCGGTAATGGTAAACTCCGGCTCCCAGCGTCCGTCTTTCGTCACTTTACCGCCAATCCCTACCTGCACATCCAGTGCTTCGTTGTGGCGGAACCACGGGTAGTTACTGCGATGCCAGTCATCCGGCGC‑TTCAACCTGTTTGCCATCG‑ACAAAAGACTGGCGGAACAAACCATATTGATAGTTCAGACCGTAACCCGTCGCAGACTGACCGACAGTTGCCATTGAGTCGAGGAAGCACGCCGCCAGACGTCCCAGACCACCGTTACCCAGCG  >  NC_000913/3551681‑3552143

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: ‑