Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,628,984 A→G D513D (GAT→GAC ettA ← energy‑dependent translational throttle protein EttA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,628,9840AG93.3% 176.2 / 1.0 60D513D (GAT→GACettAenergy‑dependent translational throttle protein EttA
Reads supporting (aligned to +/- strand):  ref base A (2/2);  new base G (34/22);  total (36/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01

CAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTACAGTTTTACAACGTATTGAAATTGCTGGTTTTGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGCGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCGGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAATCCAGAATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAGAGTTCACCAACGCGTTTACCCTGATCAACCCCTTT  >  NC_000913/4628739‑4629223
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                               
cAATGCCTGTTTCGACGCTTGCCGCGTCTTATCAGTCCTACAGTTTTACAACGTATTGAAATTGCTGGTGTTGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCATTTTACGCATTACGTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGCGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCGGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAATcc                                                                                                                                                                                                                                              <  2:474216/249‑1 (MQ=255)
           gCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTg                                                                                                                                                                                                                                  <  1:54351/250‑1 (MQ=255)
             gACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACATCGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCg                                                                                                                                                                                                                               <  1:177602/251‑1 (MQ=255)
               cGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGAt                                                                                                                                                                                                                             >  1:113038/1‑251 (MQ=255)
                gCTTACCGCGTCTTATCAGGCCTACAGTTTTACAACGTATTGAAATTGCTGGTTTTGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGCGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCGGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAATCCAGAATGTGCGTGGCGATa                                                                                                                                                                                                                            <  1:383986/251‑1 (MQ=255)
                 cTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATAc                                                                                                                                                                                                                           >  1:391172/1‑251 (MQ=255)
                  ttACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACg                                                                                                                                                                                                                          >  2:365662/1‑251 (MQ=255)
                   tACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACg                                                                                                                                                                                                                          <  2:104500/250‑1 (MQ=255)
                         gTCTTAGCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAg                                                                                                                                                                                                                    <  1:482679/250‑1 (MQ=255)
                           cTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCCTTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGa                                                                                                                                                                                                                  >  2:424237/1‑250 (MQ=255)
                             tATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGa                                                                                                                                                                                                                  <  2:272887/248‑1 (MQ=255)
                                               aCAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATa                                                                                                                                                                                              >  2:291022/1‑250 (MQ=255)
                                                 aaCGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGCGGAACCACCGGTCGTGCGAGATAAcc                                                                                                                                                                                           >  2:348484/1‑251 (MQ=255)
                                                   cGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCAt                                                                                                                                                                                         >  1:309921/1‑251 (MQ=255)
                                                              tcttGGTTTTTATGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCg                                                                                                                                                                              <  1:46511/247‑1 (MQ=255)
                                                                cttggtttGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACAATGGCACAGCCCggg                                                                                                                                                                            >  2:355752/6‑251 (MQ=255)
                                                                      ttGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTcc                                                                                                                                                                      >  1:167771/1‑251 (MQ=255)
                                                                        gTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTcc                                                                                                                                                                      >  1:30400/1‑249 (MQ=255)
                                                                        gTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTcc                                                                                                                                                                      >  2:332605/1‑249 (MQ=255)
                                                                          aGGCCGGATAAGGCGTTCGCGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGCGCGTCTGCGCCCGGCGTGCGTTTCTTGTACTCTTCGTACTCGGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAATCCAGAATGTACGTGGCGATACGGTCGAGGCACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCcagc                                                                                                                                                                   >  1:229571/1‑250 (MQ=255)
                                                                                 aaaaGGCTTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGtt                                                                                                                                                           <  2:392149/249‑1 (MQ=255)
                                                                                  tAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGATACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGGGATAACCATCGCACAGCCCGGGAACTACAGCAGGGAGttt                                                                                                                                                          >  1:482050/1‑251 (MQ=255)
                                                                                   aaGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGtttt                                                                                                                                                         >  1:343457/1‑251 (MQ=255)
                                                                                        gTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAgc                                                                                                                                                    >  1:404432/1‑251 (MQ=255)
                                                                                        gTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAgc                                                                                                                                                    <  2:184969/251‑1 (MQ=255)
                                                                                                   aTCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAgggg                                                                                                                                         >  2:256636/1‑250 (MQ=255)
                                                                                                     ccGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGttt                                                                                                                                       <  1:291022/251‑1 (MQ=255)
                                                                                                      cGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTc                                                                                                                                      <  1:443985/251‑1 (MQ=255)
                                                                                                          aTTTTACTCATTACTTCGCAATAAGCTTGTACTTGCTACGCTTCGGCTCCAGTGCGTCTGCCCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGCTACGGTCGAGGAACCACCGGTCGTGCGCGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTGCCACCGCGCACAGGGTATCggtc                                                                                                                                  >  2:388173/1‑248 (MQ=255)
                                                                                                                            cAATACACTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCACACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGGCGGtggttggt                                                                                                                >  1:381907/1‑251 (MQ=255)
                                                                                                                              attatcTTGTACGTGAGTCGCTTCGGCTCCAGTGAGTCTGCGCCCAGCGTTCGTTTCTTGTACTCTTCGTACTCAGGAAAGTTACCTTCGAAGAACGCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTTGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGTAACTCCAGCAGGGCGTTTTCCAGCGCGCGCATTGTTTCGATAGCCAGGTCGTTGGTTGGTtc                                                                                                              <  2:180784/246‑1 (MQ=255)
                                                                                                                                      gTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGCGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGagc                                                                                                       >  2:481546/1‑250 (MQ=255)
                                                                                                                                      gTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGagc                                                                                                       <  1:256636/250‑1 (MQ=255)
                                                                                                                                      gTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGACGTTGGTTGGTTCGTCGagc                                                                                                       >  2:480315/1‑250 (MQ=255)
                                                                                                                                                 gCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCCGGGTTTCGATCTCCAGGTCGTTGGTTGGTTCGTCGACCAGCAGCATGTTg                                                                                           >  2:229929/1‑251 (MQ=255)
                                                                                                                                                             gtgcgTCTGCGCCCAGCGTGAGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTAACTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCACCGCGAGCAGGGTTTCGATATCCAGGTGGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAGCCTgc                                                                               >  2:85895/3‑251 (MQ=255)
                                                                                                                                                               gcgTCTGCGCACAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGCGGAACCACCGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCCGCGCGCGCAGGGTTTCGCTATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTgc                                                                               >  2:438666/1‑249 (MQ=255)
                                                                                                                                                                 gTCTGCGCCCAGCGGGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGATGGTTGGTTCGTGGAGCAGCAGCATGTTGCCGCCAACCTgcagc                                                                            >  2:310750/1‑250 (MQ=255)
                                                                                                                                                                          cAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCtt                                                                       >  1:111755/1‑246 (MQ=255)
                                                                                                                                                                          cAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCtt                                                                       <  2:111755/246‑1 (MQ=255)
                                                                                                                                                                           aGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAg                                                                 <  1:79724/251‑1 (MQ=255)
                                                                                                                                                                             cgtgcgtTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGcc                                                                   >  1:178866/1‑247 (MQ=255)
                                                                                                                                                                               tgcgtTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCCGATGc                                                             >  2:367317/1‑251 (MQ=255)
                                                                                                                                                                                      attGTACTTTTCGTACTTAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGGCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTTGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAgacga                                                       <  2:161356/249‑1 (MQ=255)
                                                                                                                                                                                                      tCGGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAATCCAGAATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCCCCACCgg                                       >  2:349583/1‑250 (MQ=255)
                                                                                                                                                                                                       cagTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGga                                      >  1:186915/3‑250 (MQ=255)
                                                                                                                                                                                                               ttACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCAcc                             >  2:428778/1‑251 (MQ=255)
                                                                                                                                                                                                                  ccgtCGAAGAACACAATTTTACCTGCACCCTGGTAGTCCAGGATGTGCGTGGCGCTTCGGTCGAGGAACCAACGGTCGTGTGAGAGAACCATCGCACAGCCCGGTAACGCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCATCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGTTCACCCCCGGAAAGTTCACCAAc                          <  2:67165/248‑1 (MQ=255)
                                                                                                                                                                                                                  ccTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAAc                          <  1:235967/251‑1 (MQ=255)
                                                                                                                                                                                                                  ccTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAAc                          >  1:388400/1‑251 (MQ=255)
                                                                                                                                                                                                                     tCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAAcgc                        <  2:309921/250‑1 (MQ=255)
                                                                                                                                                                                                                           aaCTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTAcc                 <  2:391172/251‑1 (MQ=255)
                                                                                                                                                                                                                                aaCTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGAt            <  2:178866/251‑1 (MQ=255)
                                                                                                                                                                                                                                    ttACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAAc        >  2:81027/1‑251 (MQ=255)
                                                                                                                                                                                                                                    ttACCTTCAACCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCCGCAGCTTCGCCAGATGCAGACGACCGCGCTCCCCACCGGAACGTTCACCAACGCGTTTCCCCTGATCAAc        >  2:197030/1‑251 (MQ=255)
                                                                                                                                                                                                                                         ttCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCCTGAGACATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTCTTCAAGAGCGCGCAGGGTTTCGATATCCAGGACGTTGGTTGGTTCGTCGAGCAGCAGCATGATGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGCCCGCGCTCAACACCGGACAGTTCACCAACGCGTTTACCCTGATCAACCCCtt   >  1:214858/1‑251 (MQ=255)
                                                                                                                                                                                                                                          tCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCt    <  1:206899/249‑1 (MQ=255)
                                                                                                                                                                                                                                           cATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCttt  <  1:407893/250‑1 (MQ=255)
                                                                                                                                                                                                                                           cATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCttt  >  1:470675/1‑250 (MQ=255)
                                                                                                                                                                                                                                             tCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGtt                                 >  1:127493/1‑217 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                               
CAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTACAGTTTTACAACGTATTGAAATTGCTGGTTTTGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGCGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCGGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAATCCAGAATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAGAGTTCACCAACGCGTTTACCCTGATCAACCCCTTT  >  NC_000913/4628739‑4629223

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