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

TGGCACCGGCTCTGGTTGCTGCATTTGGTGGTAAAGAAAACATTACTAACCTCGACGCATGTATTACCCGTCTGCGCGTCAGCGTTGCTGATGTGTCTAAAGTGGATCAGGCCGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAAT‑CCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCTGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCACGGTGTCGGTATACACGCGATTTTGCCAGTTAACACCAC  >  NC_000913/1159070‑1159603
                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                           
tGGCACCGGCTCTGGTTGCTGCATTTGGTGGTAAAGAAAACATTACTAACCTCGACGCATGTATTACCCGTCTGCGCGTCAGCGTTGCTGATGTGTCTAAAGTGGATCAGGCCGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAAT‑CCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGcc                                                                                                                                                                                                                                                                     >  2:338657/1‑275 (MQ=255)
        gccctGGTTGCTGTATTTGGTGGTAAAGAAAACATTACTAACCTCGACGCATGTATTACCCGTCTGCGCGTCAGCGTTGCTGACGTGTCTCACGTGGACCAGGCCGGCCTGTATCAACCGGCCGCAGCGGGCGTATTGGTTGCTGGTTCTTGGTTTCATGGGATTTTCGGTAATAAATTCTGTTCACTGAAAACCGAGATGGATGAGTACAGTGTTACCCAATTCTCCCATAACA‑GTTGTGGTGACCAAGGCAGACAGACGGAAGCCTTTTTTAc                                                                                                                                                                                                                                                             >  4:295384/4‑275 (MQ=255)
         ctctGGTTGCTGCATTTGGTGGTAAAGAAAACATTACTAACCTCGACGCATGTATTACCCGTCTGCGCGTCAGCGTTGCTGATGTGTATAAAGTGGATCAGGCCGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAAT‑CCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTAc                                                                                                                                                                                                                                                             <  2:336955/274‑1 (MQ=255)
          tctGGTTGCTGCATTTGGTGGTAAAGAAAACATTACTAACCTCGACGCATGTATTACCCGTCTGCGCGTCAGCGTTGCTGATGTGTCTAAAGTGGATCAGGCCGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAAT‑CCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTAc                                                                                                                                                                                                                                                             <  2:417392/273‑1 (MQ=255)
                  cggcAGTTGGTGGTAAAGAAAACATGACTAACCTCGACGCATGTATTACCCGTCTGCGCGTCAGCGTTGCTGATGTGTCTAAAGGGGATGAGGCCGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGATGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAAT‑CCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTAtt                                                                                                                                                                                                                                                   <  1:154258/273‑1 (MQ=255)
                   tgcATTTGGTGGTAAAGAAAACATTACTAACCTGACCGAATGTATTACCCGGCTGCGCGTCAGCGTTGGTGATGTGTCTAAATGGGATCAGGCCGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGGTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAAT‑CCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTc                                                                                                                                                                                                                                                  <  4:25244/275‑1 (MQ=255)
                    gcATTTGGTGGTAAAGAAAACATTACTAACCGCGACGCATGTATTACCCGTCTGCGCGTCAGCGTTGCTGATGTGTCTAAAGTGGATCAGGCCGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAAT‑CCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCa                                                                                                                                                                                                                                                 <  1:114708/275‑1 (MQ=255)
                       gttGGGGGGAAATAAAACATTACTAACCTCGACGCATGTATTACCCGTCTGCGCGTCAGCGTTGCTGATGTGTCTAAAGTGGATCAGGCCGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAAT‑CCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTc                                                                                                                                                                                                                                                  <  4:338653/270‑1 (MQ=255)
                               tAAAGAAAACATTACTAACCTCGACGCATGTATTACCCGTCTGCGCGTCAGCGTTGCTGATGTGTCTAAAGTGGATCAGGCCGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCGGGTGTTCAGGCGCTTTTCGGTACTAAATCCGATACTCTGAAAACCGAGATGGATGAGTACATCCGTACACACTAAT‑CCGTAAGTCGTTGGTGAGACTAAGGCAGCCAACTGGCTGACATTCTGCCATGGTTTGTTCAGAATTGaggcg                                                                                                                                                                                                                                      >  4:381753/1‑271 (MQ=255)
                                      aaCATTACTAACCTCGACGCATGTATTACCCGTCTGCGCGTCAGCGTTGCTGATGTGTCTAAAGTGGATCAGGCCGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAAT‑CCGTAAGACGTTGGGGAGACTAAGGCAGCAAGATGGCTGCCTTTCTTTAAGGTGTTCTTCAGATTTGATCCGTGCCgg                                                                                                                                                                                                                                >  2:16904/1‑274 (MQ=255)
                                       cagttAATAACCTCGAGGAATGTATTACCCTGCTGCGCGTCAGCGTTGCTGATGTGTCTAAAGTGGATCAGGCCGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAAT‑CCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGt                                                                                                                                                                                                                               <  4:212115/271‑1 (MQ=255)
                                                aCCTCGACGCATGTATTACCCGTCTGCGCGTCAGCGTTGCTGATGTGTCCAAAGTGGATCAGGCCGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAACCCGAGATGGATGAGTACCTCCGTAACCACTAAT‑CCGTATGACGTTGGGGAGACTAAGTCAGCCAGTTGCCTGCCTTTTTTACAGGTGTTATGCAGACTCGCTACGTGCCGGTAATGCTGaa                                                                                                                                                                                                                      >  4:234518/1‑274 (MQ=255)
                                                     gACGCATGTATTACCCGTCTGCGCGTCAGCGTTGCTGATGTGTCTAAAGTGGATCAGGCCGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAAT‑CCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGc                                                                                                                                                                                                                                   >  2:212118/1‑256 (MQ=255)
                                                          aTGTATTACCCGTCTGCGCGTCAGCGTTGCTGATGTGTCTAAAGTGGATCAGGCCGGCATGAAGAAACTGGGCGCAGCTGTCGTAATGGTTGCAGGTTCTGGTGTTCCGGCGATTTTCGGTACTAAATCCGATACACTGAAAACCGAGATGGATGAGTACATAAGTAACCACTCAT‑CCGTAAGACGTTGGGGTGACTAACGCAGCCAAATGGTTGCCTTTTTTACAGGTCTACTTCAGACTTGATACGGTCCGTGAATGGTGCACATACACGgtg                                                                                                                                                                                                           >  4:58231/1‑275 (MQ=255)
                                                                tACCCGTCTGCGCGTCAGCGTTGCTGATGTGTCTAAAGTGGATCAGGCCGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAAGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAAT‑CCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGgtg                                                                                                                                                                                                                                                        <  2:372135/224‑1 (MQ=255)
                                                                 aCCCGTCTGCGCGTCAGCGTTGCTGATGTGTCTAAAGTGGATCAGGCCGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAAT‑CCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTa                                                                                                                                                                                                    <  2:381758/275‑1 (MQ=255)
                                                                      tCTGCGCGTCAGCGTTGCTGATGTGTCTAAAGTGGATCAGGCCGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAAT‑CCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGa                                                                                                                                                                                               >  2:491752/1‑275 (MQ=255)
                                                                         gcgcGTCAGCGGTGCGGATGTGGTTAAAGGGGATCAGGCCGGGCTGAAGAAACGGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAAT‑CCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAg                                                                                                                                                                                            <  2:58232/275‑1 (MQ=255)
                                                                           tcgcaaGCGGTGCTGATGTGTCTAAAGTGGATCAGGCCGTCCTGAATAAACTGGGCGCAGCGGGTGTAGTGGTTGCTGGTTCTAGTGTTCAGGCGATTTTCGGTACTAAAGCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAAT‑CCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATACAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTa                                                                                                                                                                                                    <  1:140317/260‑1 (MQ=255)
                                                                                aGCTTTGCTGATTGGTCTAAATGGGATCAGGCAGGCCGGAAGAAACTGGGCGTGACGGGCTGAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGTGATGGATGAGTACATCCGTAACCACTAAT‑CCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGACAGAACCTTc                                                                                                                                                                                     <  4:273874/275‑1 (MQ=255)
                                                                                aGCGTTGCTGATGTGTCTAAAGTGGATCAGGCCGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAAT‑CCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTc                                                                                                                                                                                     <  2:295388/275‑1 (MQ=255)
                                                                                               tCTAAAGTGGATCAGGCCGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAAT‑CCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTTGGCCGTAGACGATAGAACCTTCCACGTTGGTATCGtc                                                                                                                                                                      >  2:365338/1‑274 (MQ=255)
                                                                                                         aTCAGGCCGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAAT‑CCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGACCCTTCCACGTTGGTATCGTAGGTTTTGTCg                                                                                                                                                            >  2:463556/1‑275 (MQ=255)
                                                                                                           atgcgcgggcagtagaaACTGGGCGCCGCGAGAGTATTGGTTGTGGTTTATGGTGTTCATTCGATTTTATGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAATCACTAAT‑CCGTAAAACGTTGGGGCGACTAAGGCAGCCAGATGGCTGACGTTTTTACAGTTGGTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTACCACGTGGGTATCGTAGGTTGTATCGaa                                                                                                                                                          <  4:288954/262‑1 (MQ=255)
                                                                                                             ggccggccTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAAT‑CCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAAc                                                                                                                                                                                         >  2:291919/1‑242 (MQ=255)
                                                                                                             ggccggccTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAAT‑CCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAAc                                                                                                                                                         >  4:340248/1‑274 (MQ=255)
                                                                                                                          aaaCTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAAT‑CCGTAAGACGTTGGGGAGACTAAGGCAGCCAGAt                                                                                                                                                                                                                                                                            <  7:116214/146‑1 (MQ=255)
                                                                                                                          aaaCTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAAT‑CCGTAAGACGTTGGGGAGACTAAGGCAGCCAGAt                                                                                                                                                                                                                                                                            >  8:116214/1‑146 (MQ=255)
                                                                                                                              tGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAGTCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAAT‑CCGTATGACGTTGGGGAGACTAAGGCAGCCATCTGGCTGcc                                                                                                                                                                                                                                                                     <  5:467480/149‑1 (MQ=255)
                                                                                                                               gggCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAAT‑CCGTAAGACGTTGGGGAGACTAAGGCAGCCATCTGGCTGCCt                                                                                                                                                                                                                                                                    <  6:47435/149‑1 (MQ=255)
                                                                                                                                 gcgcAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAAT‑CCGTAAGACGTTGGGGAGACTAAGGCAGCCATCTGGCTGCCtta                                                                                                                                                                                                                                                                  <  6:313837/149‑2 (MQ=255)
                                                                                                                                                  gTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAAT‑CCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCttttt                                                                                                                                                                                                                                                                >  8:111280/1‑134 (MQ=255)
                                                                                                                                                  gTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAAT‑CCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCttttt                                                                                                                                                                                                                                                                <  7:111280/134‑1 (MQ=255)
                                                                                                                                                  gTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAAT‑CCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCa                                                                                                                                                                                                                                                 >  7:357379/1‑149 (MQ=255)
                                                                                                                                                   ttGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAAT‑CCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGt                                                                                                                   <  2:234522/274‑1 (MQ=255)
                                                                                                                                                    tGCTGGTTCTGGGGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATACGTAACCACTAAT‑CCGTAAGACGCTGGGGAGACTAAGGCATCCAGATGGCTGCCTNTTTTACAGGTGTTATTCAGAATTGACACGTGACGGTAATGCTTAAATTACGCGGTGTGCAGTAGACGATNGAACCTTCCACGCTGATATCGTAGGGTTTGTCGAACACGGTCTCGNCNCTCCCCCGGAACCAAACGttttttggca                                                                                                                 >  3:53972/1‑270 (MQ=255)
                                                                                                                                                      cTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAAT‑CCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCt                                                                                                                                                                                                                                                                    <  8:404199/126‑1 (MQ=255)
                                                                                                                                                      cTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAAT‑CCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCt                                                                                                                                                                                                                                                                    >  7:404199/1‑126 (MQ=255)
                                                                                                                                                         gTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAAT‑CCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGa                                                                                                                                                                                                                                          >  8:135008/1‑149 (MQ=255)
                                                                                                                                                              ggTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAAT‑CCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCTGGTAGc                                                                                                        <  2:340252/274‑1 (MQ=255)
                                                                                                                                                                     aGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAAT‑CCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTa                                                                                                                                                                                                                              >  8:460181/1‑149 (MQ=255)
                                                                                                                                                                      ggCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAAT‑CCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTaa                                                                                                                                                                                                                             >  8:364555/1‑149 (MQ=255)
                                                                                                                                                                        cGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAAT‑CCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATg                                                                                                                                                                                                                           <  7:135008/149‑1 (MQ=255)
                                                                                                                                                                         gATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAAT‑CCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGc                                                                                                                                                                                                                          <  6:148067/149‑1 (MQ=255)
                                                                                                                                                                                            gATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAAT‑CCGTAAGACGTTGGGGCGTCTACGCCAGCCAGATGGCTGCCTTCTTTACAGGTGCTATTCAGAATTGATACGTGCCCGTAATGCTGTATCTACGCGGTGTCCCGTAGACGCTAGTACCTTCCACGTTGGTATCCTCTGTTTTGCCCTATCGGTTATTACCTTTCCCCTGTAACACGAAGTTTTTCGTCCCCTGGTAGCgg                                                                                                      >  4:234028/1‑246 (MQ=255)
                                                                                                                                                                                            gATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAAT‑CCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGCTACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATGGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGCTATTGACGTTCCCCTGTCACGAAGAGTTTTTCTGCACCTGGTAGCGGGTGAAAGGATCAACAAGAGCGTAGCTAcc                                                                         >  4:102095/1‑275 (MQ=255)
                                                                                                                                                                                            gATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAAT‑CCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGt                                                                                                                                                              >  3:120535/1‑190 (MQ=255)
                                                                                                                                                                                                aCCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAAT‑CCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGaa                                                                                                                                                          <  1:120535/190‑1 (MQ=255)
                                                                                                                                                                                                aCCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAAT‑CCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCTGGTAGCGGGTGa                                                                                                  <  2:234032/246‑1 (MQ=255)
                                                                                                                                                                                                 ccTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAAT‑CCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTc                                                                                                                                                                                     >  4:308232/1‑162 (MQ=255)
                                                                                                                                                                                                     aaaaCCGAGATGGATGAGTACATCCGTAACCACTAAT‑CCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACg                                                                                                                                                                                 <  2:308236/162‑1 (MQ=255)
                                                                                                                                                                                                      aaaCCGAGATGGATGAGTACATCCGTAACCACTAAT‑CCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGCAACGAGAAGTTTTTCGTCACCTGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGgcg                                                               <  2:102096/275‑1 (MQ=255)
                                                                                                                                                                                                          ctagaTGGATGAGTACAACCGTAACCACTAAT‑CCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTGTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCGGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCTGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGa                                                           <  1:100479/273‑1 (MQ=255)
                                                                                                                                                                                                                   tgAGTACATCCGTAACCACTAAT‑CCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGt                                                                                                                                                                                >  6:259579/1‑149 (MQ=255)
                                                                                                                                                                                                                          aTCCGTAACCACTAAT‑CCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATc                                                                                                                                                                         >  5:3804/1‑149 (MQ=255)
                                                                                                                                                                                                                            ccGTAAACAGTAAT‑CCGTACGGATATGTGAAAACAAAGGCAGGCAGGTGGCTTGCTTTTTTACAGGTGTGATTTAGAATTGATAAGTTCCGGTAAAGCTGAAATTAAGCGGGGTGGCGTAGGCGATCGAGACGTCCACGTTGGTATCGTAGGTTTTGTCTAACAGGTTATGGACGTTCCCCTGTAACGAGAAATTTTTCGTCACCTGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGC‑GAAGGTGCCATACGGGGTCa                                         <  4:16904/274‑1 (MQ=255)
                                                                                                                                                                                                                               tAACCACTAAT‑CCGTAAGACGTTGGGGAGACTAATGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTAtt                                                                                                                                                 >  2:334880/1‑168 (MQ=255)
                                                                                                                                                                                                                               tAACCACTAAT‑CCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTACAGAGAAGTTTTTCGTCACCTGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACATTGCTCGGCGCGGAAGGTGCCATCCGTCGTCCCg                                       >  2:6564/1‑274 (MQ=255)
                                                                                                                                                                                                                                aaCCACTAAT‑CCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGc                                                                                                                                                                                                        >  3:132390/1‑112 (MQ=255)
                                                                                                                                                                                                                                   cACTAAT‑CCGTAAGACGTTGGGGAGACTAATGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACg                                                                                                                                             <  4:334876/168‑1 (MQ=255)
                                                                                                                                                                                                                                    aCTAAT‑CCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTa                                                                                                                                                                                                    <  1:132390/112‑1 (MQ=255)
                                                                                                                                                                                                                                      tAAT‑CCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGc                                                                                                                                                                                                        >  2:412460/1‑106 (MQ=255)
                                                                                                                                                                                                                                      tAAT‑CCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCTGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCACGGTGTCGGt                               >  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)
                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                           
TGGCACCGGCTCTGGTTGCTGCATTTGGTGGTAAAGAAAACATTACTAACCTCGACGCATGTATTACCCGTCTGCGCGTCAGCGTTGCTGATGTGTCTAAAGTGGATCAGGCCGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAAT‑CCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCTGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCACGGTGTCGGTATACACGCGATTTTGCCAGTTAACACCAC  >  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: ‑