Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,857,245 T→C 8.3% D387G (GAT→GGT)  yidJ ← sulfatase/phosphatase superfamily protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,857,2450TC 8.3% 186.4 / 2.6 68D387G (GAT→GGT) yidJsulfatase/phosphatase superfamily protein
Reads supporting (aligned to +/- strand):  ref base T (44/17);  new base C (3/4);  total (47/21)
Fisher's exact test for biased strand distribution p-value = 1.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGTAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAACGTCTGCAAAACGGATATCATCGATCAGGTTATGCATTTCATTTGGGTCATTACGGCGATCGTAAAGTTCATCACTGGTGAAGAGGTTGAGTACCAGTTTAAAGTCATCCGTCACCCAGCAACGCACCGGAATAAAACCGCCAAAGCTGTCATGCTCAATCTCGTAGCGGTTAAATTCCACCATCACGCCGCGTGGCTCTTTCACGGCAAGGATATTTTCCCCCGGCAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGA  >  NC_000913/3857019‑3857473
                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                    
ccccATCCAGCGGGGCCGTGCATCTTTGCGCCACGGACGCAGACTCCAGTGGTAACGGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTTCGAACGTCTGCAAAACGGATATCATCGATCAGGTTATGCATTTCATTTGGGTCATTACGGCGATCGTAAAGTTCATCACTGGTGAAGAGGTTGAGTACCAGTTTAAAGTCATc                                                                                                                                                                                                                                     <  1:540212/228‑1 (MQ=255)
     tCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGTAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAACGTCTGCAAAACGGATATCATCGCTCAGGTTATGCATTTCATTTGGTTCATTACGCCGATCGTAAAGTTCATCACTGGTGAAGAGGTTGGGTCCCAGTTTAAAGTCACCCTTcc                                                                                                                                                                                                                                >  2:712967/1‑227 (MQ=255)
     tCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGTAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAACGTCTGCAAAACGGATATCATCGATCAGGTTATGCATTTCATTTGGGTCATTACGGCGATCGTAAAGTTCATCACTGGTGAAGAGGTTTGTAACCAGTTTAAAGTCATCCGTccc                                                                                                                                                                                                                               >  1:740362/1‑227 (MQ=255)
          cgcgGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGTAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAACGTCTGCAAAACGGATATCATCGATCCGGTTATGCATTTCATTTGGGTCATTACGGCGATCGTAAAGTTCATACCTGGAGAAGAGGTTGAGTACCAGTTTTAAGTCATCCGTCACCCAGc                                                                                                                                                                                                                          >  1:725839/1‑229 (MQ=255)
              gCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGTAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAACGTCTGCAAAACGGATATCATCGATCAGGTTATGCATTTCATTTGGGTCATTACGGCGATCGTAAAGTTCATCACTGGTGAAGAGGTTGAGTACCAGTTTAAAGTCATCCGTCACCCAGCAACg                                                                                                                                                                                                                      <  1:7666/229‑1 (MQ=255)
               ccGTGCATCTTTACGCCACGGACGCAGACTCCATTGGTAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTACTGCGAACGTCTGCAAAACGGATATCATCGATCAGGTTATGCATTTCATTTGGGTCATTACCGCGATCGTAAAGTTACTCACTGGTGACGAGGTTGAGTACCAGTTTAAAGTCATCCGTCCCCAAGCAAcca                                                                                                                                                                                                                     >  2:510630/1‑227 (MQ=255)
                 cgtcgtCTTTACGCCACTGACGCAGACTCCATTGGTAACTGCGGATCGGATCGAGAATTTTGTCCATGTAATCCAATAAGGGGTCATGCATTTTGCTGCGAACGTCTGCAAAACGGATATCATCGATCAGGTTATGCATTTCATTTGGGTCATTACGGCGATCGTAAAGTTCATCACTGGTGAAGAGGTTGAGTACCAGTTTAAAGTCATCCGTCACCCAGCAACGcgc                                                                                                                                                                                                                   <  1:740624/224‑3 (MQ=255)
                            cGCCACGGACGCAGACTCCATTGGTAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAACGTCTGCAAAACGGATATCATCGATCAGGTTATGCATTTCATTTGGGTCATTACGGCGATCGTAAAGTTCATCACTGGTGAAGAGGTTGAGTACCAGTTTAAAGTCATCCGTCACCCAGCAACGCACCGGAATaa                                                                                                                                                                                                           >  1:113635/1‑226 (MQ=255)
                                 cGGACGCAGACTCCATTGGTAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAACGTCTGCAAAACGGATATCATCGATCAGGTTATGCATTTCATTTGGGTCATTACGGCGATCGTAAAGTTCATCACTGGTGAAGAGGTTTAGTACCAGTTTAAAGTCATCCGTCACCCAGCAACGCACCGGAATAAAACCGcc                                                                                                                                                                                                    >  1:166194/1‑228 (MQ=255)
                                 cGGACGCAGACTCCATTGGTAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAACGTCTGCAAAACGGATATCATCGATCAGGTTATGCATTTCATTTGGGTCATTACGGCGATCGTAAAGTTCATCACTGGTGAAGAGGTTGAGTACCAGTTTAAAGTCATCCGTCACCCAGCAACGCACCGGAATAAAACAGCCa                                                                                                                                                                                                   >  1:372949/1‑229 (MQ=255)
                                                      aCTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAACGTCTGCAAAACGGATATCATCGATCAGGTTATGCATTTCATTTGGGTCATTACGGCGATCGTAAAGTTCATCACTGGTGAAGAGGTTGAGTACCAGTTTAAAGTCATCCGTAACCCAGCAACGCACCGGAATAAAACCGCCAAAGCTGTCATGCTCAATCTCg                                                                                                                                                                              >  1:545812/1‑229 (MQ=255)
                                                        tggggAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAACGTCGGCAAAACGGATATCATCGATCAGGTTATGCATTTCATTTGGGTCATTACGGCGATCGTAAAGTTCATCACTGGTGAAGAGGTTGAGTACCAGTTTAAAGTCATCCGTCACCCAGCAACGCACCGGAATAAAACCGCCAAAGCTGTCATGCTCAATCTCGTa                                                                                                                                                                            <  1:236376/226‑1 (MQ=255)
                                                        tGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAACGTCTGCAAAACGGATATCATCGATCAGGTTATGCATTTCATTTGGGTCATTACGGCGATCGTAAAGTTCATCACTGGTGAAGAGGTTGAGTACCAGTTTAAAGTCATCCGTCACCCAGCAACGCACCGGAATAAAACCGCCAAAGCTGTCATGGCCAATCTCGt                                                                                                                                                                             >  1:686873/1‑228 (MQ=255)
                                                          cGGAACGGATCGCGATTTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAACGTCTGCAAAACGGATATCATCGATCAGGTTATGCATTTCATTTGGGTCATTACGGCGATCGTCAAGTTCATCACTGGTGAAGAGGGTGAGTACCAGTTTAAAGTCATCCGACACCCAGCAACGCACCGGAATAAAACGGCCAAAGCTGTCATGCTcctgtctctttt                                                                                                                                                                          >  1:108024/1‑218 (MQ=255)
                                                          cGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAACGTCTGCAAAACGGATATCATCGATCAGGTTATGCATTTCATTTGGGTCATTACGGCGATCGTCAAGTTCATCACTGGTGAAGAGGTTGAGTACCAGTTTAAAGTCATCCGTCACCCAGCAACGCACCGGAATAAAACAGCCAAAGCTGTCATGCTCAATCTCGTAg                                                                                                                                                                           >  1:583198/1‑228 (MQ=255)
                                                          cGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAACGTCTGCAAAACGGATATCATCGATCAGGTTATGCATTTCATTTGGGTCATTACGGCGATCGTAAAGTTCATCACTGGTGAAGAGGTTGAGTACCAGTTTAAAGTCATCCGTCACCCAGCAACGCACCGGAATAAAACCGCCAAAGCTGTCATGCTCAATCTCGTAg                                                                                                                                                                           >  1:518291/1‑228 (MQ=255)
                                                             aaCGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAACGTCTGCAAAACGGATATCATCGATCAGGTTATGCATTTCATTTGGGTCATTACGGCGATCGTAAAGTTCATCATTGGTGAAGAGGTTGAGTACCAGTTTAAAGTCATCCGTCAACCAGCAACGCACCGGAATAAAACCGCCGAAGCTGTCATGCTCAATCCCGTAGCgg                                                                                                                                                                        >  1:223493/1‑228 (MQ=255)
                                                                  aTCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAACGTCTGCAAAACGGATATCATCGATCAGGTTATGCATTTCATTTGGGTCATTACGGCGATCGGAAAGTTCATCACTGGTGAAGAGGTTGAGACCCAGTTTAAAGTAATCCGTCACCCAGCAACGCACCGAACTAAAAACGCCAAAGCTGCCACGCTCAACCTCCAAACGGTTaaa                                                                                                                                                                   >  2:270357/1‑228 (MQ=255)
                                                                  aTCGCGAATTTTGTCCATGTAATCCAATAAGGCGGCATGCATTTTGCTGCGAACGTCTGCAAAACGGATATCATCGATCAGGTTATGCATTTCATTTGGGTCATTACGGCGATCGTAAAGTTCATCACTGGTGAAGAGGTTGAGTACCAGTTTAAAGTCATCCGTCACCCAGCAACGCACCGGAATAAAACCGCCAAACCTGTAAACATCAATCTCGTAGCGGTTaaa                                                                                                                                                                   >  1:635723/1‑228 (MQ=255)
                                                                                        tCCAATAAGGCGTCCTGCATTTTGCTGCGAACGTCTGCAAAACGGATATCATCGATCAGGTTATGCATTTCATTTGGGTCATTACGGCGATCGTAAAGTTCATCACTGGTGAAGCGGTTGAGGACCAGTTTAAAACCATCAGTCACCCAGAAATGCACCGGGATAAAACCGCCAAAGTTGCCAGGCTCAATCTCGCAGAGGTTAAATTACCCCACCAccccgcagggcc                                                                                                                                            >  2:554686/1‑223 (MQ=255)
                                                                                             tAAGGCGTCATGCATTTTGCTGCGAACGTCTGCAAAACGGATATCATCGATCAGGTTATGCATTTCATTTGGGTCATTACGGCGATCGTAAAGTTCATCACTGGTGAAGAGGTTGAGTACCAGTTTAAAGTCATCCGTCACCCAGCAACGCACCGGAATAAAACCGCCAAAGCTGTCATGCTCAATCTCGTAGCGGTTAAATTCCACCATCACGCCGCGTGGCTCGTTc                                                                                                                                       >  1:272644/1‑229 (MQ=255)
                                                                                             tAAGGCGTCATGCATTTTGCTGCGAACGTCTGCAAAACGGATATCATCGATCAGGTTATGCATTTCATTTGGGTCATTACGGCGATCGTAAAGTTCATCACTGGTGAAGAGGTTGAGTACCAGTTTAAAGTCATCCGTCAACCAGCAACGCACCGGAATAAAACCGCCAAAGCTGTCATGCCCAATCTCGTACCGGTTAAATTCCACCATCACGCCCCGTGGCTATTTc                                                                                                                                       >  1:205474/1‑229 (MQ=255)
                                                                                                  gcgaaTGCTTTTTGCTGCGAACGTCTGAAAAACGGATATCATCGATCAGGTTATGCATTTCATTTGGGTAATTACGGCGATCGTAAAGTTCATCACTGGTGAAGAGGTTGAGTACCAGTTTAAAGTCCCCCGTCACCCAGCAACGCACCGGAATAAAACCGCCAAAGCTGTCATGCTCAATCTCGTAGCGGTTAAATTCCACCATCACGCCGCGTGGCTCTTTCACGGc                                                                                                                                  <  1:538403/225‑1 (MQ=255)
                                                                                                  cGTCATGCATTTTGCTGCGAACGTCTGCAAAACGGATATCATCGATCAGGTTATGCATTTCATTTGGGTCATTACGGCGATCGTAAAGTTCATCACTGGTGAAGAGGTTGAGTACCAGTTTTAAGTCACCCGTCACCCAGCAACGCACCGGAATAAAACCGCCAAAGATGGCATGCTCAATCTCCGAGCGGTTAAATTCCACCATCACGCCGCGTGGCCCTTTCACGGc                                                                                                                                  >  2:538403/1‑229 (MQ=255)
                                                                                                  cGTCATGCATTTTGCTGCGAACGTCTGCAAAACGGATATCATCGATCAGGTTATGCATTTCATTTGGGTCATTACGGCGATCGTAAAGTTCATCACTGGTGAAGAGGTTGAGTACCAGTTTAAAGTCATCCGTCACCCAGCAACGCACCGGAATAAAACCGCCAAAGCTGTCATGCTCAATCTCGTAGCGGTTAAATTCCCCCATCACGCCACGAGGCCCTTTCACGGc                                                                                                                                  >  1:725925/1‑229 (MQ=255)
                                                                                                  cGTCATGCATTTTGCTGCGAACGTCTGCAAAACGGATATCATCGATCAGGTTATGCATTTCATTTGGGTCATTACGGCGATCGTAAAGTTCATCACTGGTGAAGAGGTTGAGTACCAGTTTAAAGTCATCCGTCACCCAGCAACGCACCGGAATAAAACCGCCAAAGCTGTCATGCTCAATCTCGTAGCGGTGAAATTCCACCATCACGCCCCGTGGCTCTTTCACgca                                                                                                                                  >  1:56585/1‑227 (MQ=255)
                                                                                                  cGTCATGCATTTTGCTGCGAACGTCTGCAAAACGGATATCATCGATCAGGTTATGCATTTCATTTGGGTCATTACGGCGATCGTAAAGTTCATCACTGGTGAAGAGGTTGAGTACCAGTTTAAAGTCATCCGTCACCCAGCAACGCACCGGAATAAAACCGCCAAAGCTGTCATCCTCAATCTCGTAGCGGTTAAATTCCACCATAACCCCCCGTGGCCCTTTccaggc                                                                                                                                  >  1:644890/1‑224 (MQ=255)
                                                                                                  cGTCATGCATTTTGCTGCGAACGTCTGCAAAACGGATATCATCGATCAGGTTATGCATTTCATTTGGGTCATTACGGCGATCGTAAAGTTCATCACTGGTGAAGAGGTTGAGTACCAGTTTAAAGTCACCCGTCACCCAGCAACGCAACGGAAAAAAACCGCCAAAGCTGTCATGCTCAACTTCGAGCCGGTTAAATTCACACAACACCCCGCTGGGCCCTTTCCCgga                                                                                                                                  >  2:538389/1‑228 (MQ=255)
                                                                                                  cGTCATGCATTTTGCTGCGAACGTCTGAAAAACGGATATCATCGATCAGGTTATGCATTTCATTTGGGTCATTACGGCGATCGTAAAGTTCATCACTGGTGAAGAGGTTGAGTACCAGTTTAAAGTCACCCGTCACCCAGCAACGCACCGGAATAAAACCGCCAAAGCTGTCATGCTCAATCTCGTAGCGGTTAAATTCCACCATCACGCCGCGTGGCTCTTTCACGGc                                                                                                                                  <  1:538389/229‑1 (MQ=255)
                                                                                                   ctCATGCATTTTGCTGCGAACGTCTGCAAAACGGATATCATCGATCAGGTTATGCATTTCATTTGGGTCATTACGGCGATCGTAAAGTTCATCACTGGTGAAGAGGTTGAGTACCAGTTTAAAGTCATCCGTCCCCCAGCAACGCACCGGAATAAAACCGCCAAAGCTGTCATGCTCAATCTCGTAGCGGTTAAATTCCACCATCACGCCGCGTGGCTCTTTCACGGCa                                                                                                                                 <  1:535252/228‑1 (MQ=255)
                                                                                                    ttttgatttttGCTGCAAACGGCTGAAAAAGGGAATTCATCGATTGGTTTATGCCTTTAATTTGGGTTCTTACGGCGATCGTAAAGTTCATCACGGGTGAAGAGGTTGAGTACCAGTTTAAAGTCATCCGTCCCCCAGCAACGCACCGGAATAAAACCGCCAAAGCTGTCATGCTCAATCTCGTAGCGGTTAAATTCCACCATCACGCCGCGTGGCTCTTTCACGGCaa                                                                                                                                <  2:725838/222‑1 (MQ=255)
                                                                                                    tCATGCATTTTGCTGCGAACGTCTGTAAAACGGATATCATCGATCAGGTTATGCATTTCATTTGGGTCATTACGGCGATCGTAAAGTTCATCACTGGTGAAGAGGTTGAGTACCAGTTTAAAGTCATCCGTCACCCAGCAACGCACCGGAATAAAACCGCCAAAGCTGTCATGCTCAATCTCGTAGCGGTTAAATTCCACCATCACGCCCCGTGGCCCTTTCACGGCa                                                                                                                                 >  1:500070/1‑228 (MQ=255)
                                                                                                    tCATGCATTTTGCTGCGAACGTCTGCAAAACGGATATCATCGATCAGGTTATGCATTTCATTTGGGTCATGACGGCGATCGTAAAGTTCATCACTGGTGAAGAGGTTGAGTACCAGTTTAAAACCATCCGTCCCCCAGCAACGCACCGGAATAAAACCGCCAAAACTGCTATGCTCAAACTCGTAGCGGTTAAATTCCACCAACAACCCCCAGGGCCGTTTAACGGCaa                                                                                                                                >  2:328216/1‑229 (MQ=255)
                                                                                                            tttGCTGCGAACGTCTGCAAAACGGATATCATCGATCAGGTTATGCATTTCATTTGGGTCATTACGGCGATCGTAAAGTTCATCACTGGTGAAGAGGTTGAGTACCAGTTTAAAGTCATCCGTCACCCAGCAACGCACCGGAATAAAACCGCCAAAGCTGTCATGCTCAATCTCGTAGCGGTTAAATTCCACCATCACGCCGCGTGGCTCTTTCACGGCAAGGATAttt                                                                                                                        >  1:483651/1‑229 (MQ=255)
                                                                                                                           tGCAAAACGGATATCATCGATCAGGTTATGCATTTCATTTGGGTCATTACGGCGATCGTACAGTTCATCACTGGTGAAGCGGTTGAGTCCCAGTTTAACGTCATCCGTCACCCAGCAACGCCCCGGAATAAAACCGCCAAAGCTGTCCTGCTCAATCTCTCAGCGTGTCAATTCCACCATCACTCCGCCCGACCCTTTCACGGCAAGGATATTTTCCCCCGGccggat                                                                                                          >  1:603924/1‑223 (MQ=255)
                                                                                                                           tGCAAAACGGATATCATCGATCAGGTTATGCATTTCATTTGGGTCATTACGGCGATCGTAAAGTTCATCACTGGTGAAGAGGTTGAGTACCAGTTTAAAGTCATCCGTCACCCAGCAACGCACCGGAATAAAACCGCCAAAGCTGTCATGCTCAATCTCGTAGCGGTTAAATTCCACCATCACGCCGCGTGGCTCTTTCACGGCAAGGAAATTTTTCCCCGGCACAAtc                                                                                                         >  1:672150/1‑229 (MQ=255)
                                                                                                                             aaaaaCGGATTTCATCGATCAGGTTATGCATTTCATTTGGGTCATTACGGCGATCGTAAAGTTCATCACTGGTGAAGAGGTTGAGTACCAGTTTAAAGTCATCCGTCACCCAGCAACGCACCGGAATAAAACCGCCAAAGCTGTCATGCTCAATCTCGTAGCGGTTAAATTCCACCATCACGCCGCGTGGCTCTTTCACGGCAAGGATATTTTCCCCCGGCAGAAtctc                                                                                                       <  1:14867/228‑1 (MQ=255)
                                                                                                                                     atatCATCGATCAGGTTATGCATTTCATTTGGGTCATTACGGCGATCGTAAAGTTCATCACTGGTGAAGAGGTTGAGTACCAGTTTAAAGTCATCCGTCACCCAGCAACGCACCGGAATAAAACCGCCAAAGCTGTCATGCTCAATCTCGTAGCGGTTAAATTCCACCATCACGCCGCGTGGCTCTTTCACGGCAAGGATATTTTCCCCCGGCAGAATCTCTGGttttt                                                                                               >  1:127181/1‑229 (MQ=255)
                                                                                                                                            cGATCAGGTTATGCATTTCATTTGGGTCATTACGGCGATCGTAACGTTCATCACTGGTGAAGAGGTTGAGTACCAGTTTAAAGTCATCCGTCACCCAGCAACGCACCGGAATAAAACCGCCAAAGCTGTCATGCTCAATCTCGTAGCGGTTAAATTCCACCATCACGCCGCGTGGCTCTTTCACGGCaa                                                                                                                                >  2:40974/1‑189 (MQ=255)
                                                                                                                                            cGATCAGGTTATGCATTTCATTTGGGTCATTACGGCGATCGTAAAGTTCATCACTGGTGAAGAGGTTGAGTACCAGTTTAAAGTCATCCGTCACCCAGCAACGCACCGGAATAAAACCGCCAAAGCTGTCATGCTCAATCTCGTAGCGGTTAAATTCCACCATCACGCCGCGTGGCTCTTTCACGGCaa                                                                                                                                <  1:40974/189‑1 (MQ=255)
                                                                                                                                                      aTGCATTTCATTTGGGTCATTACGGCGATCGTAAAGTTCATCACTGGTGATGAGGTTGAGTACCAGTTTAAAGTCATCCGTCACCCAGCAACGCACCGGAATAAAACCGCCAAAGCTGTCATGCTCAATCTCGTAGCGGTTAAATTCCACCCTCACGCCGCCTGGCTCTTTCACGGTAAGGAGATTTTCCCCCAGCAGACTCTCTGGCTTTTCAATATCTGCCAgccca                                                                              >  2:219370/1‑226 (MQ=255)
                                                                                                                                                      aTGCATTTCATTTGGGTCATTACGGCGATCGTAAAGTTCATCACTGGTGAAGAGGTTGAGTACCAGTTTAAAGTCATCCGTCACCCAGCAACGCACCGGAATAAAACCGCCAAAGCTGTCATGCTCAATCTCGTAGCGGTTAAATTCCCCCATCACGCCGCGTGGCTCTTTCACGGc                                                                                                                                  >  2:66804/1‑177 (MQ=255)
                                                                                                                                                      aTGCATTTCATTTGGGTCATTACGGCGATCGTAAAGTTCATCACTGGTGAAGAGGTTGAGTACCAGTTTAAAGTCATCCGTCACCCAGCAACGCACCGGAATAAAACCGCCAAAGCTGTCATGCTCAATCTCGTAGCGGTTAAATTCCACCATCACGCCGCGTGGCTCTTTCACGGc                                                                                                                                  <  1:66804/177‑1 (MQ=255)
                                                                                                                                                         catttcatttGGGTCATTACGGCGATCGTAAAGTTCATCACTGGTGAAGAGGTTGAGTACCAGTTTAAAGTCATCCGTCACCCAGCAACGCACCGGAATAAAACCGCCAAAGCTGTCATGCTCAATCTCGTAGCGGTTAAATTCCACCATCACGCCGCGTGGCTCTTTCCCGGCAAGGGTATTTTCCCCCGGCAGAATCTCTGGTTTTTCAAAATCTCCCAGCGCcctc                                                                           >  1:524793/1‑229 (MQ=255)
                                                                                                                                                             tcatttGGGTCATTACGGCGATCGTAAAGTTCATCACTGGTGAAGAGGTTGAGTACCAGTTTAAAGTCATCCGTCACCCAGCAACGCACCGGAATAAAACCGCCAAAGCTGTCATGCTCAATCTCGTAGCGGTTAAATTCCACCATCACGCCGCGTGGCTCTTTCACGGCAAGGATATTTTCCCCCGGCAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCAttt                                                                       >  1:162072/1‑228 (MQ=255)
                                                                                                                                                             tcatttGGGTCATTACGGCGATCGTAAAGTTCATCACTGGTGAAGAGGTTGAGTACCAGTTTAAAGTCATCCGTCACCCAGCAACGCACCGGAATAAAACCGCCAAAGCTGTCATGCTCAATCTCGTAGCGGTTAAATTCCACCATCACGCCGCGGGGCTCTTTCACGGCAAGGATATTTTCCCCCGGCAAAAATCCTGGTTTTTCAATATCTGCCATAGCCATCATTg                                                                       >  2:549333/1‑229 (MQ=255)
                                                                                                                                                                    ggTCATTACGGCAAGCGTAAAGTTCCTCACTGGTGAAGAGGTTGAATACCAGTTTAAAGTCATCCGTCCCCCAGCAACGCACCGGAATAAAACCGCCAAAGCTGTCATGCTCAATCTCGTAGCGGTTAAATTCCACCATCACGCCGCGTGGCTCTTTCACGGCAAGGATATTTTCCCCCGGCAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCa                                                                 <  1:510630/228‑1 (MQ=255)
                                                                                                                                                                       cATTTCGGCAATCGTAAAGTTCCTCCCTGGAGAAGAGGTTGAGTACCAGTTTAAAGTCACCCGTCCCCCAGCAACGCCCCGGAATAAAACCGCCAAAGCTGTCATGCTCAATCTCGTAGCGGTTAAATTCCACCATCACGCCGCGTGGCTCTTTCACGGCAAGGATATTTTCCCCCGGCAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTaa                                                             <  2:272644/229‑1 (MQ=255)
                                                                                                                                                                               cGATCGTAAAGTTAAGCCCTGGTGAATAGGTTGAGTACCAGTTTAAAGTCATCCGTCACCCAGCAACCGACCGGAATAAAACCGCCAAAGCTGGCATGCTCAATTTAGTAGCGGTTAAATTCCACCATCACGCCGCGTGGCTCTTTCACGGCAAGGATATTTTCCCCCGGCAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGata                                                      <  1:270357/228‑1 (MQ=255)
                                                                                                                                                                                   cGTAAAGTTCATCACTGGTGAAGAGGTTGAGTACCAGTTTAAAGTCATCCGTCACCCAGCAACGCACCGGAATAAAACCGCCAAAGCTGTCATGCTCAATCTCGTAGCGGTTAAATTCCACCATCACGCCGCGTGGCTCTTTCACGGCCAGGATATTTTCTCCCGGCAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGTTAtgact                                                 >  1:215579/1‑229 (MQ=255)
                                                                                                                                                                                   cGTAAAGTTCATCACTGGTGAAGAGGTTGAGTACCAGTTTAAAGTCATCCGTCACCCAGCAACGCACCGGAATAAAACCGCCAAAGCTGTCATGCTCAATCTCGTAGCGGTTAAATTCCACCATCACGCCGCGTGGCTCTTTCACGGCCAGGATATTTTCCCCCGGCAGAGTCTCTGGTTTTTCAATCTCTGCCAGCGCCATCATTGTCGGCAGTAAATCGCTAtgac                                                  >  1:112097/1‑228 (MQ=255)
                                                                                                                                                                                        attcctccctgggTGAAGAGGTTGAGAACAAGTTTAAAGTCATCCGTCTCCCAGCACCGCCCGGGAATAAAACCGCCAAAGCTGTCCTGCTCAATCTCGTAGCGGTTAAATTCCCCCATCACGCCGCGTGGCTCTTTCACGGCAAGGATATTTTCCCCCGGCAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGactgactg                                            <  2:518291/218‑1 (MQ=255)
                                                                                                                                                                                             atcaCTGGTGAAGAGGTTGAGTACCAGTTTAAAGTCATCCGTCACCCAGCAACGCACCGGAATAAAACCGCCAAAGCTGTCATGCTCAATCTCGTAGCGGTTAAATTCCACCATCACGCCGCGTGGCTCTTTCACGGCAAGGATATTTTCCCCCGGCCGGATCTCTGGTTTTTCAATATCGTCAAGCGCACCCATTTTCGCGAGTAAAACGATAttactacacggccta                                       >  2:597442/1‑219 (MQ=255)
                                                                                                                                                                                              cctcgtgggAAGGGGTGGAATCCCAGTTGAAAGCCATCCTTCACACAGCAACGCCCAGAAATAAAACCGCCAAAGCTGTCATGCTCAATCTCGTAGCGGTTAAATTCCACCATCACGCCGCGTGGCTCTTTCACGGCAAGGATATTTTCCCCCGGCAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTAt                                      <  2:740361/221‑1 (MQ=255)
                                                                                                                                                                                                  tGGTGAAGAGGTTGAGTACCAGTTTAAAGTCATCCGTCACCCAGCAACGCACCGGAATAAAACCGCCAAAGCTGTCATGCTCAATCTCGTAGCGGTTAAATTCCACCATCACGCCGCGTGGCTCTTTCACGGCAAGGATATTTTCCCCCGGCAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTCTCGGCAGTAAATCGATATGACTGACTGGCGTATCGa                                   >  1:287871/1‑228 (MQ=255)
                                                                                                                                                                                                    gTGAAGAGGTTGAGTACCAGTTTAAAGTCATCCGTCACCCAGCAACGCACCGGAATAAAACCGCCAAAGCTGTCATGCTCAATCTCGTAGCGGTTAAATTCCACCATCACGCCGCGTGGCTCTTTCACGGCAAGGATATTTTCCCCCGGCAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGAcc                                 >  1:259088/1‑228 (MQ=255)
                                                                                                                                                                                                       aaGAGGTTGAGTACCAGTTTAAAGTCATCCGTCACCCAGCAACGCACCGGAATAAAACTGCCAAAGCTGTCATGCTCAATCTCGTAGCGGTTAAATTCCACCATACCGCCGCGAGGCTCTTTCAAGGCAAGGATATTTTACCCCGGCCATATCTCTTGCTTTTCAATATATGCCAAGCCCATCTATGTCGGCAGTAAAACTAAATGACTAACTGGCGTATcacactgtc                             >  2:605753/1‑221 (MQ=255)
                                                                                                                                                                                                       aaGAGGTTGAGTACCAGTTTAAAGTCATCCGTCACCCAGCAACGCACCGGAATAAAACCGCCAAAGCTGTCATGCTCAATCTCGTAGCGGTTAAATTCCACCATCACGCCGCGTGGCTCTTTCACGGCAAGGATATTTTCCCCCGGCAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTTTCGGCAGTAAATCGATATGACTTACTGGCGTATCGACCTg                               >  1:587411/1‑227 (MQ=255)
                                                                                                                                                                                                         ggggTTGAGATCCGGTTTAAAGGTCTCCGTCCACCAGCAACGCACCGGAATAAACCCGCCAAAGCTGTCATGCTCAATCTCGTAGCGGTTAAATTCCGCCATCACGCCGCGTGGCTCTTTCACGGCAAGGATATTTTCCCCCGGCAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTcgc                           <  1:554687/227‑1 (MQ=255)
                                                                                                                                                                                                         gagGTTGAGTACCAGTTTAAAGTCATCCGTCACCCAGCAACGCACCGGAATAAAACCGCCAAAGCTGTCATGCTCAATCTCGTAGCGGTTAAATTCCACCATCACGCCGCGTGGCTCTTTCACGGCAAGGATATTTTCCCCCGGCAGAATCTCTGGTTTTTCAATATCATCCAGAGCCACCATTTACGGCAGTAAATCGATATGATGGACTGTGGTACCGACCGTTcgc                           >  1:434789/1‑229 (MQ=255)
                                                                                                                                                                                                              tGAGTACCAGTTTAAAGTCATCCGTCACCCAGCAACGCACCGGAATAAAACCGCCAAAGCTGTCATGCTCAATCTCGTAGCGGTTAAATTCCACCATCACGCCGCGTGGCTCTTTCACGGCAAGGATATTTTCCCCCGGCAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGCCAGTAAATCGATATGACTGACTGGCGGCTCGACCTGTCGCCtctc                      >  1:422170/1‑229 (MQ=255)
                                                                                                                                                                                                               gAGTACCAGTTTAAAGTCATCCGTCACCCAGCAACGCACCGGAATAAAACCGCCAAAGCTGTCATGCTCAATCTCGTAGCGGTTAAATTCCACCATCACGCCGCGTGGCTCTTTCACGGCAAGGATATTTTCCCCCGGCGGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACCGGCGTATCGCCCTTTCGCCCCTcc                     >  1:602050/1‑229 (MQ=255)
                                                                                                                                                                                                                      aGTTTAAAGTCATCCGTCACCCAGCAACGCACCGGAATAAAACCGCCAAAGCTGTCATGCTCAATCTCGTAGCGGTTAAATTCCACCATCACGCCGCGTGGCTCTTTCACGGCAAGGATATTTTCCCCCGGCAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTCTCGGCAGTAAATCGATATGACTGTCTTGGGAATCGACCTATCGCCGCTCCACTTgcg              >  1:753836/1‑229 (MQ=255)
                                                                                                                                                                                                                       gTTTAAAGTCATCCGTCACCCAGCAACGCACCGGAATAAAACCGCCAAAGCTGTCATGCTCAATCTCGTAGCGGTTAAATTCCACCATCACGCCGCGTGGCTCTTTCACGGCAAGGATATTTTCCCCCGGCAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGGGTATCGACCTGTCGCCGCTCCCCttcggg             >  1:474051/1‑225 (MQ=255)
                                                                                                                                                                                                                         ttAAAGTCACCCGTCCCCCAGCACAGCCCCGGAATAAAACCGCCAAAGCTGTCATTGTCAATCTCGTAGCGGTTAAATTCCACCATCACGCCGCGTGGCTCTTTCACGGCAAGGATATTTTCCCCCGGCAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTgcggcg           <  2:644889/229‑1 (MQ=255)
                                                                                                                                                                                                                            aaGTCCTCCGTCACCCAGCACCGCCCCGGAATAAACCCGCCAAAGCTGTCATGCTCAATCTCGTAGCGGTTAAATTCCACCATCACGCCGCGTGGCTCTTTCACGGCAAGGATATTTTCCCCCGGCAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAAc        <  1:231382/229‑1 (MQ=255)
                                                                                                                                                                                                                            aaGTCATACGTCCCCGAGCAACGAACCAGAATAAAACCGCCAAAACTGTCATGCTCAATCTCGTAGCGGTTAAATTCCACCATCACGCCGCGTGGCTCTTTCACGGCAAGGATATTTTCCCCCGGCAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGaa         <  1:712968/228‑1 (MQ=255)
                                                                                                                                                                                                                              gccATCCTTCCCCCACCAGGGCGCGGGAATAAAACCGCGAAAGCGGTCATGCTCAACCCCGTAGCGGTTAAATTCCACCAGCACGCCGCGTGGCTCTTTCACGGCAAGGATATTTTCCCCCGGCAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACgg      <  2:125440/227‑1 (MQ=255)
                                                                                                                                                                                                                                  tCCGTCACCCAGCAACGCACCGGAATAAAACCGCCAAAGCTGTCATGCTCAATCTCGTAGCGGTTAAATTCCACCATCACGCCGCGTGGCTCTTTCACGGCAAGGATATTTTCCCCCGGCAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCCTTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGgatgc  >  1:142286/1‑228 (MQ=255)
                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                    
CCCCATCCAGCGCGGCCGTGCATCTTTACGCCACGGACGCAGACTCCATTGGTAACTGCGGAACGGATCGCGAATTTTGTCCATGTAATCCAATAAGGCGTCATGCATTTTGCTGCGAACGTCTGCAAAACGGATATCATCGATCAGGTTATGCATTTCATTTGGGTCATTACGGCGATCGTAAAGTTCATCACTGGTGAAGAGGTTGAGTACCAGTTTAAAGTCATCCGTCACCCAGCAACGCACCGGAATAAAACCGCCAAAGCTGTCATGCTCAATCTCGTAGCGGTTAAATTCCACCATCACGCCGCGTGGCTCTTTCACGGCAAGGATATTTTCCCCCGGCAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGA  >  NC_000913/3857019‑3857473

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: ‑