Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,159,336 A→C 7.1% intergenic (+34/+26) ptsG → / ← fhuE fused glucose‑specific PTS enzymes: IIB component/IIC component/ferric‑rhodotorulic acid outer membrane transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,159,3360AC 7.1% 236.2 / 7.0 77intergenic (+34/+26)ptsG/fhuEfused glucose‑specific PTS enzymes: IIB component/IIC component/ferric‑rhodotorulic acid outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base A (36/33);  new base C (2/4);  total (39/38)
Fisher's exact test for biased strand distribution p-value = 4.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.55e-01

TGGCACCGGCTCTGGTTGCTGCATTTGGTGGTAAAGAAAACATTACTAACCTCGACGCATGTATTACCCGTCTGCGCGTCAGCGTTGCTGATGTGTCTAAAGTGGATCAGGCCGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCTGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCACGGTGTCGGTATACACGCGATTTTGCCAGTTAACACCAC  >  NC_000913/1159070‑1159603
                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                           
tGGCACCGGCTCTGGTTGCTGCATTTGGTGGTAAAGAAAACATTACTAACCTCGACGCATGTATTACCCGTCTGCGCGTCAGCGTTGCTGATGTGTCTAAAGTGGATCAGGCCGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGcc                                                                                                                                                                                                                                                                     >  2:338657/1‑275 (MQ=255)
        gccctGGTTGCTGTATTTGGTGGTAAAGAAAACATTACTAACCTCGACGCATGTATTACCCGTCTGCGCGTCAGCGTTGCTGACGTGTCTCACGTGGACCAGGCCGGCCTGTATCAACCGGCCGCAGCGGGCGTATTGGTTGCTGGTTCTTGGTTTCATGGGATTTTCGGTAATAAATTCTGTTCACTGAAAACCGAGATGGATGAGTACAGTGTTACCCAATTCTCCCATAACAGTTGTGGTGACCAAGGCAGACAGACGGAAGCCTTTTTTAc                                                                                                                                                                                                                                                             >  4:295384/4‑275 (MQ=255)
         ctctGGTTGCTGCATTTGGTGGTAAAGAAAACATTACTAACCTCGACGCATGTATTACCCGTCTGCGCGTCAGCGTTGCTGATGTGTATAAAGTGGATCAGGCCGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTAc                                                                                                                                                                                                                                                             <  2:336955/274‑1 (MQ=255)
          tctGGTTGCTGCATTTGGTGGTAAAGAAAACATTACTAACCTCGACGCATGTATTACCCGTCTGCGCGTCAGCGTTGCTGATGTGTCTAAAGTGGATCAGGCCGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTAc                                                                                                                                                                                                                                                             <  2:417392/273‑1 (MQ=255)
                  cggcAGTTGGTGGTAAAGAAAACATGACTAACCTCGACGCATGTATTACCCGTCTGCGCGTCAGCGTTGCTGATGTGTCTAAAGGGGATGAGGCCGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGATGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTAtt                                                                                                                                                                                                                                                   <  1:154258/273‑1 (MQ=255)
                   tgcATTTGGTGGTAAAGAAAACATTACTAACCTGACCGAATGTATTACCCGGCTGCGCGTCAGCGTTGGTGATGTGTCTAAATGGGATCAGGCCGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGGTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTc                                                                                                                                                                                                                                                  <  4:25244/275‑1 (MQ=255)
                    gcATTTGGTGGTAAAGAAAACATTACTAACCGCGACGCATGTATTACCCGTCTGCGCGTCAGCGTTGCTGATGTGTCTAAAGTGGATCAGGCCGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCa                                                                                                                                                                                                                                                 <  1:114708/275‑1 (MQ=255)
                       gttGGGGGGAAATAAAACATTACTAACCTCGACGCATGTATTACCCGTCTGCGCGTCAGCGTTGCTGATGTGTCTAAAGTGGATCAGGCCGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTc                                                                                                                                                                                                                                                  <  4:338653/270‑1 (MQ=255)
                               tAAAGAAAACATTACTAACCTCGACGCATGTATTACCCGTCTGCGCGTCAGCGTTGCTGATGTGTCTAAAGTGGATCAGGCCGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCGGGTGTTCAGGCGCTTTTCGGTACTAAATCCGATACTCTGAAAACCGAGATGGATGAGTACATCCGTACACACTAATCCGTAAGTCGTTGGTGAGACTAAGGCAGCCAACTGGCTGACATTCTGCCATGGTTTGTTCAGAATTGaggcg                                                                                                                                                                                                                                      >  4:381753/1‑271 (MQ=255)
                                      aaCATTACTAACCTCGACGCATGTATTACCCGTCTGCGCGTCAGCGTTGCTGATGTGTCTAAAGTGGATCAGGCCGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCAAGATGGCTGCCTTTCTTTAAGGTGTTCTTCAGATTTGATCCGTGCCgg                                                                                                                                                                                                                                >  2:16904/1‑274 (MQ=255)
                                       cagttAATAACCTCGAGGAATGTATTACCCTGCTGCGCGTCAGCGTTGCTGATGTGTCTAAAGTGGATCAGGCCGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGt                                                                                                                                                                                                                               <  4:212115/271‑1 (MQ=255)
                                                aCCTCGACGCATGTATTACCCGTCTGCGCGTCAGCGTTGCTGATGTGTCCAAAGTGGATCAGGCCGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAACCCGAGATGGATGAGTACCTCCGTAACCACTAATCCGTATGACGTTGGGGAGACTAAGTCAGCCAGTTGCCTGCCTTTTTTACAGGTGTTATGCAGACTCGCTACGTGCCGGTAATGCTGaa                                                                                                                                                                                                                      >  4:234518/1‑274 (MQ=255)
                                                     gACGCATGTATTACCCGTCTGCGCGTCAGCGTTGCTGATGTGTCTAAAGTGGATCAGGCCGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGc                                                                                                                                                                                                                                   >  2:212118/1‑256 (MQ=255)
                                                          aTGTATTACCCGTCTGCGCGTCAGCGTTGCTGATGTGTCTAAAGTGGATCAGGCCGGCATGAAGAAACTGGGCGCAGCTGTCGTAATGGTTGCAGGTTCTGGTGTTCCGGCGATTTTCGGTACTAAATCCGATACACTGAAAACCGAGATGGATGAGTACATAAGTAACCACTCATCCGTAAGACGTTGGGGTGACTAACGCAGCCAAATGGTTGCCTTTTTTACAGGTCTACTTCAGACTTGATACGGTCCGTGAATGGTGCACATACACGgtg                                                                                                                                                                                                           >  4:58231/1‑275 (MQ=255)
                                                                tACCCGTCTGCGCGTCAGCGTTGCTGATGTGTCTAAAGTGGATCAGGCCGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAAGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGgtg                                                                                                                                                                                                                                                        <  2:372135/224‑1 (MQ=255)
                                                                 aCCCGTCTGCGCGTCAGCGTTGCTGATGTGTCTAAAGTGGATCAGGCCGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTa                                                                                                                                                                                                    <  2:381758/275‑1 (MQ=255)
                                                                      tCTGCGCGTCAGCGTTGCTGATGTGTCTAAAGTGGATCAGGCCGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGa                                                                                                                                                                                               >  2:491752/1‑275 (MQ=255)
                                                                         gcgcGTCAGCGGTGCGGATGTGGTTAAAGGGGATCAGGCCGGGCTGAAGAAACGGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAg                                                                                                                                                                                            <  2:58232/275‑1 (MQ=255)
                                                                           tcgcaaGCGGTGCTGATGTGTCTAAAGTGGATCAGGCCGTCCTGAATAAACTGGGCGCAGCGGGTGTAGTGGTTGCTGGTTCTAGTGTTCAGGCGATTTTCGGTACTAAAGCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATACAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTa                                                                                                                                                                                                    <  1:140317/260‑1 (MQ=255)
                                                                                aGCTTTGCTGATTGGTCTAAATGGGATCAGGCAGGCCGGAAGAAACTGGGCGTGACGGGCTGAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGTGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGACAGAACCTTc                                                                                                                                                                                     <  4:273874/275‑1 (MQ=255)
                                                                                aGCGTTGCTGATGTGTCTAAAGTGGATCAGGCCGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTc                                                                                                                                                                                     <  2:295388/275‑1 (MQ=255)
                                                                                               tCTAAAGTGGATCAGGCCGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTTGGCCGTAGACGATAGAACCTTCCACGTTGGTATCGtc                                                                                                                                                                      >  2:365338/1‑274 (MQ=255)
                                                                                                         aTCAGGCCGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGACCCTTCCACGTTGGTATCGTAGGTTTTGTCg                                                                                                                                                            >  2:463556/1‑275 (MQ=255)
                                                                                                           atgcgcgggcagtagaaACTGGGCGCCGCGAGAGTATTGGTTGTGGTTTATGGTGTTCATTCGATTTTATGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAATCACTAATCCGTAAAACGTTGGGGCGACTAAGGCAGCCAGATGGCTGACGTTTTTACAGTTGGTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTACCACGTGGGTATCGTAGGTTGTATCGaa                                                                                                                                                          <  4:288954/262‑1 (MQ=255)
                                                                                                             ggccggccTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAAc                                                                                                                                                                                         >  2:291919/1‑242 (MQ=255)
                                                                                                             ggccggccTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAAc                                                                                                                                                         >  4:340248/1‑274 (MQ=255)
                                                                                                                          aaaCTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGAt                                                                                                                                                                                                                                                                            <  7:116214/146‑1 (MQ=255)
                                                                                                                          aaaCTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGAt                                                                                                                                                                                                                                                                            >  8:116214/1‑146 (MQ=255)
                                                                                                                              tGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAGTCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTATGACGTTGGGGAGACTAAGGCAGCCATCTGGCTGcc                                                                                                                                                                                                                                                                     <  5:467480/149‑1 (MQ=255)
                                                                                                                               gggCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCATCTGGCTGCCt                                                                                                                                                                                                                                                                    <  6:47435/149‑1 (MQ=255)
                                                                                                                                 gcgcAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCATCTGGCTGCCtta                                                                                                                                                                                                                                                                  <  6:313837/149‑2 (MQ=255)
                                                                                                                                                  gTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCttttt                                                                                                                                                                                                                                                                >  8:111280/1‑134 (MQ=255)
                                                                                                                                                  gTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCttttt                                                                                                                                                                                                                                                                <  7:111280/134‑1 (MQ=255)
                                                                                                                                                  gTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCa                                                                                                                                                                                                                                                 >  7:357379/1‑149 (MQ=255)
                                                                                                                                                   ttGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGt                                                                                                                   <  2:234522/274‑1 (MQ=255)
                                                                                                                                                    tGCTGGTTCTGGGGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATACGTAACCACTAATCCGTAAGACGCTGGGGAGACTAAGGCATCCAGATGGCTGCCTNTTTTACAGGTGTTATTCAGAATTGACACGTGACGGTAATGCTTAAATTACGCGGTGTGCAGTAGACGATNGAACCTTCCACGCTGATATCGTAGGGTTTGTCGAACACGGTCTCGNCNCTCCCCCGGAACCAAACGttttttggca                                                                                                                 >  3:53972/1‑270 (MQ=255)
                                                                                                                                                      cTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCt                                                                                                                                                                                                                                                                    <  8:404199/126‑1 (MQ=255)
                                                                                                                                                      cTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCt                                                                                                                                                                                                                                                                    >  7:404199/1‑126 (MQ=255)
                                                                                                                                                         gTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGa                                                                                                                                                                                                                                          >  8:135008/1‑149 (MQ=255)
                                                                                                                                                              ggTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCTGGTAGc                                                                                                        <  2:340252/274‑1 (MQ=255)
                                                                                                                                                                     aGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTa                                                                                                                                                                                                                              >  8:460181/1‑149 (MQ=255)
                                                                                                                                                                      ggCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTaa                                                                                                                                                                                                                             >  8:364555/1‑149 (MQ=255)
                                                                                                                                                                        cGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATg                                                                                                                                                                                                                           <  7:135008/149‑1 (MQ=255)
                                                                                                                                                                         gATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGc                                                                                                                                                                                                                          <  6:148067/149‑1 (MQ=255)
                                                                                                                                                                                            gATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGCGTCTACGCCAGCCAGATGGCTGCCTTCTTTACAGGTGCTATTCAGAATTGATACGTGCCCGTAATGCTGTATCTACGCGGTGTCCCGTAGACGCTAGTACCTTCCACGTTGGTATCCTCTGTTTTGCCCTATCGGTTATTACCTTTCCCCTGTAACACGAAGTTTTTCGTCCCCTGGTAGCgg                                                                                                      >  4:234028/1‑246 (MQ=255)
                                                                                                                                                                                            gATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGCTACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATGGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGCTATTGACGTTCCCCTGTCACGAAGAGTTTTTCTGCACCTGGTAGCGGGTGAAAGGATCAACAAGAGCGTAGCTAcc                                                                         >  4:102095/1‑275 (MQ=255)
                                                                                                                                                                                            gATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGt                                                                                                                                                              >  3:120535/1‑190 (MQ=255)
                                                                                                                                                                                                aCCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGaa                                                                                                                                                          <  1:120535/190‑1 (MQ=255)
                                                                                                                                                                                                aCCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCTGGTAGCGGGTGa                                                                                                  <  2:234032/246‑1 (MQ=255)
                                                                                                                                                                                                 ccTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTc                                                                                                                                                                                     >  4:308232/1‑162 (MQ=255)
                                                                                                                                                                                                     aaaaCCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACg                                                                                                                                                                                 <  2:308236/162‑1 (MQ=255)
                                                                                                                                                                                                      aaaCCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGCAACGAGAAGTTTTTCGTCACCTGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGgcg                                                               <  2:102096/275‑1 (MQ=255)
                                                                                                                                                                                                          ctagaTGGATGAGTACAACCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTGTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCGGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCTGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGa                                                           <  1:100479/273‑1 (MQ=255)
                                                                                                                                                                                                                   tgAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGt                                                                                                                                                                                >  6:259579/1‑149 (MQ=255)
                                                                                                                                                                                                                          aTCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATc                                                                                                                                                                         >  5:3804/1‑149 (MQ=255)
                                                                                                                                                                                                                            ccGTAAACAGTAATCCGTACGGATATGTGAAAACAAAGGCAGGCAGGTGGCTTGCTTTTTTACAGGTGTGATTTAGAATTGATAAGTTCCGGTAAAGCTGAAATTAAGCGGGGTGGCGTAGGCGATCGAGACGTCCACGTTGGTATCGTAGGTTTTGTCTAACAGGTTATGGACGTTCCCCTGTAACGAGAAATTTTTCGTCACCTGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGAAGGTGCCATACGGGGTCa                                         <  4:16904/274‑1 (MQ=255)
                                                                                                                                                                                                                               tAACCACTAATCCGTAAGACGTTGGGGAGACTAATGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTAtt                                                                                                                                                 >  2:334880/1‑168 (MQ=255)
                                                                                                                                                                                                                               tAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTACAGAGAAGTTTTTCGTCACCTGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACATTGCTCGGCGCGGAAGGTGCCATCCGTCGTCCCg                                       >  2:6564/1‑274 (MQ=255)
                                                                                                                                                                                                                                aaCCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGc                                                                                                                                                                                                        >  3:132390/1‑112 (MQ=255)
                                                                                                                                                                                                                                   cACTAATCCGTAAGACGTTGGGGAGACTAATGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACg                                                                                                                                             <  4:334876/168‑1 (MQ=255)
                                                                                                                                                                                                                                    aCTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTa                                                                                                                                                                                                    <  1:132390/112‑1 (MQ=255)
                                                                                                                                                                                                                                      tAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGc                                                                                                                                                                                                        >  2:412460/1‑106 (MQ=255)
                                                                                                                                                                                                                                      tAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCTGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCACGGTGTCGGt                               >  1:140586/1‑275 (MQ=255)
                                                                                                                                                                                                                                           ccGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTa                                                                                                                                                                                                    <  4:412455/106‑1 (MQ=255)
                                                                                                                                                                                                                                             gTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAg                                                                                                                                                       <  7:460181/149‑1 (MQ=255)
                                                                                                                                                                                                                                              tAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGCCGTTCCCCTGTAACGAGAAGTTTTTCGTCACCTGGTAGCGGGTGAAGAGATCCACCAGCGCGTCGATACCTTGCTCGGCGCGGAAGGTGCCAGACGGCATCACGGTGCCGGTATACCcg                        >  4:43244/1‑275 (MQ=255)
                                                                                                                                                                                                                                               aaGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGt                                                                                                                                                     <  5:259579/149‑1 (MQ=255)
                                                                                                                                                                                                                                                             aCTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCTGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCACGGTGTCGGTATACACGCGATTTTGCCAGtt          >  2:353060/1‑274 (MQ=255)
                                                                                                                                                                                                                                                               tAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCTGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGCCGCGGAAGGTGCCCTACGGCGTCCCGGTGCCGGTATACACGCGATTTTGCCAGTTaa        >  2:71040/1‑274 (MQ=255)
                                                                                                                                                                                                                                                                aaGGCAGCCATCTGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACa                                                                                                                                                        >  6:35089/1‑129 (MQ=255)
                                                                                                                                                                                                                                                                aaGGCAGCCATCTGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACa                                                                                                                                                        <  5:35089/129‑1 (MQ=255)
                                                                                                                                                                                                                                                                     aGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACga                                                                                                                               >  5:232602/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                      gCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCTGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGAAGGTGc                                                     >  4:271848/1‑222 (MQ=255)
                                                                                                                                                                                                                                                                       ccAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGTTACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCNTTTAACGAGACGTTTTTCGTCAcc                                                                                                               >  3:130619/1‑163 (MQ=255)
                                                                                                                                                                                                                                                                       ccAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCTGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCACGGTGTCGGTATACACGCGATTTTGCCAGTTAAcaccac  >  2:165824/1‑272 (MQ=255)
                                                                                                                                                                                                                                                                        cAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGaa                                                                                                                            <  6:3804/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                          gATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCTGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATa                                                 <  2:271852/222‑1 (MQ=255)
                                                                                                                                                                                                                                                                           aTGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCTGGt                                                                                                           <  1:130619/163‑1 (MQ=255)
                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                           
TGGCACCGGCTCTGGTTGCTGCATTTGGTGGTAAAGAAAACATTACTAACCTCGACGCATGTATTACCCGTCTGCGCGTCAGCGTTGCTGATGTGTCTAAAGTGGATCAGGCCGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCTGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCACGGTGTCGGTATACACGCGATTTTGCCAGTTAACACCAC  >  NC_000913/1159070‑1159603

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