Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,065,570 T→A 14.8% T66S (ACC→TCC)  yihO ← putative sulphoquinovose importer

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,065,5700TA14.8% 159.3 / 2.3 77T66S (ACC→TCC) yihOputative sulphoquinovose importer
Reads supporting (aligned to +/- strand):  ref base T (43/18);  new base A (2/10);  total (48/29)
Fisher's exact test for biased strand distribution p-value = 7.61e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01

AATTGCCGGGATCATCGCACCATACGAGCAGTTCATCAGACTGTATGAAAGTCCAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGACATAATCACTCCTTATAACGCCCCGCAGGGCGTCGTTTATTATTAT  >  NC_000913/4065334‑4065808
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                   
aaTTGCCGGGATCAGCGCACCATACGAGCAGTTCATCAGACTGTATGAAAGTCCAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTaaa                                                                                                                                                                                                                                                  <  2:502888/239‑1 (MQ=255)
 aTTGCCGGGATCATCGCACCATACGAGCAGTTCATCAGACTGTATGAAAGTCCAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGCTaa                                                                                                                                                                                                                                                   >  2:482376/1‑238 (MQ=255)
   ttccGGGATCATCGCCCCATACGAGCAGTTCATCAGACTGTATGAAAGTCCAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAAcg                                                                                                                                                                                                                                               <  2:221190/237‑1 (MQ=255)
             atcGCACCATACGAGCAGTTCATCAGACTGTATGAAAGTCCAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCACTATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTaaaa                                                                                                                                                                                                                                      >  2:115847/1‑237 (MQ=255)
                   ccATACGAGCAGTTCAACAGACTGTATGAAAGTCCAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGc                                                                                                                                                                                                                              <  2:614912/239‑1 (MQ=255)
                         gAGCAGTTCATCAGACTGTATGAAAGTCCAAACATCATAAACAGCGCGGTAGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAg                                                                                                                                                                                                                       <  2:664544/239‑1 (MQ=255)
                           gCAGTTCATCAGACTGTATGAAAGTCCAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTAAGCTTATCGGTAAACCCGGTAAAAAACTTCGCCACCCGaa                                                                                                                                                                                                                    >  2:799733/1‑239 (MQ=255)
                           gCAGTTCATCAGACTGTATGAAAGTCCAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGATGGCACTAAAATGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGACCCAATATTTTTCCGCAGGACGAGTAAAAAACCGGGTAGCTAACCGGAAAAAACGGTAAAAAACTTACCCAGCAGaa                                                                                                                                                                                                                    >  1:33682/1‑239 (MQ=255)
                            cAGTTCATCAGACTGTATGAAAGTCCAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGaaa                                                                                                                                                                                                                   >  2:203950/1‑239 (MQ=255)
                              gTTCATCAGACTGTATGAGAGTCCAAACCTCATAAACAGCGCGATGGCAATCGTCGTTTTCACCGGCAAACAAAAGGTGGTGGCGATAAACTGGCGCGTTGCAATTAACGCAGCCGGAACCGCTGCATATACATTAAAGGGACGGAATTTGCCTTTAAGCCCAATCTTTTTTCTCCAGTCGAGGAAAAATCCAGTGAGGATACCTGTAACACAGGTAAAAAAATTCCCCAACACaaaaa                                                                                                                                                                                                                 >  1:828693/1‑237 (MQ=255)
                              gTTCATCAGACTGTATGAAAGTCCAAACATCATAAACAGCGCGGTGACAATCGTCGTTTTAACCGGCAAACAAAAGGAGGTGGCGCTAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCGGCATATAAAATAACGGGCCGAAATTTGCCTATAGGCCCAAAATTTTTCCACGAGTCGAGGAAAAATCCGGTGAACATTACAGAAAAAGAGGAaaaaaaatacaccaaccaaaaaa                                                                                                                                                                                                                   >  1:453070/1‑222 (MQ=255)
                                 catcaGACTGTATGAAAGTCCAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAAATTCGCCACCAaaaaaaaaa                                                                                                                                                                                                              >  1:968819/1‑234 (MQ=255)
                                         tGTATGAAAGTCCAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATCAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCTTATCGGTAAACGCGGTAAAAAACTTCCCCACCAGAAAGATAATCCAGccc                                                                                                                                                                                                      >  2:646725/1‑238 (MQ=255)
                                         tGTATGAAAGTCCAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAAAAAGAAAATCCAGCCa                                                                                                                                                                                                      >  2:303756/1‑239 (MQ=255)
                                           tATGAAGGTCCAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATa                                                                                                                                                                                                    >  2:855563/1‑239 (MQ=255)
                                           tATGAAAGTCCAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATACCGGTAAAACCGGTAAAAAACTTCCCCCCCAAAAAGAAAATCCCGCCaaa                                                                                                                                                                                                    >  1:790901/1‑237 (MQ=255)
                                           tATGAAAGGCCAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTCACCGGCAACCAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGCCGACATTTGCCTTTAGGCCCAATATTTTTCCTCTCGCCGCGGAAAAATCCGGTACGACTATCGGTCAACGCGGCAACAAACTTCGCCACCTGAAAAGAAATCCCCCCaaa                                                                                                                                                                                                    >  1:235580/1‑237 (MQ=255)
                                                   tCCAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAACATCCGGGGAGCCTATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCCTAGTAAGCAg                                                                                                                                                                                            >  2:685772/1‑239 (MQ=255)
                                                     cAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAACTAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAAAGCGGTAAAAAACTTCGCCACCAGAAAGATAATACCGCCAAAATAAgaaaaa                                                                                                                                                                                          >  1:206812/1‑234 (MQ=255)
                                                     cAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGc                                                                                                                                                                                          >  2:675078/1‑239 (MQ=255)
                                                     cAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGc                                                                                                                                                                                          >  2:617983/1‑239 (MQ=255)
                                                     cAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACAAGAAGGATAATCCCGCCATAGTAAGCAGGc                                                                                                                                                                                          >  2:144491/1‑239 (MQ=255)
                                                     cAAACATCATAAACAGCGCGGTGGCAATCGACGTTTTAAACGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTCGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGACTATCGGGAAAATCGGCAAACAACTTAGCCACAAAAAAGAAAATCCAGCCaaaaaaaagagga                                                                                                                                                                                          >  1:259288/1‑227 (MQ=255)
                                                       aaCATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGAGGCGATAAACTGGAGCGTTGCAATTAACGCAGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCAGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTGAAAAACTTCGCCCCCAGAAAGATAATCCCTCCATAGGAAGCAGGCa                                                                                                                                                                                        >  2:338499/1‑238 (MQ=255)
                                                           tcatAAACAGCGCGGTGGCAACCGTCGTTTTAACCGGCAAACAAAAGGTGGGGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGcc                                                                                                                                                                                                                             <  2:660485/200‑1 (MQ=255)
                                                                 cacgggcggggaacagcgtcctTATAACCGAAAAAAAAAGGGTGGTGGGAAAAAAAGGGAGGGTTGAAATTACCGCCGGCGGAACCGTTGCATAAAAAATAAAGGGACGAAATTGGCCTTTAGGCCCAATATTTTCCCGGGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCCCCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCa                                                                                                                                                                             <  1:203950/223‑1 (MQ=255)
                                                                    gcgcGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCg                                                                                                                                                                          >  2:844933/1‑239 (MQ=255)
                                                                     cgcgGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGAAAATCCCGCCATAGTAAGCAGCCATCCCTAATTCATCgg                                                                                                                                                                         >  2:880498/1‑239 (MQ=255)
                                                                      gcgGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCGGAAAGATAATCCCGCCATAGTAAGCAGGAATCCCTAATTCATCgga                                                                                                                                                                        >  2:401717/1‑238 (MQ=255)
                                                                      gcgGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGt                                                                                                                                                                        >  2:679296/1‑239 (MQ=255)
                                                                             cAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAAt                                                                                                                                                                 >  2:701408/1‑239 (MQ=255)
                                                                                 tcgttttttCAACGGCAAACAAAGAGGTGGTGGAAAAAAATGGGATGTTGAAAAAAACCGCCCCGGAAACAGTGCCATATAAAAAAAGGAAACGAATTGTGCCTTTAGGCACAATTTTTTACCGCAGGCCCAGAAAAAATCCGGTGAGTATACCGGTAAGCGGATAAAAAAACTTCACCACAAGAAAGAAAAACCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTg                                                                                                                                                              <  1:844933/234‑1 (MQ=255)
                                                                                  tccttttaccccGGCAAAAAGGGGGTGGTGCCTATAGACTGACGCGTGGAAATTGGCCTCGGCAGACCCGTGGCATAAAAAATAAAGAGACGTAATGTGCGTTTAGGCCCAATTTTATTCCGCGAGTCGATTAAAATTCGTGTGAGCACTTCGATAAAGGAGGTAAAAACCTCCACCAGAGCAAATATAAACCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTgag                                                                                                                                                            <  1:401717/229‑1 (MQ=255)
                                                                                    cgtTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCACTTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGGGTCGAGAAAAAATCCGGGGAGCATAACGGCAACCGCGGTAAAAAACTTCCCCAACAGAAAGCTAATCCCGCCATTTAAACAGCCAACCCCAAATTCATCGGAAAAAAATTTgaggg                                                                                                                                                          >  1:861753/1‑238 (MQ=255)
                                                                                            cccgaaaacaaaAGGGTGGGGGGAAAAACTTGGACCTTTAAAATTAACGCCGCCGGACCCGTGAAATAAAAAAAAAAGGGACAGAATTTGCTTTTAGGCCAAATATTTTTCGGGGAGCAAATAAAAAATCGGGGGAGTATATCGGAAAACGCGGAAAAAATCTCCGCCACAAGAAAGATAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACag                                                                                                                                                   <  1:126601/234‑1 (MQ=255)
                                                                                            ccGGGAAAAAAAAGGGGGTGGGGAAAAAAAGGGGTGGTGCAATTTACCGCGGCCGGAACCTATGTATAAAAAAAAAAGGGAAAAATTTTGCCTTTAGCACAAAATTTTTCCCGGGAGTCGGTAAAAGACCCGGTGAGCATTGTGAAAAACGGGGAAAAAAATCTCGCCACCAGAAAAAAAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACag                                                                                                                                                   <  1:880498/239‑1 (MQ=255)
                                                                                                caaacaaaAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTCTAAAATTTCAGGAGATCCACAGTa                                                                                                                                               >  2:346972/1‑238 (MQ=255)
                                                                                                              gtggCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGAACCGATGCAGAACAt                                                                                                                                 >  1:161833/1‑239 (MQ=255)
                                                                                                                             gAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCCTCGGCATAAAATTTGAGGAAATACCAAGGACCGAGTCAAAGAATTAAAATCAAAcag                                                                                                                    >  1:526505/1‑236 (MQ=255)
                                                                                                                                 gtgtaaaTTAACGCCGCGAGACCCGATGAATAAAAAAAAAAGGGACAAAATTGCCTTTAAGGCCAAATTTTTTTCCGCGAGGCGATAAAAAATCCGGTGAGCATTTCGGTAAACGCGGTAAAAAACTCCGCCACAAAAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGt                                                                                                                <  2:396458/233‑1 (MQ=255)
                                                                                                                                 ctttgaataaagccccgcggGACCCGTCGGATAAAAAAAAAAGGGAAGAATTTGCCTTATAGGCCAAATAGTTGCCCGGAAGAGGAGTAAAACTCGGGGGGTAATTCGGTTAAACGGGAAAAAAAATTTCGCCACCAGAAAGAAAAGCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACg                                                                                                                 <  1:701408/224‑1 (MQ=255)
                                                                                                                                  tttaaatttacccccgcggGACCCGCTGAATATAAAAAAAAGGGACGAATTTTTCCGTTAGCCCAAATTTTTTCCCGCGTGTCGAGAAAAAAGCCGGTGAGAATATCGGTAACCGCGAAAAAAAACTCCGCCACCAGAAATATACTCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGt                                                                                                                <  1:679296/225‑1 (MQ=255)
                                                                                                                                  ttGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACTCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGAATAAAATTTGCGGAGATACAGAGAACCGATGCAGAGCATTAAATTCGAAACGACGTc                                                                                                               >  1:474700/1‑238 (MQ=255)
                                                                                                                                  ttGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGCAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCTTAGTAAGCAGGCATCCCTAATTCCTCGGTATAAAATTTGAGTAGCTACAGAGTCCCGCTCCAGAGCATTAAATTAAAACCAACggg                                                                                                               >  1:781368/1‑236 (MQ=255)
                                                                                                                                    gCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAAGGCATTAAATTCGAACCGACGTAGc                                                                                                             >  2:253255/1‑238 (MQ=255)
                                                                                                                                           aaCGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGGGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATAAAAAGTACCAAGGCAAAGAATTAAATTCAAACCGACGGGGCccaaacc                                                                                                      >  1:396458/1‑239 (MQ=255)
                                                                                                                                                                  aaaaaaaaaaGGGAAAAATTTCGCCTTTAGCCAAAATTTTTCCCCGGGAGTCAAAAAAAAACCCGGTGAGCTTATCGTAAACCGCGAAAAAAACCTTCGCCACCAAAAAGATAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAGCCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGa                                                                                <  2:18276/237‑1 (MQ=255)
                                                                                                                                                                    aaaatgaaGTGACAAATTTTGTATGTAGCCCATATTTGTTGCCGAGAGTCGAGTAAAACTGCGGTGGAAATTAAGATACACGCGATAAAAAGAGTCGCAACCAAAAAGATATACCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGAtt                                                                              <  1:346972/239‑1 (MQ=255)
                                                                                                                                                                    aaaataaaGGGACAAAATTTGCTTTTAGCCCCAATTTTTTTCCGGGAGTCGAGAAAAAACCCGGTGAGCATATCGGAAAACGGGGAAAAAACCTTCCCCCCCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGAtt                                                                              <  2:790901/239‑1 (MQ=255)
                                                                                                                                                                    aaaaaaagGGGACAAATTGTTCCTTTGGCCTAAATTTTTTACCGAGAGGTGAAAAAAAATCCGGGGTGTAAAACGGTAAACGGGGAAAAAACCTCCGCACCAAGAAAGTAAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGAtt                                                                              <  1:144491/239‑1 (MQ=255)
                                                                                                                                                                     aaaaaaaGGGAAATAATTTGTTTTTGGGCCCAAATATTTTCCAGAAGTCAGGAAAAAAACCGGTGAGTAACACGGAAAAAGTGGGAAAAAAACTGCGCCACAAAAAAAAAAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGAtt                                                                              <  1:503544/239‑1 (MQ=255)
                                                                                                                                                                         aaaGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGTCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAAAGTACCGAt                                                                                                                                          >  2:510331/1‑174 (MQ=255)
                                                                                                                                                                         aaaGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCg                                                                                                                                                                          >  2:956731/1‑142 (MQ=255)
                                                                                                                                                                         aaaGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAGCCGACGTCCCCCATACCATAGGGGATTTTTTCCCGCAAAATAAGtt                                                                         >  1:18276/1‑239 (MQ=255)
                                                                                                                                                                           aGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATaaaa                                                                                                                                                                  >  2:411182/1‑148 (MQ=255)
                                                                                                                                                                               aCGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAAAGCATTAAATTCGAACCCACGTCGCCCATACCCTAGGCGATTTTTTCCCGCAGATTCCTTTTTAAt                                                                    >  1:69020/1‑238 (MQ=255)
                                                                                                                                                                                cGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCccata                                                                                                        >  2:342531/1‑201 (MQ=255)
                                                                                                                                                                                cGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGACCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAAATTCAGTTTTaaatt                                                                  >  1:329211/1‑236 (MQ=255)
                                                                                                                                                                                cGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCAGCCATAGTAAGCAGGCATCCCTAATTCATCGGGATAAAATTTGACGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGGCGCCCATACCAAAGGCGATTTTTTCCCGCACATTCAGTTTTaaag                                                                   >  1:828655/1‑236 (MQ=255)
                                                                                                                                                                                  aaaTTTGCCTTTAGGCCAAATTTTTTCCGGCAGATAGAAAAAAACCCCGGTGAGCATTTCGGAAAACGCGATAAAAACCTTCGCCACCAAAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCa                                                                <  2:968819/239‑1 (MQ=255)
                                                                                                                                                                                    tttttCCTTGAGGCCCAAATTTTTTCCGCGGGTCGGATAAAAATCCGGGGAGCATATCGGAAAACGCGGTAAAAACCTTCGCCACAAAAAAGATAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCa                                                                <  2:728637/236‑1 (MQ=255)
                                                                                                                                                                                                aGGCCAAATTTTTATCCGCGAGGCGAGAAAAAATCCGGGGACAATTGCGGAAAAGGGGGAAAAAACCTGCGCCAGCAGAAAGATAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGt                                                    <  2:526505/239‑1 (MQ=255)
                                                                                                                                                                                                 tgCCCAATTTTTTTCCCGGTGTTGCGAATAAATCCCGTTGACCATATCGGTACACGGTAAAAAAACCTTCGCCACAAAAAAGAAAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGt                                                    <  2:151944/238‑1 (MQ=255)
                                                                                                                                                                                                   cccAATATTTTGCCGGGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACAAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCa                                                  <  2:329211/238‑1 (MQ=255)
                                                                                                                                                                                                      aaTATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAACCTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTCTAAAATTTGAGGAGATACAGAGTACCGATGAATAGCATTACATTCGAACCGACGCCGCCCCTACCATAGGCGATTTTTTCCCCAACATTCAGTTTTATCTGCAGGCGATTATGGtacaaac                                              >  1:168383/1‑233 (MQ=255)
                                                                                                                                                                                                             ttcccGCGAGTCGAGTAAAACCGCGGTGAGAATAGCGAAAACCGCGGAAAAAACCTCCGCCACCAAAAAGATAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGACATAATCa                                        <  2:933724/236‑1 (MQ=255)
                                                                                                                                                                                                             ttCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCCCAGTTTCAGTTTTAATTTCAGGGGATTATGGTCaaaataaatca                                        >  1:54684/1‑228 (MQ=255)
                                                                                                                                                                                                                 cccgtgggtcgtataaaaTCCGGGGGAGTATATCGGTAAGCTCAGAAAACAGCTTCCCAATAAAAAAAATAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATTCCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGACATAATCACt                                      <  1:482376/225‑1 (MQ=255)
                                                                                                                                                                                                                   cgagtcgagtAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGCTTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGACATAATCACTCCtt                                  >  2:852460/1‑239 (MQ=255)
                                                                                                                                                                                                                   cgagtcgagtAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAAAGAACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAAATTCCGTGGGTTATGGTCAGAAAAAAACACCCCtt                                  >  1:933724/1‑239 (MQ=255)
                                                                                                                                                                                                                                       cGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAAt                                                                    >  2:338167/1‑188 (MQ=255)
                                                                                                                                                                                                                                       cGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGCTTATGGTCAGACATAATCACTCCTTTTAACGCCCCGCAGGGc                 >  2:147277/1‑239 (MQ=255)
                                                                                                                                                                                                                                         ggtgccaataaCGGTATCCGCGGTAAAAACCTGCCCCCCAGAAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGATTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGACATAATCACTCCTTATAACGCCCCGCAGGGcg                <  2:206812/232‑1 (MQ=255)
                                                                                                                                                                                                                                          gtagagatatCGGAAAACGGGGTAAAAAACGTCGCCCAGAAAAAAATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATGAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGACATAATCACTCCTTATAACGCCCCGCAGGGcg                <  2:861753/232‑1 (MQ=255)
                                                                                                                                                                                                                                               cATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGACATAATAAGTCCTTTTAAAGCCCCGCAGGACTCCGttt          >  2:130714/1‑239 (MQ=255)
                                                                                                                                                                                                                                                 acgcGGTAAACCGGGAAAAAACTTCCCCAACAAAAAAGAAAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGACATAATCACTCCTTATAACGCCCCGCAGGGCGTCGTttat        <  2:474700/236‑1 (MQ=255)
                                                                                                                                                                                                                                                  tacGGGAAACGGCGGAAAAAACCTCCGCCACAAGAAAAAAAATCCCGCCATAGTAAGCAGGCATCCCTAATTCAGCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGACATAATCACTCCTTATAACGCCCCGCAGGGCGTCGTttat        <  2:69020/237‑1 (MQ=255)
                                                                                                                                                                                                                                                        tAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGACATAATCACACCTTATAACGCCCCGCAGGGCGTCGTtttttattat  >  1:657078/1‑239 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                   
AATTGCCGGGATCATCGCACCATACGAGCAGTTCATCAGACTGTATGAAAGTCCAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGACATAATCACTCCTTATAACGCCCCGCAGGGCGTCGTTTATTATTAT  >  NC_000913/4065334‑4065808

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: