Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,065,579 T→A 16.0% T63S (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,5790TA16.0% 124.9 / 4.7 68T63S (ACC→TCC) yihOputative sulphoquinovose importer
Reads supporting (aligned to +/- strand):  ref base T (43/9);  new base A (2/9);  total (49/19)
Fisher's exact test for biased strand distribution p-value = 8.14e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TACGAGCAGTTCATCAGACTGTATGAAAGTCCAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGACATAATCACTCCTTATAACGCCCCGCAGGGCGTCGTTTAT  >  NC_000913/4065356‑4065802
                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                 
tACGAGCAGTTCATCAGACTGTATGAAAGTCCAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGGGTAAAAAACCGGTGATCATATCGTTAAACGCCGGaaaaa                                                                                                                                                                                                                              >  2:624630/1‑229 (MQ=255)
      cAGTTCATCAGACTGTATGAAAGTCCAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCg                                                                                                                                                                                                                       >  2:166368/1‑229 (MQ=255)
      cAGTTCATCAGACTGTATGAAAGTCCAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGGGAGCATATCGGTAAACGAGGTAAAAAAATTCg                                                                                                                                                                                                                       >  2:121218/1‑229 (MQ=255)
        ggtcaacaggatgttgtaaaGTCCAAACAATATAAACGGGGGGGTTGAAATCGGCGTTTTAACCGGCAAACAAAAGGTGGTGGCTATAAACTGGAGCGGTGAAATTAACTCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAAGTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGcc                                                                                                                                                                                                                    <  2:583262/212‑1 (MQ=255)
             acaGACTTTATGAAAGCCCAAACAACAATAACGGCGGGGTGGCAACCGGCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAg                                                                                                                                                                                                               <  2:139625/228‑1 (MQ=255)
                 aCTGTATGAAAGTCCAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCAAGTCGAGTAAAATTCCGGGGTGCATCTCGGTAAACGCGGTAAAAAACTTAGCCACCATaaa                                                                                                                                                                                                            >  1:294231/1‑228 (MQ=255)
                   tGTATGAAAGTCCAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCAGTGGGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGaa                                                                                                                                                                                                         >  2:473877/1‑228 (MQ=255)
                     tATGAAAGTCCAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCAACCAGAAAGATaa                                                                                                                                                                                                       >  1:376814/1‑229 (MQ=255)
                     tATGAAAGTCCAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAACATTTTTCCGCGAGTCAGGTAAAAATCCGGTGGGCATATCGGTAAACGCGGTAAAACACTTCGccacacaaaaaaaaa                                                                                                                                                                                                       >  2:46373/1‑218 (MQ=255)
                     tATGAAAGTCCAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCACATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCAAATCGGTAAACGCGGTAAAAAACTTCGCACCCGAAAGAAAAt                                                                                                                                                                                                      >  2:173910/1‑229 (MQ=255)
                     tATGAAAGTCCAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTAGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCAACCGGAAAGATaa                                                                                                                                                                                                       >  1:646987/1‑229 (MQ=255)
                       tGAAAGTCCAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATc                                                                                                                                                                                                     >  1:496971/1‑229 (MQ=255)
                           aGTCCAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGc                                                                                                                                                                                                 >  1:411094/1‑229 (MQ=255)
                               cAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGCTAAACTGAGGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCAAATATTTTTCCCCAGGTCGCGGAAAAATCCGGAGAGCAAAGCGTTAACCGAGTAAAAAAACTTCGCACACCAAAAGAAAAACCCGGCaaa                                                                                                                                                                                             >  2:520249/1‑227 (MQ=255)
                               cAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCCATATTTTTCCGCGAGTCAAGTAAAAATACGGTGAGCATATCGGTAAACCCGGGACAAAACTTCGCCACCAGAACGTTAATCCCGCCATa                                                                                                                                                                                             >  1:379959/1‑229 (MQ=255)
                               cAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGGGCATATCGGTAAACGCGGTCAAAAACTTCGCCACCAGAAAAATAATCCAGCCATa                                                                                                                                                                                             >  1:208144/1‑229 (MQ=255)
                               cAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGcccaa                                                                                                                                                                                             >  1:462618/1‑226 (MQ=255)
                               cAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATa                                                                                                                                                                                             >  1:254981/1‑229 (MQ=255)
                               cAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATa                                                                                                                                                                                             >  1:237646/1‑229 (MQ=255)
                               cAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCAGATCGGTAAACGCGGAAAAAAACTTCGCCACAAGAAAGAAAAAccccacaaa                                                                                                                                                                                             >  2:99974/1‑223 (MQ=255)
                                              gcgcGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGATTAAAAATCCGGTGAGCATATCGGGAAACTCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGAAAGAAGGCATccc                                                                                                                                                                              >  1:228056/1‑229 (MQ=255)
                                              gcgcGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCAACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATccc                                                                                                                                                                              >  1:549315/1‑229 (MQ=255)
                                              gcgcGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGGAAACGCGGTAAAAAAATTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATccc                                                                                                                                                                              >  1:689746/1‑229 (MQ=255)
                                              gcgcGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGCTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACAGGAAAGATAATCCCGCCTCAGTAAGAAGGAAcccc                                                                                                                                                                              >  1:727718/1‑229 (MQ=255)
                                               cgcgGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATACCGCCATAGAAAGCAGGCATCCCt                                                                                                                                                                             >  1:74820/1‑229 (MQ=255)
                                               cgcgGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGGGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCAACAGAAAGATAATCCGGCCATAGTAAGAAGGCAcccca                                                                                                                                                                             >  2:42123/1‑228 (MQ=255)
                                                gcgGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGCAATAAACTTCGCCACCCGAAAGATAATCCCGCCATAGCAAGCATGCATCCCTa                                                                                                                                                                            >  1:191109/1‑229 (MQ=255)
                                                       taagcgtcttTTAAACCGGAAAAAAAAAGGTGGTGGCAGAAAATTGGAGGGTTAAAATTACCGCCTCCGGAACCGTTTAATAAAAAATAAAGGGACAAATTTTGCTTTTAGTCCCAATATTTTTCCGAGAGTCGAGTAAAAACCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAAATAATCCCGCTATAGTAAGCAGGCATCCCTAATTCATc                                                                                                                                                                     <  2:197675/225‑1 (MQ=255)
                                                                       cGGAAAAAAAAAGGGGGGGGAGAAAAAATGGGGGGGTGCAATAAACCCCGCAGAAACCGCTGAAAATAAAAGAAAGAGAAAAAATTTGATTTTGGGACAAATATTTTCCGGGGAGGGGAGAAAAAAGCTGGGGAGCATGCGGGAAAACGCGGTAAAAAACTCCCCCACAAAAAAGTAAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTgagg                                                                                                                                                     <  2:599472/229‑1 (MQ=255)
                                                                          caaacaaaAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGAATAAAATTTgaggga                                                                                                                                                   >  2:516409/1‑226 (MQ=255)
                                                                                       gtggCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGAGGATACAGAGAACCGa                                                                                                                                     >  1:606870/1‑229 (MQ=255)
                                                                                        tggCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCATAAGTAAGCAGGCATACCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGAt                                                                                                                                    >  1:399985/1‑229 (MQ=255)
                                                                                                       gCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATAACGAGTACCGATGCAGAGCATTAAAtt                                                                                                                     >  1:348087/1‑229 (MQ=255)
                                                                                                          ttGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGATGAGATACAGAGTACCGATGCAGAGCATTAAATTCg                                                                                                                   >  1:438582/1‑228 (MQ=255)
                                                                                                           tGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCACCCCTAATTCATCGGTATAAAATTTGAGGGAATACAGAGTACCCATCCAAAGAATAAAATTCaaa                                                                                                                 >  2:299121/1‑229 (MQ=255)
                                                                                                             cAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATCCATCGGTATAAAATTTTAGGAAAAACAGAGGACCGATCAAAAGCATTAAATTACAAcc                                                                                                               >  2:100180/1‑229 (MQ=255)
                                                                                                                  gcacgcaggagcactggttgataaAAAAAAAAGGGACAAAGATTGGCCTTAAGGCCAAATTTTTTCCCGAGAGCCGAAAAAAACGCGGGGGAGAGTACGGTGAAACAGACAAAAGAACCACCCCACCAGAAACAAAACACCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGAc                                                                                                            <  2:191109/208‑1 (MQ=255)
                                                                                                                  accgccgcCGGAACCGATGTATATAAAATAAAGGGACGAATTTTGCTTTTAGGTCAAATTTTTTACCGAGAGTCGAGTAAAAACCGGGAGACCATACGGGAAAAGGAGGTAAAAAACTTCGCCACCAAAAATAAAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACg                                                                                                           <  2:720422/226‑1 (MQ=255)
                                                                                                                        gccgGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCGGGTGAGCATATCCGTAAACGCGGTAAAAAACATCGCCACCAGAAAGATAATCCCGACATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGCGGAGAAACACTGTCCCCATTCAGAGCATTAAATTCAAACCGACGTGGCcca                                                                                                    >  2:71198/1‑229 (MQ=255)
                                                                                                                         ccgGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCCTTAAATTCGAACCGACGTCGCccat                                                                                                   >  1:396688/1‑229 (MQ=255)
                                                                                                                                   tGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCACAATATTTTTCCGCGAGTTGAGTAACAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGACTTCCCTAATTCATCGGTCTAAAATTTGAGTCGATACAGAGTACCGATCCAGATCATTAAATTCGTACCGACGCAGTCCATACCTTCGGc                                                                                          >  1:575784/1‑228 (MQ=255)
                                                                                                                                            aataaaGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTcc                                                                                >  1:457959/1‑229 (MQ=255)
                                                                                                                                             ataaaGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCACTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGAGTTTTTccc                                                                               >  1:676284/1‑229 (MQ=255)
                                                                                                                                               aaaGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATAACAAAGGCGATTTTTTccccc                                                                             >  2:461292/1‑227 (MQ=255)
                                                                                                                                               aaaGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCAAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGAAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATCCGAACCGAAGTCGCCCATACCCTAGCGGATTTTTTCACGc                                                                             >  1:612182/1‑229 (MQ=255)
                                                                                                                                                 aGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTgag                                                                                                                                                      >  1:18679/1‑154 (MQ=255)
                                                                                                                                                 aGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAAATTACCGATGCAAAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTccccccg                                                                           >  1:462331/1‑225 (MQ=255)
                                                                                                                                                 aGGGAAGAATTTCGCCTTTAGGCCAAATATTTTGCCGCGAGCCGAGAAAAACTCCGGGGAGCATAGCGATAAACGCGATAAAAAACTCCACCACCAGAAAAAAAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTgag                                                                                                                                                      <  2:18679/154‑1 (MQ=255)
                                                                                                                                                    gACAAAATTTGCCTTGAGCCCAAATTTTTTCCCGAGAGTCGAAAAAAAACCCGGGTACCATTCCGGTAAACGCGATAAAAAACTCCCCCCCCAAAAAGATAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGAtt                                                                        <  1:121218/229‑1 (MQ=255)
                                                                                                                                                    gACAAAATTGGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAATTCCGGTGAGCATATCGGTAAACGGGGAAAAAAACTCCGCCACAAGAAAAATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGAtt                                                                        <  1:624631/229‑1 (MQ=255)
                                                                                                                                                      acgaaTTTTCCGTTAAGGCCAATTTTTTTCCGGCAGACGCAAAAAAAACCCGGTGATAATAGGGGAAAGCGCGAAAAAAACCTTCGCCACCAAAAAGAAAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTc                                                                       <  1:473877/226‑1 (MQ=255)
                                                                                                                                                               gCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGAGGTCGCCCAAACCACAGGGGATTTTTTCCCCCAGATCAAGTTTTAAt                                                              >  2:514638/1‑229 (MQ=255)
                                                                                                                                                                    tAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGAGGTCGCCCATACCATAGGCGATTTTTTCCCCCAGATTCACTTTTAATGTCAg                                                         >  1:10415/1‑229 (MQ=255)
                                                                                                                                                                               tttctCGGCGGGCCGAATAAAAACCCGAGGACACATCTGGAAAAGCGGTAAAAAACCTTCCACACAAAAAAGAAAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTAGAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGt                                              <  2:208144/229‑1 (MQ=255)
                                                                                                                                                                                 tttCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCACTACCATAGGCGCTTTTTTCCCGAAGATTCAGTTTTAATTTCCATGGATTATGGTc                                             >  2:475705/1‑228 (MQ=255)
                                                                                                                                                                                  ttCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCa                                            >  1:663621/1‑228 (MQ=255)
                                                                                                                                                                                       cgagtcgagtAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGcc                                                                                                      >  1:69591/1‑165 (MQ=255)
                                                                                                                                                                                       cgagtcgagtAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGAATAAAATTTGAGGAGATACAGAGTACCGATGCAGCGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGAGTATggcgggctct                                       >  1:742409/1‑220 (MQ=255)
                                                                                                                                                                                         agtcgagtAAAAATCCGGTGAGCATATCGGTAAACGCGGCAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGACGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGCCATAAt                                    >  1:661028/1‑229 (MQ=255)
                                                                                                                                                                                           ttgtagaaaaaaTCCGGTGAACATATCGGAAACCTGTGAAAAAAACTTCGCACCAAAAAGATAACCCCGGCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGACATAATCa                                  <  1:173910/222‑1 (MQ=255)
                                                                                                                                                                                           aaaaaaaaaaaaGCCGGGGCGCATATGGTAAACGGGAGAAAAAAACTCCGACACAAGAAAGAAAATGCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGACATAAt                                    <  1:669453/221‑1 (MQ=255)
                                                                                                                                                                                            tcgaaaaaaaccccGGTGACCATAGCGGTAAACGCAGAAAAAAACCTCCCCACCAAAAAGATAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGACATAATCa                                  <  1:520249/217‑1 (MQ=255)
                                                                                                                                                                                               aaaaaaaaaCCGGGGAGACATAGCGGAAAACGGGAAAAAAAACTCCGCCACAAAAAAGATAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGACATAATCACt                                <  1:100180/226‑1 (MQ=255)
                                                                                                                                                                                                aataaaaaaacgggtaggaagtCGGAAAACGCGGGTAAAAACCCCGCACCACAAAAAGAAAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGACATAATCACt                                <  2:74820/208‑1 (MQ=255)
                                                                                                                                                                                                      tcccGGTGATAATCTGGATAACCGGGGTAAAAACCTGCGCCACAAAAAATATATCCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGACATAATCACTCCTtata                         <  2:228056/226‑1 (MQ=255)
                                                                                                                                                                                                        ccGGGGATCATATCGGAAAACGGTGAAAAAAACTCCCCCACCAAAAAGAAAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGACATAATCACTCCTTATAACg                      <  1:150850/229‑1 (MQ=255)
                                                                                                                                                                                                                 cATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGAATAAAATTTGAGGAGGTACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATTGCAGCGGATTATGGTCAGTCATACTCACCCCTtatattgtctcttt               >  2:669452/1‑218 (MQ=255)
                                                                                                                                                                                                                              cgcgGAAAAAAACTTCGCCACCAGAAAGATAATGCCGACATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGACATAATCACTCCTTATAACGCCCCGCAGGGCGTCGTttat  <  1:299121/229‑1 (MQ=255)
                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                 
TACGAGCAGTTCATCAGACTGTATGAAAGTCCAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGACATAATCACTCCTTATAACGCCCCGCAGGGCGTCGTTTAT  >  NC_000913/4065356‑4065802

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: