Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,065,580 T→A 13.4% 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,5800TA13.4% 125.2 / 2.1 64T63S (ACC→TCC) yihOputative sulphoquinovose importer
Reads supporting (aligned to +/- strand):  ref base T (43/8);  new base A (2/7);  total (49/15)
Fisher's exact test for biased strand distribution p-value = 4.51e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TACGAGCAGTTCATCAGACTGTATGAAAGTCCAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCG‑AGTAAAAATCCGGTGAG‑CATATCGGTAAACGCG‑GTAAAAAA‑CTTCGC‑CACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGACATAATCACTCCTTATAACGCCCCGCAGGGCGTCGTTTAT  >  NC_000913/4065357‑4065803
                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                 
tACGAGCAGTTCATCAGACTGTATGAAAGTCCAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCG‑GGTAAAAAACCGGTGAT‑CATATCGTTAAACGCC‑GGaaaaa                                                                                                                                                                                                                              >  2:624630/1‑229 (MQ=255)
      cAGTTCATCAGACTGTATGAAAGTCCAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCG‑AGTAAAAATCCGGTGAG‑CATATCGGTAAACGCG‑GTAAAAAA‑CTTCg                                                                                                                                                                                                                       >  2:166368/1‑229 (MQ=255)
      cAGTTCATCAGACTGTATGAAAGTCCAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCG‑AGTAAAAATCCGGGGAG‑CATATCGGTAAACGAG‑GTAAAAAA‑ATTCg                                                                                                                                                                                                                       >  2:121218/1‑229 (MQ=255)
             acaGACTTTATGAAAGCCCAAACAACAATAACGGCGGGGTGGCAACCGGCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCG‑AGTAAAAATCCGGTGAG‑CATATCGGTAAACGCG‑GTAAAAAA‑CTTCGC‑CACCAg                                                                                                                                                                                                               <  2:139625/228‑1 (MQ=255)
                 aCTGTATGAAAGTCCAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCAAGTCG‑AGTAAAATTCCGGGGTG‑CATCTCGGTAAACGCG‑GTAAAAAA‑CTTAGC‑CACCATaaa                                                                                                                                                                                                            >  1:294231/1‑228 (MQ=255)
                   tGTATGAAAGTCCAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCG‑AGTAAAAATCCAGTGGG‑CATATCGGTAAACGCG‑GTAAAAAA‑CTTCGC‑CACCAGAAAGaa                                                                                                                                                                                                         >  2:473877/1‑228 (MQ=255)
                     tATGAAAGTCCAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCG‑AGTAAAAATCCGGTGAG‑CATATCGGTAAACGCG‑GTAAAAAA‑CTTCGC‑AACCAGAAAGATaa                                                                                                                                                                                                       >  1:376814/1‑229 (MQ=255)
                     tATGAAAGTCCAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAACATTTTTCCGCGAGTCA‑GGTAAAAATCCGGTGGG‑CATATCGGTAAACGCG‑GTAAAACA‑CTTCGc‑cacacaaaaaaaaa                                                                                                                                                                                                       >  2:46373/1‑218 (MQ=255)
                     tATGAAAGTCCAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCACATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCG‑AGTAAAAATCCGGTGAG‑CAAATCGGTAAACGCG‑GTAAAAAA‑CTTCG‑‑CACCCGAAAGAAAAt                                                                                                                                                                                                      >  2:173910/1‑229 (MQ=255)
                     tATGAAAGTCCAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTAGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCG‑AGTAAAAATCCGGTGAG‑CATATCGGTAAACGCG‑GTAAAAAA‑CTTCGC‑AACCGGAAAGATaa                                                                                                                                                                                                       >  1:646987/1‑229 (MQ=255)
                       tGAAAGTCCAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCG‑AGTAAAAATCCGGTGAG‑CATATCGGTAAACGCG‑GTAAAAAA‑CTTCGC‑CACCAGAAAGATAATc                                                                                                                                                                                                     >  1:496971/1‑229 (MQ=255)
                           aGTCCAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCG‑AGTAAAAATCCGGTGAG‑CATATCGGTAAACGCG‑GTAAAAAA‑CTTCGC‑CACCAGAAAGATAATCCCGc                                                                                                                                                                                                 >  1:411094/1‑229 (MQ=255)
                               cAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGCTAAACTGAGGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCAAATATTTTTCCCCAGGTCG‑CGGAAAAATCCGGAGAG‑CAAAGCGTTAACCGAG‑TAAAAAAA‑CTTCGCACACCA‑AAAGAAAAACCCGGCaaa                                                                                                                                                                                             >  2:520249/1‑227 (MQ=255)
                               cAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCCATATTTTTCCGCGAGTCA‑AGTAAAAATACGGTGAG‑CATATCGGTAAACCCG‑GGACAAAA‑CTTCGC‑CACCAGAACGTTAATCCCGCCATa                                                                                                                                                                                             >  1:379959/1‑229 (MQ=255)
                               cAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCG‑AGTAAAAATCCGGTGGG‑CATATCGGTAAACGCG‑GTCAAAAA‑CTTCGC‑CACCAGAAAAATAATCCAGCCATa                                                                                                                                                                                             >  1:208144/1‑229 (MQ=255)
                               cAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCG‑AGTAAAAATCCGGTGAG‑CATATCGGTAAACGCG‑GTAAAAAA‑CTTCGC‑CACCAGAAAGATAATCCCGcccaa                                                                                                                                                                                             >  1:462618/1‑226 (MQ=255)
                               cAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCG‑AGTAAAAATCCGGTGAG‑CATATCGGTAAACGCG‑GTAAAAAA‑CTTCGC‑CACCAGAAAGATAATCCCGCCATa                                                                                                                                                                                             >  1:254981/1‑229 (MQ=255)
                               cAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCG‑AGTAAAAATCCGGTGAG‑CATATCGGTAAACGCG‑GTAAAAAA‑CTTCGC‑CACCAGAAAGATAATCCCGCCATa                                                                                                                                                                                             >  1:237646/1‑229 (MQ=255)
                               cAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCG‑AGTAAAAATCCGGTGAG‑CAGATCGGTAAACGCG‑GAAAAAAA‑CTTCGC‑CACAAGAAAGAAAAAccccacaaa                                                                                                                                                                                             >  2:99974/1‑223 (MQ=255)
                                              gcgcGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCG‑ATTAAAAATCCGGTGAG‑CATATCGGGAAACTCG‑GTAAAAAA‑CTTCGC‑CACCAGAAAGATAATCCCGCCATAGAAAGAAGGCATccc                                                                                                                                                                              >  1:228056/1‑229 (MQ=255)
                                              gcgcGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCG‑AGTAAAAATCCGGTGAG‑CATATCGGTAAACGCG‑GTAAAAAA‑CTTCGC‑AACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATccc                                                                                                                                                                              >  1:549315/1‑229 (MQ=255)
                                              gcgcGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCG‑AGTAAAAATCCGGTGAG‑CATATCGGGAAACGCG‑GTAAAAAA‑ATTCGC‑CACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATccc                                                                                                                                                                              >  1:689746/1‑229 (MQ=255)
                                              gcgcGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGCTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCG‑AGTAAAAATCCGGTGAG‑CATATCGGTAAACGCG‑GTAAAAAA‑CTTCGC‑CACAGGAAAGATAATCCCGCCTCAGTAAGAAGGAAcccc                                                                                                                                                                              >  1:727718/1‑229 (MQ=255)
                                               cgcgGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCG‑AGTAAAAATCCGGTGAG‑CATATCGGTAAACGCG‑GTAAAAAA‑CTTCGC‑CACCAGAAAGATAATACCGCCATAGAAAGCAGGCATCCCt                                                                                                                                                                             >  1:74820/1‑229 (MQ=255)
                                               cgcgGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCG‑AGTAAAAATCCGGGGAG‑CATATCGGTAAACGCG‑GTAAAAAA‑CTTCGC‑CAACAGAAAGATAATCCGGCCATAGTAAGAAGGCAcccca                                                                                                                                                                             >  2:42123/1‑228 (MQ=255)
                                                gcgGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCG‑AGTAAAAATCCGGTGAG‑CATATCGGTAAACGCG‑GCAATAAA‑CTTCGC‑CACCCGAAAGATAATCCCGCCATAGCAAGCATGCATCCCTa                                                                                                                                                                            >  1:191109/1‑229 (MQ=255)
                                                       taagcgtcttTTAAACCGGAAAAAAAAAGGTGGTGGCAGAAAATTGGAGGGTTAAAATTACCGCCTCCGGAACCGTTTAATAAAAAATAAAGGGACAAATTTTGCTTTTAGTCCCAATATTTTTCCGAGAGTCG‑AGTAAAAACCCGGTGAG‑CATATCGGTAAACGCG‑GTAAAAAA‑CTTCGC‑CACCAGAAAAATAATCCCGCTATAGTAAGCAGGCATCCCTAATTCATc                                                                                                                                                                     <  2:197675/225‑1 (MQ=255)
                                                                       cGGAAAAAAAAAGGGGGGGGAGAAAAAATGGGGGGGTGCAATAAACCCCGCAGAAACCGCTGAAAATAAAAGAAAGAGAAAAAATTTGATTTTGGGACAAATATTTTCCGGGGAGGGG‑AGAAAAAAGCTGGGGAG‑CATGCGGGAAAACGCG‑GTAAAAAA‑CTCCCC‑CACAAAAAAGTAAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTgagg                                                                                                                                                     <  2:599472/229‑1 (MQ=255)
                                                                          caaacaaaAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCG‑AGTAAAAATCCGGTGAG‑CATATCGGTAAACGCG‑GTAAAAAA‑CTTCGC‑CACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGAATAAAATTTgaggga                                                                                                                                                   >  2:516409/1‑226 (MQ=255)
                                                                                       gtggCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCG‑AGTAAAAATCCGGTGAG‑CATATCGGTAAACGCG‑GTAAAAAA‑CTTCGC‑CACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGAGGATACAGAGAACCGa                                                                                                                                     >  1:606870/1‑229 (MQ=255)
                                                                                        tggCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCG‑AGTAAAAATCCGGTGAG‑CATATCGGTAAACGCG‑GTAAAAAA‑CTTCGC‑CACCAGAAAGATAATCCCGCATAAGTAAGCAGGCATACCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGAt                                                                                                                                    >  1:399985/1‑229 (MQ=255)
                                                                                                       gCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCG‑AGTAAAAATCCGGTGAG‑CATATCGGTAAACGCG‑GTAAAAAA‑CTTCGC‑CACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATAACGAGTACCGATGCAGAGCATTAAAtt                                                                                                                     >  1:348087/1‑229 (MQ=255)
                                                                                                          ttGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCG‑AGTAAAAATCCGGTGAG‑CATATCGGTAAACGCG‑GTAAAAAA‑CTTCGC‑CACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGATGAGATACAGAGTACCGATGCAGAGCATTAAATTCg                                                                                                                   >  1:438582/1‑228 (MQ=255)
                                                                                                           tGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCG‑AGTAAAAATCCGGTGAG‑CATATCGGTAAACGCG‑GTAAAAAA‑CTTCGC‑CACCAGAAAGATAATCCCGCCATAGTAAGCAGGCACCCCTAATTCATCGGTATAAAATTTGAGGGAATACAGAGTACCCATCCAAAGAATAAAATTCaaa                                                                                                                 >  2:299121/1‑229 (MQ=255)
                                                                                                             cAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCG‑AGTAAAAATCCGGTGAG‑CATATCGGTAAACGCG‑GTAAAAAA‑CTTCGC‑CACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATCCATCGGTATAAAATTTTAGGAAAAACAGAGGACCGATCAAAAGCATTAAATTACAAcc                                                                                                               >  2:100180/1‑229 (MQ=255)
                                                                                                                  accgccgcCGGAACCGATGTATATAAAATAAAGGGACGAATTTTGCTTTTAGGTCAAATTTTTTACCGAGAGTCG‑AGTAAAAACCGGGAGAC‑CATACGGGAAAAGGAG‑GTAAAAAA‑CTTCGC‑CACCAAAAATAAAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACg                                                                                                           <  2:720422/226‑1 (MQ=255)
                                                                                                                        gccgGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCG‑AGTAAAAATCGGGTGAG‑CATATCCGTAAACGCG‑GTAAAAAA‑CATCGC‑CACCAGAAAGATAATCCCGACATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGCGGAGAAACACTGTCCCCATTCAGAGCATTAAATTCAAACCGACGTGGCcca                                                                                                    >  2:71198/1‑229 (MQ=255)
                                                                                                                         ccgGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCG‑AGTAAAAATCCGGTGAG‑CATATCGGTAAACGCG‑GTAAAAAA‑CTTCGC‑CACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCCTTAAATTCGAACCGACGTCGCccat                                                                                                   >  1:396688/1‑229 (MQ=255)
                                                                                                                                   tGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCACAATATTTTTCCGCGAGTTG‑AGTAACAATCCGGTGAG‑CATATCGGTAAACGCG‑GTAAAAAA‑CTTCGC‑CACCAGAAAGATAATCCCGCCATAGTAAGCAGACTTCCCTAATTCATCGGTCTAAAATTTGAGTCGATACAGAGTACCGATCCAGATCATTAAATTCGTACCGACGCAGTCCATACCTTCGGc                                                                                          >  1:575784/1‑228 (MQ=255)
                                                                                                                                            aataaaGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCG‑AGTAAAAATCCGGTGAG‑CATATCGGTAAACGCG‑GTAAAAAA‑CTTCGC‑CACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTcc                                                                                >  1:457959/1‑229 (MQ=255)
                                                                                                                                             ataaaGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCG‑AGTAAAAATCCGGTGAG‑CATATCGGTAAACGCG‑GTAAAAAA‑CTTCGC‑CACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCACTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGAGTTTTTccc                                                                               >  1:676284/1‑229 (MQ=255)
                                                                                                                                               aaaGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCG‑AGTAAAAATCCGGTGAG‑CATATCGGTAAACGCG‑GTAAAAAA‑CTTCGC‑CACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATAACAAAGGCGATTTTTTccccc                                                                             >  2:461292/1‑227 (MQ=255)
                                                                                                                                               aaaGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCA‑AGTAAAAATCCGGTGAG‑CATATCGGTAAACGCG‑GTAAAAAA‑CTTCGC‑CACCAGAAAGATAATCCCGCCATAGTAAGAAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATCCGAACCGAAGTCGCCCATACCCTAGCGGATTTTTTCACGc                                                                             >  1:612182/1‑229 (MQ=255)
                                                                                                                                                 aGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCG‑AGTAAAAATCCGGTGAG‑CATATCGGTAAACGCG‑GTAAAAAA‑CTTCGC‑CACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTgag                                                                                                                                                      >  1:18679/1‑154 (MQ=255)
                                                                                                                                                 aGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCG‑AGTAAAAATCCGGTGAG‑CATATCGGTAAACGCG‑GTAAAAAA‑CTTCGC‑CACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAAATTACCGATGCAAAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTccccccg                                                                           >  1:462331/1‑225 (MQ=255)
                                                                                                                                                 aGGGAAGAATTTCGCCTTTAGGCCAAATATTTTGCCGCGAGCCG‑AGAAAAACTCCGGGGAG‑CATAGCGATAAACGCG‑ATAAAAAA‑CTCCAC‑CACCAGAAAAAAAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTgag                                                                                                                                                      <  2:18679/154‑1 (MQ=255)
                                                                                                                                                    gACAAAATTTGCCTTGAGCCCAAATTTTTTCCCGAGAGTCG‑AAAAAAAACCCGGGTAC‑CATTCCGGTAAACGCG‑ATAAAAAA‑CTCCCC‑CCCCAAAAAGATAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGAtt                                                                        <  1:121218/229‑1 (MQ=255)
                                                                                                                                                    gACAAAATTGGCCTTTAGGCCCAATATTTTTCCGCGAGTCG‑AGTAAAATTCCGGTGAG‑CATATCGGTAAACGGG‑GAAAAAAA‑CTCCGC‑CACAAGAAAAATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGAtt                                                                        <  1:624631/229‑1 (MQ=255)
                                                                                                                                                      acgaaTTTTCCGTTAAGGCCAATTTTTTTCCGGCAGACGCAAAAAAAACCCGGTGAT‑AATAGGGGAAAGCGCG‑AAAAAAAC‑CTTCGC‑CACCAAAAAGAAAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTc                                                                       <  1:473877/226‑1 (MQ=255)
                                                                                                                                                              gCCTTTAGGCCCAATATTTTTCCGCGAGTCG‑AGTAAAAATCCGGTGAG‑CATATCGGTAAACGCG‑GTAAAAAA‑CTTCGC‑CACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGAGGTCGCCCAAACCACAGGGGATTTTTTCCCCCAGATCAAGTTTTAAt                                                              >  2:514638/1‑229 (MQ=255)
                                                                                                                                                                   tAGGCCCAATATTTTTCCGCGAGTCG‑AGTAAAAATCCGGTGAG‑CATATCGGTAAACGCG‑GTAAAAAA‑CTTCGC‑CACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGAGGTCGCCCATACCATAGGCGATTTTTTCCCCCAGATTCACTTTTAATGTCAg                                                         >  1:10415/1‑229 (MQ=255)
                                                                                                                                                                              tttctCGGCGGGCCG‑AATAAAAACCCGAGGA‑‑CACATCTGGAAAAGCG‑GTAAAAAACCTTCCA‑CACAAAAAAGAAAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTAGAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGt                                              <  2:208144/229‑1 (MQ=255)
                                                                                                                                                                                tttCCGCGAGTCG‑AGTAAAAATCCGGTGAG‑CATATCGGTAAACGCG‑GTAAAAAA‑CTTCGC‑CACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCACTACCATAGGCGCTTTTTTCCCGAAGATTCAGTTTTAATTTCCATGGATTATGGTc                                             >  2:475705/1‑228 (MQ=255)
                                                                                                                                                                                 ttCCGCGAGTCG‑AGTAAAAATCCGGTGAG‑CATATCGGTAAACGCG‑GTAAAAAA‑CTTCGC‑CACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCa                                            >  1:663621/1‑228 (MQ=255)
                                                                                                                                                                                      cgagtcg‑agtAAAAATCCGGTGAG‑CATATCGGTAAACGCG‑GTAAAAAA‑CTTCGC‑CACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGcc                                                                                                      >  1:69591/1‑165 (MQ=255)
                                                                                                                                                                                      cgagtcg‑agtAAAAATCCGGTGAG‑CATATCGGTAAACGCG‑GTAAAAAA‑CTTCGC‑CACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGAATAAAATTTGAGGAGATACAGAGTACCGATGCAGCGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGAGTATggcgggctct                                       >  1:742409/1‑220 (MQ=255)
                                                                                                                                                                                        agtcg‑agtAAAAATCCGGTGAG‑CATATCGGTAAACGCG‑GCAAAAAA‑CTTCGC‑CACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGACGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGCCATAAt                                    >  1:661028/1‑229 (MQ=255)
                                                                                                                                                                                          ttgtagaaaaaaTCCGGTGAA‑CATATCGGAAACCTGT‑GAAAAAAA‑CTTCG‑‑CACCAAAAAGATAACCCCGGCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGACATAATCa                                  <  1:173910/222‑1 (MQ=255)
                                                                                                                                                                                          aaaaaaaaaaaaGCCGGGGCG‑CATAT‑GGTAAACGGGAGAAAAAAA‑CTCCGA‑CACAAGAAAGAAAATGCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGACATAAt                                    <  1:669453/221‑1 (MQ=255)
                                                                                                                                                                                              aaaaaaaaaCCGGGGAGACATAGCGGAAAACGGG‑AAAAAAAA‑CTCCGC‑CACAAAAAAGATAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGACATAATCACt                                <  1:100180/226‑1 (MQ=255)
                                                                                                                                                                                                     tcccGGTGAT‑AATCTGGATAACCGGG‑GTAAAAAC‑CTGCGC‑CACAAAAAATATATCCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGACATAATCACTCCTtata                         <  2:228056/226‑1 (MQ=255)
                                                                                                                                                                                                       ccGGGGAT‑CATATCGGAAAACGGT‑GAAAAAAA‑CTCCCC‑CACCAAAAAGAAAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGACATAATCACTCCTTATAACg                      <  1:150850/229‑1 (MQ=255)
                                                                                                                                                                                                                cATATCGGTAAACGCG‑GTAAAAAA‑CTTCGC‑CACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGAATAAAATTTGAGGAGGTACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATTGCAGCGGATTATGGTCAGTCATACTCACCCCTtatattgtctcttt               >  2:669452/1‑218 (MQ=255)
                                                                                                                                                                                                                            cgcg‑GAAAAAAA‑CTTCGC‑CACCAGAAAGATAATGCCGACATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGACATAATCACTCCTTATAACGCCCCGCAGGGCGTCGTttat  <  1:299121/229‑1 (MQ=255)
                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                 
TACGAGCAGTTCATCAGACTGTATGAAAGTCCAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCG‑AGTAAAAATCCGGTGAG‑CATATCGGTAAACGCG‑GTAAAAAA‑CTTCGC‑CACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGACATAATCACTCCTTATAACGCCCCGCAGGGCGTCGTTTAT  >  NC_000913/4065357‑4065803

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