Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,104,337 T→C 7.7% E58G (GAG→GGG)  wbbK ← lipopolysaccharide biosynthesis protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,104,3370TC 7.7% 188.8 / 2.2 67E58G (GAG→GGG) wbbKlipopolysaccharide biosynthesis protein
Reads supporting (aligned to +/- strand):  ref base T (36/23);  new base C (3/3);  total (40/27)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.04e-01

GTATAGCATTTTAAATAAGAAAAAGCTAGGCTCCATAAGAATTTCACGGAATAAAATTCCTTTATAAAAAGGTGCAGGGTTATGACAATACACATATCTTTTTTTAGTGACGACATTGGCCGTAATATCATGCAGACAAATCCAATGCGTAGCATTCAGCTCTTTTGAAAGTTTTTTACAAACTACATATTCAAAGTGCAAACGTTTTAGCCACGACCCTTTAACCTCAGGAAACTCAATGAATTTAACCCATGGATAACTTTCTTTTAACTCTTTAGCAGAATGGACTAATGCGATAAAACTGACATTTTCTTTATTATTAGTTGCTGCCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGACTATGCTTTTTCCCATAATTTGGTCTCATGATTGTATTTCTTTATGATTTTTGCTGGTACTCCCGCAATGACAGT  >  NC_000913/2104102‑2104568
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                       
gTATAGCATTTTAAATAAGAAAAAGCTAGGCTCCATAAGAATTTCACGGAATAAAATTCCTTTATAAAAAGGTGCAGGGTTATGACAATACACATATCTTTTTTTAGTGACGACATTGGCCGTAATATCATGCAGACAAATCCAATGCGTAGCATTCAGCTCTTTTGAAAGTTTTTTACAAACTACATATTCAAAGTGCAAACGTTTTAGCCACGACCCTTTAACCTCAGGAAACTCaa                                                                                                                                                                                                                                      <  1:950928/239‑1 (MQ=255)
 tataGCATTTTAAATAAGAAAAAGCTAGGCTCCATAAGAATTTCACGGAATAAAATTCCTTTATAAAAAGGTGCAGGGTTATGACAATACACATATCTTTTTTTAGTGACGACATTGGCCGTAATATCATGCAGACAAATCCAATGCGTAGCATTCAGCTCTTTTGAAAGTTTTTTACAAACTACATATTCAAAGTGCAAACGTTTTAGCCACGACCCTTTAACCTCAGGAAACTCAAt                                                                                                                                                                                                                                     >  2:800223/1‑239 (MQ=255)
 tataGCATTTTAAATAAGAAAAAGCTAGGCTCCATAAGAATTTCACGGAATAAAATTCCTTTATAAAAAGGTGCAGGGTTATGACAATACACATATCTTTTTTTAGTGACGACATTGGCCGTAATATCATGCAGACAAATCCAATGCGTAGCATTCAGCTCTTTTGAAAGTTTTTTACAAACTACATATTCAAAGTGCAAACGTTTTAGCCACGACCCTTTAACCTCAGGAAACTCAAt                                                                                                                                                                                                                                     >  2:262167/1‑239 (MQ=255)
 tataGCATTTTAAATAAGAAAAAGCTAGGCTCCATAAGAATTTCACGGAATAAAATTCCTTTATAAAAAGGTGCAGGGTTATGACAATACACATATCTTTTTTTAGTGACGACATTGGCCGTAATATCATGCAGACAAATCCAATGCGTAGCATTCAGCTCTTTTGAAAGTTTTTTACAAACTACATATTCAAAGTGCAAACGTTTTAGCCACGACCCTTTAACCTCAGGAAACTCAAt                                                                                                                                                                                                                                     >  1:516233/1‑239 (MQ=255)
       aTTTTAAATAAGAAAAAGCTAGGCTCCATAAGAATTTCACGGAATAAAATTCCTTTATAAAAAGGTGCAGGGTTATGACAATACACATATCTTTTTTTAGTGACGACATTGGCCGTAATATCATGCAGACAAATCCAATGCGTAGCATTCAGCTCTTTTGAAAGTTTTTTACAAACTACATATTCAAAGTGCAAACGTTTTAGCCACGACCCTTTAACCTCAGGAAACTCAATGAAttt                                                                                                                                                                                                                               <  1:800223/239‑1 (MQ=255)
        ttttAAATAAGAAAAAGCTAGGCTCCATAAGAATTTCACGGAATAAAATTCCTTTATAAAAAGGTGCAGGGTTATGACAATACACATATCTTTTTTTAGTGACGACATTGGCCGTAATATCATGCAGACAAATCCAATGCGTAGCATTCAGCTCTTTTGAAAGTTTTTTACAAACTACATATTCAAAGTGCAAACGTTTTAGCCACGACCCTTTAACCCACGGAAACACAAAGAAttaa                                                                                                                                                                                                                              >  2:950928/1‑237 (MQ=255)
           tAAATAAGATAAAGCTAGGCTCCATAAGAATTTCACGGAATAAAATTCCTTTATAAAACGGGGCAGGGTTATGACAATACACATATCTTTTTTTAGTGACGACATTGGCCGTAATATCATGCAGACAAATCCAATGCGTAACATTCAGATCTTTTGAAAGTTTTTTCCAAACTACCTATTCAAAGGTCCACAGTTTTAGCCACGCCCCTTTAACCTCCGGAAACCCAATGAATTTAAcc                                                                                                                                                                                                                           >  2:785559/1‑239 (MQ=255)
                             gCTCCATAAGAATTTCACGGAATAAAATTCCTTTATAAAAAGGTGCAGGGTTATGACAATACACATATCTTTTTTTAGTGACGACATTGGCCGTAATATCATGCAGACAAATCCAATGCGTAGCATTCAGCTCTTTTGAAAGTTTTTTCCAAACTACATATTCAAAGTGCAAACGTTTTAGCCACGACCCTTTAACCTCAGGAAAATCAAGGAATTTAACCCATGGAAAAAtttttttt                                                                                                                                                                                                         >  2:777973/1‑239 (MQ=255)
                                ccATAAGTATTTCACGGAATAAAATTCCTTTAAAAAAAGGGTCAGGGTTAGGACAATACATATATCTTTTTTTAGGGACGACAGTGGCCGTAATATCATGCAGACAAATCCAATGCGTAGCATTCAGCTCTTTTGAAAGTTTTTTACAAACTACATATTCAAAGTGCAAACGTTTTAGCCACGACCCTTTAACCTCAGGAAACTCAATGAATTTAACCCATGGATAACTTTCTTTTaa                                                                                                                                                                                                       <  1:785559/238‑1 (MQ=255)
                                       aaTTTCACGGAATAAAATTCCTTTATAAAAAGGTGCAGGGTTATGACAATACACATATCTTTTTTTAGTGACGACATTGGCCGTAATATCATGCAGACAAATCCAATGCGTAGCATTCAGCTCTTTTGAAAGTTTTTTACAAACTACATATTCAAAGTGCAAACGTTTTAGCCACGACCCTTTAACCTCAGGAAACTCAATGAATTTAACCCATGGATAACTTTCTTTTAACTCTTTAg                                                                                                                                                                                               <  1:262167/239‑1 (MQ=255)
                                                 aaTAAAATTCCTTTATAAAAAGGTGCAGGGTTATGACAATACACATATCTTTTTTTAGTGACGACATTGGCCGTAATATCATGCAGACAAATCCAATGCGTAGCATTCAGCTCTTTTGAAAGTTTTTTACAAACTACATATTCAAAGTGCAAACGTTTTAGCCACGACCCTTTAACCTCAGGAAACTCAATGAATTTAACCCATGGATAACTTTCTTTTaa                                                                                                                                                                                                       <  2:87829/221‑1 (MQ=255)
                                                 aaTAAAATTCCTTTATAAAAAGGTGCAGGGTTATGACAATACACATATCTTTTTTTAGTGACGACATTGGCCGTAATATCATGCAGACAAATCCAATGCGTAGCATTCAGCTCTTTTGAAAGTTTTTTACAAACTACATATTCAAAGTGCAAACGTTTTAGCCACGACCCTTTAACCTCAGGAAACTCAATGAATTTAACCCATGGATAACTTTCTTTTaa                                                                                                                                                                                                       >  1:87829/1‑221 (MQ=255)
                                                           cTTTATAAAAAGGTGCAGGGTTATGACAATACACATATCTTTTTTTAGTGACGACATTGGCCGTAATATCATGCAGACAAATCCAATGCGTAGCATTCAGCTCTTTTGAAAGTTTTTTACAAACTACATATTCAAAGTGCAAACGTTTTAGCCACGACCCTTTAACCTCAGGAAACTCAATGAATTTAACCCATGGATAAc                                                                                                                                                                                                                 >  1:176327/1‑201 (MQ=255)
                                                           cTTTATAAAAAGGTGCAGGGTTATGACAATACACATATCTTTTTTTAGTGACGACATTGGCCGTAATATCATGCAGACAAATCCAATGCGTAGCATTCAGCTCTTTTGAAAGTTTTTTACAAACTACATATTCAAAGTGCAAACGTTTTAGCCACGACCCTTTAACCTCAGGAAACTCAATGAATTTAACCCATGGATAAc                                                                                                                                                                                                                 <  2:176327/201‑1 (MQ=255)
                                                           cTTTATAAAAAGGTGCAGGGTTATGACAATACACATATCTTTTTTTAGTGACGACATTGGCCGTAATATCATGCAGACAAATCCAATGCGTAGCATTCAGCTCTTTTGAAAGTTTTTTACAAACTACATATTCAAAGTGCAAACGTTTTAGCCACGACCCTTTAACCTCAGGAAACTCAATGAATTTAACCCATGGATAACTTTCTTTTaa                                                                                                                                                                                                       >  1:118206/1‑211 (MQ=255)
                                                           cTTTATAAAAAGGTGCAGGGTTATGACAATACACATATCTTTTTTTAGTGACGACATTGGCCGTAATATCATGCAGACAAATCCAATGCGTAGCATTCAGCTCTTTTGAAAGTTTTTTACAAACTACATATTCAAAGTGCAAACGTTTTAGCCACGACCCTTTAACCTCAGGAAACTCAATGAATTTAACCCATGGATAACTTTCTTTTaa                                                                                                                                                                                                       <  2:118206/211‑1 (MQ=255)
                                                           cTTTATAAAAAGGTGCAGGGTTATGACAATACACATATCTTTTTTTAGTGACGACATTGGCCGTAATATCATGCAGACAAATCCAATGCGTAGCATTCAGCTCTTTTGAAAGTTTTTTACAAACTACATATTCAAAGTGCAAACGTTTTAGCCACGACCCTTTAACCTCAGGAAACTCAATGAATTTAACCCATGGATAACTTTCTTTTaa                                                                                                                                                                                                       >  1:147029/1‑211 (MQ=255)
                                                                taAAAAGGTGCAGGGTTATGAAAATACACATATCTTTTTTTAGTGACGACATTGGCCGTAATATCATGCAGAAAAATCCAATGCGTAGCATTCAGCTCTTTTGAAAGTTTTTTACAAACTACATATTCAAAGTGCAAACGTTTTAGCCACGACCCTTTAACCTCAGGAAACTCAATGAATTTAACCCATGGATAACTTTCTTTTaa                                                                                                                                                                                                       <  1:920867/206‑1 (MQ=255)
                                                                     aGGTGCAGGGTTATGACAATACACATATCTTTTTTTAGCGACGACATTGGCCGTAATATCATGCAGACAAATCTAAAGCGTAGCATTCAGCTCTTTTGAACGTTTTTTAAAAACTACATATTCAAATTGCAAACGTTTTAGACACGACCCTTTAACCTCAGAAAACCACATGAATTTAACCCAAGGATAACTTTATTTTAACCTTTCAGCAAACATGAAAAATGCGAAAAAAcaaaaat                                                                                                                                                                 >  2:845518/1‑233 (MQ=255)
                                                                          cAGGGTTATGACAATACACATATCTTTTTTTAGTGACGACATTGGCCGTAATATCATGCAGACAAATCCAATGCGTAGCATTCAGCTCTTTTGAAAGTTTTTTACAAACTACATATTCAAAGTGCAAACGTTTTAGCCACGACCCTTTAACCTCAGGAAACTCAATGAATTTAACCCATGGATAACTTTCTTTTAACTCTTTAg                                                                                                                                                                                               >  2:233810/1‑204 (MQ=255)
                                                                          cAGGGTTATGACAATACACATATCTTTTTTTAGTGACGACATTGGCCGTAATATCATGCAGACAAATCCAATGCGTAGCATTCAGCTCTTTTGAAAGTTTTTTACAAACTACATATTCAAAGTGCAAACGTTTTAGCCACGACCCTTTAACCTCAGGAAACTCAATGAATTTAACCCATGGATAACTTTCTTTTAACTCTTTAg                                                                                                                                                                                               <  1:233810/204‑1 (MQ=255)
                                                                                 aTGACAATCCACATATCTTTTTTTAGTGACGACATTGGCCGTAATATCATGCAGACAAATCCAATGCGTAGCATTCAGCTCTTTTGAAAGTTTTTTACAAACTACATATTCAAAGTGCAAACGTTTTAGCCACGACCCTTTAACCTCAGGAAACTCaa                                                                                                                                                                                                                                      <  2:472593/158‑1 (MQ=255)
                                                                                 aTGACAATACACATATCTTTTTTTAGTGACGACATTGGCCGTAATATCATGCAGACAAATCCAATGCGTAGCATTCAGCTCTTTTGAAAGTTTTTTACAAACTACATATTCAAAGTGCAAACGTTTTAGCCACGACCCTTTAACCTCAGGAAACTCaa                                                                                                                                                                                                                                      >  1:472593/1‑158 (MQ=255)
                                                                                       aTACACATATCTTTTTTTAGTGACGACATTGTCCGTAATATCATGCAGACAAATCCAATGCGTAGCATTCGGCTCTTTTGAAAGTTTTTTACAAACTACATATTCAAAGTGCAAACGTTTTAGCCACGACCCTTTAACCTCAGGAAACTCAATGAATTTAACCCATGGATAACTTTCTTTTAACCCTTTAg                                                                                                                                                                                               >  2:382605/1‑191 (MQ=255)
                                                                                       aTACACATATCTTTTTTTAGTGACGACATTGTCCGTAATATCATGCAGACAAATCCAATGCGTAGCATTCAGCTCTTTTGAAAGTTTTTTACAAACTACATATTCAAAGTGCAAACGTTTTAGCCACGACCCTTTAACCTCAGGAAACTCAATGAATTTAACCCATGGATAACTTTCTTTTAACTCTTTAg                                                                                                                                                                                               <  1:382605/191‑1 (MQ=255)
                                                                                       aTACACATATCTTTTTTTAGTGACGACATTGGCCGTAATATCATGCAGACAAATCCAATGCGTAGCATTCAGCTCTTTTGAAAGTTTTTTACAAACTACATATTCAAAGTGCAAACGTTTTAGCCACGACCCTTTAACCTCAGGAAACTCAATGAATTTAACCCa                                                                                                                                                                                                                         >  2:567565/1‑165 (MQ=255)
                                                                                       aTACACATATCTTTTTTTAGTGACGACATTGGCCGTAATATCATGCAGACAAATCCAATGCGTAGCATTCAGCTCTTTTGAAAGTTTTTTACAAACTACATATTCAAAGTGCAAACGTTTTAGCCACGACCCTTTAACCTCAGGAAACTCAATGAATTTAACCCa                                                                                                                                                                                                                         <  1:567565/165‑1 (MQ=255)
                                                                                       aTACACATATCTTTTTTTAGTGACGACATTGGCCGTAATATCATGCAGACAAATCCAATGCGTAGCATTCAGCTCTTTTGAAAGTTTTTTACAAACTACATATTCAAAGTGCAAACGTTTTAGCCACGACCCTTTAACCTCAGGAAACTCAATGAATTTAACCCATGGATAACTTTCTTTTAACTCTTTAGCAGAATGGACTAATGCGATAAAACTGACATTTTTTTTTTTATTAttt                                                                                                                                                >  2:26185/1‑238 (MQ=255)
                                                                                       aTACACATATCTTTTTTTAGTGACGACATTGGCCGTAATATCATGCAGACAAATCCAATGCGTAGCATTCAGCTCTTTTGAAAGTTTTTTACAAACTACATATTCAAAGTGCAAACGTTTTAGCCACGACCCTTTAACCTCAGGAAACTCAATGAATTTAACCCATGGATAACTTTCTTTTAACTCTTTACAAGAATGGACTAAATCCAAAAAAATAAAATTTTATTAATTATAAGAtg                                                                                                                                               >  2:474501/1‑239 (MQ=255)
                                                                                              gatCTTTTTTGAGTGACGACCTTGGCCGTAATATCATGCAGACAAATCCAACGCGTAGCATTCAGCTCTTTTGTAAGTTTTTTACACACTACACATTCAACGTGCAAACGTTTTAGCCACGCCCCTTTAACCTCCGGAAACTCACTGACTTTAACCCATGGATAACTTTCTTTTAACTCTTTTCCCGAATGACCACATCCGTTACAACTGCACTTTTTTTTTTTCTTGTTTGCTTCCaa                                                                                                                                        >  2:110588/2‑239 (MQ=255)
                                                                                                tCTTTTTTTAGTGACGACATTGGCCGTAATATCATGCAGACAAATCCAATGCGTAGCATTCAGCTCTTTTGAAAGTTTTTTACAAACTACATATTCAAAGTGCAAACGTTTTAGCCACGACCCTTTAACCTCAGGAAACTCAATGAATTTAACCCATGGATAACTTTCTTTTAACTCTTTAGCAGAATGGACTAATGCGATAAAACTGACATTTTCTTTATTATTAGTTGCTGCCaaaa                                                                                                                                      >  1:711632/1‑239 (MQ=255)
                                                                                                             acgacATTGGCCTTTATTTCATGCAGACAAATCCAATGCGTAGCATTCAGCTCTTTTGAAAGTTTTTTACAAACTACAGATTCAAAGTGCAAACGTTTTAGCCACGACCCTTTAACCTCAGGAAACTCAATGAATTTAACCCATGGATAACTTTCTTTTAACTCTTTAGCAGAATGGACTAATGCGATAAAACTGACATTTTCTTTATTATTAGTTGCTGCCAAAAATTTTTTCAAAAt                                                                                                                         <  2:516233/239‑1 (MQ=255)
                                                                                                                     ggCCGTAATATCATGCAGACAAATCCAATGCGTAGCATTCAGCTCTTTTGAAAGTTTTTTACAAACTACATATTCAAAGTGCAAACGTTTTAGCCACGACCCTTTAACCTCAGGAAACTCAATGAATTTAACCCATGGATAACTTTCTTTTAACTCTTTAg                                                                                                                                                                                               >  1:235924/1‑161 (MQ=255)
                                                                                                                     ggCCGTAATATCATGCAGACAAATCCAATGCGTAGCATTCAGCTCTTTTGAAAGTTTTTTACAAACTACATATTCAAAGTGCAAACGTTTTAGCCACGACCCTTTAACCTCAGGAAACTCAATGAATTTAACCCATGGATAACTTTCTTTTAACTCTTTAg                                                                                                                                                                                               <  2:235924/161‑1 (MQ=255)
                                                                                                                           aaTATCATGCAGACAAATCCAATGCGTAGCATTCAGCTCTTTTGAAAGTTTTTTACAAACTACATATTCAAAGTGCAAACGTTTTAGCTACGACCCTTTAACCTCAGGAAACTCAATGAATTTAACCCATGGATAACTTTCTTTTAACTCTTTAGCAGa                                                                                                                                                                                           <  2:410255/159‑1 (MQ=255)
                                                                                                                           aaTATCATGCAGACAAATCCAATGCGTAGCATTCAGCTCTTTTGAAAGTTTTTTACAAACTACATATTCAAAGTGCAAACGTTTTAGCCACGACCCTTTAACCTCAGGAAACTCAATGAATTTAACCCATGGATAACTTTCTTTTAACTCTTTAGCAGa                                                                                                                                                                                           >  1:410255/1‑159 (MQ=255)
                                                                                                                           aaTATCATGCAGACAAATCCAATGCGTAGCATTCAGCTCTTTTGAAAGTTTTTTACAAACTACATATTCAAAGTGCAAACGTTTTAGCCACGACCCTTTAACCTCAGGAAACTCAATGAATTTAACCCATGGATAACTTTCTTTTAACTCTTTAGCAGATTGGACTAATGCGATAAAACTGACATTTTCTTTATTATTAGTTGCTGCCAAAAATTTTTTCAAAATGGAAAAGGGACCg                                                                                                            >  2:676270/1‑238 (MQ=255)
                                                                                                                           aaTATCATGCAGAAAAATCCAATGCGTAGCATTCAGCTCTTTTGAAAGTTTTTTACAAACTACATATTCAAAGTGCAAACGTTTTAGCCACGACCCTTTAACCTCAGGAAACTCAATGAATTTAACCCATGGATCACTTCCTTTTAAATCTTTAGCAGAATGGACTACTGCTATAAAACTGCCATTTTCTTTATTATTAGTTGCTGCCAAAAATTTTTTCAAAATGGTAAATGGACCg                                                                                                            >  1:846322/1‑238 (MQ=255)
                                                                                                                                cATGCAGACAAATCCAATGCGTAGCATTCAGCTCTTTTGAAAGTTTTTTACAAACTACATATTCAAAGTGCAAACGTTTTAGCCACGACCCTTTAACCTCAGGAAACTCAATGAATTTAATCCTTGGATAACTGTCTTTTAACTCTTTATCAAAATGGACTAATGCGATAAAACTGCACTTTTCTTTATTATTAGTTGCTGCCCAAAATTTTTTCAAAATGGAAAATAGAACGCCCGTg                                                                                                      >  2:421082/1‑239 (MQ=255)
                                                                                                                                       aCAAAGCCAATGGGTAGAATTCAGCTCTTTTGAAAGTTTTTTACAAACTACATATTCAAAGTGCAAACGTTTTAGCCACGACCCTTTAACCTCAGGAAACTCAATGAATTTAACCCATGGATAACTTTCTTTTAACTCTTTAGCAGAATGGACTAATGCGATAAAACTGACATTTTCTTTATTATTAGTTGCTGCCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAg                                                                                                        <  1:645750/230‑1 (MQ=255)
                                                                                                                                         aaaTCCAATGCGTAGCATTCAGCTCTTTTGAAAGTTTTTTACAAACTACATATTCAAAGTGCAAACGTTTTAGCCAAGACCCTTTAACCTCAGGAAACTCAATGAATTTAACCCATGGATAACTTTCTTTTAACTCTTTAGCAGAATGGACTAATGCGATAAAACTGACATTTTCTTTATAATTAGTTGCTGCCAAAAATTTTTTCAAAATGGTAAATGGACTGCAAGTGGTAAAATTg                                                                                             >  1:296461/1‑239 (MQ=255)
                                                                                                                                                 tGCGTAGCATTCAGCTCTTTTGAAAGTTTTTTACAAACTACATATTCAAAGTGCAAACGTTTTAGCAACGACCCTTTAACCTCAGGAAACTCAATGAATTTAACCCATGGATAACTTTCTTTTAACTCTTTAGCAGAATGGACTAATGCGATAAAACTGACATTTTCTTTTTTTTTAGTTGCTGCCAAAAATTTTTTCAAAATGGTAAATGAACCGCCAGTGGCAAAAtttctgtcttt                                                                                     >  2:645750/1‑230 (MQ=255)
                                                                                                                                                     tAGCATTCAGCTCTTTTGAAAGTTTTTTACAAACTACATATTCAAAGTGCAAACGTTTTAGCCACGACCCTTTAACCTCAGGAAACTCAATGAATTTAACCCATGGATAACTTTCTTTTAACTCTTTAGCAGAATGGACTAATGCGATAAAACTGACATTTTCTTTATTATTAGTTGCTGCCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAaac                                                                                  >  2:618165/1‑238 (MQ=255)
                                                                                                                                                     tAGCATTCAGCTCTTTTGAAAGTTTTTTACAAACTACATATTCAAAGTGCAAACGTTTTAGCCACGACCCTTTAACCTCAGGAAACTCAATGAATTTAACCCATGGATAACTTTCTTTTAACTCTTTAGCAGAATGGACTAATGCGATAAAACTGACATTTTCTTTATTATTAGTTGCTGCCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAaac                                                                                  <  1:618165/238‑1 (MQ=255)
                                                                                                                                                      aGCATTCAGCTCTTTTGAAAGTTTTTTACAAACTACATATTCAAAGTGCAAACGTTTTAGCCACGACCCTTTAACCTCAGGAAACTCAATGAATTTAACCCATGGATAACTTTCTTTTAACTCTTTAGCAGAATGGACTAATGCGATAAAACTGACCTTTTCTTTATTATTAGTTGCTGCCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCCGAacca                                                                                 >  1:308934/1‑238 (MQ=255)
                                                                                                                                                       gCATTCAGCTCTTTTGAAAGTTTTTTACACACTACATATTCAAAGTGCAAACGTTTTAGCCACGACCCTTTAACCTCAGGAAACTCAATGAATTTAACCCATGGATAACTTTCTTTTAACTCTTTAGCAGAATGGACTAATGCGATAAAACTGACATTTTCTTTATTATTAGTTGCTGCCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGCGTTAAAATTGACCGCAGCaaaaac                                                                               >  1:421093/1‑239 (MQ=255)
                                                                                                                                                               ctctTTTGAAAGTTTTTTACAAACTACATATTCAAAGTGCAAACGTTTTAGCCACGACACTTTAACCTCAGGAAACTCAATGAATTTAACCCATGGATCACTTTCTTTTAACTCTTTAGCAGAATGGACTAATGCGATAAAACTGACATTTTCTTTATTATTAGTTGCTGCCAAAAATTTTTTCAAAATGGTAAATGGACCGCACGTGGTAAAATTGACCGCAGAAACAACGAatat                                                                         >  2:151762/1‑237 (MQ=255)
                                                                                                                                                                tctTTTGAAAGTTTTTTACAAACTACAGATACAAAGTGCAACCGTTTTAGCCACGACCCTTTAACCTAAGGAAACACAATGAATTTAACCCTTGGATAGCTTTCTTTTAACCCTTTAGCAGAATGCATTAATGCGATAAAACTGACATTTTCTTTATTATTAGTTGCTCCCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGACTATGCt                                                                      <  1:165866/239‑1 (MQ=255)
                                                                                                                                                                  ttttTAAAGTTTTTTACAAACTAAATATTCAAAGTGCAAACGTTTTAGCCACGACCCTTTAACCTCAGGAAACCCAATGAATTTAACCCATGGATAACTTTCTTTTAACTCTTTAGCAGAATGGACTAATGCGATAAAACTGACATTTTCTTTATTATTAGTTGCTGCCAAAAATTTTTTCAAAATGGTaaa                                                                                                                   <  2:805038/192‑1 (MQ=255)
                                                                                                                                                                  ttttGAAAGTTTTTTACAAACTACATATTCAAAGTGCAAACGTTTTAGCCACGACCCTTTAACCTCAGGAAACCCAATGAATTTAACCCATGGATAACTTTCTTTTAACTCTTTAGCAGAATGGACTAATGCGATAAAACTGACATTTTCTTTATTATTAGTTGCTGCCAAAAATTTTTTCAAAATGGTaaa                                                                                                                   >  1:805038/1‑192 (MQ=255)
                                                                                                                                                                         aGTTTTTTACAAACTACATATTCAAAGTGCAAACGTTTTAGCCACGACCCTTTAACCTCAGGAAACTCAATGAATTTAACCCATGGATAACTTTCTTTTAACTCTTTAGCAGCATGGACTAATGCGATAATACTGACATTTTCTTTATTATTCGTTGCTGCCAAAACTTTTTTCAAAATGGTAAATGGACCGCCAGTGGTCAAATTGACCGCATCAAAAACGCCTAAGCTTTTTCCCAt                                                             >  2:248495/1‑239 (MQ=255)
                                                                                                                                                                                     aCTACATATTCAAAGTGCAAACGTTTTAGCCACGACCCTTTAACCTCAGGAAACTCAATGAATTTAACCCATGGATAACTTTCTTTTAACTCTTTAGCAGAATGGACTAATGCGATAAAACTGACATTTTCTTTATTATTAGTTGCTGCCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGCCTATGCTTTTTCCCATACTTTGGtct                                                   >  1:850644/1‑237 (MQ=255)
                                                                                                                                                                                          atatTCAAAGTGCAAACGTTTTAGCCACGACCCTTTAACCTCAGGAAACTCAATGAATTTAACCCATGGATAACTTTCTTTTAACTCTTTAGCAGAATGGACTAATGCGATAAAACTGACATTTTCTTTATTATTAGTTGCTGCCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGACTATGCTTTTTCCCATAATTTGGTCTCATGAtt                                            <  1:777973/239‑1 (MQ=255)
                                                                                                                                                                                             ttCAAAGTGGAAATGTTTTCGCACGGGCCCTTTACCCCAATGAACCCAAATGAACTTAACCCATGGTAAACTTTGTTTTAACCTTTTAGAAGAATGGACTAATGAGATAAAACTGACTTTTTCTTTATTATTATTTTCTGTCAAAAATTGTTTCAAATTGTTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGACTATGCTTTTTCCCATAATTTGGTCTCATGATTGTa                                         <  2:549708/239‑1 (MQ=255)
                                                                                                                                                                                              tCAAAGTTAAACCGTTTTAGCCCCGACCTTTTAACCCCAGGAAACTCAATGAATTTAACCCATGAATAACTTTTTTTTAACTCTTTAGCAGAATGGACTAATTCGATAAAACTGACATTTTCTTTATTATTAGTTGCTGCCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGACTATGCTTTTTCCCATAATTTGGTCTCATGATTGTAt                                        <  2:993770/239‑1 (MQ=255)
                                                                                                                                                                                                 aaGTGCAAACGTTTTAGCCACGACCCTTTAACCTCAGGAAACTCAATGAATTTAACCCATGGa                                                                                                                                                                                                                     >  1:842558/1‑63 (MQ=255)
                                                                                                                                                                                                 aaGTGCAAACGTTTTAGCCACGACCCTTTAACCTCAGGAAACTCAATGAATTTAACCCATGGa                                                                                                                                                                                                                     <  2:842558/63‑1 (MQ=255)
                                                                                                                                                                                                 aaGTGCAAACGTTTTAGCCACGACCCTTTAACCTCAGGAAACTCAATGAATTTAACCCATGGATAACTTTCTTTTAACTCTTTAGCAGAATGGACTAATGCGATAAAACTGACATTTTCTTTATTATTAGTTGCTGCCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGACTATGCTTTTTCCCATAATTTGGTCTCATGATTGTAtt                                       >  1:993770/1‑237 (MQ=255)
                                                                                                                                                                                                 aaGTGCAAACGTTTTAACCACGACCCTTTAACCTCAGGTAACTCAATGAATTTAACCCATGGATAACTTTCTTTTAACTCTTTAGCAGAATGGACTAATGCGATAAAACTGACATTTTCTTTATTATTAGTTGCTGCCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACAACTATGCTTTTTCCCATAATTTGGTCTCATGATTGTAtt                                       >  1:549708/1‑237 (MQ=255)
                                                                                                                                                                                                             tttAGCCACGACCCTTTAACCTCAGGAACCTCAATGAATTTAACCCATGGATAACTTTCTTTTAACTCTTTAGCAGAATGGACTAATGCGATAAAACTGACATTTTCTTTATTATTAGTTGCTGCCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGACTATGCTTTTTCCCATAATTTGGTCTCATGATTGTATTTCTTTATGAttttt                         <  1:26185/239‑1 (MQ=255)
                                                                                                                                                                                                                    accaCCCTTTAACCTAAGAAAACCCATAGAATTTAACCGATGTATAACTTTCTTTTAAGTCTTGAGCAAGATGGACTAAAGCAATAAAACTGACATTTTCTTTATTTTAGGTTGCTGCCAAAAATTTTTTAAAAATGGTAAATGGGCCGCCAGTTGTAAAATTGACCGCAGAAACAACGACTATGCTTTTTCCCATAATTTGGTCTCATGATTGTATTTCTTTATGATTTTTGCTGGTa                  <  1:845518/236‑1 (MQ=255)
                                                                                                                                                                                                                        cccTTTAACCTAAGGAAACTAAATGAATTTAACCCATGGATAACTTTCTTTTAACTCTTTAGCAGAATGGACTAATGCGATAAAACTGACATTTTCTTTATTATTAGTTGCTGCCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGACTATGCTTTTTCCCATAATTTGGTCTCATGATTGTATTTCTTTATGATTTTTGCTGGTACTcc              <  1:474501/239‑1 (MQ=255)
                                                                                                                                                                                                                             tAACCTCAGGAAACTCAATGAATTTAACCCATGGATAACTTTCTTTTAACTCTTTAGCAGAATGGACTAATGCGATAAAACTGACATTTTCTTTATTATTAGTTGCTGCCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGACTATGCTTTTTCCCATAATTTGGTCTCATGATTGTATTTCTTTTTGATTTTTGCTGGTACTCCCGCaa         >  1:1013199/1‑239 (MQ=255)
                                                                                                                                                                                                                                   cAGGAAACTCAATGAATTTAACCCATGGATAACTTTCTTTTAACTCTTTAGCAGAATGGACTAATGCGATAAAACTGACATTTTCTTTATTATTAGTTGCTGCCAAAAATTTTTTCAAAATGGTAAATGGACCGcc                                                                                                          >  1:353557/1‑136 (MQ=255)
                                                                                                                                                                                                                                   cAGGAAACTCAATGAATTTAACCCATGGATAACTTTCTTTTAACTCTTTAGCAGAATGGACTAATGCGATAAAACTGACATTTTCTTTATTATTAGTTGCTGCCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGACTATGCTTTTTCCCATAATTTGGTCTCATGATTGTATTTCTTTATGATTTTTGCTGGTACTCCCGCAATGACAg   >  1:857193/1‑239 (MQ=255)
                                                                                                                                                                                                                                   aaGGAAACTCAATGAATTTAACCCATGGATAACTTTCTTTTAACTCTTTAGCAGAATGGACTAATGCGATAAAACTGACATTTTCTTTATTATTAGTTGCTGCCAAAAATTTTTTCAAAATGGTAAATGGACCGcc                                                                                                          <  2:353557/135‑1 (MQ=255)
                                                                                                                                                                                                                                    aGGAAACTCAATGAATTTAACCCATGGATAACTTTCTTTTAACTCTTTAGCAGAATGGACTAATGCGATAAAACTGACATTTTCTTTATTATTAGTTGCTGCCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGACTATGCTTTTTCCCATAATTTGGTCTCATGATTGTATTTCTTTATGATTTTTGCTGGTACTCCCGCAATGACAGt  >  1:494429/1‑239 (MQ=255)
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                       
GTATAGCATTTTAAATAAGAAAAAGCTAGGCTCCATAAGAATTTCACGGAATAAAATTCCTTTATAAAAAGGTGCAGGGTTATGACAATACACATATCTTTTTTTAGTGACGACATTGGCCGTAATATCATGCAGACAAATCCAATGCGTAGCATTCAGCTCTTTTGAAAGTTTTTTACAAACTACATATTCAAAGTGCAAACGTTTTAGCCACGACCCTTTAACCTCAGGAAACTCAATGAATTTAACCCATGGATAACTTTCTTTTAACTCTTTAGCAGAATGGACTAATGCGATAAAACTGACATTTTCTTTATTATTAGTTGCTGCCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGACTATGCTTTTTCCCATAATTTGGTCTCATGATTGTATTTCTTTATGATTTTTGCTGGTACTCCCGCAATGACAGT  >  NC_000913/2104102‑2104568

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