Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,509,880 A→C 14.2% D84A (GAC→GCC)  yfeO → putative ion channel protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,509,8800AC14.2% 143.4 / 4.1 65D84A (GAC→GCC) yfeOputative ion channel protein
Reads supporting (aligned to +/- strand):  ref base A (46/7);  new base C (2/10);  total (48/17)
Fisher's exact test for biased strand distribution p-value = 5.60e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGTTGTTATTATCGCTCCCCGCCGTGGCAATTGGGATTGCGTCCAGTCTTATTCTGATTGTGGTGATGAAAATCGCCTCGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATT‑CACCCCTCTGGATCATCGGTGTATTAAC‑GCTAACGGGTATTGCGGTGGGGTT‑GGTTATCCGTTTCAGCCAGGGTCATGCC‑GGACCAGACC‑CCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGACCGTCA‑ATATCGCCCTTGCGGTTGCCATTGGC‑GCTCGTCTGTTACCGCGCGTCAACCGAATGGAGTGGACTATTTTAGCCTCTGCCGGAACCATCGGCGCACTGTTTGG  >  NC_000913/2509655‑2510111
                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                          
tgttgtTATTATCGCTCCCCGCCGTGGCAATTGGGATTGCGTCCAGTCTTATTCTGATTGTGGTGATGAAAATCGCCTCGGTATTACAGACTTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATT‑CACCCCTCTGGATCATCGGGTTTTTAAC‑GCAAACGGGTATTGCGGTGGGGTT‑GGTTATCCGTTTCAGCCAGGGTCATGCA‑GGACCAGCAC‑CCGCCTGTGa                                                                                                                                                                                                                               >  1:354332/1‑238 (MQ=255)
tgttgtTATTATCGCTCCCCGCCGTGGCAATTGGGATTGCGTCCAGTCTTATTCTGATTGTGGTGATGAAAATCGCCTCGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATT‑CACCCCTCTGGATCATCGGTGTATTAAC‑GCTAACGGGCATTGCGGTGGGGTT‑GGTTATCCGTTTCAGCCAGGGTCATGCC‑GGACCAGACC‑CCGCCTGTGa                                                                                                                                                                                                                               >  1:377090/1‑238 (MQ=255)
tgttgtTATTATCGCTCCCCGCCGTGGCAATTGGGATTGCGTCCAGTCTTATTCTGATTGTGGTGATGAAAATCGCCTCGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTATGGGGATAGCCCAGGATT‑CACCCCACCGGATCATCGGTGTATTAAC‑GCTAACGGGAATTGCGGGGGGTTT‑GGTTTACCGTTCTAGCAAGGGCAATGCC‑AGACCAACCC‑CCGCCTGGGa                                                                                                                                                                                                                               >  2:276679/1‑238 (MQ=255)
     ttattaTCGCTCCCCGCCGTGGCAATTGGGATTGCGTCCAGTCTTATTCTGATTGTGGTGATGAAAATCGCCTCGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATT‑CACCCCTCTGGATCATCGGTGTATTAAC‑GCTAACGGGTATTGCGGTGGGGTT‑GGTTATCCGTTTCAGCCAGGGTCATGCC‑GGACCAGACC‑CCGCCTGTGAAcc                                                                                                                                                                                                                            >  1:860410/1‑236 (MQ=255)
       attatCGCTCCCCGCCGTGGCAATTGGGATTGCGTCCAGTCTTATTCTGATTGTGGTGATGAAAATCGCCTCGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATT‑CACCCCTCTGCATCATCGGTGTATTAAC‑GCTAACGGGTATTGCGGCGGGGTT‑GGTTATCCGTTTCAGCCAGGGGCATGCT‑GCCCCCGACC‑CCGCCGTGGAACCGCTg                                                                                                                                                                                                                        >  1:833805/1‑238 (MQ=255)
               tCCCCGCCGTGGCAATTGGGATTGCGTCCAGTCTTATTCTGATTGTGGTGATGAAAATCGCCTCGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATT‑CACCCCTCTGGATCATCGGTGTATTAAC‑GCTAACGGGTATTGCGGTGGGGTT‑GGTTATCCGTTTCAGCCAGGGTCTTGCC‑GGACCAGACC‑CCGCCTGTGAACCGTTGATCGgcgc                                                                                                                                                                                                                >  1:151668/1‑238 (MQ=255)
                     ccgTGGCAATTGGGATTGCGTCCAGTCTTATTCTGATTGTGGTGATGAAAATCGCCTCGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATT‑CACCCCTCTGGATCATCGGTGTATTAAC‑GCTAACGGGTATTGCGGTGGGGTT‑GGTTATCCGTTTAAGCCAGGGTCATGCC‑GGACCAGACC‑CCGCCTGGGAACCGCCGATCGGCGCCCCGGt                                                                                                                                                                                                          >  2:866728/1‑238 (MQ=255)
                      cgTGGCAATTGGGATTGCGTCCAGTCTTATTCTGATTGTGGTGATGAAAATCGCCTCGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATT‑CACCCCTCTGGATCATCGGGGTATTAAC‑GCTAACGGGTATTGCGGTGGGGTT‑GGTTATCCGTTTCAGCCAGGGTCATGCC‑GGGCCAGACC‑CCGCCTGTGAACCCCCGACCGGCGAACCGGtt                                                                                                                                                                                                         >  1:613433/1‑238 (MQ=255)
                         ggCAATTGGGATTGCGTCCAGTCTTATTCTGATTGTGGTGATGAAAATCGCCTCGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATT‑CAACCCTCTGGATCATCGGGAGATTAAC‑GCTAACGGGTATTCTGGTGGGGCT‑GGTTATCCGTTTCAGCAAGGGACCTTCG‑GGCTCAGACC‑CCCCCTGGGGACAGTGGAACGGCGCAcgcggtacg                                                                                                                                                                                                      >  2:691489/1‑230 (MQ=255)
                            aaTTGGGATGTTGCCCCGTTTTGTTTTGTTTGTGGTGAACAAAATCCCCTGGGTATTACAGAATTTGCTCTGGCAACGACCGCCGGGACCTCTGGGGATAGCCCAGGATT‑CCCCCCGCTGGATCCTCGGTGTATTAAC‑GCTACCGGGTATTGCGTTGGGTTT‑GGTTATCCGTTTCAGCCAGGGTCATGCC‑GGCCCAGCCC‑CCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCct                                                                                                                                                                                                  <  2:799033/239‑1 (MQ=255)
                                           cAGTCTTATTCTGATTGTGGTGATGAAAATCGCCTCGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATT‑CACCCCTCTGGATCATCGGTGTATTAAC‑GCTAACGGGTATTGCGGTGGGGTT‑GGTTATCCGTTTCAGCCAGGGTCATGCC‑GGCCCAGACC‑CCGCCTGTGAACCGCGGACCGGCGCACCGGTTTCGCCCTCTGCGCTACCTgg                                                                                                                                                                                     >  1:487390/1‑237 (MQ=255)
                                                        aTTGTGGTGATGAAAATCGCCTCGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATT‑CACCCCTCTGGATCATCGGTGTATTAAC‑GCTAACGGGTATTGCGGTGGGGTT‑GGTTATCCGTTTCAGCCAGGGTCATGCC‑GGACCAGACC‑CCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCCtt                                                                                                                                                                      >  1:890954/1‑239 (MQ=255)
                                                        aTTGTGGTGATGAAAATCGCCTCGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATT‑CACCCCTCTGGATCATCGGTGTATTAAC‑GCTAACGGGTATTGCGGTGGGGTT‑GGTTATCCGTTTCAGCCAGGGTCATGCC‑GGACCAGACC‑CCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCAtt                                                                                                                                                                      >  1:556668/1‑239 (MQ=255)
                                                        aTTGTGGTGATGAAAATCGCCTCGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATT‑CACCCCTCTGGATCATCGGTGTATTAAC‑GCTAACGGGTATTGCGGTGGGGTT‑GGTTATCCGTTTCAGCCAGGGTCATGCC‑GGACCAGACC‑CCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTATGCGCTACCTACACTTATCTTAGCGtt                                                                                                                                                                      >  1:552572/1‑239 (MQ=255)
                                                        aTTGTGGTGATGAAAATCGCCTCGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATT‑CACCCCTCTGGATCATCGGTGTATTAAC‑GCTAACGGGTATTGCGGTGGGGTT‑GGTTATCCGTTTCAGCCAGGGTCATGCC‑GGACCAGACC‑CCGCCTGGGAACCGCTGATCGGCGCACCGGTTCCGCCCCCGTCGCCACCTGGACTTTTGGGAGCAtt                                                                                                                                                                      >  2:695095/1‑239 (MQ=255)
                                                        aTTGTGGTGATGAAAATCGCCTCGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATT‑CACCCCTCTGGATCATCGGTGTATTAAC‑GCTAACGGGTATTGCGGTGGGGTT‑GGTTATCCGTTTCAGCCAGGGGCATGCC‑GGACCAGACC‑CCGCCGGTGAACCGCTGATCGGCGCACCGGTTCCGCCCCCTGCGCCACCTGGCCTTATCGGACCAtt                                                                                                                                                                      >  2:478089/1‑239 (MQ=255)
                                                             ggtgATGAAAATCGCCTCGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATT‑CACCCCTCTGGATCATCGGTGTATTAAC‑GCTAACGGGTATTGCGGTGGGGTT‑GGTTATCCGTTTCAGCCAGGGTCATGCC‑GGACCAGACC‑CCGCCTGTGAACCGCTGATCGGCGCCCCGGTTCCGCCCTCTCCCCTACCTGGACTTATCGTAGCATTaaa                                                                                                                                                                   >  2:460169/1‑236 (MQ=255)
                                                              gtgATGAAAATCGCCTCGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATT‑CACCCCTCTGGATCATCGGTGTATTAAC‑GCTAACGGGTATTGCGGTGGGGTT‑GGTTATCCGTTTCAGCCAGGGTCATGCC‑GGACCAGACC‑CGACCTGTGAACCGCTGACCGGGGCCCCGGGTCCGCCCTCCGCCCCACCTGGGCTTTACGTacacttatttt                                                                                                                                                                 >  2:27546/1‑228 (MQ=255)
                                                                gatgaAAATCGCCTCGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATT‑CACCCCTCTGGATCATCGGTGTATTAAC‑GCTAACGGGTATTGCGGTGGGGTT‑GGTTATCCGTTTCAGCCAGGGTCATGCC‑GGACCAGACC‑CCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTTTCGTAGCATTAATTCTCg                                                                                                                                                              >  1:224810/1‑239 (MQ=255)
                                                                                         aaTTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATT‑CACCCCTCTGGATCATCGGTGTATTAAC‑GCTAACGGGTATTGCGGTGGGGTT‑GGTTATCCGTTTCAGCCAGGGTCATGCC‑GGACCAGACC‑CCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTCCCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTggg                                                                                                                                     >  1:797454/1‑239 (MQ=255)
                                                                                           tttGCTCTGGCAACGACTGCCGGGAGCTCTGGGGATAGCCCAGGATT‑CACCCCTCTGGATCATCGGTGTATTAAC‑GCTAACGGGAATTGCGGTGGGGTT‑GGTTATCCGTTTCAGCCAGGGTCATGCC‑GGACCAGACC‑CCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTCCCTGGGCTTATCGTACACTTAATTTCCGGTCTTGCTGCGGGCGCCAcccggtgcc                                                                                                                                   >  2:401714/1‑233 (MQ=255)
                                                                                           tttGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATT‑CACCCCTCTGGATCATCGGTGTATTAAC‑GCTAACGGGTATTGCGGTGGGGTT‑GGTTATCCGTTTCAGCCAGGGTCATGCC‑GGACCAGACC‑CTGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCCTTAATTCTCGGCCTTCCGTGCGGCGTACCCCTgggcc                                                                                                                                   >  2:536187/1‑237 (MQ=255)
                                                                                           tttGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATT‑CACCCCTCTGGATCATCGGTGTATTAAC‑GCTAACGGGTATTGCGGTGGGGTT‑GGTTATCCGTTTCAGCCAGGGTCATGCC‑GGACCAGACC‑CCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGcc                                                                                                                                                                                                    >  1:668708/1‑174 (MQ=255)
                                                                                           tttGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATT‑CACCCCTCTGGATCATCGGTGTATTAAC‑GCTAACGGGTATTGCGGTGGGGTT‑GGTTATCCGTTTCAGCCAGGGTCATGCC‑GGACCAGACC‑CCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTgggg                                                                                                                                    >  1:247245/1‑238 (MQ=255)
                                                                                           tttGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATT‑CACCCCTCTGGATCATAGGTGTATTAAC‑GCTAACGGGTATTGCGGTGGGGTT‑GGTTATCCGTTTCAGCCAGGGTCATGCC‑GGACCCGACC‑CCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCACTCTGCGCTACCTGGCCTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTAGGCCTtggg                                                                                                                                    >  1:438080/1‑238 (MQ=255)
                                                                                           tttGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGACT‑CACCCCTCTGGATCATCGGTGTATTAAC‑GCTAACGGGTATTGCGGTGGGGTT‑GGTTATCCGTTTCAGCCAGGGTCATGCC‑GGACCAGACC‑CCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCACCCTCTGCGCTACCTGGACTTATCGGAGCCTTACTTCTCGGTCTTTCTGCGGGCGTAAGCCTGGGGc                                                                                                                                   >  2:313245/1‑239 (MQ=255)
                                                                                           tttGCTCTGGAAACGACTGCAGGGAACTCTGGGGATAGCCCAGGATT‑CACCCCTCTGGATCATCGGTGTATTAAC‑GCTAAGGGGTATTGCGGGGGGGTT‑GGTTATCCGTTTCAGCCAGGGTCATGCC‑GGACCAGACC‑CCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGc                                                                                                                                   <  1:478089/239‑1 (MQ=255)
                                                                                              gCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATT‑CACCCCTCTGGATCATCGGTGTATTAAC‑GCTAACGGGTATTGCGGTGGGGTT‑GGTTATCCGTTTCAGCCAGGGTCATGCC‑GGACCAGACC‑CCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATGAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCgg                                                                                                                                >  1:156488/1‑239 (MQ=255)
                                                                                                   ggCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATT‑CACCCCTCTGGATCATCGGTGTATTAAC‑GCTAACGGGTATTGCGGTGGGGTT‑GGTTATCCGTTTCAGCCAGGGTCATGCC‑GGACCAGACC‑CCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCTTCTGCGCTACCTGGACTTATCGTAg                                                                                                                                                                          >  1:60805/1‑192 (MQ=255)
                                                                                                   ggCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATT‑CACCCCTCTGGATCATCGGTGTATTAAC‑GCTAACGGGTATTGCGGTGGGGTT‑GGTTATCCGTTTCAGCCAGGGTCATGCC‑GGACCAGACC‑CCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCTTAGCATTAATTTTCGGTCTTGCTGTCGGCGTCAGCCTGTGCCCGGacaa                                                                                                                            >  1:524838/1‑235 (MQ=255)
                                                                                                   ggCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATT‑CACCCCTCTGGATCATCGGTGTATTAAC‑GCTAACGGGTATTGCGGTGGGGTT‑GGTTATCCGTTTCAGCCAGGGTCATGCC‑GGACCAGACC‑CCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACa                                                                                                                            >  1:10561/1‑238 (MQ=255)
                                                                                                    gcaaggcggtccGGGAACTGTGGGGATACCCCGGGTTC‑CCCCCCTTGGGCTCAGCGGT‑TATCAACGGTTAAGGGGTTGGGGGGGGGGGTT‑GGTTAACCGTCTCGGCAGGGGTCATGCC‑GGCCCAGCCC‑CCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACAt                                                                                                                           <  2:151668/228‑1 (MQ=255)
                                                                                                            ttccGGGACCTTTGGGG‑TAGCCCCGGGTTCCCCCCCTCTGGTTCATGGTTGTATAAAC‑CTAAAGGGGTTTGGGGGTGGGGTT‑GGTTATCCGTTCCAGCCGGGGTCATGCC‑GGCCCAGCCC‑CCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATg                                                                                                                  <  2:10561/237‑1 (MQ=255)
                                                                                                                                gCCCAGGATT‑CACCCCTCTGGATCATCGGTGTATTAAC‑GCTAACGGGTATTGCGGTGGGGTT‑GGTTATCCGTTTCAGCCAGGGTCATGCC‑GGACCAGACC‑CCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCTTCAGCCTGGGGCCGGAACATCCGATCATGACCGTCA‑ATATTGCCCtctg                                                                                             >  1:342087/1‑236 (MQ=255)
                                                                                                                                 cccAGGATT‑CACCCCTCTGGATCATCGGTGTATTAAC‑GCTAACGGGTATTGCGGTGGGGTT‑GGTTATCCGTTTCAGCCAGGGTCATGCC‑GGACCAGACC‑CCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTTCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGACCGTCA‑ATATCGCCCTTGCg                                                                                            >  1:408116/1‑239 (MQ=255)
                                                                                                                                    aGGATT‑CACCCCTCTGGATCATCGGTGTATTAAC‑GCTAACGGGTATTGCGGTGGGGTT‑GGTTATCCGTTTCAGCCAGGGTCATGCC‑GGACCAGACC‑CCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCCCGGCCTTGCTGGCGGCGTCAGCCTGCGGCCGGAACATCCGTTCATGACCGTCA‑ATATCGCCCTTGCGGtt                                                                                         >  1:481285/1‑239 (MQ=255)
                                                                                                                                            aCCCCTCTGGATCATCGGTGTATTAAC‑GCTAACGGGTATTGCGGTGGGGTT‑GGTTATCCGTTTCAGCCAGGGTCATGCC‑GGACCAGACC‑CCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCTTGACCTCCA‑ATATCGCCCTTGCGGTTGcacttg                                                                                  >  1:367785/1‑234 (MQ=255)
                                                                                                                                            aCCCCTCTGGATCATCGGTGTATTAAC‑GCTAACGGGTATTGCGGTGGGGTT‑GGTTATCCGTTTCAGCCAGGGTCATGCC‑GGACCAGACC‑CCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGACCGTCA‑ATATCGCCCTTGCGGTTGCCATTg                                                                                  >  1:143763/1‑239 (MQ=255)
                                                                                                                                            aCCCCTCTGGATCATCGGTGTATTAAC‑GCTAACGGGTATTGCGGTGGGGTT‑GGTTATCCGTTTCAGCCAGGGTCATGCC‑GGACCAGACC‑CCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAACATTAATTCTCGGTCTTGCCGGCGGCCTCAGCCTGGGGCCGGAACCTCCGGCCAGAACCGTCACATGT‑GGCCTTGCGGTTGGCATTg                                                                                  >  2:84922/1‑239 (MQ=255)
                                                                                                                                                  ctGGATCATGGGTGTATTAGC‑GCTAACGGGTATTGGGGTGGGTTT‑GGTTATCCGTTTCAGCCGGGGTAATGCG‑GGCCCAGACC‑CCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGACCGTCA‑ATATCGCCCTTGCGGTTGCCATTGGC‑GCTc                                                                           <  1:474340/239‑1 (MQ=255)
                                                                                                                                                          tcGGTGTATTAAC‑GCTAACGGGTATTGCGGTGGGGTT‑GGTTATCCGTTTCAGCCAGGGTCATGCC‑GGACCAGACC‑CCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGCTCCTGACCGTCA‑ATATCGCCCTTGCGGTTGCCATTGGC‑GCCCGCCTGtt                                                                    >  2:790481/1‑238 (MQ=255)
                                                                                                                                                          tcGGTGTATTAAC‑GCTAACGGGTATTGCGGTGGGGTT‑GGTTATCCGTTTCAGCCAGGGTCATGCC‑GGACCAGACC‑CCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTAGGGCCGGAACATCCGATCATGACCGGCA‑ATATCGCCCTTGCGGTTGCCATTGGC‑GCTCGTCTGtt                                                                    >  1:552726/1‑238 (MQ=255)
                                                                                                                                                               gtATTAAC‑GCTAACGGGTATTGTGGTGGGGTT‑GGTTATCCGTTTCAGCCAGGGTAATGCC‑GGACCAGACC‑CCGCATGGGAACCGCTGATCGGCGCATCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGGAACATTCATTCTCGTGCTTGCTTGCGGCGCCAGGCGCGGGCCGGAACATCCGGTCATGACCGCAA‑CTATGGCCCCTGCTGTTGCATATGGC‑GCTCGTCTGTTAcccgc                                                              >  2:718568/1‑236 (MQ=255)
                                                                                                                                                               gtATTAAC‑GCTAACGGGTATTGCGGTGGGGTT‑GGTTATCCGTTTCAGCCAGGGTCATGCC‑GGACCAGACC‑CCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGACCATCCGATCATGACCGTCA‑ATATCGCCCTTGCGGTTGCCATTTGC‑GCTCGTCTGTTACcgc                                                               >  1:366533/1‑238 (MQ=255)
                                                                                                                                                                     aC‑GCTAACGGGTATTGCGGTGGGGTT‑GGTTATCCGTTTCAGCCAGGGTCATGCC‑GGACCAGACC‑CCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGACCGGCA‑ATATCGCCCTTGCGGTTGCCATTGGC‑GCTCGTCTGTTACCGCCCGTCa                                                         >  1:395044/1‑238 (MQ=255)
                                                                                                                                                                           aaCGGGTATTGCGGTGGGGTT‑GGTTATCCGTTTCAGCCAGGGTCATGCC‑GGACCAGACC‑CCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGACCGTCA‑ATATCGCCCTTGCGGTTGc                                                                                       >  1:466861/1‑204 (MQ=255)
                                                                                                                                                                             cGGGTATTGCGGTGGGGTT‑GGTTATCCGTTTCAGCCAGGGTCATGCC‑GGACCAGACC‑CCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGACCGTCA‑ATATCGCCCTTGCGGTTGCCATTGGC‑GCTCGTCTGTTACCGCGCGTCAACCGAATg                                                 >  1:449902/1‑239 (MQ=255)
                                                                                                                                                                                ttATGGCGGGGGGGTG‑TGTTTCCCGTTTCTGCC‑GGGTCCTTCCGGGCCCAGCCC‑CCTCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGACCGTCA‑ATATCGCCCTTGCGGTTGCCATTGGC‑GCTCGTCTGTTACCGCGCGTCAACCGAATGGa                                               <  1:105557/237‑1 (MQ=255)
                                                                                                                                                                                gttttGGGGTGGGGTT‑TGTTTCCCTTTTCGGCCGGGGTCCTGCC‑CGACCAGCCC‑CCGCCTGTGAACCGCTGATCGGCGCACCGGTTTCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCCGCCTGGGGCCGGAACATCCGATCATGACCGTCA‑ATATCGCCCTTGCGGTTGCCATTGGC‑GCTCGTCTGTTACCGCGCGTCAACCGAATGGa                                               <  2:889275/235‑1 (MQ=255)
                                                                                                                                                                                    tggggTGGGGTT‑TGTTATCCGTTCACGCCAGGTTCAGTCC‑GGCCCAGACC‑CCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGACCGTCA‑ATATCGCCCTTGCGGTTGCCATTGGC‑GCTCGTCTGTTACCGCGCGTCAACCGAATGGAGTGGa                                          <  1:313245/236‑1 (MQ=255)
                                                                                                                                                                                    ggggggggggGT‑GTTTACCGGTTCCAGCCAGGGTCAGGCC‑GGACCAGACC‑CCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGACCGTCA‑ATATCGCCCTTGCGGTTGCCATTGGC‑GCTCGTCTGTTACCGCGCGTCAACCGAATGGAGTGGa                                          <  1:679222/233‑1 (MQ=255)
                                                                                                                                                                                     gggggggggTTGGGTTTCCCGTTTCAGCCGGGTTCGGGCC‑GCACCAGACC‑CCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGACCGTCA‑ATATCGCCCTTGCGGTTGCCATTGGC‑GCTCGTCTGTTACCGCGCGTCAACCGAATGGAGTGGa                                          <  1:401714/234‑1 (MQ=255)
                                                                                                                                                                                           gggTT‑GGTTATCCGTTTCAGCCAGGGTCATGCC‑GGACCAGACC‑CCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGGTCATGACCGTCA‑ATATCGCCCTTGCGGTTGCCTTTGTG‑GCCCGCCTTTTACCCCCCGTCAACCGAATGGGGGGGAAAAtttt                                   >  2:511829/1‑239 (MQ=255)
                                                                                                                                                                                           gggTT‑GGTTATCCGTTTCAGCCAGGGTCATGCC‑GGACCAGACC‑CCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGACCGTCA‑AAATCGCCCTTGCGGTTGCCCTTGAG‑CCTCGCCTGTTCCCGCGCGTCACCCACATGGAGTgggaaatttt                                   >  2:839146/1‑231 (MQ=255)
                                                                                                                                                                                                           ttCAGCCAGGGTCATGCC‑GGACCAGACC‑CCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGACCGTCA‑ATATCGCCCTTGCGGTTGCCATTGGC‑GCTCGTCTGTTACCGCGCGTCAACCGAATGGAGTGGCCTATTTTAGCCTCTGCCGGaa                     >  1:636645/1‑238 (MQ=255)
                                                                                                                                                                                                                   gggTCATGCC‑GGACCAGCCC‑CCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGACCGTCA‑ATATCGCCCTTGCGGTTGCCATTGGC‑GCTCGTCTGTTACCGCGCGTCAACCGAATGGAGTGGACTATTTTAGCCTCTGCCGGAACCATCGgcg            <  1:536187/239‑1 (MQ=255)
                                                                                                                                                                                                                    ggTCATGCC‑GGACCAGACC‑CCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGCGGCCGGATAATTCGCTCCTGCACGTCA‑CTGTCGCCCTTGCGGGTGCCATTGGCGGCACTTGT‑TTACCCGTCTGCAAACGCCGGGAGGGGACCATTTTTGCCTCTTCCGGAACCAACggtga           >  2:263442/1‑236 (MQ=255)
                                                                                                                                                                                                                      tgaTGCC‑GGCCGAGCCC‑CCGCCTGGGACCCGCTGATCGGCCCACCGGGTCCGCCCTCTGCGCTCCCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGACCGTCA‑ATATCGCCCTTGCGGTTGCCATTGGC‑GCTCGTCTGTTACCGCGCGTCAACCGAATGGAGTGGACTATTTTAGCCTCTGCCGGAACCATCGGCGCa          <  2:524838/236‑1 (MQ=255)
                                                                                                                                                                                                                      tccggCC‑GGACCAGCCC‑CCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGACCGTCA‑ATATCGCCCTTGCGGTTGCCATTGGC‑GCTCGTCTGTTACCGCGCGTCAACCGAATGGAGTGGACTATTTTAGCCTCTGCCGGAACCATCGGCGCa          <  1:27546/234‑1 (MQ=255)
                                                                                                                                                                                                                      tcagtcc‑GGACC‑GACCGCCGCCTGTGAACCGCTGATCGGCGCTCCGGTTCCGCCCTCTGCGCTACCTGGTCTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGACCGTCA‑ATATCGCCCTTGCGGTTGCCATTGGC‑GCTCGTCTGTTACCGCGCGTCAACCGAATGGAGTGGACTATTTTAGCCTCTGCCGGAACCATCGGCGCa          <  1:26623/233‑1 (MQ=255)
                                                                                                                                                                                                                      ccagcgc‑GGACCAGCCC‑CCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGACCGTCA‑ATATCGCCCTTGCGGTTGCCATTGGC‑GCTCGTCTGTTACCGCGCGTCAACCGAATGGAGTGGACTATTTTAGCCTCTGCCGGAACCATCGGCGCa          <  1:596035/232‑1 (MQ=255)
                                                                                                                                                                                                                         tGCC‑GGACCAGACC‑CCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGACCGTCA‑ATATCGCCCTTGCGGTTGCCCTTGGC‑GCTCGTCTGTTACCGCGCGTCACCCGAATGGAGTGGACTATTTTAGCCTCTGCCGGaa                     >  1:395883/1‑224 (MQ=255)
                                                                                                                                                                                                                              ggcccAGCCC‑CCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGACCGTCA‑ATATCGCCCTTGCGGTTGCCATTGGC‑GCTCGTCTGTTACCGCGCGTCAACCGAATGGAGTGGACTATTTTAGCCTCTGCCGGAACCATCGGCGCACTGTTTgg  <  1:460169/236‑1 (MQ=255)
                                                                                                                                                                                                                               gACCAGACC‑CCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGACCGTCA‑ATATCGCCCTTGCGGTTGCCATTGGC‑GCTCGTCTGTTACCGCGCGTCa                                                         <  1:663067/183‑1 (MQ=255)
                                                                                                                                                                                                                               gACCAGACC‑CCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGACCGTCA‑ATATCGCCCTTGCGGTTGCCATTGGC‑GCTCGTCTGTTACCGCGCGTCa                                                         >  2:663067/1‑183 (MQ=255)
                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                          
TGTTGTTATTATCGCTCCCCGCCGTGGCAATTGGGATTGCGTCCAGTCTTATTCTGATTGTGGTGATGAAAATCGCCTCGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATT‑CACCCCTCTGGATCATCGGTGTATTAAC‑GCTAACGGGTATTGCGGTGGGGTT‑GGTTATCCGTTTCAGCCAGGGTCATGCC‑GGACCAGACC‑CCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGACCGTCA‑ATATCGCCCTTGCGGTTGCCATTGGC‑GCTCGTCTGTTACCGCGCGTCAACCGAATGGAGTGGACTATTTTAGCCTCTGCCGGAACCATCGGCGCACTGTTTGG  >  NC_000913/2509655‑2510111

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