Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,320,007 G→A H219H (CAC→CAT phnI ← carbon‑phosphorus lyase core complex subunit PhnI

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,320,0070GA98.0% 152.4 / NA 51H219H (CAC→CATphnIcarbon‑phosphorus lyase core complex subunit PhnI
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (21/29);  total (22/29)
Fisher's exact test for biased strand distribution p-value = 4.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTACTCCGGGGCTTGCAGAGCGCGGTCGACCAGCGCCATCGCCATCGCTTTGCGCTCGCTCATGCCGAACACCAGCCCGTAGCCGCGCGTGAAGTGCGGCGGCTCACCCGGCGGGTCGATAAAACCGTTAACCATTTCACACTCAGTCATCAGCAGTTCGCCGACGTTCACCGCAAATCCCAGCTCTTCCGGCACAATCGACACGTCGATATAGCCGCTGCGGATCTCGCCTGCGAACGGGTGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGTCGCCGCGCATCAACTGCTGCAAACGGGAGGAGCGTGAGCAGGGGTAAACCGGCGGCGTGCGGGTGATGTCATCCGGCTGTGCGCCGCTATCCTCTTCAAACTTCGCCAGCCCCTGACGCGCCAGCAGGCTGAAAACGTGCGGCGACGGCTGCTGTTCGCTGTCGGCGGTGGTCAGCGTCGGCGCTTC  >  NC_000913/4319766‑4320252
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                     
cGTACTCCGGGGCTTGCAGAGCGCGGTCGACCAGCGCCATCGCCATCGCTTTGCGCTCGCTCATGCCGAACACCAGCCCGTAGCCGCGCGTGAAGTGCGGCGGCTCACCCGGCGGGTCGATAAAACCGTTAACCATTTCACACTCAGTCATCAGCAGTTCGCCGACGTTCACCGCAAATCCCAGCTCTTCCGGCACAATCTACACGTCGATATAGCCGCTGCGGATCTCGCCTGCGAACGGATGATTgcgc                                                                                                                                                                                                                                              <  1:178751/251‑1 (MQ=255)
 gTACTCCGGGGCTTGCAGAGCGCGGTCGACCAGCGCCATCGCCATCGCTTTGCGCTCGCTCATGCCGAACACCAGCCCGTAGCCGCGCGTGAAGTGCGGCGGCTCACCCGGCGGGTCGATAAAACCGTTAACCATTTCACACTCAGTCATCAGCAGTTCGCCGACGTTCACCGCAAATCCCAGCTCTTCCGGCACAATCGACACGTCGATATAGCCGCTGCGGATCTCGCCTGCGAACGGATGATTGCGcc                                                                                                                                                                                                                                             <  1:101137/251‑1 (MQ=255)
               gCAGAGCGCGGTCGACCAGCGCCATCGCCATCGCTTTGCGCTCGCTCATGCCGAACACCAGCCCGTAGCCGCGCGTGAAGTGCGGCGGCTCACCCGGCGGGTCGATAAAACCGTTAACCATTTCACACTCAGTCATCAGCAGTTCGCCGACGTTCACCGCAAATCCCAGCTCTTCCGGCACAATCGACACGTCGATATAGCCGCTGCGGATCTCGCCTGCGAACGGATGATTGCGCCCGTAACCGCGTTgg                                                                                                                                                                                                                               <  1:21748/251‑1 (MQ=255)
                  gagCGCGGTCGACCAGCGCCATCGCCATCGCTTTGCGCTCGCTCATGCCGAACACCAGCCCGTAGCCGCGCGTGAAGTGCGGCGGCTCACCCGGCGGGTCGATAAAACCGTTAACCATTTCACACTCAGTCATCAGCAGTTCGCCGACGTTCACCGCAAATCCCAGCTCTTCCGGCACAATCGACACGTCGATATAGCCGCTGCGGATCTCGCCTGCGAACGGATGATTGCGCCCGTAACCGCGTTGGGTg                                                                                                                                                                                                                            >  1:87266/1‑251 (MQ=255)
                    gcgcGGTCGACCAGCGCCATCGCCATCGCTTTGCGCTCGCTCATGCCGAACACCAGCCCGTAGCCGCGCGTGAAGTGCGGCGGCTCACCCGGCGGGTCGATAAAACCGTTAACCATTTCACACTCAGTCATCAGCAGTTCGCCGACGTTCACCGCAAATCCCAGCTCTTCCGGCACAATCGACACGTCGATATAGCCGCTGCGGATCTCGCCTGCGAACGGATGATTGCGCCCGTAACCGCGTTGGGTGGa                                                                                                                                                                                                                          <  1:252756/251‑1 (MQ=255)
                       cgGTCGACCAGCGCTATCGCCCTCGCTTTGCGCTCGCTCATGCCGAACACCAGCCCGTAGACGCGCGTGAAGTGCGGCGGCTCACCCGGCGGCTCGATAAAACCGTTAACCATTTCACACTCAGTCATCAGCAGGTCGCCGACGTTCACCGCAAATCCCAGCTCTTCCGGCACAATCGCCACGTCGCTACAGCCGCTGCGGATCTCGCCTGCGAACGGATGATTGCGCCCGTCACCGCGCTGGGTGGAggg                                                                                                                                                                                                                       >  2:170265/1‑249 (MQ=255)
                            gACCAGCGCCATCGCCATCGCTTTGCGCTCGCTCATGCCGAACACCAGCCCGTAGCCGCGCGTGAAGTGCGGCGGCTCACCCGGCGGGTCGATAAAACCGTTAACCATTTCACACTCAGTCATCAGCAGTTCGCCGACGTTCACCGCAAATCCCAGCTCTTCCGGCACAATCGACACGTCGATATAGCCGCTGCGGATCTCGCCTGCGAACGGATGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGgcca                                                                                                                                                                                                                  <  1:61462/251‑1 (MQ=255)
                            gACCAGCGACATCGCCATCGCTTTGCGCTCGCTCATGCCGAACACCAGCCCGTAGCCGCGCGTGAAGTGCGGCGGCTCACCCGGCGGGTCGATAAAACCGTTAACCATTTCACACTCAGTCATCAGCAGTTCGCCGACGTTCACCGCAAATCCCAGCTCTTCCGGCACAATCGACACGTCGATATAGCCGCTGCGGATCTCGCCTGCGGACGGATGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGgcc                                                                                                                                                                                                                   >  1:87383/1‑250 (MQ=255)
                             cccatcgccatcgccatCGCTTTGCGCTCGCTCATGCCGAACACCAGACCGTAGCCGCGCGGGAAGTGCGGCGGCTCACCCGGCGGGTCGATAAAACCGTTAACCATTTCACACTCAGTCATCAGCAGTGCGCCGACGTTCGCCGCAAATCCCAGCTCTACCGGCACAATCGACACGTCGATATAGCCGCTGCGGATCTCGCCTGCGAACGGATGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGgcca                                                                                                                                                                                                                  <  2:55075/245‑1 (MQ=255)
                              ccAGCGCCATCGCCATCGCTTTGCGCTCGCTCATGCCGAACACCAGCCCGTAGCCGCGCGTGAAGTGCGGCGGCTCACCCGGCGGGTCGATAAAACCGTTAACCATTTCACACTCAGTCATCAGCAGTTCGCCGACGTTCACCGCAAATCCCAGCTCTTCCGGCACAATCGACACGTCGATATAGCCGCTGCGGATCTCGCCTGCGAACGGATGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGgccagc                                                                                                                                                                                                                >  1:90868/1‑251 (MQ=255)
                                    ccatcgccatcgCTTTGCGCTCGCTCATGCCGAACACCAGCCCGTAGCCGCGCGTGAAGTGCGGCGGCTCACCCGGCGGGTCGATAAAACCGTTAACCATTTCACACTCAGTCATCAGCAGTTCGCCGACGTTCACCGCAAATCCCAGCTCTTCCGGCACAATCGACACGTCGCTATAGCCGCCGCGGATCTCGCCTGCGAACGGATGATTGCGCCCGTAACCGCGTTGGGGGGAGTAGgccagcgcaagc                                                                                                                                                                                                          >  1:210401/1‑251 (MQ=255)
                                        cgccatcgcTTTGAGCTCGCTCATGCCGAACACCAGCCCGTAGCCGCGCGTGAAGTGCGGCGGCTCACCCGGCGGGTCGATAAAACCGTTAACCATTTCACACTCAGTCATCAGCAGTTCGCCGACGTTCACCGCAAATCCCAGCTCTTCCGGCACAATCGACACGTCGATATAGCCGCTGCGGATCTCGCCTGCGAACGGATGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAAt                                                                                                                                                                                                      >  1:217095/1‑251 (MQ=255)
                                                   tGCGCTCGCTCATGCCGAGCACCAGCCCGTAGCCGCGCGTGAAGGGCGGCGGCTCGCCCGGCGGGTCGACAAAACCGTTAACCCTTTCTCACTCAGTCATCAGCAGTTCGCCGACGTTCACCGCAAACCCCAGCTCTTCCGGCACAATCGACACGCCGATATAGCCGCGGCGGATCTCGCCTGCGAACGGATGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTtcgtc                                                                                                                                                                                            <  2:124471/250‑1 (MQ=255)
                                                       ctcgctcATGCCGAACACCAGCCCGTAGCCGCGCGTGAAGTGCGGCGGCTCACCCGGCGGGTCGATAAAACCGTTAACCATTTCACACTCAGTCATCAGCAGTTCGCCGACGTTCACCGCAAATCCCAGCTCTTCCGGCACAATCGACACGTCGATATAGCCGCTGCGGATCTCGCCTGCGAACGGATGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGTcgccgc                                                                                                                                                                                       >  2:91011/1‑251 (MQ=255)
                                                             cATGCCGAACACCAGCCCGTAGCCGCGCGTGAAGTGCGGCGGCTCACCCGGCGGGTCGATAAAACCGTTAACCATTTCACACTCAGTCATCAGCAGTTCGCCGACGTTCACCGCAAATCCCAGCTCTTCCGGCACAATCGACACGTCGATATAGCCGCTGCGGATCTCGCCTGCGAACGGATGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGTCGCCGCGCATCa                                                                                                                                                                                 <  2:171215/251‑1 (MQ=255)
                                                              aTGCCGAACACCAGCCCGTAGCCGCGCGTGAAGTGCGGCGGCTCACCCGGCGGGTCGATAAAACCGTTAACCATTTCACACTCAGTCATCAGCAGTTCGCCGACGTTCACCGCAAATCCCAGCTCTTCCGGCACAATCGACACGTCGATATAGCCGCTGCGGATCTCGCCTGCGAACGGATGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGTCGCCGCGCATCaa                                                                                                                                                                                <  2:122211/251‑1 (MQ=255)
                                                                              cGTAGCCGCGCGTGAAGTGCGGCGGCTCACCCGGCGGGTCGATAAAACCGTTAACCATTTCACACTCAGTCATCAGCAGTTCGCCGACGTTCACCGCAAATCCCAGCTCTTCCGGCACAATCGACACGTCGATATAGCCGCTGCGGATCTCGCCTGCGAACGGATGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGTCGCCGCGCATCAACTGCTGCAAACGggag                                                                                                                                                                <  1:231875/251‑1 (MQ=255)
                                                                              cGTAGACGCGCGTGCAGTGCGGCGGCTCACCCGGCGGGTCGATAAAACCGTTAACCATTTCACACTCAGTCATCAGCAGTTCGCCGACGTTCACCGCACATCTCAGCTCTTCCGGCACAATCGACACGTCGATATAACCGCTGCGGATCTCGCCTGCGAACGGGTAAGTGCGCCCGTAACCGCGGTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGTCGCCGCGCATTAACTGCTGCCAACggg                                                                                                                                                                  >  2:122907/1‑249 (MQ=255)
                                                                                                  ggcggcTCACCCGGCGGGTCGATAAAACCGTTAACCATTTCACACTCAGTCATCAGCAGTGCGCCGACGTTCACCGCAAATCCCAGCTCTTCCGGCACAATCGACACGTCGATATAGCCGCTGCGGATCTCGCCTGCGAACGGATGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCGAGCAAATAGCCTGCGTCGCCGCGCATCAACTGCTGCAAACGGGAGGAGCGTGAGCAGGGGTAAAc                                                                                                                                            <  1:253624/251‑1 (MQ=255)
                                                                                                     ggcTCACCCGGCGGGTCGATAAAACCGTTAACCATTTCACACTCAGTCATCAGCAGTTCGCCGACGTTCACCGCAAATCCCAGCTCTTCCGGCACAATCGACACGTCGATATAGCCGCTGCGGATCTCGCCTGCGAACGGATGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCGTCGTCGCCGCGCATCAACTGCTGCAAACGGGAGGAGCGTGAGCAGGGGTAAACcgg                                                                                                                                         <  2:58343/251‑1 (MQ=255)
                                                                                                         cACCCAGCGGGTCGATCAAACCGTTAACCATTTCACACTCAGTCATCAGCAGTTCGCCGCCGTTCACCGCAAATCCCAGCTCTTCCGGCACAATCGACACGTCGATATAACCGCTGCGGGTCTCGCCTGCGAACGGATGATTGCGCCCGTAACCGCGTTGGGGGGAGTAGGACAGCGCCCGCAAATAGCCTTCGTCGCCGCGCATCAACTGCTGCCAACGGGAGGGGCGTGAGCAGGGGTAAACCggcggc                                                                                                                                     >  2:94603/1‑251 (MQ=255)
                                                                                                            ccGGCGGGTCGATAAAACCGTTAACCATTTCACACTCAGTCATCAGCAGTTCGCCGACGTTCACCGCAAATCCCAGCTCTTCCGGCACAATCGACACGTCGATATAGCCGCTGCGGATCTCGCCTGCGAACGGATGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGTCGCCGCGCATCAACTGCTGCAAACGGGAGGAGCGTGAGCAGGGGTAAACCGGCGGCGTg                                                                                                                                  <  1:68235/251‑1 (MQ=255)
                                                                                                             cggcggGTCGATAAAACCGTTAACCATTTCACACTCAGTCATCAGCAGTTCGCCGACGTTCACCGCAAATCCCAGCTCTTCCGGCACAATCGACACGTCGATATAGCCGCTGCGGATCTCGCCTGCGAACGGATGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGTCGCCGCGCATCAACTGCTGCAAACGGGAGGAGCGTGAGCAGGGGTAAACCGGCGGCGTGc                                                                                                                                 <  1:67158/251‑1 (MQ=255)
                                                                                                              ggcggGTCGATAAAACCGTTAACCATTTCACACTCAGTCATCAGCAGTTCGCCGACGTTCACCGCAAATCCCAGCTCTTCCGGCACAATCGACACGTCGATATAGCCGCTGCGGATCTCGCCTGCGAACGGATGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGTCGCCGCGCATCAACTGCTGCAAACGGGAGGAGCGTGAGCAGGGGGAAACCGGCGGCGTGCg                                                                                                                                >  2:3609/1‑251 (MQ=255)
                                                                                                                cggGTCGATAAAACCGTTAACCATTTCACACTCAGTCATCAGCAGTTCGCCGACGTTCACCGCAAATCCCAGCTCTTCCGGCACAATCGACACGTCGATATAGCCGCTGCGGATCTCGCCTGCGAACGGATGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGTCGCCGCGCATCAACTGCTGCAAACGGGAGGAGCGTGAGCAGGGGTAAACCGGCGGCGTGCggg                                                                                                                              >  1:62410/1‑251 (MQ=255)
                                                                                                                          aaaCCGTTAACCATTTCACACTCAGCCATCAGCAGTTCGCCGGCGTTCAACGCAAATCCCAGCTCTTCCGGCACAATCGACACGTCGATATAGCCGCTGCGGATCTCGCCTGCGAACGGATGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGTCGCCGCGCATCAACTGCTGCAAACGGGAGGAGCGTGAGCAGGGGTAAACCGGCGGCGTGCGGGTGATGTCATc                                                                                                                    <  2:210401/251‑1 (MQ=255)
                                                                                                                              cGTTAACCATTTCACACTCAGTCATCAGCAGTTCGCCGACGTTCACCGCAAATCCCAGCTCTTCCGGCACAATCGACACGTCGATATAGCCGCTGCGGATCTCGCCTGCGAACGGATGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGTCGCCGCGCATCAACTGCTGCAAACGGGAGGAGCGTGAGCAGGGGGAAACCGGCGGCGTGCGGGTGATGTCATCCGGc                                                                                                                >  2:196469/1‑251 (MQ=255)
                                                                                                                                  aaCCATTTCACACTTAGTCATCAGCAGTTCGCCGACGGTCACCGCAAATCCCAGCTCTTCCGGCACAATGGACACTTCGATATAGCCGCTGCGGATGGCGACTGCGAATGGATGATTGCGCCCGTAACCGCGTGGGGTGGAGTAGGCCAGTGCCAGCTAATTGCCTCCGTCGCCGCGTATCAACTGCTGCAAACGGGATGAGGGGGAGCAGGGGTAAACCGGCGGCGTGCGGGTGATGTCATCCGGCtgt                                                                                                             <  1:158801/250‑1 (MQ=255)
                                                                                                                                  aaCCATTTCACACTCAGTCATCAGCAGTTCGCCGACGTTCACCGCAAATCCCAGCTCTTCCGGCACAATCGACACGTCGATATAGCCGCTGCGGATCTCGCCTGCGAACGGATGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGTCGCCGCGCATCAACTGCTGCAAACGGGAGGAGCGTGAGCAGGGGTAAACCGGCGGCGTGCGGGTGATGTCATCCGGCtgt                                                                                                             <  2:37614/250‑1 (MQ=255)
                                                                                                                                        ttCACACTCAGTCATCAGCAGTTCGCCGACGTTCACCGCAAATCCCAGCTCTTCCGGCACAATCGACACGTCGATATAGCCGCTGCGGATCTCGCCTGCGAACGGATGATTGCGCCCGTAACCGCGTTGGGTGGCGTACGCCCGCGCCAGCAAATAGCCTTCGTCGCCGCGCATCAACTGCTGCTAACGGGCGGAGCGTGAGCCGGGGTACACCGGCGGCGTGCGGGTGATGTCATCCGGCTGTGcgccgc                                                                                                      >  2:64748/1‑251 (MQ=255)
                                                                                                                                          cacaCTCAGTCATCAGCAGTTCGCCGACGTTCACCGCAAATCCCAGCTCTTCCGGCACAATCGACACGTCGATATAGCCGCTGCGGATCTCGCCTGCGAACGGATGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGTCGCCGCGCATCAACTGCTGCAAACGGGAGGAGCGTGAGAAGGGGTAAACCGGCGGCGTGCGGGTGATGTCATCCGGCTGTGCGCCGCtg                                                                                                    >  1:145821/1‑250 (MQ=255)
                                                                                                                                                 aGTCATCAGCAGTTCGCCGACGTTCACCGCAAATCCCAGCTCTTCCGGCACAATCGACACGTCGATATAGCCGCTGCGGATCTCGCCTGCGAACGGATGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGTCGCCGCGCATCAACTGCTGCAAACGGGAGGAGCGTGAGCAGGGGTAAACCGGCGGCGTGCGGGTGATGTCATCCGGCTGTGCGCCGCTATCCTCtt                                                                                             <  2:215893/251‑1 (MQ=255)
                                                                                                                                                 aGTCATCAGCAGTTCGACGACGTTCAACGCACATCCCAGCTCTTCCGGCACAATCGACCCGTCGATATAGCCGCTGCGGATCTCGCCTGCGAACGGATGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGTCGCCGCGCATCAACTGCTGCACACGGGAGGAGCGTGAGCAGGGGTAACCCGGCGGCGTGCGGGTGCTGTCATCCGGCTGAGCGCCGCTATCCTCtt                                                                                             >  2:169815/1‑251 (MQ=255)
                                                                                                                                                   ccatcaGCAGTTCGCCGACGTTCACCGCAAATCCCAGCTCTTCCGGCACAATCGACACGTCGATATAGCCGCTGCGGATCTCGCCTGCGAACGGATGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGTCGCCGCGCATCAACTGCTGCAAACGGGAGGAGCGTGAGCAGGGGTAAACCGGCGGCGTGCGGGTGATGTCATCCGGCTGTGCGCCGCTATCCTCTTCa                                                                                           <  2:251209/250‑1 (MQ=255)
                                                                                                                                                        agcagTTCGCCGACGTTCACCGCAAATCCCAGCTCTTCCGGCACAATCGACACGTCGATATAGCCGCTGCGGATCTCGCCTGCGAACGGATGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGTCGCCGCGCATCAACTGCTGCAAACGGGAGGAGCGTGAGCAGGGGTAAACCGGCGGCGTGCGGGTGATGTCATCCGGCTGTGCGCCGCTATCCTCTTCAAACtt                                                                                      >  1:47663/1‑251 (MQ=255)
                                                                                                                                                         gcagTTCGCCGACGTTCACCGCAAATCTCAGCTCTTCCGGCACAATCGACACGGCGATATAGCCGCTGCGGATCTCGCCTGCGAACGGATGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGTCGCCGCGCATCAACTGCTGCAAACGGGAGGAGCGTGAGCAGGGGTAAACCGGCGGCGTGCGGGTGATGTCATCCGGCTGTGCGCCGCTATCCTCTTCAAACTTc                                                                                     <  1:122907/251‑1 (MQ=255)
                                                                                                                                                                     cGTTCCCCGCAAATCCCAGCGCTCCCGGCACAATCGACACGTCGATATAGCCGCTGCGGATCTCGCCTGCGAACGGATGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGTCGCCGCGCATCAACTGCTGCAAACGGGAGGAGCGTTAGCAGGGGTAAACCGGCGGCGTGCGGGTGATGTCATCCGGCTGTGCGCCGCTATCCTCTTCAAACTTCGCCAGCCCCTGa                                                                         <  2:34993/251‑1 (MQ=255)
                                                                                                                                                                      gttccccccAAACCCCAGCTCTGACGGCGCAATCTACACGTCGATATCGCCGCTGCGGATCGCGCCTGCGAACGGATGATTGCGCCCGTAGCCGCGTTGGGTGGAGTGGGCCAGCGCCAGCAAAAAGCCTTCGGCGCCGCGCATCAACTGCTGCAAACGGGAGGAGCGTGAGCAGGGGTAAACCGGCGCCGTGCGGGTGATGTCATCCGGCTGTGCGCCGCTAGCCTTTGCAAACCTCGCCAGCCCCTGAc                                                                        <  2:94919/243‑1 (MQ=255)
                                                                                                                                                                      gtcccccGCAAATCCCAGCTCTTCCGGCACAAGCGACACGTCGATATAGCCGCTGCGGATCTCGCCTGCGAACGGATGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGTCGCCGCGCATCAACTGCGGCAAACGGGAGGAGCGTGAGCAGGGGTAAACCTGCGGCGTGCGGGTGATGTCATCCGGCTGTGCGCCGCTTTCCTCTTCAAACTTCGCCAGCCCCTGAc                                                                        <  2:65170/246‑1 (MQ=255)
                                                                                                                                                                      gTTCACCGCAAATCACAGCTCTTCTGGCACAATCGACACGTCGCTATAGCCGCTGCGGATCTCGCCTGCGAACGGATGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGTCGCCGCGCATCAACTGCTGCAAACGGGAGGAGCGTGAGCAGGGGTAAACCGGCGGCGTGCGGGTGATGTCAACCGGCTGTGCGCCGCTATCCTCTTCAAACTTCGCCAGCCCCTGAc                                                                        >  1:208722/1‑251 (MQ=255)
                                                                                                                                                                       ttCACCGCAAATCCCAGCTCTTCCGGCACAATCGACACGTCGATATAGCCGCTGCGGATCTCGCCTGCGAACGGATGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGTCGCCGCGCATCAACTGCTGCAAACGGGAGGAGCGTGAGCAGGGGTAAACCGGCGGCGTGCGGGTGATGTCATCCGGCTGTGCGCCGCTATCCTCTTCAAACTTCGCCAGCCCCTGAcg                                                                       >  1:13206/1‑251 (MQ=255)
                                                                                                                                                                           ccGCAAATCCCCGCTCTTCCGGCACAATCGACACGTCGATATAGCCGCTGCGGATCTCGCCTGCGAACGGATGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGTCGCCGCTCATCAACTGCTGCAAACGGGAGGAGCGTGAGCAGGGGTAAACCGGCGGCGTGCGGGTGATGTCATCCGGCTGTGCGCCGCTATCCTCTTCAAACTTCGCCAGCCCCTGACGCGcc                                                                   <  2:174112/251‑1 (MQ=255)
                                                                                                                                                                           ccGCAAATCCCAGCTCTTCCGGCACAATCGACACGTCGATATAGCCGCTGCGGATCTCGCCTGCGAACGGATGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGTCGCCGCGCATCAACTGCTGCAAACGGGAGGAGCGTGAGCAGGGGTAAACCGGCGGCGTGCGGGTGATGTCATCCGGCTGTGCGCCGCTATCCTCTTCAAACTTCGCCAGCCCCTGACGCGcc                                                                   >  1:74634/1‑251 (MQ=255)
                                                                                                                                                                                             ccGGCACAATCGACACGTCGATATAGCCGCTGCGGATCTCGCCTGCGAACGGATGATTGCGCCCGTACCCGCGTTGGGTGGAGTAGGCCAGCGCCAGCCAATAGCCTTCGTCGCCGCGCATCCACTGCTGCAAACGGGAGGAGCGTGAGCAGGGGTAAACCGGCGGCGTGCGGGTGCTGTCCTCCGGATGTGCGCCGCTATCCTCTTCAAACTTCGCCAGCCCCTGACGCGCCAGCAGGCTGAAAACGTgc                                                 >  2:241882/1‑251 (MQ=255)
                                                                                                                                                                                                   caATCGACACGTCGATATAGCCGCTGCGGATCTCGCCTGCGAACGGATGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGTCGCCGCGCATCAACTGCTGCAAACGGGAGGAGCGTGAGCAGGGGTAAACCGGCGGCGTGCGGGTGATGTCATCCGGCTGTGCGCCGCTATCCTCTTCAAACTTCGCCAGCCCCTGACGCGCCAGCAGGCTGAAAACGTGCGGCGAc                                           <  1:19952/251‑1 (MQ=255)
                                                                                                                                                                                                                gATATAGCCGCTGCGGATCTCGCCTGCGAACGGATGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGTCGCCGCGCATCAACTGCTGCAAACGGGAGGAGCGTGAGCAGGGGTAAACCGGCGGCGTGCGGGTGATGTCCTCCGGCTGTGCGCCGCTATCCTCTTCAAACTTCGCCAGCCCCTGACGCGCCAGCAGGCTGAACACGTGCGGCGCCGGCTGCTGTTCGc                              >  2:238396/1‑251 (MQ=255)
                                                                                                                                                                                                                gATATAGCCGCTGCGGATCTCGCCTGCGAACGGATGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGTCGCCGCGCATCAACTGCTGCAAACGGGAGGAGCGTGAGCAGGGGTAAACCGGCGGCGTGCGGGTGATGTCATCCGGCTGTGCGCCGCTATCCTCTTCAAACTTCGCCAGCCCCTGACGCGCCAGCAGGCTGAAAACGTGCGGCGACGGCTGCTGGTCGc                              >  2:49547/1‑251 (MQ=255)
                                                                                                                                                                                                                              ggATCTCGCCTGCGAACGGATGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGTCGCCGCGCATCAACTGCTGCAAACGGGAGGAGCGTGAGCAGGGGTAAACCGGCGGCGTGCGGGTGATGTCATCCGGCTGTGCGCCGCTATCCTCTTCAAACTTCGCCAGCCCCTGACGCGCCAGCAGGCTGAAAACGTGCGGCGACGGCTGCTGTTCGCTGTCGGCGGTGGTc                <  1:91011/251‑1 (MQ=255)
                                                                                                                                                                                                                                   ccGCCTGCGAACGGATGATTGCGCCCGTAACCGCGATGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGTCGCCGCGCAGCAACTGCTGCAAACGGGAGGAGCGTGAGCAGGGGTAAACCGGCGGCGTGCGGGTGATGTCATCCGGCTGTGCGCCGCTATCCTCTTCAAACTTCGCCAGCCCCTGACGCGCCAGCAGGCTGAAAACGTGCGGCGACGGCTGCTGTTCGCTGTCGGCGGTGGTCAGCGt           <  1:94603/250‑1 (MQ=255)
                                                                                                                                                                                                                                        ggCGAACGGATGAGTGCGCCCTTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGTCGCCGCGCATCAACTGTTGCAAACGGGAGGGGCGTGAGCAGGGGTAAACCGGCGGCGTGCGGGTGATGTCAGCCGGCTGTGCGCCGGTATCTTCTTCAAACTTCGCCAGCGCCTTACGCGCCAGCGGGCTGAAATCGTTCGGCGACGGCTTCTGTTCGCTGTCGGCGGTGGTCATCGTCGgcg      <  2:136200/250‑1 (MQ=255)
                                                                                                                                                                                                                                            aaCGGATGATTGCGCCCGTAACCGCGTTGGGGGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGTCGCCGCGCATCAACTGCGGCAAACGGGAGGAGCGTGAGCAGGGGTAAACCGGCGGCGTGCGGGTGATGTCATCCGGCTGTGCGCCGCTATCCTCTTCAAACGTCGCCAGCCCCTGACGCGCCAGCAGGCTGAAAACGTGCGGCGACGGCTGCTGTTCGCTGTCGGCGGTGGTCAGCGTCGGCGCTTc  <  2:74634/251‑1 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                     
CGTACTCCGGGGCTTGCAGAGCGCGGTCGACCAGCGCCATCGCCATCGCTTTGCGCTCGCTCATGCCGAACACCAGCCCGTAGCCGCGCGTGAAGTGCGGCGGCTCACCCGGCGGGTCGATAAAACCGTTAACCATTTCACACTCAGTCATCAGCAGTTCGCCGACGTTCACCGCAAATCCCAGCTCTTCCGGCACAATCGACACGTCGATATAGCCGCTGCGGATCTCGCCTGCGAACGGGTGATTGCGCCCGTAACCGCGTTGGGTGGAGTAGGCCAGCGCCAGCAAATAGCCTTCGTCGCCGCGCATCAACTGCTGCAAACGGGAGGAGCGTGAGCAGGGGTAAACCGGCGGCGTGCGGGTGATGTCATCCGGCTGTGCGCCGCTATCCTCTTCAAACTTCGCCAGCCCCTGACGCGCCAGCAGGCTGAAAACGTGCGGCGACGGCTGCTGTTCGCTGTCGGCGGTGGTCAGCGTCGGCGCTTC  >  NC_000913/4319766‑4320252

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