Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,509,878 A→C 15.1% 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,8780AC15.1% 151.1 / 4.3 73D84A (GAC→GCC) yfeOputative ion channel protein
Reads supporting (aligned to +/- strand):  ref base A (50/7);  new base C (2/14);  total (52/21)
Fisher's exact test for biased strand distribution p-value = 2.98e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGTTGTTATTATCGCTCCCCGCCGTGGCAATTGGGATTGCGTCCAGTCTTATTCTGATTGTGGTGATGAAAATCGCCTCGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAGCCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGACCGTCAATATCGCCCTTGCGGTTGCCATTGGCGCTCGTCTGTTACCGCGCGTCAACCGAATGGAGTGGACTATTTTAGCCTCTGCCGGAACCATCGGCGCACTGTTTGG  >  NC_000913/2509653‑2510109
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                          
tgttgtTATTATCGCTCCCCGCCGTGGCAATTGGGATTGCGTCCAGTCTTATTCTGATTGTGGTGATGAAAATCGCCTCGGTATTACAGACTTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCTCTGGATCATCGGGTTTTTAACGCAAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAGCCAGGGTCATGCAGGACCAGCACCCGCCTGTGa                                                                                                                                                                                                                               >  2:354332/1‑238 (MQ=255)
tgttgtTATTATCGCTCCCCGCCGTGGCAATTGGGATTGCGTCCAGTCTTATTCTGATTGTGGTGATGAAAATCGCCTCGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCTCTGGATCATCGGTGTATTAACGCTAACGGGCATTGCGGTGGGGTTGGTTATCCGTTTCAGCCAGGGTCATGCCGGACCAGACCCCGCCTGTGa                                                                                                                                                                                                                               >  2:377090/1‑238 (MQ=255)
tgttgtTATTATCGCTCCCCGCCGTGGCAATTGGGATTGCGTCCAGTCTTATTCTGATTGTGGTGATGAAAATCGCCTCGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTATGGGGATAGCCCAGGATTCACCCCACCGGATCATCGGTGTATTAACGCTAACGGGAATTGCGGGGGGTTTGGTTTACCGTTCTAGCAAGGGCAATGCCAGACCAACCCCCGCCTGGGa                                                                                                                                                                                                                               >  1:276679/1‑238 (MQ=255)
     ttattaTCGCTCCCCGCCGTGGCAATTGGGATTGCGTCCAGTCTTATTCTGATTGTGGTGATGAAAATCGCCTCGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAGCCAGGGTCATGCCGGACCAGACCCCGCCTGTGAAcc                                                                                                                                                                                                                            >  2:860410/1‑236 (MQ=255)
       attatCGCTCCCCGCCGTGGCAATTGGGATTGCGTCCAGTCTTATTCTGATTGTGGTGATGAAAATCGCCTCGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCTCTGCATCATCGGTGTATTAACGCTAACGGGTATTGCGGCGGGGTTGGTTATCCGTTTCAGCCAGGGGCATGCTGCCCCCGACCCCGCCGTGGAACCGCTg                                                                                                                                                                                                                        >  2:833805/1‑238 (MQ=255)
               tCCCCGCCGTGGCAATTGGGATTGCGTCCAGTCTTATTCTGATTGTGGTGATGAAAATCGCCTCGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAGCCAGGGTCTTGCCGGACCAGACCCCGCCTGTGAACCGTTGATCGgcgc                                                                                                                                                                                                                >  2:151668/1‑238 (MQ=255)
                     ccgTGGCAATTGGGATTGCGTCCAGTCTTATTCTGATTGTGGTGATGAAAATCGCCTCGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTAAGCCAGGGTCATGCCGGACCAGACCCCGCCTGGGAACCGCCGATCGGCGCCCCGGt                                                                                                                                                                                                          >  1:866728/1‑238 (MQ=255)
                      cgTGGCAATTGGGATTGCGTCCAGTCTTATTCTGATTGTGGTGATGAAAATCGCCTCGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCTCTGGATCATCGGGGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAGCCAGGGTCATGCCGGGCCAGACCCCGCCTGTGAACCCCCGACCGGCGAACCGGtt                                                                                                                                                                                                         >  2:613433/1‑238 (MQ=255)
                         ggCAATTGGGATTGCGTCCAGTCTTATTCTGATTGTGGTGATGAAAATCGCCTCGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCAACCCTCTGGATCATCGGGAGATTAACGCTAACGGGTATTCTGGTGGGGCTGGTTATCCGTTTCAGCAAGGGACCTTCGGGCTCAGACCCCCCCTGGGGACAGTGGAACGGCGCAcgcggtacg                                                                                                                                                                                                      >  1:691489/1‑230 (MQ=255)
                            aaTTGGGATGTTGCCCCGTTTTGTTTTGTTTGTGGTGAACAAAATCCCCTGGGTATTACAGAATTTGCTCTGGCAACGACCGCCGGGACCTCTGGGGATAGCCCAGGATTCCCCCCGCTGGATCCTCGGTGTATTAACGCTACCGGGTATTGCGTTGGGTTTGGTTATCCGTTTCAGCCAGGGTCATGCCGGCCCAGCCCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCct                                                                                                                                                                                                  <  1:799033/239‑1 (MQ=255)
                                           cAGTCTTATTCTGATTGTGGTGATGAAAATCGCCTCGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAGCCAGGGTCATGCCGGCCCAGACCCCGCCTGTGAACCGCGGACCGGCGCACCGGTTTCGCCCTCTGCGCTACCTgg                                                                                                                                                                                     >  2:487390/1‑237 (MQ=255)
                                                        aTTGTGGTGATGAAAATCGCCTCGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAGCCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCCtt                                                                                                                                                                      >  2:890954/1‑239 (MQ=255)
                                                        aTTGTGGTGATGAAAATCGCCTCGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAGCCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCAtt                                                                                                                                                                      >  2:556668/1‑239 (MQ=255)
                                                        aTTGTGGTGATGAAAATCGCCTCGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAGCCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTATGCGCTACCTACACTTATCTTAGCGtt                                                                                                                                                                      >  2:552572/1‑239 (MQ=255)
                                                        aTTGTGGTGATGAAAATCGCCTCGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAGCCAGGGTCATGCCGGACCAGACCCCGCCTGGGAACCGCTGATCGGGGCCCCGTGTCCGCCCTCCGCGCCCCCTGGAccttacggaacatt                                                                                                                                                                      >  1:474340/1‑226 (MQ=255)
                                                        aTTGTGGTGATGAAAATCGCCTCGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAGCCAGGGTCATGCCGGACCAGACCCCGCCTGGGAACCGCTGATCGGCGCACCGGTTCCGCCCCCGTCGCCACCTGGACTTTTGGGAGCAtt                                                                                                                                                                      >  1:695095/1‑239 (MQ=255)
                                                        aTTGTGGTGATGAAAATCGCCTCGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAGCCAGGGGCATGCCGGACCAGACCCCGCCGGTGAACCGCTGATCGGCGCACCGGTTCCGCCCCCTGCGCCACCTGGCCTTATCGGACCAtt                                                                                                                                                                      >  1:478089/1‑239 (MQ=255)
                                                          tgtgGTGATGAAAATCGCCTCGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAGCCAGGGTCTTGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGGGCACCCGTTCCGCCCTCTCGCCAACCTGGCCtttctggtgcattaa                                                                                                                                                                    >  1:596035/1‑226 (MQ=255)
                                                             ggtgATGAAAATCGCCTCGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAGCCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCCCCGGTTCCGCCCTCTCCCCTACCTGGACTTATCGTAGCATTaaa                                                                                                                                                                   >  1:460169/1‑236 (MQ=255)
                                                              gtgATGAAAATCGCCTCGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAGCCAGGGTCATGCCGGACCAGACCCGACCTGTGAACCGCTGACCGGGGCCCCGGGTCCGCCCTCCGCCCCACCTGGGCTTTACGTacacttatttt                                                                                                                                                                 >  1:27546/1‑228 (MQ=255)
                                                                gatgaAAATCGCCTCGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAGCCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTTTCGTAGCATTAATTCTCg                                                                                                                                                              >  2:224810/1‑239 (MQ=255)
                                                                                         aaTTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAGCCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTCCCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTggg                                                                                                                                     >  2:797454/1‑239 (MQ=255)
                                                                                           tttGCTCTGGCAACGACTGCCGGGAGCTCTGGGGATAGCCCAGGATTCACCCCTCTGGATCATCGGTGTATTAACGCTAACGGGAATTGCGGTGGGGTTGGTTATCCGTTTCAGCCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTCCCTGGGCTTATCGTACACTTAATTTCCGGTCTTGCTGCGGGCGCCAcccggtgcc                                                                                                                                   >  1:401714/1‑233 (MQ=255)
                                                                                           tttGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAGCCAGGGTCATGCCGGACCAGACCCTGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCCTTAATTCTCGGCCTTCCGTGCGGCGTACCCCTgggcc                                                                                                                                   >  1:536187/1‑237 (MQ=255)
                                                                                           tttGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAGCCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGcc                                                                                                                                                                                                    >  2:668708/1‑174 (MQ=255)
                                                                                           tttGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAGCCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTgggg                                                                                                                                    >  2:247245/1‑238 (MQ=255)
                                                                                           tttGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCTCTGGATCATAGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAGCCAGGGTCATGCCGGACCCGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCACTCTGCGCTACCTGGCCTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTAGGCCTtggg                                                                                                                                    >  2:438080/1‑238 (MQ=255)
                                                                                           tttGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGACTCACCCCTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAGCCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCACCCTCTGCGCTACCTGGACTTATCGGAGCCTTACTTCTCGGTCTTTCTGCGGGCGTAAGCCTGGGGc                                                                                                                                   >  1:313245/1‑239 (MQ=255)
                                                                                           tttGCTCTGGAAACGACTGCAGGGAACTCTGGGGATAGCCCAGGATTCACCCCTCTGGATCATCGGTGTATTAACGCTAAGGGGTATTGCGGGGGGGTTGGTTATCCGTTTCAGCCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGc                                                                                                                                   <  2:478089/239‑1 (MQ=255)
                                                                                              gctgtgccgaggcttGCCGGGACCTGGGGGGAGTGCCCATGTTTCCCCCCCGGGGTGCAGTGTTGTAGTAGCGCAAGCGGGTTGGGGGGTGGGTTGGGTTCGCCGTTTCGGCAAGGGTAATGCCGGACCAGCCCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGTCTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCgg                                                                                                                                <  1:506569/225‑1 (MQ=255)
                                                                                              gCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAGCCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATGAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCgg                                                                                                                                >  2:156488/1‑239 (MQ=255)
                                                                                                   ggCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAGCCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCTTCTGCGCTACCTGGACTTATCGTAg                                                                                                                                                                          >  2:60805/1‑192 (MQ=255)
                                                                                                   ggCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAGCCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCTTAGCATTAATTTTCGGTCTTGCTGTCGGCGTCAGCCTGTGCCCGGacaa                                                                                                                            >  2:524838/1‑235 (MQ=255)
                                                                                                   ggCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAGCCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACa                                                                                                                            >  2:10561/1‑238 (MQ=255)
                                                                                                    gcaaggcggtccGGGAACTGTGGGGATACCCCGGGTTCCCCCCCTTGGGCTCAGCGGTTATCAACGGTTAAGGGGTTGGGGGGGGGGGTTGGTTAACCGTCTCGGCAGGGGTCATGCCGGCCCAGCCCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACAt                                                                                                                           <  1:151668/228‑1 (MQ=255)
                                                                                                            ttccGGGACCTTTGGGGTAGCCCCGGGTTCCCCCCCTCTGGTTCATGGTTGTATAAACCTAAAGGGGTTTGGGGGTGGGGTTGGTTATCCGTTCCAGCCGGGGTCATGCCGGCCCAGCCCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATg                                                                                                                  <  1:10561/237‑1 (MQ=255)
                                                                                                                               tcccccgggttcccccccTCTGGATCTGCGGTTGTTAACCGTTACGGGTTTTTGGGGGGGGTTTGTTTTCCCGTTCCAGCAGGGGCCAGCCCGGCCCAGCCCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGACCGTCAATATCGCCCt                                                                                                <  1:556668/224‑1 (MQ=255)
                                                                                                                                gCCCAGGATTCACCCCTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAGCCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCTTCAGCCTGGGGCCGGAACATCCGATCATGACCGTCAATATTGCCCtctg                                                                                             >  2:342087/1‑236 (MQ=255)
                                                                                                                                gCCCAGGATTCACCCCTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAGCCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGCCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGGCTTGCTGGGGGCGTCACCCCTGGGCCGGGCCACCCCATCactaccctacaaaacgccctttg                                                                                              >  1:790977/1‑217 (MQ=255)
                                                                                                                                 cccAGGATTCACCCCTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAGCCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTTCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGACCGTCAATATCGCCCTTGCg                                                                                            >  2:408116/1‑239 (MQ=255)
                                                                                                                                    aGGATTCACCCCTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAGCCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCCCGGCCTTGCTGGCGGCGTCAGCCTGCGGCCGGAACATCCGTTCATGACCGTCAATATCGCCCTTGCGGtt                                                                                         >  2:481285/1‑239 (MQ=255)
                                                                                                                                            aCCCCTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAGCCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCTTGACCTCCAATATCGCCCTTGCGGTTGcacttg                                                                                  >  2:367785/1‑234 (MQ=255)
                                                                                                                                            aCCCCTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAGCCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGACCGTCAATATCGCCCTTGCGGTTGCCATTg                                                                                  >  2:143763/1‑239 (MQ=255)
                                                                                                                                            aCCCCTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAGCCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAACATTAATTCTCGGTCTTGCCGGCGGCCTCAGCCTGGGGCCGGAACCTCCGGCCAGAACCGTCACATGTGGCCTTGCGGTTGGCATTg                                                                                  >  1:84922/1‑239 (MQ=255)
                                                                                                                                                  ctGGATCATGGGTGTATTAGCGCTAACGGGTATTGGGGTGGGTTTGGTTATCCGTTTCAGCCGGGGTAATGCGGGCCCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGACCGTCAATATCGCCCTTGCGGTTGCCATTGGCGCTc                                                                           <  2:474340/239‑1 (MQ=255)
                                                                                                                                                     ggtccgcgtgtttatcgcGCTAGGGGTTATTGGGGTGGGGTGTGTTCCCGTCTCCGCCCAGGGCTATGCGGGACCCGCCCCCGCCTTTGCACCGGTGATCGGAGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCTTTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGACCGTCAATATCGCCCTTGCGGTTGCCATTGGCGCTCGt                                                                         <  1:833805/221‑1 (MQ=255)
                                                                                                                                                          tcGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAGCCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGCTCCTGACCGTCAATATCGCCCTTGCGGTTGCCATTGGCGCCCGCCTGtt                                                                    >  1:790481/1‑238 (MQ=255)
                                                                                                                                                          tcGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAGCCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTAGGGCCGGAACATCCGATCATGACCGGCAATATCGCCCTTGCGGTTGCCATTGGCGCTCGTCTGtt                                                                    >  2:552726/1‑238 (MQ=255)
                                                                                                                                                               gtATTAACGCTAACGGGTATTGTGGTGGGGTTGGTTATCCGTTTCAGCCAGGGTAATGCCGGACCAGACCCCGCATGGGAACCGCTGATCGGCGCATCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGGAACATTCATTCTCGTGCTTGCTTGCGGCGCCAGGCGCGGGCCGGAACATCCGGTCATGACCGCAACTATGGCCCCTGCTGTTGCATATGGCGCTCGTCTGTTAcccgc                                                              >  1:718568/1‑236 (MQ=255)
                                                                                                                                                               gtATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAGCCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGACCATCCGATCATGACCGTCAATATCGCCCTTGCGGTTGCCATTTGCGCTCGTCTGTTACcgc                                                               >  2:366533/1‑238 (MQ=255)
                                                                                                                                                                     aCGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAGCCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGACCGGCAATATCGCCCTTGCGGTTGCCATTGGCGCTCGTCTGTTACCGCCCGTCa                                                         >  2:395044/1‑238 (MQ=255)
                                                                                                                                                                           aaCGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAGCCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAGCATCCGATCATGACCGGCACTATCGCCCTTGCGGTGTCCCTGGTCGCCCGTCTGTTAcgccccgcacaccgca                                                   >  1:557912/1‑224 (MQ=255)
                                                                                                                                                                           aaCGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAGCCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGACCGTCAATATCGCCCTTGCGGTTGc                                                                                       >  2:466861/1‑204 (MQ=255)
                                                                                                                                                                             cGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAGCCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGACCGTCAATATCGCCCTTGCGGTTGCCATTGGCGCTCGTCTGTTACCGCGCGTCAACCGAATg                                                 >  2:449902/1‑239 (MQ=255)
                                                                                                                                                                                ttATGGCGGGGGGGTGTGTTTCCCGTTTCTGCCGGGTCCTTCCGGGCCCAGCCCCCTCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGACCGTCAATATCGCCCTTGCGGTTGCCATTGGCGCTCGTCTGTTACCGCGCGTCAACCGAATGGa                                               <  2:105557/237‑1 (MQ=255)
                                                                                                                                                                                gttttGGGGTGGGGTTTGTTTCCCTTTTCGGCCGGGGTCCTGCCCGACCAGCCCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTTCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCCGCCTGGGGCCGGAACATCCGATCATGACCGTCAATATCGCCCTTGCGGTTGCCATTGGCGCTCGTCTGTTACCGCGCGTCAACCGAATGGa                                               <  1:889275/235‑1 (MQ=255)
                                                                                                                                                                                    tggggTGGGGTTTGTTATCCGTTCACGCCAGGTTCAGTCCGGCCCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGACCGTCAATATCGCCCTTGCGGTTGCCATTGGCGCTCGTCTGTTACCGCGCGTCAACCGAATGGAGTGGa                                          <  2:313245/236‑1 (MQ=255)
                                                                                                                                                                                    ggggggggggGTGTTTACCGGTTCCAGCCAGGGTCAGGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGACCGTCAATATCGCCCTTGCGGTTGCCATTGGCGCTCGTCTGTTACCGCGCGTCAACCGAATGGAGTGGa                                          <  2:679222/233‑1 (MQ=255)
                                                                                                                                                                                     gggggggggTTGGGTTTCCCGTTTCAGCCGGGTTCGGGCCGCACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGACCGTCAATATCGCCCTTGCGGTTGCCATTGGCGCTCGTCTGTTACCGCGCGTCAACCGAATGGAGTGGa                                          <  2:401714/234‑1 (MQ=255)
                                                                                                                                                                                       gtgggggggtgttgttcccctttCCGCCCGGGTCGTTGCCGGCCCGCCCCCCGCCGTGGAACCGGTGTACGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGTAACATCCGAGCATGACCGTCAATATCGCCCTTGCGGTTGCCATTGGCGCTCGTCTGTTAGCGTGCGACAACCGAATGGAGTGGa                                          <  1:366533/219‑1 (MQ=255)
                                                                                                                                                                                           gggTTGGTTATCCGTTTCAGCCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGGTCATGACCGTCAATATCGCCCTTGCGGTTGCCTTTGTGGCCCGCCTTTTACCCCCCGTCAACCGAATGGGGGGGAAAAtttt                                   >  1:511829/1‑239 (MQ=255)
                                                                                                                                                                                           gggTTGGTTATCCGTTTCAGCCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGACCGTCAAAATCGCCCTTGCGGTTGCCCTTGAGCCTCGCCTGTTCCCGCGCGTCACCCACATGGAGTgggaaatttt                                   >  1:839146/1‑231 (MQ=255)
                                                                                                                                                                                                            ttCAGCCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGACCGTCAATATCGCCCTTGCGGTTGCCATTGGCGCTCGTCTGTTACCGCGCGTCAACCGAATGGAGTGGCCTATTTTAGCCTCTGCCGGaa                     >  2:636645/1‑238 (MQ=255)
                                                                                                                                                                                                                     gggTCATGCCGGACCAGCCCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGACCGTCAATATCGCCCTTGCGGTTGCCATTGGCGCTCGTCTGTTACCGCGCGTCAACCGAATGGAGTGGACTATTTTAGCCTCTGCCGGAACCATCGgcg            <  2:536187/239‑1 (MQ=255)
                                                                                                                                                                                                                      ggTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGCGGCCGGATAATTCGCTCCTGCACGTCACTGTCGCCCTTGCGGGTGCCATTGGCGGCACTTGTTTACCCGTCTGCAAACGCCGGGAGGGGACCATTTTTGCCTCTTCCGGAACCAACggtga           >  1:263442/1‑236 (MQ=255)
                                                                                                                                                                                                                         tgaTGCCGGCCGAGCCCCCGCCTGGGACCCGCTGATCGGCCCACCGGGTCCGCCCTCTGCGCTCCCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGACCGTCAATATCGCCCTTGCGGTTGCCATTGGCGCTCGTCTGTTACCGCGCGTCAACCGAATGGAGTGGACTATTTTAGCCTCTGCCGGAACCATCGGCGCa          <  1:524838/236‑1 (MQ=255)
                                                                                                                                                                                                                         tccggCCGGACCAGCCCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGACCGTCAATATCGCCCTTGCGGTTGCCATTGGCGCTCGTCTGTTACCGCGCGTCAACCGAATGGAGTGGACTATTTTAGCCTCTGCCGGAACCATCGGCGCa          <  2:27546/234‑1 (MQ=255)
                                                                                                                                                                                                                         tcagtccGGACCGACCGCCGCCTGTGAACCGCTGATCGGCGCTCCGGTTCCGCCCTCTGCGCTACCTGGTCTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGACCGTCAATATCGCCCTTGCGGTTGCCATTGGCGCTCGTCTGTTACCGCGCGTCAACCGAATGGAGTGGACTATTTTAGCCTCTGCCGGAACCATCGGCGCa          <  2:26623/233‑1 (MQ=255)
                                                                                                                                                                                                                         ccagcgcGGACCAGCCCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGACCGTCAATATCGCCCTTGCGGTTGCCATTGGCGCTCGTCTGTTACCGCGCGTCAACCGAATGGAGTGGACTATTTTAGCCTCTGCCGGAACCATCGGCGCa          <  2:596035/232‑1 (MQ=255)
                                                                                                                                                                                                                            tGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGACCGTCAATATCGCCCTTGCGGTTGCCCTTGGCGCTCGTCTGTTACCGCGCGTCACCCGAATGGAGTGGACTATTTTAGCCTCTGCCGGaa                     >  2:395883/1‑224 (MQ=255)
                                                                                                                                                                                                                                 ggcccAGCCCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGACCGTCAATATCGCCCTTGCGGTTGCCATTGGCGCTCGTCTGTTACCGCGCGTCAACCGAATGGAGTGGACTATTTTAGCCTCTGCCGGAACCATCGGCGCACTGTTTgg  <  2:460169/236‑1 (MQ=255)
                                                                                                                                                                                                                                  gACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGACCGTCAATATCGCCCTTGCGGTTGCCATTGGCGCTCGTCTGTTACCGCGCGTCa                                                         <  2:663067/183‑1 (MQ=255)
                                                                                                                                                                                                                                  gACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGACCGTCAATATCGCCCTTGCGGTTGCCATTGGCGCTCGTCTGTTACCGCGCGTCa                                                         >  1:663067/1‑183 (MQ=255)
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                          
TGTTGTTATTATCGCTCCCCGCCGTGGCAATTGGGATTGCGTCCAGTCTTATTCTGATTGTGGTGATGAAAATCGCCTCGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCTCTGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAGCCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGACCGTCAATATCGCCCTTGCGGTTGCCATTGGCGCTCGTCTGTTACCGCGCGTCAACCGAATGGAGTGGACTATTTTAGCCTCTGCCGGAACCATCGGCGCACTGTTTGG  >  NC_000913/2509653‑2510109

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: