Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,998,523 T→A 6.3% N208I (AAC→ATC)  ycaC‑IV ← putative hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,998,5230TA 6.3% 184.3 / 2.1 63N208I (AAC→ATC) ycaC‑IVputative hydrolase
Reads supporting (aligned to +/- strand):  ref base T (31/28);  new base A (2/2);  total (33/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 = 9.00e-01

TTGTCATGGAGTGCGGTGTTCATCGGCGTTCCTTGGGATGAGCTTGCCCGCATCCACTGTGCAACGGCGGACGGGCGCTGGATTGGATAGGTGTGCTGCCCTGGGCGGGCACCGATTGGCGGTGCCCGCGGCGATTGGGCTTACTTGTTGGCGTTGAAGGCGTTGTAGCTGGTCATCAGGTTGCGGTAGTCCGGGATGTGGTTGGAGCACAGCGCGG‑CCAGGCCTTCGATGTCGTTGCGCCAGTCGCGGTGCAGCTCACAGGCCACGCCGAACCAGGTCATCAGTTGCGCA  >  NC_002947/5998376‑5998666
                                                                                                                                                   |                                                                                                                                                
ttGTCATGGAGTGCGGTGTTCATCGGCGTTCCTTGGGATGAGCTTGCCCGCATCCACTTTGCAACGGCGGACGGGCGCTGGATTGGATAGGTGTGCTGCCCTGGGCGGGCACCGATTGGCGGTGCCCGCGGCGATTGGGCTTACttgtt                                                                                                                                                 <  8:7059/149‑1 (MQ=255)
 tGTCATGGAGTGCGGTGTTCATCGTGGTTCCTTGGGATGAGCTTGCCCGCAGCCACTGTGCAACGGCGGCCGGGCGCTGTATTGGATAGTTGTGCTGCCCTGGGCGGGCACCGATTGGGGTTGCCCGCGGCTATTGGTCTTACttgttg                                                                                                                                                <  2:70037/149‑1 (MQ=255)
    cATGGAGTGCGGTGTTCATCGGCGTTCCTTGGGATGAGCTTGCCCGCATCCACTGTGCAACGGCGGACGGGCGCTGGATTGGATAGGTGTGCTGCCCTGGGCGGGCACCGATTGGCGGTGCCCGCGGCGATTGGGCTTACTTGTTGGCg                                                                                                                                             >  1:376534/1‑149 (MQ=255)
     aTGGAGTGCGGTGTTCATCGGCGTTCCTTGGGATGAGCTTGCCCGCATCCACTGTGCAACGGCGGACGGGCGCTGGATTGTATAGGTGTCCTGCCCTGGGCGGGCACCTATTGGCGGTGCCCGCGGCGATTGGTCTTACTTGATGGctt                                                                                                                                            >  2:311891/1‑147 (MQ=255)
     aTGGAGTGCGGTGTTCATCGGCGTTCCTTGGGATGAGCTTGCCCGCATCCACTGTGCAACGGCGGACGGGCGCTGGATTGGATAGGTGTGCTGCCCTGGGCGGGCACCTATTGGAGGTGCCCGCGGCGATTGGGCTTACTTGATGGCGt                                                                                                                                            >  2:316456/1‑149 (MQ=255)
             cgttgtTCTTCGGCGTTCTTTGGGATGAGCTTGCCCGCATCCACTGTGCAACTGCGTACGGGCGGTGGATTGGTTAGGTGTGCTGCCCTGGGCGGGCACCGATTGGCGGTGCCCGCGGCGATTGGGCTTACTTGATGGCGTTGAAGGCg                                                                                                                                    <  1:311891/146‑1 (MQ=255)
             cGGTGGTCATCGGCGTTCCTTGGGATGAGCTTGCCCGCATCCACTGTGCAACGGCGGACGGGCGCTGGATTGGATAGGTGTGCTGCCCTGGGCGGGCACCGATTGGCGGTGCCCGCGGCGATTGGGCTTACTTGTTGGCGTTGAAgg                                                                                                                                      <  6:325718/147‑1 (MQ=255)
             cGGTGGTCATCGGCGTGCCTTGGGATGAGCGTGCCCGCATCCACTGGGCAACGGCGGACGGGCGCTGGATTGGATAGGTGCGCTGCCCTGGGCGGGCACCGATTGGCGGGGCCCGCGGCGATTGGGCTTACTTGATGGCGTTGAAGGCg                                                                                                                                    <  1:316456/149‑1 (MQ=255)
              ggTGTTCATCGGCGTTCCTTGGGATGAGCTTGCCCGCATCCACTGTGCAACGGCGGACGGGCGCTGGATTGGATAGGTGTGCTGCCCTGGGCGGGCACCGATTGGCGGTGCCCGCGGCGATTGGGCTTACTTGTTGGCGTTGAAGGCGt                                                                                                                                   <  8:295347/149‑1 (MQ=255)
               gtgtTCATCGGCGTTCCTTGGGATGAGCTTGCCCGCATCCACTGTGCAACGGCGGACGGGCGCTGGATTGGATAGGTGTGCTGCCCTGGGCGGGCACCGATTGGCGGTGCCCGCGGCGATTGGGCTTACTTGTTGGCGTTGAAGGCGtt                                                                                                                                  <  5:278607/149‑1 (MQ=255)
                     aTCGGCGTTCCTTGGGATGAGCTTGCCCGCATCCACTGTGCAACGGCGGACGGGCGCTGGATTGGATAGGTGTGCTGCCCTGGGCGGGCACCGATTGGCGGTGCCCGCGGCGATTGGGCTTACTTGTTGGCGTTGAAGGCGTTGTAGCt                                                                                                                            >  3:114311/1‑149 (MQ=255)
                      tCGGGGTTCATTGTGTTTAGCTTGCCCGCATCCACTGTGAAACGGCGTCCGGGCGCTGGATTGGATAGGTGGGCTGCCCTGGGCGGGCACCGATTGGCGGTGCCCGCGGCGATTGGGCTTACTTGTTGGCGTTGAAGGCGTTGTAGCTg                                                                                                                           <  8:339374/149‑1 (MQ=255)
                      tCGGCGTTCCTTGGGATGAGCTTGCCCGCATCCACTGTGCAACGGCGGACGGGCGCTGGATTGGATAGGTGTGCTGCCCTGGGCGGGCACCGATTGGCGGTGCCCGCGGCGATTGGGCTTACTTGTTGGCGTTGAAGGCGTTGTAGCTg                                                                                                                           <  7:183665/149‑1 (MQ=255)
                              ccTTGGGATGAGCTTGCCCGCATCCACTGTGCAACGGCGGACGGGCGCTGGATTGGATAGGTGTGCTGCCCTGGGCGGGCACCGATTGGCGGTGCCCGCGGCGATTGGGCTTACTTGTTGGCGTTGAAGGCGTTGTAGCTGGtcatca                                                                                                                    <  5:89263/148‑1 (MQ=255)
                                   ggATGAGCTTGCCCGCATCCACTGTGCAACGGCGGACGGGCGCTGGATTGGATAGGTGTGCTGCCCTGGGCGGGCACCGATTGGCGGTGCCCGCGGCGATTGGGCTTACTTGTTGGCGTTGAAGGCGTTGTAGCTGGTCATCAGGTTGc                                                                                                              >  2:17235/1‑149 (MQ=255)
                                     aTGAGCTTGCCCGCATCCACTGTGCAACGGCGGACGGGCGCTGGATTGGATAGGTGTGCTGCCCTGGGCGGGCACCGATTGGCGGTGCCCGCGGCGATTGGGCTTACTTGTTGGCGTTGAAGGCGTTGTAGCTGGTCATCAGGTTGCgg                                                                                                            <  4:114311/149‑1 (MQ=255)
                                           ttGCCCGCATCCACTGTGCAACGGCGGACGGGCGCTGGATTGGATAGGTGTGCTGCCCTGGGCGGGCACCGATTGGCGGTGCCCGCGGCGATTGGGCTTACTTGTTGGCGTTGAAGGCGTTGTAGCTGATCATCAGGTTGCGGTAGTcc                                                                                                      >  1:128562/1‑149 (MQ=255)
                                               ccGCATCCACTGTGCAACGGCGGACGGGCGCTGTATTGTATAGGTGTGCTGCCCTGGGCGGGCACCGATTGGCGGTGCCCGCGGCGATTGGGCTTACTTGTTGGCGTTGAAGGCGTTGTAGCTGGTCATCAGGTTGCGGTAGTCCGGGa                                                                                                  <  1:169317/149‑1 (MQ=255)
                                               ccGCATCCACTGTGCAACGGCGGACGGGCGCTGGATTGGATAGGTGTGCTGCCCTGGGCGGGCACCGATTGGCGGTGCCCGCGGCGATTGGGCTTACTTGTTGGCGTTGAAGGCGTTGTAGCTGGTCATCAGGTTGCGGTAGTCCGGGa                                                                                                  <  1:217788/149‑1 (MQ=255)
                                               ccGCATCCACTGTGCAACGGCGGACGGGCGCTGGATTGGATAGGTGTGCTGCCCTGGGCGGGCACCGATTGGCGGTGCCCGCGGCGATTGGGCTTACTTGTTGGCGTTGAAGGCGTTGTAGCTGGTCATCAGGTTGCGGTAGTCCGGGa                                                                                                  >  1:339513/1‑149 (MQ=255)
                                                     ccACTGTGCAACGGCGGACGGGCGCTGGATTGGATAGGTGTGCTGCCCTGGGCGGGCACCGATTGGCGGTGCCCGCGGCGATTGGGCTTACTTGTTGGCGTTGAAGGCGTTGTAGCTGGTCATCAGGTTGCGGTAGTCCGGGATGTGGt                                                                                            >  5:389450/1‑149 (MQ=255)
                                                       aCTGTGCAACGGCGGACGGGCGCTGGATTGGATAGGTGTGCTGCCCTGGGCGGGCACCGATTGGCGGTGCCCGCGGCGATTGGGCTTACttgt                                                                                                                                                  <  5:159640/93‑1 (MQ=255)
                                                       aCTGTGCAACGGCGGACGGGCGCTGGATTGGATAGGTGTGCTGCCCTGGGCGGGCACCGATTGGCGGTGCCCGCGGCGATTGGGCTTACttgt                                                                                                                                                  >  6:159640/1‑93 (MQ=255)
                                                          gttcAAGGGGGGACGGGCGCTGGATTGTTTTGGTGTGCTTCCCTGGGCGGGCACCGATTGGCGGTGCCCGCGGCGATTGGGCTTACTTGTTGGCGTTGAAGGCGTTGTAGCTGGTCATCAGGTTGCGGTAGTCCGGGATGTGGTTGGa                                                                                        <  3:371068/145‑1 (MQ=38)
                                                             cAACGGCGGACGGGCGCTGGATTGGATAGGTGTGCTGCCCTGGGCGGGCACCGATTGGCGGTGCCCGCGGCGATTGGGCTTACTTGTTGGCGTTGAAGGCGTTGTAGCTGGTCATCAGGTTGCGGTAGTCCGGGATGTGGTTGGAGca                                                                                     >  1:101247/1‑148 (MQ=255)
                                                              aaCGGCGGACGGGCGCTGGTTTGGTTTGGTGTGGTGCCCTGGGGGGGCACGGTTTGGCGGTGCCCGCGGCGATTGGGCTTACTTGTTGGCGTTGAAGGCGTTGTAGCTGGTCATCAGGTTGCGGTAGTCCGGGATGTGGTTGGAGcaca                                                                                   <  2:376534/149‑1 (MQ=38)
                                                              aaCGGCGGACGGGCGCTGGATTGGATAGGTGTGCTGCCCTGGGCGGGCACCGATTGGCGGTGCCCGCGGCGATTGGGCTTACTTGTTGGCGTTGAAGGCGTTGTAGCTGGTCATCAGGTTGCGGTAGTCCGGGATGTGGTTGGAGca                                                                                     >  4:72771/1‑147 (MQ=255)
                                                               aCGGCGGACGGGCGCTGGATTGGATAGGTGTGCTGCCCTGGGCGGGCACCGATTGGCGGTGCCCGCGGCGATTGGGCTTACTTGTTGGCGTTGAAGGCGTTGTAGCTGGTCATCAGGTTGCGGTAGTCCGGGATGTGGTTGGAGCACAg                                                                                  <  7:197692/149‑1 (MQ=255)
                                                                        gggCGCTGGATTGGATAGGTGTGCTGCCCTGGGCGGGCACCGATTGGCGGTGCCCGCGGCGATTGGGCTTACTTGTTGGCGTTGAAGGCGTTGTAGCTGATCATCAGGTTGCGGTAGTCCGGGATGTGGTTTGAGCACAGCGCGG‑CCAg                                                                        <  2:128562/149‑1 (MQ=255)
                                                                         ggCGCTGGATTGGATAGGTGTGCTGCCCTGGGCGGGCACCGATTGGCGGTGCCCGCGGCGATTGGGCTTACTTGTTGGCGTTGAAGGCGTTGTAGCTGGTCATCAGGTTGCGGTAGTCCGGGATGTGGTTGGAGCACAGCGCGG‑CCAgg                                                                       <  1:17235/149‑1 (MQ=255)
                                                                          gcgcTGGATTGGATAGGTGTGCTGCCCTGGGCGGGCACCGATTGGCGGTGCCCGCGGCGATTGGGCTTACTTGTTGGCGTTGAAGGCGTTGTAGCTGGTCATCAGGTTGCGGTAGTCCGGGATGTGGTTg                                                                                          >  4:310052/1‑130 (MQ=255)
                                                                          gcgcTGGATTGGATAGGTGTGCTGCCCTGGGCGGGCACCGATTGGCGGTGCCCGCGGCGATTGGGCTTACTTGTTGGCGTTGAAGGCGTTGTAGCTGGTCATCAGGTTGCGGTAGTCCGGGATGTGGTTg                                                                                          <  3:310052/130‑1 (MQ=255)
                                                                                gattggatAGGTGTGCTGCCCTGGGCGGGCACCGATTGGCGGTGCCCGCGGCGATTGGGCTTACTTGTTGGCGTTGAAGGCGTTGTAGCTGGTCATCAGGTTGCg                                                                                                             <  8:68718/105‑1 (MQ=255)
                                                                                gattggatAGGTGTGCTGCCCTGGGCGGGCACCGATTGGCGGTGCCCGCGGCGATTGGGCTTACTTGTTGGCGTTGAAGGCGTTGTAGCTGGTCATCAGGTTGCg                                                                                                             >  7:68718/1‑105 (MQ=255)
                                                                                gattggatAGGTGTGCTGCCCTGGGCGGGCACCGATCGGCGGTGCCCGCGGCGATTGGGCTTACTTGTTGGCGTTGAAGGCGTTGTAGCTGGTCATCAGGTTGCGGTAGTCCGGGATGTGGTTGGAGCACAGCGCGG‑CCAGGCCTTCGa                                                                >  2:238363/1‑149 (MQ=255)
                                                                                 attggatAGGTGTGCTGCCCTGGGCGGGCACCGATTGGCGGTGCCCGCGGCGATTGGGCTTACTTGTTGGCGTTGAAGGCGTTGTAGCTGGTCATCAGGTTGCGGTAGTCCGGGATGTGGTTGGAGCACAGCGCGG‑CCAGGCCTTCGAt                                                               >  1:70874/1‑149 (MQ=255)
                                                                                    ggatAGGTGTGCTGCCCTGGGCGGGCACCGATTGGCGGTGCCCGCGGCGATTGGGCTTACTTGTTGGCGTTGAAGGCGTTGTAGCTGGTCATCAGGTTGCGGTAGTCCGGGATGTGGTTGGAGCACAGCGCGG‑CCAGGCCTTCGATGt                                                             <  2:19025/148‑1 (MQ=255)
                                                                                     gttAGGTGTGCTGCCCTGGGCGGGCACCGATTGGCGGTGCCCGCGGCGATTGGGCTTACTTGTTGGCGTTGAAGGCGTTGTAGCTGGTCATCAGGTTGCGGTAGTCCGGGATGTGGTTGGAGCACAGCGCGG‑CCAGGCCTTCGATGTCg                                                           >  2:78989/3‑149 (MQ=255)
                                                                                     gatAGGTGTGCTGCCCTTGGCGGGCACCGAATGGCGGTGCCAGCGGCGATTGGGCTTACTTGTTGTCGTTGAAGGCGTTGTAGCTGGTCATCAGCTTGCGGTAGTCCGGGATGTGGTTGGAGCACAGCGCGG‑CCAGGCCTTCGATGTCg                                                           >  5:75010/1‑149 (MQ=255)
                                                                                      atAGGTGTGCTGCCCTGGGCGGGCACCGATTGGCGGTGCCCGCGGCGATTGGGCTTACTTGTTGGCGTTGAAGGCGTTGTAGCTGGTCATCAGGTTGCGGTAGTCCGGGATGTGGTTGGAGCACAGCGCGG‑CCAGGCCTTCGATGTCGt                                                          >  2:349913/1‑149 (MQ=255)
                                                                                         ggggtgCTGCCCGTGGCGGGCACCGTTTGGCGGTGCCCGCGGCGATTGGGCTTACTTGTTGGCGTTGAAGGCGTTGTAGCTGGTCATCAGGTTGCGGTAGTCCGGGATGTGGTTGGAGCACAGCGCGG‑CCAGGCCTTCGATGTCGTTg                                                        <  6:389450/145‑1 (MQ=255)
                                                                                         ggTGTGCTGCCCTGGGCGGGCACCGATTGGCGGTGCCCGCGGCGATTGGGCTTACTTGTTGGCGTTGAAGGCGTTGTAGCTGGTCATCAGGTTGCGGTAGTCCGGGATGTGGTTGGAGCACAGCGCGG‑CCAGGCCTTCGATGTCGTTgc                                                       <  3:72771/149‑1 (MQ=255)
                                                                                          gtgtGCTGCCCTGGGCGGGCACCGATTGGCGGTGCCCGCGGCGATTGGGCTTACTTGTTGGCGTTGAAGGCGTTGTAGCTGGTCATCAGGTTGCGGTAGTCCGGGATGTGGTTGGAGCACAGCGCGG‑CCAGGCCTTCGATGTCGTTgcg                                                      >  6:331499/1‑149 (MQ=255)
                                                                                           tgtgCTGCCCTGGGCGGGCACCGATTGGCGGTGCCCGCGGCGATTGGGCTTACTTGTTGGCGTTGAAGGCGTTGTAGCTGGTCATCAGGTTGCGGTAGTCCGGGATGTGGTTGGAGCACAGCGCGG‑CCAGGCCTTCGATGTCGTTgcgc                                                     >  5:342943/1‑149 (MQ=255)
                                                                                                 gcCCTGGGCGGGCACCGATTGGCGGTGCCCGCGGCGATTGGGCTTACTTGTTGGCGTTGAAGGCGTTGTAGCTGGTCATCAGGTTGCGGTAGTCCGGGATGTGGTTGGAGCACAGCGCGG‑CCAGGCCTTCGATGTCGTTGCGCCAGTcg                                               >  5:375174/1‑149 (MQ=255)
                                                                                                   ccTGGGCGGGCACCGATTGGCGGTGCCCGCGGCGATTGGGCTTACTTGTTGGCGTTGAAGGCGTTGTAGCTGGTCATCAGGTTGCGGTAGTCCGGGATGTGGTTGGAGCACAGCGCGG‑CCAGGCCTTCGATGTCGTTGCGCCAGTcgcg                                             <  1:349913/149‑1 (MQ=255)
                                                                                                   ccTGGGCGGGCACCGATTGGCGGTGCCCGCGGCGATTGGGCTTACTTGTTGGCGTTGAAGGCGTTGTAGCTGGTCATCAGGTTGCGGTAGTCCGGGATGTGGTTGGAGCACAGCGCGG‑CCAGGCCTTCGATGTCGTTGCGCCAGTcgcg                                             >  1:161516/1‑149 (MQ=255)
                                                                                                       ggcgggcACCGATCGGCGGTGCCCGCGGCGATTGGGCTTACTTGTTGGCGTTGCAGGCGTTGTAGCTGGTCATGAGGTTGCGGTAGTCCGGGATGTGGTTGGAGCACAGGGCGGCCCAGGCCTTCGATGTCGTTGCGCCAGTCGCgg                                            <  1:238363/147‑1 (MQ=25)
                                                                                                             cACCGATTGGCGGTGCCCGCGGCGATTGGGCTTACTTGTTGGCGTTGAAGGCGTTGTAGCTGGTCATCAGGTTGCGGTAGTCCGGGATGTGGTTGGAGCACAGCGCGG‑CCAGGCCTTCGAt                                                               >  7:86178/1‑121 (MQ=255)
                                                                                                             cACCGATTGGCGGTGCCCGCGGCGATTGGGCTTACTTGTTGGCGTTGAAGGCGTTGTAGCTGGTCATCAGGTTGCGGTAGTCCGGGATGTGGTTGGAGCACAGCGCGG‑CCAGGCCTTCGAt                                                               <  8:86178/121‑1 (MQ=255)
                                                                                                              aCCGATTGGCGGTGCCCGCGGCGATTGGGCTTACTTGTTGGCGTTGAAGGCGTTGTAGCTGGTCATCAGGTTGCGGTAGTCCGGGATGTGGTTGGAGCACAGCGCGG‑CCAGGCCTTCGATGTCGTTGCGCCAGTCGCGGTGCAGCTca                                   >  2:360767/1‑148 (MQ=255)
                                                                                                              aCCGATTGGCGATGCCCGCGGCGATTGGGCTTACTTGTTGGCGTTGAAGGCGTTGTAGCTGGTCATCAGGTTGCGGTAGTCCGGGATGTGGTTGGAGCACAGCGCGG‑CCAGGCCTTCGATGTCGTTGCGCCAGTCGCGGTGCAGCTcac                                  >  2:345122/1‑149 (MQ=255)
                                                                                                                cGATTGGCGGTGCCCGCGGCGATTGGGCTTACTTGTTGGCGTTGAAGGCGTTGTAGCTGGTCATCAGGTTGCGGTAGTCCGGGATGTGGTTGGAGCACAGCGCGG‑CCAGGCCTTCGATGTCGTTGCGCCAGTCGCGGTGCAGCTCACAg                                <  2:101247/149‑1 (MQ=255)
                                                                                                                cGATGGGCGGTGCCCGCGGCGATTGGGCTTACTTGTTGGCGTTGAAGGCGTTGTAGCTGGTCATCAGGTTGCGGTAGTCCGGGATGTGGTTGGAGCACAGCGCGG‑CCAGGCCTTCGATGTCGTTGCGCCAGTCGCGGTGCAGCTCACAg                                <  3:1310/149‑1 (MQ=255)
                                                                                                                  aTTGGCGGTGCCCGCGGCGATTGGGCTTACTTGTTGGCGTTGAAGGCGTTGTAGCTGGTCATCAGGTTGCGGTAGTCCGGGATGTGGTTGGAGCACAGCGCGG‑CCAGGCCTTCGATGTCGTTGCGCCAGTCGCGGTg                                          >  4:42571/1‑137 (MQ=255)
                                                                                                                  aTTGGCGGTGCCCGCGGCGATTGGGCTTACTTGTTGGCGTTGAAGGCGTTGTAGCTGGTCATCAGGTTGCGGTAGTCCGGGATGTGGTTGGAGCACAGCGCGG‑CCAGGCCTTCGATGTCGTTGCGCCAGTCGCGGTg                                          <  3:42571/137‑1 (MQ=255)
                                                                                                                      gCGGTGCCCGCGGCGATTGGGCTTACTTGTTGGCGTTGAAGGCGTTGTAGCTGGTCATCAGGTTGCGGTAGTCCGGGATGTGGTTGGAGCACAGCGCGG‑CCAGGCCTTCGATGTCGTTGCGCCAGTCGCGGTGCAGCTCACAGGCCACg                          >  4:54639/1‑149 (MQ=255)
                                                                                                                      gCGGTGCCCGCGGCGATTGGGCTTACTTGTTGGCGTTGAAGGCGTTGTAGCTGGTCATCAGGTTGCGGTAGTCCGGGATGTGGTTGGAGCACAGCGCGG‑CCAGGCCTTCGATGTCGTTGCGCCAGTCGCGGTGCAGCTCACAGGCCACg                          >  8:146022/1‑149 (MQ=255)
                                                                                                                      gCGGTGCCCGCGGCGATTGGGCTTACTTGTTGGCGTTGAAGGCGTTGTAGCTGGTCATCAGGTTGCGGTAGTCCGGGATGTGGTTGGAGCACAGCGCGG‑CCAGGCCTTCGATGTCGTTGCGCCAGTCGCGGTGCAGCTCACAGGCCACg                          >  5:42345/1‑149 (MQ=255)
                                                                                                                       cGGTGCCCGCGGCGATTGGGCTTACTTGTTGGCGTTGAAGGCGTTGTAGCTGGTCATCAGGTTGCGGTAGTCCGGGATGTGGTTGGAGCACAGCGCGG‑CCAGGCCTTCGATGTCGTTGCGCCAGTCGCGGTGCAGCTCACAGGCCACGc                         <  3:54639/149‑1 (MQ=255)
                                                                                                                       cGGTGCCCGCGGCGATTGGGCTTACTTGTTGGCGTTGAAGGCGTTGTAGCTGGTCATCAGGTTGCGGTAGTCCGGGATGTGGTTGGAGCACAGCGCGG‑CCAGGCCTTCGATGTCGTTGCGCCAGTCGCGGTGCAGCTCACAGGCCACGc                         <  3:357564/149‑1 (MQ=255)
                                                                                                                       cGGTGCCCGCGGCGATTGGGCTTACTTGTTGGCGTTGAAGGCGTTGTAGCTGGTCATCAGGTTGCGGTAGTCCGGGATGTGGTTGGAGCACAGCGCGG‑CCAGGCCTTCGATGTCGTTGCGCCAGTCGCGGTGCAGCTCACAGGCCACGc                         <  6:342943/149‑1 (MQ=255)
                                                                                                                           gCCCGCGGCGATTGGGCTTACTTGTTGGCGTTGAAGGCGTTGTAGCTGGTCATCAGGTTGCGGTAGTCCGGGATGTGGTTGGAGCACAGCGCGG‑CCAGGCCTTCGATGTCGTTGCGCCAGTCGCGGTGCAGCTCACAGGCCACGCCGaa                     >  5:107932/1‑149 (MQ=255)
                                                                                                                                     aTTGGGCTTACTTGTTGGCGTTGAAGGCGTTGTAGCTGGTCATCAGGTTGCGGTAGTCCGGGATGTGGTTGGAGCACAGCGCGG‑CCAGGCCTTCGATGTCGTTGCGCCAGTCGCGGTGCAGCTCACAGGCCACGCCGAACCAggcttga           >  5:354982/1‑144 (MQ=255)
                                                                                                                                      ttGGGCTTACTTGTTGGCGTTGAAGGCGTTGTAGCTGGTCATCAGGTTGCGGTAGTCCGGGATGTGGTTGGAGCACAGCGCGG‑CCAGGCCTTCGATGTCGTTGCGCCAGTCGCGGTGCAGCTCACAGGCCACGCCGAACCAGGtcat            >  4:255754/1‑147 (MQ=255)
                                                                                                                                            ttACTTGTTGGCGTTGAACGCGTTGTAGCTGGTCATCAGGTTGCGGTAGTCCGGGATGTGGTTCGAGCACAGCGCGG‑CCAGGCCTTCGATGTCGTTGCGCCAGTCGCGGTGCAGCTCACAGGCCACGCCGACCCAGTTCATCAGTTtcg    >  2:123812/1‑149 (MQ=21)
                                                                                                                                              aCTTGTTGGCGTTGAAGGCGTTGTAGCTGGTCATCAGGTTGCGGTAGTCCGGGATGTGGTTGGAGCACAGCGCGG‑CCAGGCCTTCGATGTCGTTGCGCCAGTCGCGGTGCAGCTCACAGGCCACGCCGAACCAGGTCATCAGTTGCGCa  >  4:33655/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
TTGTCATGGAGTGCGGTGTTCATCGGCGTTCCTTGGGATGAGCTTGCCCGCATCCACTGTGCAACGGCGGACGGGCGCTGGATTGGATAGGTGTGCTGCCCTGGGCGGGCACCGATTGGCGGTGCCCGCGGCGATTGGGCTTACTTGTTGGCGTTGAAGGCGTTGTAGCTGGTCATCAGGTTGCGGTAGTCCGGGATGTGGTTGGAGCACAGCGCGG‑CCAGGCCTTCGATGTCGTTGCGCCAGTCGCGGTGCAGCTCACAGGCCACGCCGAACCAGGTCATCAGTTGCGCA  >  NC_002947/5998376‑5998666

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: ‑