Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,132,855 A→G 5.8% L102L (CTT→CTC)  PP_3636 ← sulfonate ABC transporter substrate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,132,8550AG 5.8% 200.5 / 2.2 69L102L (CTT→CTC) PP_3636sulfonate ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (30/35);  new base G (2/2);  total (32/37)
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.40e-01

GGAATGGCCACTGTCTTGCCACCGAGCTGTTTCACGGCGGATATATCCGGCGCAACAGTCAGGCCCGAACCGCCCACGTGGTTCCAGGCCACCACCTTGGCTGGCACCTTGCTGCCATAGCGCGCCCATACCGTCATCGGCGAAAGCAGGTGGATGACATTGACCTGGCCGGAGATGAATGCCTCGATCACCTGCGCCCAGCTGCGCAGCAGCACCGGGCGCTCGGCCTTGATGCCTTCGGCCTCGAACAGGCCGTTGTTATGTGCCACCAGCAGCGGCGTGGCATCGGTGA  >  NC_002947/4132712‑4133003
                                                                                                                                               |                                                                                                                                                    
ggAATGGCCACTGTCTTGCCACCGAGCTGTTTCACGGCGGATATATCCGGCGCAACAGTCAGGCCCGAACCGCCCACGTGGTTCCAGGCCACCACCTTGGCTGGCACCTTGCTGCCATAGCGCGCCCATACCGTCATCGGCGAAAGCAg                                                                                                                                                 >  4:309498/1‑149 (MQ=255)
  aaTGGCCACTGTCTTGCCACCGAGCTGTTTCACGGCGGATATATCCGGCGCAACAGTCAGGCCCGAACCGCCCACGTGGTTCCAGGCCACCACCTTGGCTGGCACCTTGCTGCCATAGCGCGCCCATACCGTCATCGGCGAAAGCAGGt                                                                                                                                               <  8:114515/149‑1 (MQ=255)
   aTGGCCACTGTCTTGCCACCGTGCTGTTTCACGGCGGATATATCCGGCGCAACAGTCAGGCCCGAACCGCCCACGTGGTTCCAGTCCACCACCTTGGCTGGCACCTTGCTGCCATAGCGCGCCCATACCGTCATCGGCGAAAGCAGGTg                                                                                                                                              <  3:309498/149‑1 (MQ=255)
   aTGGCCACTGTCTTGCCACCGAGCTGTTTCACGGCGGATATATCCGGCGCAACAGTCAGGCCCGAACCGCCCACGTGGTTCCAGGCCACCACCTTGGCTGGCACCTTGCTGCCATAGCGCGCCCATACCGTCATCGGCGaaa                                                                                                                                                     >  6:283861/1‑142 (MQ=255)
   aTGGCCACTGTCTTGCCACCGAGCTGTTTCACGGCGGATATATCCGGCGCAACAGTCAGGCCCGAACCGCCCACGTGGTTCCAGGCCACCACCTTGGCTGGCACCTTGCTGCCATAGCGCGCCCATACCGTCATCGGCGaaa                                                                                                                                                     <  5:283861/142‑1 (MQ=255)
   aTGGCCACTGTCTTGCCACCGAGCTGTTTCACGGCGGATATATCCGGCGCAACAGTCAGGCCCGAACCGCCCACGTGGTTCCAGGCCACCACCTTGGCTGGCACCTTGCTGCCATAGCGCGCCCATACCGTCATCGGCGAAAGCAGGTg                                                                                                                                              <  2:268698/149‑1 (MQ=255)
    tGGCCACTGTCTTGCCACCGAGCTGTTTCACGGCGGATATATCCGGCGCAACAGTCAGGCCCGAACCGCCCACGTGGTTCCAGGCCACCACCTTGGCTGGCACCTTGCTGCCATAGCGCGCCCATACCGTCATCGGCGAAAGCAGGTgg                                                                                                                                             >  3:233672/1‑149 (MQ=255)
         aCTGTCTTGCCACCGAGCTGTTTCACGGCGGATATATCCGGCGCAACAGTCAGGCCCGAACCGCCCACGTGGTTCCAGGCCACCACCTTGGCTGGCACCTTGCTGCCATAGCGCGCCCATACCGTCATCGGCGAAAGCAGGTGGATGa                                                                                                                                         >  4:291157/1‑148 (MQ=255)
         aCTGTCTTGCCACCGAGCTGTTTCACGGCGGATATATCCGGCGCAACAGTCAGGCCCGAACCGCCCACGTGGTTCCAGGCCACCACCTTGGCTGGCACCTTGCTGCCATAGCGCGCCCATACCGTCATCGGCGAAAGCAGGTGGATGAc                                                                                                                                        <  5:70417/149‑1 (MQ=255)
           tGTCTTGCCACCGAGCTGTTTCACGGCGGATATATCCGGCGCAACAGTCAGGCCCGAACCGCCCACGTGGTTCCAGGCCACCACCTTGGCTGGCACCTTGCTGCCATAGCGCGCCCATACCGTCATCGGCGAAAGCAGGTGGATGACAt                                                                                                                                      >  6:140095/1‑149 (MQ=255)
                        aGCTGTTTCACGGCGGATATATCCGGCGCAACAGTCAGGCCCGAACCGCCCACGTGGTTCCAGGCCACCACCTTGGCTGGCACCTTGCTGCCATAGCGCGCCCATACCGTCATCGGCGAAAGCAGGTGGATGACATTGACCTGGCCGga                                                                                                                         <  5:138526/149‑1 (MQ=255)
                        aGCTGTTTCACGGCGGATATATCCGGCGCAACAGTCAGCCCCGAACCGCCCACGTGGTTCCAGGCCACCACCTTGGCTGGCACCTTGCTGCCATAGCGCGCCCATACCGTCATCGGCGAAAGCAGGTGGATGACATTGACCTGGCCGga                                                                                                                         <  4:29933/149‑1 (MQ=255)
                           tGTTTCACGGCGGATATATCCGGCGCAACAGTCAGGCCCGAACCGCCCACGTGGTTCCAGGCCACCACCTTGGCTGGCACCTTGCTGCCATAGCGCGCCCATACCGTCATCGGCGAAAGCAGGTGGATGACATTGACCTGGCCGGAGAt                                                                                                                      >  2:8975/1‑149 (MQ=255)
                                 aCGGCGGATATATCCGGCGCAACAGTCAGGCCCGAACCGCCCACGTGGTTCCAGGCCACCACCTTGGCTGGCACCTTGCTGCCATAGCGCGCCCATACCGTCATCGGCGAAAGCAGGTGGATGACATTGACCTGGCCGGAGATGAATGc                                                                                                                >  1:184008/1‑149 (MQ=255)
                                   ggcggATATATCTGGCGCAACAGTCAGGCCCGAACCGCCCACGTGGTTCCAGGCCACCACCTTGGCTGGCACCTTGCTGCCATAGCGCGCCCATACCGTCATCGGCGAAAGCAGGTGGATGACATTGACCTGGCCGGAGATGAATGCCt                                                                                                              <  7:165406/149‑1 (MQ=255)
                                   ggcggATATATCCGGCGCAACAGTCAGGCCCGAACCGCCCACGTGGTTCCAGGCCACCACCTTGGCTGGCACCTTGCTGCCATAGCGCGCCCATACCGTCATCGGCGAAAGCAGGTGGATGACATTGACCTGGCCGGAGATGAATGCCt                                                                                                              <  2:318492/149‑1 (MQ=255)
                                     cggATATATCCGGCGCAACAGTCAGGCCCGAACAGCCCACGTGGTTCCAGGCCACCACCTTGGCTGGCACCTTGCTGCCATAGCGCGCCCATACCGTCATCGGCGAAAGCAGGTGGATTACATTGACCTGGCCGGAGATGAATGCCTc                                                                                                             >  1:262937/1‑148 (MQ=255)
                                      ggATATATCCGGCGCAACAGTCAGTCCCGAACCGCCCACGTGGTTCCAGGCCACCACCTTGGCTGGCACCTTGCTGCCATAGCGCGCCCATACCGTCATCGGCGAAAGCAGGTGGATGACATTGACCTGGCCGGAGATGAATGCCTCGa                                                                                                           <  4:239161/149‑1 (MQ=255)
                                      ggATATATCCGGCGCAACAGTCAGGCCCGAACCGCCCACGTGGTTCCAGGCCACCACCTTGGCTGGCACCTTGCTGCCATAGCGCGCCCATACCGTCATCGGCGAAAGCAGGTGGATGACATTGACCTGGCCGGAGATGAATGCCTCGa                                                                                                           <  3:121017/149‑1 (MQ=255)
                                      ggATATATCCGGCGCAACAGTCAGGCCCGAACCGCCCACGTGGTTCCAGGCCACCACCTTGGCTGGCACCTTGCTGCCATAGCGCGCCCATACCGTCATCGGCGAAAGCAGGTGGATGACATTGACCTGGCCGGAGATGAATGCCTCGa                                                                                                           >  5:95434/1‑149 (MQ=255)
                                      ggATATATCCGGCGCAACAGTCAGGCCCGAACCGCCCACGTGGTTCCAGGCCACCACCTTGGCTGGCACCTTGCTGCCATAGCGCGCCCATACCGTCATCGGCGAAAGCAGGTGGATGACATTGACCTGGCCGGAGATGAATGCCTCGa                                                                                                           <  5:140095/149‑1 (MQ=255)
                                            atCCGGCGCAACAGTCAGGCCCGAACCGCCCACGTGGTTCCAGGCCACCACCTTGGCTGGCACCTTGCTGCCATAGCGCGCCCATACCGTCATCGGCGAAAGCAGGTGGATGACATTGACCTGGCCGGAGATGAATGCCTCGATCACCt                                                                                                     <  1:142294/149‑1 (MQ=255)
                                            atCCGGCGCAACAGTCAGGCCCGAACCGCCCACGTGGTTCCAGGCCACCACCTTGGCTGGCACCTTGCTGCCATAGCGCGCCCATACCGTCATCGGCGAAAGCAGGTGGATGACATTGACCTGGCCGGAGATGAATGCCTCGATCACCt                                                                                                     >  3:68004/1‑149 (MQ=255)
                                            aaccGGCGCAACAGACAGGCACCAACCGCCCCCGTGGTTCCAGGACCCCAACTTGGCTCGACCCTTGCTGCCATAGCGCGCCCATAACGTCATCGGCGAAAGCAGGTGGATGACCTAGACCTGGCAGGAGATGAATGCGTCGCTAACCt                                                                                                     >  4:296506/3‑149 (MQ=255)
                                             tCCGGCGCAACAGTCAGGCCCGAACCGCCCACGTGGTTCCAGGCCACCACCTTGGCTGGCACCTTGCTGCCATAGCGCGCCCATACCGTCATCGGCGAGAGCAGGTGGATGACATTGACCTGGCCGGAGATGAATGCCTCGATCACCTg                                                                                                    >  3:352103/1‑149 (MQ=255)
                                             tCCGGCGCAACAGTCAGGCCCGAACCGCCCACGTGGTTCCAGGCCACCACCTTGGCTGGCACCTTGCTGCCATAGCGCGCCCATACCGTCATCGGCGAGAGCAGGTGGATGACATTGACCTGGCCGGAGATGAATGCCTCGATCACCTg                                                                                                    >  3:357929/1‑149 (MQ=255)
                                             tCCGGCGCAACAGTCAGGCCCGAACCGCCCACGTGGTTCCAGGCCACCACCTTGGCTGGCACCTTGCTGCCATAGCGCGCCCATACCGTCATCGGCGAAAGCAGGTGGATGACATTGACCTGGCCGgag                                                                                                                        >  4:221666/1‑129 (MQ=255)
                                             tCCGGCGCAACAGTCAGGCCCGAACCGCCCACGTGGTTCCAGGCCACCACCTTGGCTGGCACCTTGCTGCCATAGCGCGCCCATACCGTCATCGGCGAAAGCAGGTGGATGACATTGACCTGGCCGgag                                                                                                                        <  3:221666/129‑1 (MQ=255)
                                             tCCGGCGCAACAGTCAGGCCCGAACCGCCCACGTGGTTCCAGGCCACCACCTTGGCTGGCACCTTGCTGCCATAGCGCGCCCATACCGTCATCGGCGAAAGCAGGTGGAGGACATTGACCTGGCCGGAGATGAATGCCTCGATCACCTg                                                                                                    >  1:375929/1‑149 (MQ=255)
                                               cGGCGCAACAGTCAGGCCCGAACCGCCCACGTGGTTCCAGGCCACCACCTTGGCTGGCACCTTGCTGCCATAGCGCGCCCATACCGTCATCGGCGAAAGCAGGTGGATGACATTGACCTGGCCGGAGATGAATGCCTCGATCACCTgc                                                                                                   <  2:375929/148‑1 (MQ=255)
                                                   gcAACAGTCAGGCCCGAACCGCCCACGTGGTTCCAGGCCACCACCTTGGCTGGCACCTTGCTGCCATAGCGCGCCCATACCGTCATCGGCGAAAGCAGGTGGATGACATTGACCTGGCCGGAGATGAATGCCTCGATCACCTGCGcc                                                                                                >  5:244126/1‑147 (MQ=255)
                                                    cAACAGTCAGGCCCGAACCGCCCACGTGGTTCCAGGCCACCACCTTGGCTGGCACCTTGCAGCCATAGCGCGCCCATACCGTCATCGGCGAAAGCAGCTGGATGACATTCACCTGGCCGCAGATGAAAGCCTCGATCACCTGCCCCCAg                                                                                             >  6:321288/1‑149 (MQ=255)
                                                        aGTCAGGCCCGAACCGCCCACGTGGTTCCAGGCCACCACCTTGGCTGGCACCTTGCTGCCATAGCGCGCCCATACCGTCATCGGCGAAAGCAGGTGGATGACATTGACCTGGCCGGAGATGAATGCCTCGATCACCTGCGCCCAGCTgc                                                                                         <  3:207360/149‑1 (MQ=255)
                                                              gCCCGAACCGCGCACGTGGTTCCAGGCCACCACCTTGGCTGGCACCTTGCTGCCATAGCGCGCCCATACCGTCATCGGCGAAAGCAGGTGGATGACATTGACCTGGCCGGAGATGAATGCCTCGATCACCTGCGCCCAGCTGCtcagca                                                                                   <  5:321288/149‑1 (MQ=255)
                                                                      cGCCCACGTGGTTCCAGGCCACCACCTTGGCTGGCACCTTGCTGCCATAGCGCGCCCATACCGTCATCGGCGAAAGCAGGTGGATGACATTGACCTGGCCGGAGATGAATGCCTCGATCACCTGCGCCCAGCTGCGCAGCAGCACCggg                                                                           <  4:233672/149‑1 (MQ=255)
                                                                        cccACGTGGTTCCAGGCCACCACCTTGGCTGGCACCTTGCTGCCATAGCGCGCCCATACCGTCATCGGCGAAAGCAGGTGGATGACATTGACCTGGCCGGAGATGAATGCCTCGATCACCTGCGCCCAGCTGCGCAGCAGCACCGGgcg                                                                         >  1:178060/1‑149 (MQ=255)
                                                                         ccACGTGGTTCCAGGCCACCACCTTGGCTGGCACCTTGCTGCCATAGCGCGCCCATACCGTCATCGGCGAAAGCAGGTGGATGACATTGACCTGGCCGGAGATGAATGCCTCGATCACCTGCGCCCAGCTGCGCAGCAGCACCGGgcgc                                                                        <  7:52340/149‑1 (MQ=255)
                                                                          cACGTGGTTCCAGGCCACCACCTTGGCTGGCACCTTGCTGCCATAGCGCGCCCATACCGTCATCGGCGAAAGCa                                                                                                                                                  >  6:107564/1‑74 (MQ=255)
                                                                          cACGTGGTTCCAGGCCACCACCTTGGCTGGCACCTTGCTGCCATAGCGCGCCCATACCGTCATCGGCGAAAGCa                                                                                                                                                  <  5:107564/74‑1 (MQ=255)
                                                                             gTGGTTCCAGGCCACCACCTTGGCTGGCACCTTGCTGCCATAGCGCGCCCATACCGTCATCGGCGAAAGCAGGTGGATGACATTGACCTGGCCGGAGATGAATGCCTCGATCACCTGCGCCCAGCTGCGCAGCAGCACCGGGCGCTCgg                                                                    <  2:356399/149‑1 (MQ=255)
                                                                                 ttCCAGGCCACCACCTTGGCTGGCACCTTGCTGCCATAGCGCGCCCATACCGTCATCGGCGAAAGCAGGTGGATGACATTGACCTGGCCGGAGATGAATGCCTCGATCACCTGCGCCCAGCTGCGCAGCAGCACCGGGCGCTCGGCCtt                                                                >  5:277799/1‑149 (MQ=255)
                                                                                    cAGGCCACCACCTTGGCTGGCACCTTGCTGCCATAGCGCGCCCATACCGTCATCGGCGAAAGCAGGTGGATGACATTGACCTGGCCGGAGATGAATGCCTCGATCACCTGCGCCCAGCTGCGCAGCAGCACCGGGCGCTCGGCCTTGAt                                                             <  6:244126/149‑1 (MQ=255)
                                                                                    cAGGCCACCACCTTGGCTGGCACCTTGCTGCCATAGCGCGCCCATACCGTCATCGGCGAAAGCAGGTGGATGACATTGACCTGGCCGGAGATGAATGCCTCGATCACCTGCGCCCAGCTGCGCAGCAGCACCGGGCGCTCGGCCTTGAt                                                             >  5:2009/1‑149 (MQ=255)
                                                                                    cAGGCCACCACCTTGGCTGGCACCTTGCTGCCATAGCGCGCCCATACCGTCATCGGCGAAAGCAGGTGGATGACATTGACCTGGCCGGAGATGAATGCCTCGATCACCTGCGCCCAGCTGCGCAGCAGCACCGGGCGCTCGGCCTTGAt                                                             >  5:193798/1‑149 (MQ=255)
                                                                                    cAGGCCACCACCTTGGCTGGCACCTTGCTGCCATAGCGCGCCCATACCGTCATCGGCGAAAGCAGGTGGATGACATTGACCTGGCCGGAGATGAATGCCTCGATCACCTGCGCCCAGCTGCGCAGCAGCACCGGGCGCTCGGCCTTGAt                                                             <  4:68004/149‑1 (MQ=255)
                                                                                     aGGCCACCACCTTGGCTGGCACCTTGCTGCCATAGCGCGCCCATACCGTCATCGGCGAAAGCAGGTGGATGACATTGACCTGGCCGGAGATGAATGCCTCGATCACCTGCGCCCAGCTGCGCAGCAGCACCGGGCGCTCGGCCTTGATg                                                            >  6:372976/1‑149 (MQ=255)
                                                                                     aGGCCACCACCTTGGCTGGCACCTTGCTGCCATAGCGCGCCCATACCGTCATCGGCGAAAGCAGGTGGATGACATTGACCTGGCCGGAGATGAATGCCTCGATCACCTGCGCCCAGCTGCGCAGCAGCACCGGGCGCTCGGCCTTGATg                                                            >  2:329804/1‑149 (MQ=255)
                                                                                       gCCACCACCTTGGCTGGCACCTTGCTGCCATAGCGCGCCCATACCGTCATCGGCGAAAGCAGGTGGATGACATTGACCTGGCCGGAGATGAATGCCTCGATCACCTGCGCCCAGCTGCGCAGCAGCACCgg                                                                            <  3:75287/131‑1 (MQ=255)
                                                                                       gCCACCACCTTGGCTGGCACCTTGCTGCCATAGCGCGCCCATACCGTCATCGGCGAAAGCAGGTGGATGACATTGACCTGGCCGGAGATGAATGCCTCGATCACCTGCGCCCAGCTGCGCAGCAGCACCgg                                                                            >  4:75287/1‑131 (MQ=255)
                                                                                             accTTGGCTGGCACCTTGCTGCCATAGCGCGCCCATACCGTCATCGGCGAAAGCAGGTGGATGACATTGACCTGGCCGGAGATGAATGCCTCGATCACCTGCGCCCAGCTGCGCAGCAGCACCGGGCGCTCGGCCTTGATGCCTTCGGc                                                    <  3:291157/149‑1 (MQ=255)
                                                                                               cTTGGCTGGCACCTTGCTGCCATAGCGCGCCCATACCGTCATCGGCGAAAGCAGGTGGATGACATTGACCTGGCCGGAGATGAATGCCTCGATCACCTGCGCCCAGCTGCGCAGCAGCACCGGGCGCTCGGCCTTGATGCCTTCGGCCt                                                  >  7:35320/1‑149 (MQ=255)
                                                                                                  ggctggcACCTTGCTGCCATAGCGCGCCCATACCGTCATCGGCGAAAGCAGGTGGATGACATTGACCTGGCCGGAGATGAATGCCTCGATCACCTGCGCCCAGCTGCGCAGCAGCACCGGGCGCTCGGCCTTGATGCCTTCGGCCTCg                                                >  5:381786/1‑148 (MQ=255)
                                                                                                          ccTTGCTGCCATAGCGCGCCCATACCGTCATCGGCGAAAGCAGGTGGATGACATTGACCTGGCCGGAGATGAATGCCTCGATCACCTGCGCCCAGCTGCGCAGCAGCACCGGGCGCTCGGCCTTGATGCCTTCGGCCTCGAACAGGCCg                                       <  1:8975/149‑1 (MQ=255)
                                                                                                              gctgcCATAGCGCGCCCATACCGTCATCGGCGAAAGCAGGTGGATGACATTGACCTGGCCGGAGATGAATGCCTCGATCACCTGCGCCCAGCTGCGCAGCAGCACCGGGCGCTCGGCCTTGATGCCTTCGGCCTCGAACAGGCCgttg                                    <  6:2009/148‑1 (MQ=255)
                                                                                                              gctgcCATAGCGCGCCCATACCGTCATCGGCGAAAGCAGGTGGATGACATTGACCTGGCCGGAGATGAATGCCTCGATCACCTGCGCCCAGCTGCGCAGCAGCACCGGGCGCTCGGCCTTGATGCCTTCGGCCTCGAACAGGCCgtt                                     <  6:78393/147‑1 (MQ=255)
                                                                                                              gctgcCATAGCGCGCCCATACCGTCATCGGCGAAAGCAGGTGGATGACATTGACCTGGCCGGAGATGAATGCCTCGATCACCTGCGCCCAGCTGCGCAGCAGCACCGGGCGCTCGGCCTTGATGCCTTCGGCCTCGAACAGGCCgtt                                     <  6:193798/147‑1 (MQ=255)
                                                                                                               ttgacttAGCGCGCCGATAGCGTGAGCTTCGAATGCATGTGGATGACATTGACCTTGCCGGTGATGAATGGCTCGATCACCTGCGCCCAGCTGCGCAGCAGCACCGGGCGCTCTGCCTTGATGCCTTCGGCCTCGAAAAGGCCgttgtt                                  <  8:187672/143‑1 (MQ=255)
                                                                                                                tgcCATAGCGCGCCCATACCGTCATCGGCGAAAGCAGGTGGATGACATTGACCTGGCCGGAGATGAATGCCTCGATCACCTGCGCCCAGCTGCGCAGCAGCACCGGGCGCTCGGCCTTGATGCCTTCGGCCTCGAACAGGCCGTTGTTa                                 >  6:168225/1‑149 (MQ=255)
                                                                                                                  ccATAGCGCGCCCATACCGTCATTGGCGAAAGCAGGTGGATGACATTGACCTGGCCGGAGATGAATGCCTCGATCACCTGCGCCCAGCTGCGCAGCAGCACCGGGCGCTCGGCCTTGATGCCTTCGGCCTCGAACAGGCCGTTGTTAtg                               >  8:363733/1‑149 (MQ=255)
                                                                                                                   cATAGCGCGCCCATACCGTCATCGGCGAAAGCAGGTGGATGACATTGACCTGGCCGGAGATGAATGCCTCGATCACCTGCGCCCAGCTGCGCAGCAGCACCGGGCGCTCGGCCTTGATGCCTTCGGCCTCGAACAGGCCGTTGTTAtgt                              <  8:180935/149‑1 (MQ=255)
                                                                                                                   cATAGCGCGCCCATACCGTCATCGGCGAAAGCAGGTGGATGACATTGACCTGGCCGGAGATGAATGCCTCGATCACCTGCGCCCAGCTGCGCAGCAGCACCGGGCGCTCGGCCTTGATGCCTTCGGCCTCGAACAGGCCGTTGTTAtgt                              <  1:156511/149‑1 (MQ=255)
                                                                                                                        cgcgcCCATACCGTCATCGGCGAGAGCAGGTGGATGACATTGACCTGGCCGGAGATGAATGCCTCGATCACCTGCGCCCAGCTGCGCAGCAGCACCGGGCGCTCGGCCTTGATGCCTTCGGCCTCGAACAGGCCGTTGTTATGTGccac                         <  4:357929/149‑1 (MQ=255)
                                                                                                                        cgcgcCCATACCGTCATCGGCGAGAGCAGGTGGATGACATTGACCTGGCCGGAGATGAATGCCTCGATCACCTGCGCCCAGCTGCGCAGCAGCACCGGGCGCTCGGCCTTGATGCCTTCGGCCTCGAACAGGCCGTTGTTATGTGccac                         <  4:352103/149‑1 (MQ=255)
                                                                                                                           gcCCATACCGTCATCGGCGAAAGCAGGTGGATGACATTGACCTGGCCGGAGATGAATGCCTCGATGACCTGCGCCCAGCTGCGCAGCAGCACCGGGCGCTCGGCCTTGATGCCTTCGGCCTCGAACAGGCCGTTGTTATGTGCCACcag                      <  5:262027/149‑1 (MQ=255)
                                                                                                                           gcCCATACCGTCATCGGCGAAAGCAGGTGGATGACATTGACCTGGCCGGAGATGAATGCCTCGATCACCTGCGCCCAGCTGCGCAGCAGCACCGGGCGCTCGGCCTTGATGCCTTCGGCCTCGAACAGGCCGTTGTTATGTGCCACcag                      >  4:114129/1‑149 (MQ=255)
                                                                                                                             ccATACCGTCATCGGCGAAAGCAGGTGGATGACATTGACCTGGCCGGAGATGAATGCCTCGATCACCTGCGCCCAGCTGCGCAGCAGCACCGGGCGCTCGGCCTTGATGCCTTCGGCCTCGAACAGGCCGTTGTTATGTGCCACcag                      >  6:282762/1‑147 (MQ=255)
                                                                                                                                   cGTCATCGGCGAAAGCAGGTGGATGACATTGACCTGGCCGGAGATGAATGCCTCGATCACCTGCGCCCAGCTGCGCAGCAGCACCGGGCGCTCGGCCTTGATGCCTTCGGCCTCGAACAGGCCGTTGTTATGTGCCACCAGCAgcggcg              >  6:228675/1‑149 (MQ=255)
                                                                                                                                     tCATCGGCGAAAGCAGGTGGATGACATTGACCTGGCCGGAGATGAATGCCTCGATCACCTGCGCCCAGCTGCGCAGCAGCACCGGGCGCTCGGCCTTGATGCCTTCGGCCTCGAACAGGCCGTTGTTATGTGCCACCAGCAGCGGCGTg            <  5:168225/149‑1 (MQ=255)
                                                                                                                                       aTCGGCGAAAGCAGGTGGATGACATTGACCTGGCCGGAGATGAATGCCTCGATCACCTGCGCCCAGCTGCGCAGCAGCACCGGGCGCTCGGCCTTGATGCCTTCGGCCTCGAACAGGCCGTTGTTATGTGCCACCAGCAGCGGCGTgg           <  5:372976/148‑1 (MQ=255)
                                                                                                                                        tCGGCGAAAGCAGGTGGATGACATTGACCTGGCCGGAGATGAATGCCTCGATCACCTGCGCCCAGCTGCGCAGCAGCACCGGGCGCTCGGCCTTGATGCCTTCGGCCTCGAACAGGCCGTTGTTATGTGccacca                       <  8:206807/135‑1 (MQ=255)
                                                                                                                                           gCGAAAGCAGGTGGATGACATTGACCTGGCCGGAGATGAATGCCTCGATCACCTGCGCCCAGCTGCGCAGCAGCACCGGGCGCTCGGCCTTGATGCCTTCGGCCTCGAACAGGCCGTTGTTATGTGCCACCAGCAGCGGCGTGGCATCg      >  6:306449/1‑149 (MQ=255)
                                                                                                                                               aaGCAGGTGGATGACATTGACCTGGCCGGAGATGAATGCCTCGATCACCTGCGCCCAGCTGCGCAGCAGCACCGGGCGCTCGGCCTTGATGCCTTCGGCCTCGAACAGGCCGTTGTTATGTGCCACCAGCAGCGGCGTGGCATCGGTGa  <  5:228675/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
GGAATGGCCACTGTCTTGCCACCGAGCTGTTTCACGGCGGATATATCCGGCGCAACAGTCAGGCCCGAACCGCCCACGTGGTTCCAGGCCACCACCTTGGCTGGCACCTTGCTGCCATAGCGCGCCCATACCGTCATCGGCGAAAGCAGGTGGATGACATTGACCTGGCCGGAGATGAATGCCTCGATCACCTGCGCCCAGCTGCGCAGCAGCACCGGGCGCTCGGCCTTGATGCCTTCGGCCTCGAACAGGCCGTTGTTATGTGCCACCAGCAGCGGCGTGGCATCGGTGA  >  NC_002947/4132712‑4133003

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