Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,789,632 C→A 5.4% P150P (CCG→CCT waaH ← LPS(HepIII)‑glucuronic acid glycosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,789,6320CA 5.4% 221.5 / 2.0 63P150P (CCG→CCTwaaHLPS(HepIII)‑glucuronic acid glycosyltransferase
Reads supporting (aligned to +/- strand):  ref base C (37/21);  new base A (2/2);  total (40/23)
Fisher's exact test for biased strand distribution p-value = 6.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGACGTTGATAATTAAGGTTTTTATTCCCTTGTCTATGCAACCGACTCACTGACGTATTATGCAGATAATATTTATATAATGATTGCTCGGTATATCGCGCACGCAGCGCGTTAAACATGAATTCTGTTGTCCAGACAATATCCTGATGATGTAATCCGGCAATAAATTTAATGTTATTTTTAACAATAACATCACGGCGATAAACCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGCCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGTCAGCGTTGCACTGCGCCACGTCGAGGTCGTCCTCTAACGCCATGGTCATCAGCGTTTCGTACATGGTGGGATAGACTTCATCGTCAGCATCGACAAAAGCGACATATTTGCCCGTTGCCACTTCAATCCCACGATTACGCGCCACCGATGCGCCAGCATTCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTGCGTA  >  NC_000913/3789369‑3789916
                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                            
tGACGTTGATAATTAAGGTTTTTATTCCCTTGTCTATGCAACCGACTCACTGACGTATTATGCAGATAATATTTATATAATGATTGCTCGGTATATCGCGCACGCAGCGCGTTAAACATGAATTCTGTTGTCCAGACAATATCCTGATGATGTAATCCGGCAATAAATTTAATGTTATTTTTAACAATAACATCACGGCGATAAACCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCATTCAGGGCCGGTTAATACCCCGGTTGAGc                                                                                                                                                                                                                                                                     >  2:337021/1‑289 (MQ=255)
tGACGTTGATAATTAAGGTTTTTATTCCCTTGTCTATGCAACCGACTCACTGACGTATTATGCAGATAATATTTATATAATGATTGCTCGGTATATCGCGCACGCAGCGCGTTAAACATGAATTCTGTTGTCCAGACAATATCCTGATGATGTAATCCGGCAATAAATTTAATGTTATTTTTAACAATAACATCACGGCGATAAACCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGCCGGTTAATACGCCGGTTGAGc                                                                                                                                                                                                                                                                     >  1:159940/1‑289 (MQ=255)
tGACGTTGATAATTAAGGTTTTTATTCCCTTGTCTATGCAACCGACTCACTGACGTATTATGCAGATAATATTTATATAATGATTGCTCGGTATATCGCGCACGCAGCGCGTTAAACATGAATTCTGTTGTCCAGACAATATCCTGATGATGTAATCCGGCAATAAATTTAATGTTATTTTTAACAATAACATCACGGCGATAAACCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGCCGGTTAATACGCCGGTTGAGc                                                                                                                                                                                                                                                                     >  1:69524/1‑289 (MQ=255)
tGACGTTGATAATTAAGGTTTTTATTCCCTTGTCTATGCAACCGACTCACTGACGTATTATGCAGATAATATTTATATAATGATTGCTCGGTATATCGCGCACGCAGCGCGTTAAACATGAATTCTGTTGTCCAGACAATATCCTGATGATGTAATCCGGCAATAAATTTAATGTTATTTTTAACAATAACATCACGGCGATAAACCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCAATCCGCAGCCAGTCTGGGCCGGTTAATACGCCGGTTGAGc                                                                                                                                                                                                                                                                     >  1:73342/1‑289 (MQ=255)
         tAATTAAGGTTTTTATTCCCTTGTCTATGCAACCGACTCACTGACGTATTATGCAGATAATATTTATATAATGATTGCTCGGTATATCGCGCACGCAGCGCGTTAAACATGAATTCTGTTGTCCAGACAATATCCTGATGATGTAATCCGGCAATAAATTTAATGTTATTTTTAACAATAACATCACGGCGATAAACCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGCCGGTTAATACGCCGGTTGAGCGAAGGCGAt                                                                                                                                                                                                                                                            <  1:678165/289‑1 (MQ=255)
                                         ccGACTCACTGACGTATTATGCAGATAATATTTATATAATGATTGCTCGGTATATCGCGCACGCAGCGCGTTAAACATGAATTCTGTTGTCCAGACAATATCCTGATGATGTAATCCGGCAATAAATTTAATGTTATTTTTAACAATAACATCACGGCGATAAACCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGCCGGTTAATACg                                                                                                                                                                                                                                                                               >  2:707099/1‑238 (MQ=255)
                                         ccGACTCACTGACGTATTATGCAGATAATATTTATATAATGATTGCTCGGTATATCGCGCACGCAGCGCGTTAAACATGAATTCTGTTGTCCAGACAATATCCTGATGATGTAATCCGGCAATAAATTTAATGTTATTTTTAACAATAACATCACGGCGATAAACCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGCCGGTTAATACg                                                                                                                                                                                                                                                                               <  1:707099/238‑1 (MQ=255)
                                         ccGACTCACTGACGTATTATGCAGATAATATTTATATAATGATTGCTCGGTATATCGCGCACGCAGCGCGTTAAACATGAATTCTGTTGTCCAGACAATATCCTGATGATGTAATCCGGCAATAAATTTAATGTTATTTTTAACAATAACATCACGGCGATAAACCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGCCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTTCTCCCGtt                                                                                                                                                                                                                            >  1:236810/1‑289 (MQ=255)
                                         ccGACTCACTGACGTATTATGCAGATAATATTTATATAATGATTGCTCGGTATATCGCGCACGCAGCGCGTTAAACATGAATTCTGTTGTCCAGACAATATCCTGATGATGTAATCCGGCAATAAATTTAATGTTATTTTTAACAATAACATCACGGCGATAAACCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGCCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTTCTCCCGtt                                                                                                                                                                                                                            >  1:298160/1‑289 (MQ=255)
                                            aCTCACTGACGTATTATGCAGATAATATTTATATAATGATTGCTCGGTATATCGCGCACGCAGCGCGTTAAACATGAATTCTGTTGTCCAGACAATATCCTGATGATGTAATCCGGCAATAAATTTAATGTTATTTTTAACAATAACATCACGGCGATAAACCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGCCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTTCTCCCGTTTc                                                                                                                                                                                                                          >  1:555548/1‑288 (MQ=255)
                                            aCTCACTGACGTATTATGCAGATAATATTTATATAATGATTGCTCGGTATATCGCGCACGCAGCGCGTTAAACATGAATTCTGTTGTCCAGACAATATCCTGATGATGTAATCCGGCAATAAATTTAATGTTATTTTTAACAATAACATCACGGCGATAAACCCCCATCCAGACAACGTAACTCCAACGGCGCGAAGCAAGCCCCATCCGCAGCCACTCCGGGCCGGTTACTACGCCGGTTGAGCGAAGGCACTCGGTGGGGATGGATTGCCCGGTTCCTCCCGTTTc                                                                                                                                                                                                                          >  2:624449/1‑288 (MQ=255)
                                                         ttatGCAGATAATATTTATATAATGATTGCTCGGTATATCGCGCACGCAGCGCGTTAAACATGAATTCTGTTGTCCAGACAATATCCTGATGATGTAATCCGGCAATAAATTTAATGTTATTTTTAACAATAACATCACGGCGATAAACCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGCCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGt                                                                                                                                                                                                            <  1:553764/289‑1 (MQ=255)
                                                         ttatGCAGATAATATTTATATAATGATTGCTCGGTATATCGCGCACGCAGCGCGTTAAACATGAATTCTGTTGTCCAGACAATATCCTGATGATGTAATCCGGCAATAAATTTAATGTTATTTTTAACAATAACATCACGGCGATAAACCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGCCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGt                                                                                                                                                                                                            <  1:520918/289‑1 (MQ=255)
                                                         ttatGCAGATAATATTTATATAATGATTGCTCGGTATATCGCGCACGCAGCGCGTTAAACATGAATTCTGTTGTCCAGACAATATCCTGATGATGTAATCCGGCAATAAATTTAATGTTATTTTTAACAATAACATCACGGCGATAAACCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCAGGGCCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGt                                                                                                                                                                                                            <  2:298160/289‑1 (MQ=255)
                                                         ttatGCAAGTAATATTTATATAATGATTACTCGGTATATCGCGCACGCAGCGCGTTAAACATGAATTCTGTTGTCCAGACAATATCCTGATGATGTAATCCGGCAATAAATTTAATGTTATTTTTAACAATAACATCACGGCGATAAACCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGCCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGt                                                                                                                                                                                                            <  2:616356/289‑1 (MQ=255)
                                                                 ataataTTTATATAATGATTGCTCGGTATATCGCGCACGCAGCGCGTTAAACATGAATTCTGTTGTCCAGACAATATCCTGATGATGTAATCCGGCAATAAATTTAATGTTATTTTTAACAATAACATCACGGCGATAAACCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGCCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGTCAGCGTTg                                                                                                                                                                                                    >  1:379158/1‑289 (MQ=255)
                                                                        ttATATAATGATTGCTCGGTATATCGCGCACGCAGCGCGTTAAACATGAATTCTGTTGTCCAGACAATATCCTGATGATGTAATCCGGCAATAAATTTAATGTTATTTTTAACAATAACATCACGGCGATAAACCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGCCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGTCAGCGTTGCACTgcg                                                                                                                                                                                             <  2:555548/289‑1 (MQ=255)
                                                                                tGATTGCTCGGTATATCGCGCACGCAGCGCGTTAAACATGAATTCTGTTGTCCAGACAATATCCTGATGATGTAATCCGGCAATAAATTTAATGTTATTTTTAACAATAACATCACGGCGATAAACCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGCCGGt                                                                                                                                                                                                                                                                                      >  2:111054/1‑192 (MQ=255)
                                                                                tGATTGCTCGGTATATCGCGCACGCAGCGCGTTAAACATGAATTCTGTTGTCCAGACAATATCCTGATGATGTAATCCGGCAATAAATTTAATGTTATTTTTAACAATAACATCACGGCGATAAACCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGCCGGt                                                                                                                                                                                                                                                                                      <  1:111054/192‑1 (MQ=255)
                                                                                                     acgcaGCGCGTTAAACATGAATTCTGTTGTCCAGACAATATCCTGATGATGTAATCCGGCAATAAATTTAATGTTATTTTTAACAATAACATCACGGCGATAAACCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGCCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGTCAGCGTTGCACTGCGCCACGTCGAGGTCGTCCTCTAACGCCATg                                                                                                                                                                >  1:453109/1‑289 (MQ=255)
                                                                                                       gcaGCGCGTTAAACATGAATTCTGTTGTCCAGACAATATCCTGATGATGTAATCCGGCAATAAATTTAATGTTATTTTTAACAATAACATCACGGCGATAAACCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGCCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGTCAGCGTTGCACTGCGCCACGTCGAGGTCGTCCTCTAACGCCATGGt                                                                                                                                                              <  2:159940/289‑1 (MQ=255)
                                                                                                            gcgTTAAACATGAATTCTGTTGTCCAGACAATATCCTGATGATGTAATCCGGCAATAAATTTAATGTTATTTTTAACAATAACATCACGGCGATAAACCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGCCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGTCAGCGCTGCACTGCGCCACGTCGAGGTCGTCCTCTAACGCCATGGtcatca                                                                                                                                                         <  1:479834/289‑1 (MQ=255)
                                                                                                                           tCTGTTGTCCAGACAATATCCTGATGAGTTAATCCGGCAATAAATTTAATGTTATTTTTAACAATAACATCACGGCGATAAACCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGCCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGTCAGCGTTGCACTGCGCCACGTCGAGGTCGTCCTCTAACGCCATGGTCATCAGCGTTTCGTACAtgg                                                                                                                                          <  1:337021/289‑1 (MQ=255)
                                                                                                                               ttgtCCAGACAATATCCTGATGATGTAATCCGGCAATAAATTTAATGTTATTTTTAACAATAACATCACGGCGATAAACCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGCCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGTCAGCGTTGCACTGCGCCACGTCGAGGTCGTCCTCTCACGCCATGGTCATCAGCGTTTCGTACATGGTggg                                                                                                                                      >  1:238989/1‑289 (MQ=255)
                                                                                                                                            atCCTGATGATGTAATCCGGCAATAAATTTAATGTTATTTTTAACAATAACATCACGGCGATAAACCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCACATCCGCAGCCAGTCCGGGCCGGTTAATACGCCGGTTGAGCGAAGGCGAACGGTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACAAGTCAGCGTTGCACTGCGCTACGTCGAGGTCGTCCTCTAACGCCATGGCAATCAGCGTTTCGTACATGGTGGGATTGACttaaacg                                                                                                                         >  2:719792/1‑284 (MQ=255)
                                                                                                                                              ccTGATGATGTAATCCGGCAATAAATTTAATGTTATTTTTAACAATAACATCACGGCGATAAACCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGCCGGTTAATACGCCGGTTGAGCGAAGGCGCTCGGAGGGGATGGATCGCCAGGTTTCCCCCGTTTCACGAAAACACCAGTCAGCGTTGCACTGCGCCACGTCGAGGTCGTCCTCTAACGCCATGGTCACCAGCGTTTCGTACATGGTCGGATAGACTTCATCGTc                                                                                                                       >  1:641902/1‑289 (MQ=255)
                                                                                                                                                 gatgatGTAATCCGGCAATAAATTTAATGTTATTTTTAACAATAACATCACGGCGATAAACCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGCCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGTCAGCGTTGCACTGCGCCACGTCGAGGTCGTCCTCTAACGCCATGGTCATCAGCGTTTCGTACAtgg                                                                                                                                          <  1:485156/267‑1 (MQ=255)
                                                                                                                                                 gatgatGTAATCCGGCAATAAATTTAATGTTATTTTTAACAATAACATCACGGCGATAAACCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGCCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGTCAGCGTTGCACTGCGCCACGTCGAGGTCGTCCTCTAACGCCATGGTCATCAGCGTTTCGTACAtgg                                                                                                                                          >  2:485156/1‑267 (MQ=255)
                                                                                                                                                 gatgatGTAATCCGGCAATAAATTTAATGTTATTTTTAACAATAACATCACGGCGATAAACCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGCCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGTCAGCGTTGCACTGCGCCACGTCGAGGTCGTCCTCTAACGCCATGGTCATCAGCGTTTCGTACATGGCGGGACAGACTTCATCGTccg                                                                                                                     >  2:262567/1‑286 (MQ=255)
                                                                                                                                                     atgTAATCCGGCAATAAATTTAATGTTATTTTTAACAATAACATCACGCCGATAAACCACCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCAATCCGCAGCCAGTCCGGGCCGGTTAATACGCCGGTTGAGCGAAGGCGACCGGTGGGGATGGATTGCCAGGTTTCTCCCGTTCCACGAAAACACCAGTCAGCGTTGCACTTCGCCACGCCGAGGTCTTCCTCTAACGCCACGGTCATCCGCGGTTCGTACATGGCGGGACCGCCTCCATCGTCAGCTTCg                                                                                                                >  2:679643/1‑289 (MQ=255)
                                                                                                                                                       gTAATCCGGCAATAAATTTAATGTTATTTTTAACAATAACATCACGGCGATAAACCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGCCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGTCAGCGTTGCACTGCGCCACGTCGAGGTCGTCCTCTAACGCCATGGTCATCAGCGTTTCGTACATGGTGGGATAGACTTCATCGTCAGCATCGAc                                                                                                              >  2:388810/1‑289 (MQ=255)
                                                                                                                                                       gTAATCCGGCAATAAATTTAATGTTATTTTTAACAATAACATCACGGCGATAAACCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGCCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGTCAGCGTTGCACTGCGCCACGTCGAGGTCGTCCTCTAACGCCATGGTCATCAGCGTTTCGTACATGGTGGGATAGACTTCATCGTCAGCATCGAc                                                                                                              >  1:316581/1‑289 (MQ=255)
                                                                                                                                                               gCAATAAATTTAATGTTATTTTTAACAATAACATCACGGCGATAAACCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGCCGGTTAATACGCCGGTTGAGCGAAGGCGATCgg                                                                                                                                                                                                                                                         >  2:637573/1‑142 (MQ=255)
                                                                                                                                                               gCAATAAATTTAATGTTATTTTTAACAATAACATCACGGCGATAAACCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGCCGGTTAATACGCCGGTTGAGCGAAGGCGATCgg                                                                                                                                                                                                                                                         <  1:637573/142‑1 (MQ=255)
                                                                                                                                                                                      aaCAATAACATCACGGCGATAAACCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGCCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGTCAGCGTTGCACTGCGCCACGTCGAGGTCGTCCTCTAACGCCATGGTCATCAGCGTTTCGTACATGGTGGGATAGACTTCATCGTCAGCATCGACAAAAGCGACATATTTTCGCGTTGCCACCTCa                                                                               >  2:643213/1‑289 (MQ=255)
                                                                                                                                                                                           tAACATCACGGCGATAAACCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGCCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGTCAGCGTTGCACTGCGCCACGTCGAGGTCGTCCTCTAACGCCATGGTCATCAGCGTTTCGTACATGGTGGGATAGACTTCATCGTCAGCATCGACAAAAGCGACATATTTGCCCGTTGCCACTTCAATcc                                                                           >  2:123596/1‑288 (MQ=255)
                                                                                                                                                                                                  aCGGCGATAAACCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGCCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGTCAGCGTTGCACTGCGCCACGTCGAGGTCGTCCTCTAACGCCATGGTCATCAGCGTTTCGTACATGGTGGGATAGACTTCATCGTCAGCATCGACAAAAGCGACATATTTGCCCGTTGCCACTTCAATCCCACGAtt                                                                    >  1:447320/1‑288 (MQ=255)
                                                                                                                                                                                                   cGGCGATAAACCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGCCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGTCAGCGTTGCACTgcgc                                                                                                                                                                                            >  1:625351/1‑167 (MQ=255)
                                                                                                                                                                                                   cGGCGATAAACCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGCCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGTCAGCGTTGCACTgcgc                                                                                                                                                                                            <  2:625351/167‑1 (MQ=255)
                                                                                                                                                                                                        aTAAACCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGCCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGTCAGCGTTGCACTGCGCCACGTCGAGGTCGTCCTCTAACGCCATGGTCATCAGCGTTTCGTACATGGTGGGATAGACTTCATCGTCAGCATCGACAAAAGCGACATATTTGCCCGTTGc                                                                                      >  2:349019/1‑264 (MQ=255)
                                                                                                                                                                                                        aTAAACCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCAGGGCCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGTCAGCGTTGCACTGCGCCACGTCGAGGTCGTCCTCTAACGCCATGGTCATCAGCGTTTCGTACATGGTGGGATAGACTTCATCGTCAGCATCGACaaa                                                                                                           >  2:516968/1‑243 (MQ=255)
                                                                                                                                                                                                        aTAAACCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCAGGGCCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGTCAGCGTTGCACTGCGCCACGTCGAGGTCGTCCTCTAACGCCATGGTCATCAGCGTTTCGTACATGGTGGGATAGACTTCATCGTCAGCATCGACaaa                                                                                                           <  1:516968/243‑1 (MQ=255)
                                                                                                                                                                                                        aTAAACCCCCATCCACACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGCCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGTCAGCGTTGCACTGCGCCACGTCGAGGTCGTCCTCTAACGCCATGGTCATCAGCGTTTCGTACATGGTGGGATAGACTTCATCGTCAGCATCGACAAAAGCGACATATTTGCCCGTTGc                                                                                      <  1:349019/264‑1 (MQ=255)
                                                                                                                                                                                                         tAAACCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGCCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGTCAGCGTTGCACTGCGCCACGTCGAGGTCGTCCTCTAACGCCATGGTCATCAGCGTTTCGTACATGGTGGGATAGACTTCATCGTCAGCATCGACAAAAGCGACATATTTGCCCGTTGCCACTTCAATCCCACGATTACGCGcc                                                             >  2:520910/1‑288 (MQ=255)
                                                                                                                                                                                                                       gACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGCCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGTCAGCGTTGCACTGCGCCACGTCGAGGTCGTCCTCTAACGCCATGGTCATCAGCGTTTCGTACATGGTGGGATAGACTTCATCGTCAGCATCGACAAAAGCGACATATTTGCCCGTTGc                                                                                      <  2:90239/249‑1 (MQ=255)
                                                                                                                                                                                                                       gACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGCCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGTCAGCGTTGCACTGCGCCACGTCGAGGTCGTCCTCTAACGCCATGGTCATCAGCGTTTCGTACATGGTGGGATAGACTTCATCGTCAGCATCGACAAAAGCGACATATTTGCCCGTTGc                                                                                      >  1:90239/1‑249 (MQ=255)
                                                                                                                                                                                                                        aCAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGCCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGTCAGCGTTGCACTGCGCCACGTCGAGGTCGTCCTCTAACGCCATGGTCATCAGCGTTTCGTACATGGTGGGATAGACTTCATCGTCAGCATCGACAAAAGCGACATATTTGCCCGTTGCCACTTCCATCCCACGATTACCCGCCCCCGATGCGCCAGCAt                                             >  2:195380/1‑289 (MQ=255)
                                                                                                                                                                                                                        aCAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGCCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGTCAGCGTTGCACTGCGCCACGTCGAGGTCGTCCTCTAACGCCATGGTCATCAGCGTTTCGTACATGGTGGGATAGACTTCATCGTCAGCATCGACAAAAGCGACATATTTGCCCGTTGCCACTTCAATCCCACGATTACGCGCCACCGATGCGCCAGCAt                                             >  2:66293/1‑289 (MQ=255)
                                                                                                                                                                                                                            cGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGCCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGTCAGCGTTGCACTGCGCCACGTCGAGGTCGTCCTCTAACGCCATGGTCATCAGCGTTTCGTACATGGTGGGATAGACTTCATCGTCAGCATCGACAAAAGCGACATATTTGCCCGTTGCCACTTCAATCCCACGATTACGCGCCACCGATGCGCCAGCATTCg                                          >  2:365460/1‑288 (MQ=255)
                                                                                                                                                                                                                            cGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGCCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGTCAGCGTTGCACTGCGCCACGTCGAGGTCGTCCTCTAACGCCATGGTCATCAGCGTTTCGTACATGGTGGGATAGACTTCATCGTCAGCATCGACAAAAGCGACATATTTGCCCGTTGCCACTTCAATCCCACGATTACGCGCCACCGATGCGCCAGCATTCg                                          >  1:359599/1‑288 (MQ=255)
                                                                                                                                                                                                                            cGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGCCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGTCAGCGTTGCACTGCGCCACGTCGAGGTCGTCCTCTAACGCCATGGTCATCAGCGTTTCGTACATGGTGGGATAGACTTCATCGTCAGCATCGACAAAAGCGACATATTTGCCCGTTGCCACTTCAATCCCACGATTACGCGCCACCGATGCGCCAGCATTCg                                          >  1:701711/1‑288 (MQ=255)
                                                                                                                                                                                                                            cGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGCCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGTCAGCGTTGCACTGCGCCACGTCGAGGTCGTCCTCTAACGCCATGGTCATCAGCGTTTCGTACATGGTGGGATAGACTTCATCGTCAGCATCGACAAAAGCGACATATTTGCCCGTTGCCACTTCAATCCCACGATTACGCGCCACCGATGCGCCAGCATTCg                                          >  2:199678/1‑288 (MQ=255)
                                                                                                                                                                                                                             gTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGCCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGTCAGCGTTGCACTGCGCCACGTCGAGGTCGTCCTCTAACGCCATGGTCATCAGCGTTTCGTACATGGTGGGATAGACTTCATCGTCAGCATCGACAAAAGCGACATATTTGCCCGTTGCCACTTCAATCCCACGATTACGCGCCACCGATGCGCCAGCATTCGcc                                        <  1:123596/289‑1 (MQ=255)
                                                                                                                                                                                                                                 gTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGCCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGTCAGCGTTGCACTGCGCCACGTCGAGGTCGTCCTCTAACGCCATGGTCATCAGCGTTTCGTACATGGTGGGATAGACTTCATCGTCAGCATCGACAAAAGCGACATATTTGCCCGTTGCCACTTCAATCCCACGATTACGCGCCACCGATGCGCCAGCATTCGCCTGa                                     <  1:195380/288‑1 (MQ=255)
                                                                                                                                                                                                                                    cAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGCCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGTCAGCGTTGCACTGCGCCACGTCGAGGTCGTCCTCTAACGCCATGGTCATCAGCGTTTCGTACATGGTGGGATAGACTTCATCGTCAGCATCGACAAAAGCGACATATTTGCCCGTTGCCACTTCAATCCCACGATTACGCGCCACCGATGCGCCAGCATTCGCCTGATGca                                 <  1:365460/289‑1 (MQ=255)
                                                                                                                                                                                                                                       cGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGCCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGTCAGCGTTGCACTGCGCCACGTCGAGGTCGTCCTCTAACGCCATGGTCATCAGCGTTTCGTACATGGTGGGATAGACTTCATCGTCAGCATCGACAAAAGCGACATATTTGCCCGTTGCCACTTCAATCCCACGATTACGCGCCACCGa                                                        <  2:526793/263‑1 (MQ=255)
                                                                                                                                                                                                                                       cGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGCCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGTCAGCGTTGCACTGCGCCACGTCGAGGTCGTCCTCTAACGCCATGGTCATCAGCGTTTCGTACATGGTGGGATAGACTTCATCGTCAGCATCGACAAAAGCGACATATTTGCCCGTTGCCACTTCAATCCCACGATTACGCGCCACCGa                                                        >  1:526793/1‑263 (MQ=255)
                                                                                                                                                                                                                                            cgAAGAAAGCCCCATCCGCAGCCAGTCCGGGCCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGTCAGCGTTGCACTGCGCCACGTCGAGGTCGTCCTCTAACGCCATGGTCATCAGCGTTTCGTACATGGTGGGATAGACTTCATCGTCAGCATCGACAAAAGCGACATATTTGCCCGTTGCCACTTCAATCCCACGATTACGCGCCACCGATGCGCCAGCATTCGCCTGATGCAACAaacga                         >  2:501972/1‑289 (MQ=255)
                                                                                                                                                                                                                                             gaagaaAGCCCCATCCGCAGCCAGTCCGGGCCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGTCAGCGTTGCACTGCGCCACGTCGAGGTCGTCCTCTAACGCCATGGTCATCAGCGTTTCGTACATGGTGGGATAGACTTCATCGTCAGCATCGACAAAAGCGACATATTTGCCCGTTGCCACTTCAATCCCACGATTACGCGCCACCGATGCGCCAGCATTCGCCTGATGCAACAaacgaa                        <  1:388810/289‑1 (MQ=255)
                                                                                                                                                                                                                                                         aTCCGCAGCCAGTCCGGGCCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGTCAGCGTTGCACTGCGCCACGTCGAGGTCGTCCTCTAACGCCATGGTCATCAGCGTTTCGTACATGGTGGGATAGACTTCATCGTCAGCATCGACAAAAGCGACATATTTGCCCGTTGCCACTTCAATCCCACGATTACGCGCCACCGATGCGCCAGCATTCGCCTGATGCAACAAACGAACGTGCGGATAGt            >  1:94000/1‑289 (MQ=255)
                                                                                                                                                                                                                                                          tCCGCAGCCAGTCCGGGCCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGTCAGCGTTGCACTGCGCCACGTCGAGGTCGTCCTCTAACGCCATGGTCATCAGCGTTTCGTACATGGTGGGATAGACTTCATCGTCAGCATCGACAAAAGCGACATATTTGCCCGTTGCCACTTCAATCCCACGATTACGCGCCACCGATGCGCCAGCATTCGCCTGATGCAACAAACGAACGTGCGGATAGtt           >  1:386861/1‑289 (MQ=255)
                                                                                                                                                                                                                                                            cGCAGCCAGTCCGGGCCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGTCAGCGTTGCACTGCGCCACGTCGAGGTCGTCCTCTAACGCCATGGTCATCAGCGTTTCGTACATGGTGGGATAGACTTCATCGTCAGCATCGACAAAAGCGACATATTTGCCCGTTGCCACTTCAATCCCACGATTACGCGCCACCGATGCGCCAGCATTCGCCTGATGCAACAAACGAACGTGCGGATAGtttt         <  2:447320/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                   aGTCCGGGCCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGTCAGCGTTGCACTGCGCCACGTCGAGGTCGTCCTCTAACGCCATGGTCATCAGCGTTTCGTACATGGTGGGATAGACTTCATCGTCAGCATCGACAAAAGCGACATATTTGCCCGTTGCCACTTCAATCCCACGATTACGCGCCACCGATGCGCCAGCCTTCGCCTGATGCAACAAACGAACGTTCGGATAGTTTTATGCGTa  >  2:593461/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                            
TGACGTTGATAATTAAGGTTTTTATTCCCTTGTCTATGCAACCGACTCACTGACGTATTATGCAGATAATATTTATATAATGATTGCTCGGTATATCGCGCACGCAGCGCGTTAAACATGAATTCTGTTGTCCAGACAATATCCTGATGATGTAATCCGGCAATAAATTTAATGTTATTTTTAACAATAACATCACGGCGATAAACCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGCCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGTCAGCGTTGCACTGCGCCACGTCGAGGTCGTCCTCTAACGCCATGGTCATCAGCGTTTCGTACATGGTGGGATAGACTTCATCGTCAGCATCGACAAAAGCGACATATTTGCCCGTTGCCACTTCAATCCCACGATTACGCGCCACCGATGCGCCAGCATTCGCCTGATGCAACAAACGAACGTGCGGATAGTTTTCTGCGTA  >  NC_000913/3789369‑3789916

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