Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,222,274 T→C 7.0% intergenic (+26/+31) ycgH → / ← ymgD pseudogene; putative ATP‑binding component of a transport system/uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,222,2740TC 7.0% 165.6 / 2.9 59intergenic (+26/+31)ycgH/ymgDpseudogene; putative ATP‑binding component of a transport system/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base T (27/26);  new base C (2/2);  total (29/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.77e-01

TAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAACCTGATTCAGTTTTTCTACCACTTTGGCGCTCAGTTCCGGAGAT  >  NC_000913/1222126‑1222422
                                                                                                                                                    |                                                                                                                                                    
tAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGcgt                                                                                                                                                      >  7:224142/1‑149 (MQ=255)
 aaaCTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGcgtc                                                                                                                                                     >  5:129808/1‑149 (MQ=255)
  aaCTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGcgtc                                                                                                                                                     <  8:72522/148‑1 (MQ=255)
         cGTAAATTTCAACCTGAATCCAGCGGCCAGCGTCTGGTGAATGTGGGCGTTCAGCTGGGTTATAATGGCTACAATGACACCGCAGTGATGGTGGGACTGAAATATAAGTTCTGATGCCGACTTTAGCTAAAAAACCTCGTCGTATTCAt                                                                                                                                             <  6:77107/149‑1 (MQ=255)
                      cTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATACCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACa                                                                                                                                 >  8:23554/1‑148 (MQ=255)
                              aGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAata                                                                                                                        <  5:12208/149‑1 (MQ=255)
                                 ggCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAAtatta                                                                                                                     >  5:134308/1‑149 (MQ=255)
                                      gCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTAttttt                                                                                                                <  1:154126/149‑1 (MQ=255)
                                       cGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTAtttttt                                                                                                               >  8:7149/1‑149 (MQ=255)
                                        gTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTg                                                                                                              <  5:225905/149‑1 (MQ=255)
                                          cTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGcc                                                                                                            >  4:170655/1‑149 (MQ=255)
                                             gggAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGa                                                                                                          >  6:24744/1‑148 (MQ=255)
                                                    tgGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGc                                                                                                  <  5:24744/149‑1 (MQ=255)
                                                       gCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATg                                                                                               <  1:37785/149‑1 (MQ=255)
                                                        cGTGCAGCTTTGTGATAATTGCTACAATGTCTCCGCAGTGATGTTGGTCCTGAAATTTAAGTTCTGATGCCGCCGTTAGCTAAAAAACCGCGTCGTATTAATCGAGGCGGTACATATGAAATATTATTTTTTTTCTATAGGACTCAttg                                                                                              <  4:78523/149‑3 (MQ=255)
                                                           gCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGgtgt                                                                                           >  3:172164/1‑149 (MQ=255)
                                                              gccgtgTGATAATGGCTACAATGACACCGCAGTGATGGGGGGCCTGAAATATACGTTCCGCTCCCGCCGTTAGCTAAACAACCGAGGCGTTATCCTCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACTCATGCTGtcat                                                                                        <  8:177031/144‑1 (MQ=255)
                                                              gCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGTGCCTGAAATATATTTTCTGATCCCGTCGTTAGCTAATTAACCGCGTCGTATTTATCGACGCGGTACATATTAAATATTATTTTTTGCCGATAGCTCGCATGGGTtcat                                                                                        <  2:181894/149‑1 (MQ=255)
                                                                    tGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGtc                                                                                          >  6:151182/1‑141 (MQ=255)
                                                                    tGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGtc                                                                                          <  5:151182/141‑1 (MQ=255)
                                                                          tGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGa                                                                            >  4:214426/1‑149 (MQ=255)
                                                                              tACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCaa                                                                        >  2:81435/1‑149 (MQ=255)
                                                                              tACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCaa                                                                        >  4:223699/1‑149 (MQ=255)
                                                                              tACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGAGTCATAAATTGACGTGATCaa                                                                        >  4:111759/1‑149 (MQ=255)
                                                                               aCAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGa                                                                            >  1:150228/1‑144 (MQ=255)
                                                                               aCAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGa                                                                            <  2:150228/144‑1 (MQ=255)
                                                                                cAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAAc                                                                       <  8:97707/148‑1 (MQ=255)
                                                                                cAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATACCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACa                                                                      <  7:23554/149‑1 (MQ=255)
                                                                                 aaTGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAg                                                                     >  5:137551/1‑149 (MQ=255)
                                                                                   tGACACCGCTGTGATGGTGGGACTGAAATATAAGTTCTGATCGCGCAGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGcc                                                                                <  2:199236/136‑1 (MQ=255)
                                                                                   tGACACCGCAGTGATGGTGGGCCTGAAATATAACTTCTGAACCCGCAGTGAGCTAAAAAACCCCGTCGTATTCATCGACGCGGAACATATGAAATAAAATATTATGCCGATAGCACGCATTGTGTCATCAATTGcc                                                                                >  1:199236/1‑136 (MQ=255)
                                                                                       acCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTg                                                               >  2:232811/1‑149 (MQ=255)
                                                                                         cGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGc                                                                                                             <  6:4865/101‑1 (MQ=255)
                                                                                         cGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGc                                                                                                             >  5:4865/1‑101 (MQ=255)
                                                                                         cGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCg                                                             <  3:5259/149‑1 (MQ=255)
                                                                                          gCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCgg                                                            <  7:7149/149‑1 (MQ=255)
                                                                                           cAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACAGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCggg                                                           >  7:40568/1‑149 (MQ=255)
                                                                                             gTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTc                                                         <  1:81435/149‑1 (MQ=255)
                                                                                              tGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCt                                                        <  8:224142/149‑1 (MQ=255)
                                                                                              tGATGGTGGGCCTGAAATATAAGTTCAGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCt                                                        >  5:12292/1‑149 (MQ=255)
                                                                                                aTGGGGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGGGTCGGATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCtt                                                       <  2:206444/148‑1 (MQ=255)
                                                                                                   gtggGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAgc                                                   <  6:129808/149‑1 (MQ=255)
                                                                                                   gtggGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAg                                                    <  3:214426/148‑1 (MQ=255)
                                                                                                         cTGAAATATAAGTTCTGATCCCGACGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATTAAATATTATTTTTTGCAGATAGCACGCATGGTGACATCACTTTCCTTGATCAACAGCATTTGCTGGTCTTTACCTCAAAc                                             >  6:130462/1‑149 (MQ=255)
                                                                                                         cTGAAATATAAGTCCTGAGCCCGCAGTTAGCTAAAACCCCGCGGCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGAAAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAAc                                             <  6:43015/149‑1 (MQ=255)
                                                                                                                     ttCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAACCTGATTCAGttt                                 >  7:199028/1‑149 (MQ=255)
                                                                                                                          aTCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAACCTGATTCAGTTTTTCTa                            >  6:60618/1‑149 (MQ=255)
                                                                                                                           tCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCtt                                                       <  6:224767/121‑1 (MQ=255)
                                                                                                                           tCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCtt                                                       >  5:224767/1‑121 (MQ=255)
                                                                                                                            cccGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAACCTGATTCAGTTTTTCTAc                           >  7:38252/1‑148 (MQ=255)
                                                                                                                             ccgccgTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAACCTGATTCAGTTTTTCTAcc                          >  1:162281/1‑148 (MQ=255)
                                                                                                                                   ttAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAACCTGATTCAGTTTTTCTACCACTTTgg                   <  3:170655/149‑1 (MQ=255)
                                                                                                                                   ttAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAACCTGATTCAGTTTTTCTACCACTTTgg                   >  8:42197/1‑149 (MQ=255)
                                                                                                                                   ttAGCTAAAAAACCGCGCCGTATTCATCGTCGCGGTAGATATGAAATATTATTTTTTGCCGATTGCTCGCATGGTGTCATCAATTGCCGTGATCAACAGCTTTTGCGGGTCTTTAGCGCTAACCTGATTCAGTTTTTCTACCAc                        <  1:66195/144‑1 (MQ=255)
                                                                                                                                   ttAGCTAAAAAACCGCGCCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAACCTGATTCAGTTTTTCTACCAc                        <  1:69545/144‑1 (MQ=255)
                                                                                                                                   ttAGCTAAAAAACCGCGCCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAACCTGATTCAGTTTTTCTACCAc                        >  2:69545/1‑144 (MQ=255)
                                                                                                                                   ttAGCTAAAAAACCGCGCCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAACCTGATTAAGTTTTTCTACCAc                        >  2:66195/1‑144 (MQ=255)
                                                                                                                                    tAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGGACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCTCAAACCTGATTCAGTTTTTCTACCACTTTg                    <  6:12292/147‑1 (MQ=255)
                                                                                                                                      gCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAACCTGATTCAGTTTTTCTACCACTTTGgcgc                >  3:106394/1‑149 (MQ=255)
                                                                                                                                        tAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATg                                                                                               >  6:150581/1‑68 (MQ=255)
                                                                                                                                               ccGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAACCTGATTCAGTTTTTCTACCACTTTGGCGCTCAGTTCCg       <  6:137551/149‑1 (MQ=255)
                                                                                                                                                cgcgTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAACCTGATTCAGTTTTTCTACCACTTTGGCGCTCAGTTCCgg      <  1:232811/149‑1 (MQ=255)
                                                                                                                                                    tcgtATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAAGAGCATTTGCGGGTCTTTAGCGCAAACCTGATTCAGTTTTTCTACCACTTTGGCGCTCAGTTCCGGAGAt  <  4:172164/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAACCTGATTCAGTTTTTCTACCACTTTGGCGCTCAGTTCCGGAGAT  >  NC_000913/1222126‑1222422

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