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

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