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

AGTAAGGCCTGATAATTGATTAAAGTCCCGAGGTTAAAAGGCATGACACCGACCAGATTTTCTCTGCTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGAAAAACCTCACCGGAAGAAATATGCGCGAAAGTATTAATGTGGGCGAAAATGGGCAAGATTAA  >  NC_000913/2628421‑2628869
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                 
aGTAAGGCCTGATAATTGATTAAAGTCCCGAGGTTAAAAGGCATGACACCGACCAGATTTTCTCTGCTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGaa                                                                                                                                                                                                                              >  1:361884/1‑229 (MQ=255)
                tGATTAAAGTCCCGAGGTTAAAAGGCATGACACCGACCAGATTTTCTCTGCTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATaaa                                                                                                                                                                                                              >  1:264830/1‑229 (MQ=255)
                              aGGTTAAAAGGCATGACACCGACCAGATTTTCTCTGCTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTgaga                                                                                                                                                                                                                               >  1:163488/1‑198 (MQ=255)
                              aGGATAAAGGGAATGACCCCAACAAGATTTTCTCGCCTCGCAAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAAAAACAGCGTGGCAGAAAACACAATCTGCCAGAATATACGCTGGAAAATTAAACGGGTATCACCATGAGAAACGTTTTTACGCCGGGGAGTAAAGCCCCTGTAACCGCCCCAACTgaga                                                                                                                                                                                                                               <  2:163214/198‑1 (MQ=255)
                                   aaaaGGCATGACACCGACCAGATTTTCTCTGCTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGc                                                                                                                                                                                           >  1:1492967/1‑229 (MQ=255)
                                    aaaGGCATGACACCGACCAGATTTTCTCTGCTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGAATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCg                                                                                                                                                                                          >  1:793462/1‑229 (MQ=255)
                                           taacagcaaCCAAATTTTCTCTCCGCGCAAGTAATCCAAAAAGAGCAGAAAAAGGAAAAGGTATGAAAAAAAGCGTGGCAGAAAACAAAATCGGCAAGAATACGCGCTGGAAAATAAAAGGGGTACCCCCAGGGAAAAGGTTTTTCCGCGGGGGATAAAAGCCCTTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGttt                                                                                                                                                                                       <  2:1969567/217‑1 (MQ=255)
                                           tGACACCGACCAGATTTTCTCTGCTCGCCAGTAATCCTACCAATCCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGCAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGt                                                                                                                                                                                                     >  2:480930/1‑211 (MQ=255)
                                           tGACACCGACCAGATTTTCTCTGCTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTgaga                                                                                                                                                                                                                               >  1:765913/1‑185 (MQ=255)
                                           tGACACCGACCAGATTTTCTCTGCTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGttt                                                                                                                                                                                       >  1:1972963/1‑225 (MQ=255)
                                               acCGACCAGATTTTCTCTGCTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTa                                                                                                                                                                                   >  1:926470/1‑225 (MQ=255)
                                                 cGACCAGATTTTCTCTGCTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAACCGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATGAAATTAGCCGTTAATGATAGCGTTTCGTAAAAAcc                                                                                                                                                                             >  1:1001800/1‑229 (MQ=255)
                                                  gACCAGATTTTCTCTGCTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGAATCACCATGTGAAACGTTGTTAAGCCGGGGAGTAAAGACACTGTAACCGCCCCAACAGAGAATGATAATGATAATCAAATTaa                                                                                                                                                                                                         >  1:232739/1‑199 (MQ=255)
                                                      aGATTTCCCCCGCCCACCGGTAATCCTAAAAATGCAGCAACCGGAAAAAGTACAAAAACCAGCGTTGAAGAAAACACAATCTGCCAGAATATCCGCTGGAAAATAAAACGGTTATCCCCATGTAAAACTTTGTTACGCCGGGGAAAAAACCCCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGct                                                                                                                                                                        <  2:1490393/229‑1 (MQ=255)
                                                                            aTCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAAt                                                                                                                                                  >  2:1725080/1‑229 (MQ=255)
                                                                                tACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATccc                                                                                                                                               >  2:236082/1‑228 (MQ=255)
                                                                                     aTGCAGAAAACTAAAAAGTTACAAAAAACAGCGTTGCAGAAACCACAATCTGCAAGAATATACGGTGGAAAATAAACCGGTTATCCCCATGTAAAACTTTGTTCCGCCGGGGAGTAAACACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAgg                                                                                                                                          <  2:792073/228‑1 (MQ=255)
                                                                                        cagcaAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGa                                                                                                                                       >  1:844343/1‑228 (MQ=255)
                                                                                           caAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATccc                                                                                                                                   >  1:1892493/1‑229 (MQ=255)
                                                                                           caAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATcc                                                                                                                                    >  1:1329421/1‑228 (MQ=255)
                                                                                                        tATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATgccgcc                                                                                                                      >  2:1227740/1‑229 (MQ=255)
                                                                                                            aGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCACAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAAACCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATTCCGCCCAGc                                                                                                                  >  2:1961488/1‑229 (MQ=255)
                                                                                                                              aaCACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAAcatcatc                                                                                                    >  2:1798998/1‑225 (MQ=255)
                                                                                                                                       tgcaataattgaCGCTGGGAAATAAAACGGTAACACCCATAGGAAACGTTGTTACGCGGGGAAGAAAACACCCGTTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGatcat                                                                                         <  2:264371/216‑1 (MQ=255)
                                                                                                                                             aaTATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTaa                                                                                                                                                                                  >  1:1509225/1‑132 (MQ=255)
                                                                                                                                             aaTATACGCGGGAAAATAAAAGGGGTACCCCCTTGTAAACGTTTTTTCCGCGGGGGAGAAAAGCCCCTGAAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTaa                                                                                                                                                                                  <  2:1506612/132‑1 (MQ=255)
                                                                                                                                                ataCGTGGGAAAATTAAACGGTAATCCCCAGGGAAAACTTTGTTACGCGGGGAAGTAAAGCCCGTTTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTg                                                                                              <  1:1019154/213‑1 (MQ=255)
                                                                                                                                                 taCGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAgg                                                                             >  2:1032803/1‑229 (MQ=255)
                                                                                                                                                  aCGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTAACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATACTAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAAcc                                                                                                                                                                             >  1:1120635/1‑132 (MQ=255)
                                                                                                                                                       gggAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATa                                                                       >  1:1219840/1‑229 (MQ=255)
                                                                                                                                                       gggAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATa                                                                       >  2:1807942/1‑229 (MQ=255)
                                                                                                                                                       gggAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATa                                                                       >  2:185664/1‑229 (MQ=255)
                                                                                                                                                        ggAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTgaga                                                                                                                                                                                                                               >  1:821988/1‑76 (MQ=255)
                                                                                                                                                                aaCGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGaaa                                                              >  1:978581/1‑229 (MQ=255)
                                                                                                                                                                    ggTACCACCAAGTGAATCGTGTTGCCGCGGGGGATAAAAGACCGTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGACAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGAAAAAc                                                           <  2:1613385/228‑1 (MQ=255)
                                                                                                                                                                                gtgCAACGCTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCATAGGCAGATAAAAGAGACAAACCTCACCGgaagc                                               >  2:1285272/1‑228 (MQ=255)
                                                                                                                                                                                gtgAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACtga                                                                                                                                                                                                                                 >  1:1006321/1‑51 (MQ=255)
                                                                                                                                                                                   aaaCGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGAAAAACCTCACCGGAAGAAat                                            >  1:203626/1‑229 (MQ=255)
                                                                                                                                                                                    aaCGTTTTTACGCGGGGGAGAAAAGACCCTTTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGct                                                                                                                                                                        <  2:609568/104‑1 (MQ=255)
                                                                                                                                                                                    aaCGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGct                                                                                                                                                                        >  1:610638/1‑104 (MQ=255)
                                                                                                                                                                                       gttgttACGCGGGGGAAAAAAGCCCTTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCa                                                                                                                                                                           <  2:140837/98‑1 (MQ=255)
                                                                                                                                                                                       gttgttACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCa                                                                                                                                                                           >  1:141074/1‑98 (MQ=255)
                                                                                                                                                                                       gttgttACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTCAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGAAAAACCTCACCGGAAGAAATATgc                                        >  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)
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                 
AGTAAGGCCTGATAATTGATTAAAGTCCCGAGGTTAAAAGGCATGACACCGACCAGATTTTCTCTGCTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGAAAAACCTCACCGGAAGAAATATGCGCGAAAGTATTAATGTGGGCGAAAATGGGCAAGATTAA  >  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: