Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,329,209 A→G W197R (TGG→CGG)  yfaQ ← tandem DUF2300 domain‑containing protein YfaQ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,329,2090AG98.2% 189.4 / NA 54W197R (TGG→CGG) yfaQtandem DUF2300 domain‑containing protein YfaQ
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (28/25);  minor base T (0/1);  total (28/26)
Fisher's exact test for biased strand distribution p-value = 4.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACTGCCGTCGTACTCTTTTCCGCTGTGATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCACATCCAGTAAGACTTCCTGAGCTTTGTTTTGCGCGGGTAACGTCGCTAATGAATCAAGCAGCGATTTAACCGTTACCGAGGTTTCTGGCTTCAGCGTAGTCAGAGACGCTAACCACGCAGGTGCCTGACGCTGTTGCCAGTATTGCCCCCACACATCTGCGCCAATAT  >  NC_000913/2328962‑2329442
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                         
aCTGCCGTCGTACTCTTTTCCGCTGTGATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGt                                                                                                                                                                                                                                        <  1:258677/251‑1 (MQ=255)
        cgtACTCTTTTCCGCTGTGATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCACAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCCCTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCAGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCAccc                                                                                                                                                                                                                                >  2:305738/1‑251 (MQ=255)
           aCTCTTTTCCGCTGTGATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGt                                                                                                                                                                                                                              >  1:441619/1‑250 (MQ=255)
              ctTTTCCGCTGTGATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTAc                                                                                                                                                                                                                            <  2:150984/249‑1 (MQ=255)
                  tCCGCTGTGATTTTCTTCAGCGGATAGCGGGCACACAGCTCCACCTCTACACATTGCCCACTGGCAAACTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACt                                                                                                                                                                                                                      >  1:424704/1‑251 (MQ=255)
                   ccGCTGTGATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACt                                                                                                                                                                                                                      >  2:226070/1‑250 (MQ=255)
                   ccGCTGTGATTTGCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCTCCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACTa                                                                                                                                                                                                                     <  2:295840/251‑1 (MQ=255)
                          gATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGGTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACTCCCGAGc                                                                                                                                                                                                               >  2:399346/1‑250 (MQ=255)
                           aTTTTCTTCAGCGGATAGCGGGCAAACAGCTCTACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCTCCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACTACCGAGCAt                                                                                                                                                                                                             <  2:463309/251‑1 (MQ=255)
                                  tcAGCGGATAGCGGGCCAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCg                                                                                                                                                                                                      >  1:608813/1‑251 (MQ=255)
                                           aGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGAt                                                                                                                                                                                              <  1:226070/250‑1 (MQ=255)
                                            gCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGAttt                                                                                                                                                                                            <  1:568052/251‑1 (MQ=255)
                                              gggCAAACAGCTCCACCTCTACACATGGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTc                                                                                                                                                                                          <  2:441619/251‑1 (MQ=255)
                                               ggCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCg                                                                                                                                                                                         <  2:189320/251‑1 (MQ=255)
                                                   aaCAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCCCCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTc                                                                                                                                                                                     <  2:134058/251‑1 (MQ=255)
                                                                     caTTGCCCACTGGCAACCTGCGTCGGCACTGGCAAATCCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAACCCACCTCCGGGTTTTCACCCGTAccc                                                                                                                                                                                                                          <  1:584579/196‑1 (MQ=255)
                                                                     caTTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGAACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAACACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTAccc                                                                                                                                                                                                                          >  2:584579/1‑196 (MQ=255)
                                                                      aTTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGCCACCGATTTTCGCCTCGTCCAGCACCACATCCAGt                                                                                                                                                                  >  2:518864/1‑251 (MQ=255)
                                                                                ggCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCACATCCAGTAAGACTTCCt                                                                                                                                                        >  2:188954/1‑251 (MQ=255)
                                                                                ggCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCACATCCAGTAAGACTTCCt                                                                                                                                                        <  1:619548/251‑1 (MQ=255)
                                                                                  cAACCTGCGTCGACACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATAGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGCCACCGATTTTCGCCTCGTCCAGCAACACCTCCCGTAGAACTTCCTg                                                                                                                                                       >  2:73672/1‑250 (MQ=255)
                                                                                      cTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCACATCCAGAAAGACTTCCTGAGCtt                                                                                                                                                  >  1:410150/1‑251 (MQ=255)
                                                                                          gTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCACATCCAGTAAGACTTCCTGAGCTTTGtt                                                                                                                                              <  1:479158/251‑1 (MQ=255)
                                                                                                   ggCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCACATCCAGTAAGACTTCCTGAGCTTTGTTTTGCGCgg                                                                                                                                      <  1:581453/250‑1 (MQ=255)
                                                                                                        aaCCCGATTAAGAACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCAACTGCCAGTAACCCACAGCGGTGTTCCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCACATCCAGTAAGACTTCCTGAGCTTTGTTTTGCGCGGGTAACg                                                                                                                                <  2:294601/251‑1 (MQ=255)
                                                                                                          cccGATTAAGCACCGTCGCATAACGGGGTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCCCCCGGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCACATCCAGTAAGACTTCCTGAGCTTTGTTTTGCGCGGGTAACGTc                                                                                                                              <  2:608813/251‑1 (MQ=255)
                                                                                                           ccGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCACATCCAGTAAGACTTCCTGAGCTTTGTTTTGCGCCGGTAACGTCg                                                                                                                             <  1:73672/251‑1 (MQ=255)
                                                                                                                 aaGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCACATCCAGTAAGACTTCCTGa                                                                                                                                                      <  2:352251/220‑1 (MQ=255)
                                                                                                                 aaGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCACATCCAGTAAGACTTCCTGa                                                                                                                                                      >  1:352251/1‑220 (MQ=255)
                                                                                                                           gCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCACATCCAGTAAGACTTCCTGAGCTTTGTTTTGCGCGGGTAACGTCGCTAATGAATCAagc                                                                                                               >  2:292091/1‑249 (MQ=255)
                                                                                                                                  gggatAGTACCGTTTTGCTGGCGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCACATCCAGTAAGACTTCCTGGGCTTTGTTTTGCGCGGGTAACGTCGCTAATGAATCAAGCAGCGAttt                                                                                                       >  2:170371/5‑250 (MQ=255)
                                                                                                                                        gTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCACATCCAGTAAGACTTCCTGAGCTTTGTTTTGCGCGGGTAACGTCGCTAATGAATCAAGCAGCGATTTAACCGtt                                                                                                >  2:202187/1‑251 (MQ=255)
                                                                                                                                              ttttGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCACATCCAGTAAGACTTCCTGAGCTTTGTTTTGCGCGGGTAACGTCGCTAATGAATCAAGCAGCGATTTAACCGTTACCg                                                                                            >  2:163079/1‑249 (MQ=255)
                                                                                                                                                    tGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCACATCCAGTACGACTTCCTGAGCTTTGTTTTGCGCGGGTAACGTCGCTAATGAATCAAGCAGCGATTTAACCGTTACCGAGGTTTCt                                                                                    >  1:549496/1‑251 (MQ=255)
                                                                                                                                                      gTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCACATCCAGTAAGACTTCCTGAGCTTTGTTTTGCGCGGGTAACGTCGCTAATGAATCAAGCAGCGATTTAACCGTTACCGAGGTTTCTgg                                                                                  <  2:224643/251‑1 (MQ=255)
                                                                                                                                                         ccACTGACAGTAACCCACAGCGGTGTACCATCGGTTACCCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCCTCGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCCGCACCACATCCAGTAAGACTTCCTGAGCTTTGTTTTGCGCGGGTAACGTCGCTAATGAATCAAGCAGCGATTTAACCGTTACCGAGGTTTCTGGCtt                                                                               >  2:318744/1‑251 (MQ=255)
                                                                                                                                                          cACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCACATCCAGTAAGACTTCCTGAGCTTTGTTTTGCGCGGGTAACGTCGCTAATGAATCAAGCAGCGATTTAACCGTTACCGAGGTTTCTGGCTTc                                                                              >  1:414327/1‑251 (MQ=255)
                                                                                                                                                             tGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCACATCCAGTAAGACTTCCTGAGCTTTGTTTTGCGCGGGTAACGTCGCTAATGAATCAAGCAGCGATTTAACCGTTACCGAGGTTTCTGGCTTCAGc                                                                           >  1:731809/1‑251 (MQ=255)
                                                                                                                                                                 aGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCACATCCAGTAAGACTTCCTGAGCTTTGTTTTGCGCGGGTAACGTCGCTAATGAATCAAGCAGCGATTTAACCGTTACCGAGGTTTCTGGCTTCAGCGTAg                                                                       <  2:455460/251‑1 (MQ=255)
                                                                                                                                                                    aaCCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCACATCCAGTAAGACTTCCTGAGCTTTGTTTTGCGCGGGTAACGTCGCTAATGAATCAAGCAGCGATTTAACCGTTACCGAGGTTTCTGGCTTCAGCGTAGTCa                                                                    <  1:126760/251‑1 (MQ=255)
                                                                                                                                                                         acaGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCACACCAGCTCCGGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCACATCCAGTAAGACTTCCTGAGCTTTGTTTTGCGCGGGTAACGTCGCTAATGAATCAAGCAGCGATTTAACCGTTACCGAGGTTTCTGGCTTCAGCGTAGTCAGAGAc                                                               >  2:363571/1‑251 (MQ=255)
                                                                                                                                                                           aGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCACATCCAGTAAGACTTCCTGAGCTTTGTTTTGCGCGGGTAACGTCGCTAATGAATCAAGCAGCGATTTAACCGTTACCGAGGTTTCTGGCTTCAGCGCAGTCAGAGACGc                                                             <  1:664194/251‑1 (MQ=255)
                                                                                                                                                                               gtgtACCATCGGTTAACCAGCCAGCAAAACCTCCCAGACGGATTTCCTGTTTATCATCGGCAAACCAGCGCCGGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCACATCCAGTAAGACTTCCTGAGCGTTGTTTTGCGCGGGTAACGTCGCTAATGAATCAAGCAGCGATTTAACCGTTACCGAGGTTTCTGGCTTCAGCGTAGTCAGAGACGCTAAc                                                         <  1:305738/251‑1 (MQ=255)
                                                                                                                                                                                 gtACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCACATCCAGTAAGACTTCCTGAGCTTTGTTTTGCGCGGGTAACGTCGCTAATGAATCAAGCAGCGATTTAACCGTTACCGAGGTTTCTGGCTTCAGCGTAGTCAGAGACGCTAAcc                                                        >  1:124049/1‑250 (MQ=255)
                                                                                                                                                                                             aaCCCGCCAGCAAATCCACCCTGACGGATTTCCTGTTTATCTTCGGCAAACCAGCTCCGGGTTTTCGCCCGTACCCGACTGCCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCACATCCAGTAATACTTCCTGAGCTTTTTTGTTCGCGGGTAACGTCTCTAATGATTCAAGCAGCGATGTAACCGTTACCGAGTGGTCTGGTTTCTGCGTAGTCATATTCGCTAACCACGCAGGTGCCt                                            <  2:366258/250‑1 (MQ=255)
                                                                                                                                                                                                         aaaCCGCCCTGACGGCTTTCCGGTTTCTCATCGGCAAACCAGCTCCTGGTTTTCACCCGTCCCCGACTACCGATCATGGATGCAACACCGATTTTCGCCTCGTCAAGCCCCACATCCAGTAAGACTTCCGGAGCTTTGTTTTGCGCGGGTAACGTCGCTAATGAATCAAGCAGCGATGTAACCGTTCCCGAGGTTTCTGGCTTCAGCGTAGTCAGAGACGCTAGCCACCCAGGTGCCTGACGCTGTTGCCa                               <  2:18015/251‑1 (MQ=255)
                                                                                                                                                                                                              accCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCACATCCAGTAAGACTTCCTGAGCTTTGTTTTGCGCGGGTAACGTCGCTAATGAATCAAGCAGCGATTTAACCGTTACCGAGGTTTCTGGCTTCAGCGTAGTCAGAGACGCTAACCACGCAGGTGCCTGACGCTGTTGCCAGTAtt                          >  1:485739/1‑251 (MQ=255)
                                                                                                                                                                                                               cccTGCCGGTTTTCATGTTTATCATCGGCAAACCAGCTCCGGGCTTTCACCCGTACCCGCCTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCCCCACATCCCGTAAGACTTCCTGAGCTTTGTTTTGCGCGGGTAACGTCGCTAATGAATCAAGCAGCGATTTAACCGTTACCGAGGTTTCTGGCTTCAGCGTAGTCAGAGACGCTAACCACGCAGGTGCCTGACGCTGTTGCCAGTCTTg                         >  2:310271/1‑251 (MQ=255)
                                                                                                                                                                                                                   gACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCACATCCAGTAAGACTTCCTGAGCTTTGTTTTGCGCGGGTAACGTCGCTAATGAATCAAGCAGCGATTTAACCGTTACCGAGGTTTCTGGCTTCAGCGTAGTCAGAGACGCTAACCACGCAGGTGCCTGACGCTGTTGCCAGTATTGcccc                     >  2:471483/1‑251 (MQ=255)
                                                                                                                                                                                                                   gACGGATTTACTGTTTATCATCGGCAAACCAGCTCCGGGTCTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCACATCCAGTAAGACTTCCTGCGCTTTGTTTTGCGCGGGTCACGCCGCTAATGCATCAAGCAGCGATTTAACCGTTACCGAGGTTTCTGGCTTCAGCGTAGTCAGAGACGCTAACCACGCAGGTGCCTGACGCTGTTGCCAGTATTGcccc                     >  2:58678/1‑251 (MQ=255)
                                                                                                                                                                                                                     cGGATTTCCTGTTTATCATCGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCACATCCAGTAAGACTTCCTGAGCTTTGTTTTGCGCGGGTAACGTCGCTAATGAATCAAGCAGCGATTTAACCGTTACCGAGGTTTCTGGCTTCAGCGTAGTCAGAGACGCTAACCACGCAGGTGCCTGACGCTGTTGCCAGTATTGCCCCca                   <  2:410150/251‑1 (MQ=255)
                                                                                                                                                                                                                                    tcatcGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCACATCCAGTAAGACTTCCTGAGCTTTGTTTTGCGCGGGTAACGTCGCTAATGAATCAAGCAGCGATTTAACCGTTACCGAGGGTTCTGGCTTAAGCGTAGTCAGAGACGCTAACCACGCAGGTGCCTGACGCTGTTGCCAGTATTGCCCCCACACCTCTGCGCCaa     >  1:25600/1‑250 (MQ=255)
                                                                                                                                                                                                                                     catcGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCACATCCAGTAAGACTTCCTGAGCTTTGTTTTGCGCGGGTAACGTCGCTAATGAATCAAGCAGCGATTTAACCGTTACCGAGGTTTCTGGCTTCAGCGTAGTCAGAGACGCTAACCACGCAGGTGCCTGACGCTGTTGCCAGTATTGCCCCCACACATCTGCGCCAat    >  2:557148/1‑250 (MQ=255)
                                                                                                                                                                                                                                      atcGGCAAACCAGCTCCGGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCACAGCCAGTAAGACTCCCTGAGCTTTGTTTTGCGCGGGTAACGTCGCTAATGAATCAAGCAGCGATTTAACCGTTACCGAGGTTTCTTGCTTCAGCGTATTCAGAGACGCTAACCACGCAGGTGCCTGATGCTGTTGCCAGTATTGCCCCCACACATCTGCGCCAatat  <  2:733016/251‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                         
ACTGCCGTCGTACTCTTTTCCGCTGTGATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCACATCCAGTAAGACTTCCTGAGCTTTGTTTTGCGCGGGTAACGTCGCTAATGAATCAAGCAGCGATTTAACCGTTACCGAGGTTTCTGGCTTCAGCGTAGTCAGAGACGCTAACCACGCAGGTGCCTGACGCTGTTGCCAGTATTGCCCCCACACATCTGCGCCAATAT  >  NC_000913/2328962‑2329442

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: