Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,628,644 A→T C99S (TGT→AGT)  yfgF ← cyclic‑di‑GMP phosphodiesterase, anaerobic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,628,6440AT98.7% 307.6 / NA 76C99S (TGT→AGT) yfgFcyclic‑di‑GMP phosphodiesterase, anaerobic
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base T (53/22);  total (54/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGTAAGGCCTGATAATTGATTAAAGTCCCGAGGTTAAAAGGCATGACACCGACCAGATTTTCTCTGCTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCA‑CCATGTGAAACG‑TTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGAAAAACCTCACCGGAAGAAATATGCGCGAAAGTATTAATGTGGGCGAAAATGGGCAAGATTAA  >  NC_000913/2628421‑2628869
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                 
aGTAAGGCCTGATAATTGATTAAAGTCCCGAGGTTAAAAGGCATGACACCGACCAGATTTTCTCTGCTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCA‑CCATGTGAAACG‑TTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGaa                                                                                                                                                                                                                              >  1:361884/1‑229 (MQ=255)
                tGATTAAAGTCCCGAGGTTAAAAGGCATGACACCGACCAGATTTTCTCTGCTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCA‑CCATGTGAAACG‑TTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATaaa                                                                                                                                                                                                              >  1:264830/1‑229 (MQ=255)
                              aGGTTAAAAGGCATGACACCGACCAGATTTTCTCTGCTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCA‑CCATGTGAAACG‑TTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTgaga                                                                                                                                                                                                                               >  1:163488/1‑198 (MQ=255)
                              aGGATAAAGGGAATGACCCCAACAAGATTTTCTCGCCTCGCAAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAAAAACAGCGTGGCAGAAAACACAATCTGCCAGAATATACGCTGGAAAATTAAACGGGTATCA‑CCATGAGAAACG‑TTTTTACGCCGGGGAGTAAAGCCCCTGTAACCGCCCCAACTgaga                                                                                                                                                                                                                               <  2:163214/198‑1 (MQ=255)
                                   aaaaGGCATGACACCGACCAGATTTTCTCTGCTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCA‑CCATGTGAAACG‑TTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGc                                                                                                                                                                                           >  1:1492967/1‑229 (MQ=255)
                                    aaaGGCATGACACCGACCAGATTTTCTCTGCTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGAATCA‑CCATGTGAAACG‑TTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCg                                                                                                                                                                                          >  1:793462/1‑229 (MQ=255)
                                           taacagcaaCCAAATTTTCTCTCCGCGCAAGTAATCCAAAAAGAGCAGAAAAAGGAAAAGGTATGAAAAAAAGCGTGGCAGAAAACAAAATCGGCAAGAATACGCGCTGGAAAATAAAAGGGGTACCC‑CCAGGGAAAAGG‑TTTTTCCGCGGGGGATAAAAGCCCTTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGttt                                                                                                                                                                                       <  2:1969567/217‑1 (MQ=255)
                                           tGACACCGACCAGATTTTCTCTGCTCGCCAGTAATCCTACCAATCCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCA‑CCATGTGCAACG‑TTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGt                                                                                                                                                                                                     >  2:480930/1‑211 (MQ=255)
                                           tGACACCGACCAGATTTTCTCTGCTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCA‑CCATGTGAAACG‑TTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTgaga                                                                                                                                                                                                                               >  1:765913/1‑185 (MQ=255)
                                           tGACACCGACCAGATTTTCTCTGCTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCA‑CCATGTGAAACG‑TTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGttt                                                                                                                                                                                       >  1:1972963/1‑225 (MQ=255)
                                               acCGACCAGATTTTCTCTGCTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCA‑CCATGTGAAACG‑TTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTa                                                                                                                                                                                   >  1:926470/1‑225 (MQ=255)
                                                 cGACCAGATTTTCTCTGCTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCA‑CCATGTGAACCG‑TTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATGAAATTAGCCGTTAATGATAGCGTTTCGTAAAAAcc                                                                                                                                                                             >  1:1001800/1‑229 (MQ=255)
                                                  gACCAGATTTTCTCTGCTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGAATCA‑CCATGTGAAACG‑TTGTTAAGCCGGGGAGTAAAGACACTGTAACCGCCCCAACAGAGAATGATAATGATAATCAAATTaa                                                                                                                                                                                                         >  1:232739/1‑199 (MQ=255)
                                                      aGATTTCCCCCGCCCACCGGTAATCCTAAAAATGCAGCAACCGGAAAAAGTACAAAAACCAGCGTTGAAGAAAACACAATCTGCCAGAATATCCGCTGGAAAATAAAACGGTTATCC‑CCATGTAAAACT‑TTGTTACGCCGGGGAAAAAACCCCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGct                                                                                                                                                                        <  2:1490393/229‑1 (MQ=255)
                                                                            aTCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCA‑CCATGTGAAACG‑TTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAAt                                                                                                                                                  >  2:1725080/1‑229 (MQ=255)
                                                                                tACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCA‑CCATGTGAAACG‑TTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATccc                                                                                                                                               >  2:236082/1‑228 (MQ=255)
                                                                                     aTGCAGAAAACTAAAAAGTTACAAAAAACAGCGTTGCAGAAACCACAATCTGCAAGAATATACGGTGGAAAATAAACCGGTTATCC‑CCATGTAAAACT‑TTGTTCCGCCGGGGAGTAAACACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAgg                                                                                                                                          <  2:792073/228‑1 (MQ=255)
                                                                                        cagcaAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCA‑CCATGTGAAACG‑TTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGa                                                                                                                                       >  1:844343/1‑228 (MQ=255)
                                                                                           caAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCA‑CCATGTGAAACG‑TTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATccc                                                                                                                                   >  1:1892493/1‑229 (MQ=255)
                                                                                           caAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCA‑CCATGTGAAACG‑TTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATcc                                                                                                                                    >  1:1329421/1‑228 (MQ=255)
                                                                                                        tATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCA‑CCATGTGAAACG‑TTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATgccgcc                                                                                                                      >  2:1227740/1‑229 (MQ=255)
                                                                                                            aGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCA‑CCATGTGAAACG‑TTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCACAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAAACCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATTCCGCCCAGc                                                                                                                  >  2:1961488/1‑229 (MQ=255)
                                                                                                                              aaCACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCA‑CCATGTGAAACG‑TTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAAcatcatc                                                                                                    >  2:1798998/1‑225 (MQ=255)
                                                                                                                                       tgcaataattgaCGCTGGGAAATAAAAC‑GGTAACACCCATAGGAAACG‑TTGTTACGCGGGGAAGAAAACACCCGTTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGatcat                                                                                         <  2:264371/216‑1 (MQ=255)
                                                                                                                                             aaTATACGCTGGGAAATTAAACGGGTATCA‑CCATGTGAAACG‑TTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTaa                                                                                                                                                                                  >  1:1509225/1‑132 (MQ=255)
                                                                                                                                             aaTATACGCGGGAAAATAAAAGGGGTACCC‑CCTTGT‑AAACGTTTTTTCCGCGGGGGAGAAAAGCCCCTGAAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTaa                                                                                                                                                                                  <  2:1506612/132‑1 (MQ=255)
                                                                                                                                                ataCGTGGGAAAATTAAACGGTAATCC‑CCAGGGAAAACT‑TTGTTACGCGGGGAAGTAAAGCCCGTTTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTg                                                                                              <  1:1019154/213‑1 (MQ=255)
                                                                                                                                                 taCGCTGGGAAATTAAACGGGTATCA‑CCATGTGAAACG‑TTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAgg                                                                             >  2:1032803/1‑229 (MQ=255)
                                                                                                                                                  aCGCTGGGAAATTAAACGGGTATCA‑CCATGTGAAACG‑TTGTAACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATACTAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAAcc                                                                                                                                                                             >  1:1120635/1‑132 (MQ=255)
                                                                                                                                                       gggAAATTAAACGGGTATCA‑CCATGTGAAACG‑TTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATa                                                                       >  1:1219840/1‑229 (MQ=255)
                                                                                                                                                       gggAAATTAAACGGGTATCA‑CCATGTGAAACG‑TTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATa                                                                       >  2:1807942/1‑229 (MQ=255)
                                                                                                                                                       gggAAATTAAACGGGTATCA‑CCATGTGAAACG‑TTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATa                                                                       >  2:185664/1‑229 (MQ=255)
                                                                                                                                                        ggAAATTAAACGGGTATCA‑CCATGTGAAACG‑TTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTgaga                                                                                                                                                                                                                               >  1:821988/1‑76 (MQ=255)
                                                                                                                                                                aaCGGGTATCA‑CCATGTGAAACG‑TTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGaaa                                                              >  1:978581/1‑229 (MQ=255)
                                                                                                                                                                    ggTACCA‑CCAAGTGAATCG‑TGTTGCCGCGGGGGATAAAAGACCGTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGACAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGAAAAAc                                                           <  2:1613385/228‑1 (MQ=255)
                                                                                                                                                                                gtgCAACG‑CTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCATAGGCAGATAAAAGAGACAAACCTCACCGgaagc                                               >  2:1285272/1‑228 (MQ=255)
                                                                                                                                                                                gtgAAACG‑TTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACtga                                                                                                                                                                                                                                 >  1:1006321/1‑51 (MQ=255)
                                                                                                                                                                                   aaaCG‑TTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGAAAAACCTCACCGGAAGAAat                                            >  1:203626/1‑229 (MQ=255)
                                                                                                                                                                                    aaCG‑TTTTTACGCGGGGGAGAAAAGACCCTTTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGct                                                                                                                                                                        <  2:609568/104‑1 (MQ=255)
                                                                                                                                                                                    aaCG‑TTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGct                                                                                                                                                                        >  1:610638/1‑104 (MQ=255)
                                                                                                                                                                                       g‑ttgttACGCGGGGGAAAAAAGCCCTTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCa                                                                                                                                                                           <  2:140837/98‑1 (MQ=255)
                                                                                                                                                                                       g‑ttgttACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCa                                                                                                                                                                           >  1:141074/1‑98 (MQ=255)
                                                                                                                                                                                       g‑ttgttACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTCAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGAAAAACCTCACCGGAAGAAATATgc                                        >  1:1392066/1‑229 (MQ=255)
                                                                                                                                                                                          tgttACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGAAAAACCTCACCGGAAGAAataatcgc                                      >  2:1581146/1‑224 (MQ=255)
                                                                                                                                                                                             tACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGAAAAACCTCACCGGAAGAAATATGCGCGa                                    <  1:1728026/228‑1 (MQ=255)
                                                                                                                                                                                                   ggggAGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGAAAAACCTCACCGGAAGAAATATGCGCGAAAGTAtt                             >  1:207304/1‑229 (MQ=255)
                                                                                                                                                                                                            aGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGcc                                                                                                                                                                                                       >  1:1264325/1‑50 (MQ=255)
                                                                                                                                                                                                              aCCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGttt                                                                                                                                                                                       >  1:278684/1‑64 (MQ=255)
                                                                                                                                                                                                              aCCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGttt                                                                                                                                                                                       <  2:278192/64‑1 (MQ=255)
                                                                                                                                                                                                                ccTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGAAAAACCTCACCGGAAGAAATATGCGCGAAAGTATTAATGTGGGCGaaa                >  2:446868/1‑229 (MQ=255)
                                                                                                                                                                                                                 cTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGAAAAACCTCACCGGAAGAAATATGCGCGa                                    <  1:293149/208‑1 (MQ=255)
                                                                                                                                                                                                                 cTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGAAAAACCTCACCGGAAGAAATATGCGCGa                                    >  2:292625/1‑208 (MQ=255)
                                                                                                                                                                                                                 cTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGAAAAACCTCACCGGAAGAAATATGCGCGAAAGTATTAATGTGGGCGaaaa               >  1:1177790/1‑229 (MQ=255)
                                                                                                                                                                                                                   gTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCg                                                                                                                                        <  1:33486/106‑1 (MQ=255)
                                                                                                                                                                                                                   gTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCg                                                                                                                                        >  2:33432/1‑106 (MQ=255)
                                                                                                                                                                                                                   gTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGAAAAAc                                                           <  2:384924/183‑1 (MQ=255)
                                                                                                                                                                                                                   gTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGAAAAAc                                                           >  1:385617/1‑183 (MQ=255)
                                                                                                                                                                                                                   gTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGAAAAACCTCACCGGACGAAAAATACGCGAAAGTATTACAGTGGTCAAAAATg             >  2:1184033/1‑229 (MQ=255)
                                                                                                                                                                                                                   gTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGAAAAACCTCACCGGAAGAAATATGCGCGAAAGTATTAATGTGGGCGAAAATg             >  2:1348181/1‑229 (MQ=255)
                                                                                                                                                                                                                   gTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGAAAAACCTCACCGGAAGAAATATGCGCGAAAGTATTAATGTGGGCGAAAATg             >  2:777271/1‑229 (MQ=255)
                                                                                                                                                                                                                   gTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGAAAAACCTCACCGGAAGAAATATGCGCGAAAGTATTAATGTGGGCGAAAATg             >  2:699940/1‑229 (MQ=255)
                                                                                                                                                                                                                   gTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGAAAAACCTCACCGGAAGAAATATGCGCGAAAGTATTAATGTGGGCGAAAATg             >  2:556776/1‑229 (MQ=255)
                                                                                                                                                                                                                      aCCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGAAAAACCTCACCGGAAGAAATATGCGCGAAAGTATTAATGTGGGCGAAAATGGGc          <  2:1175743/229‑1 (MQ=255)
                                                                                                                                                                                                                          ccccAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGaaa                                                              <  2:1840610/173‑1 (MQ=255)
                                                                                                                                                                                                                          ccccAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGAAAAACCTCACCGGAAGAAATATGCGCGa                                    >  1:1405958/1‑199 (MQ=255)
                                                                                                                                                                                                                          ccccAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGAAAAACCTCACCGGAAGAAATATGCGCGa                                    <  2:1403539/199‑1 (MQ=255)
                                                                                                                                                                                                                          ccccAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGAAAAACCTCACCGGAAGAAATATGCGCGAAAGTATTAATGTGGGCGAAAATGGGCAAg       >  1:291662/1‑228 (MQ=255)
                                                                                                                                                                                                                          ccccAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGAAAAACCTCACCGGAAGAAATATGCGCGAAAGTATTAATGTGGGCGAAAATGGGCAAg       >  2:1605102/1‑228 (MQ=255)
                                                                                                                                                                                                                          ccccAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCACGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGaaa                                                              >  1:1843743/1‑173 (MQ=255)
                                                                                                                                                                                                                           cccAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAAGCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGAAAAACCTCACCGGAAGAAATATGCGCGAAAGTATTAATGTGGGCGAAAATGGGCAAGAt     <  2:842865/229‑1 (MQ=255)
                                                                                                                                                                                                                            ccAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGAAAAACCTCACCGGAAGAAATATGCGCGAAAGTATTAATGTGGGCGAAAATggg           <  1:1661770/222‑1 (MQ=255)
                                                                                                                                                                                                                            ccAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGAAAAACCTCACCGGAAGAAATATGCGCGAAAGTATTAATGTGGGCGAAAATggg           >  2:1658939/1‑222 (MQ=255)
                                                                                                                                                                                                                              aaCTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAAcatcatca                                                                                                   >  1:838848/1‑132 (MQ=255)
                                                                                                                                                                                                                              aaCTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAAcatcatca                                                                                                   <  2:837380/132‑1 (MQ=255)
                                                                                                                                                                                                                              aaCTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGAAAAACCTCACCGGAAGAAATATGCGCGa                                    >  1:493110/1‑195 (MQ=255)
                                                                                                                                                                                                                              aaCTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGAAAAACCTCACCGGAAGAAATATGCGCGa                                    <  2:492241/195‑1 (MQ=255)
                                                                                                                                                                                                                              aaCTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGAAAAACCTCACCGGAAGAAATATGCGCGAAAGTATTAATGTGGGCGAAAATGGGCAAGATTaa  >  2:978271/1‑229 (MQ=255)
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                 
AGTAAGGCCTGATAATTGATTAAAGTCCCGAGGTTAAAAGGCATGACACCGACCAGATTTTCTCTGCTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCA‑CCATGTGAAACG‑TTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGAAAAACCTCACCGGAAGAAATATGCGCGAAAGTATTAATGTGGGCGAAAATGGGCAAGATTAA  >  NC_000913/2628421‑2628869

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