Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,330,968 G→A 6.8% R138R (CGG→CGA PP_3802 → cation ABC transporter ATP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,330,9680GA 6.8% 174.1 / 2.1 59R138R (CGG→CGAPP_3802cation ABC transporter ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (28/27);  new base A (2/2);  total (30/29)
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.40e-01

TTCGACCGCGGCCTGCCCATGACCGTGGACGATTTCATGGCCGCCATGTGCCAGCGCCGCCCGGCCTTCCTCGGCCTGTCGCGGCGCGTGAAGCCGGCCATTGATGCAGCCTTGGTTCAGGTTGGCATGCTCGACAAGCGCACGCGGCGCATGGGGGCACTGTCCGGCGGTGAGCGCCAGCGCGTGCTGCTGGCCCAGGGCCTGATCCCCGAGCCGCAGTTGCTGGTGCTGGACGAACCGATGTCGGCCCTCGATGAAGCCGGAATCCAGGTGTTCGAACAGTTATTG  >  NC_002947/4330822‑4331109
                                                                                                                                                   |                                                                                                                                             
ttCGACCGCGGCCTGCCCATGACCGTGGACGATTTCATGGCCGCCATGTGCCAGCGCCGCCCGGCCTTCCTCGGCCTGTCGCGGCGCGTGAAGCCGGCCATTGATGCAGCCTTGGTTCAGGTTGGCATGCTCGACAAGCGCACgcggcg                                                                                                                                             >  1:3068/1‑149 (MQ=255)
ttCGACCGCGGCCTGCCCATGACCGTGGACGATTTCATGGCCGCCATGTGCCAGCGCCGCCCGGCCTTCCTCGGCCTGTCGCGGCGCGTGAAGCCGGCCATTGATGCAGCCTTGGTTCAGGTTGGCATGCTCGACAAGCGCACgcggcg                                                                                                                                             >  8:264129/1‑149 (MQ=255)
   gACCGCGGCCTGCCCATGACCGTGGACGATTTCGATGGCCGCCATGTGCCAGCGCCGCCCGGCCTTCCTCGGCCTGTCGCGGCGCGTGAAGCCGGCCATTGATGCAGCCTTGGTTCAGGTTGGCATGCTCGACAAGCGCACGCGGCGCa                                                                                                                                           >  1:41010/1‑149 (MQ=255)
   gACCGCGGCCTGCCCATGACCGTGGACGATTTCATGGCCGCCATGTGCCAGCGCCGCCCGGCCTTCCTCGGCCTGTCGCGGCGCGTGAAGCCGGCCATTGATGCAGCCTTGGTTCAGGTTGGCATGCTCGACAAGCGCACGCGGCGCAt                                                                                                                                          >  3:184599/1‑149 (MQ=255)
       gcgTCCTGCCCATGACCGTGGACGATTTCATGGCCGCCATGTGCCAGCGCCGCCCGGCCTTCCTCGGCCTGTCGCGGCGCGTGAAGCCGGCCATTGATGCAGCCTTGGTTCAGGTTGGCATGCTCGACAAGCGCACGCGGCGCATgggg                                                                                                                                      <  3:224923/149‑1 (MQ=255)
          gCCTGCCCATGCCCGTGGACGATGTCATGGCCGCCAGGTCCCAGCGCCGCCCGGCCTTCCTCGGCCTGTCGCGGCGCGTGAAGCCGGCCATTGATGCAGCCTTGGTTCAGGTTGGCATGCTCGACAAGCGCACGCGGCGCATGGGGGCa                                                                                                                                   <  4:2456/149‑1 (MQ=255)
          gCCTGCCCATGACCGTGGACGATTTCATGGCCGCCATGTGCCAGCGCCGCCCGGCCTTCCTCGGCCTGTCGCGGCGCGTGAAGCCGGCCATTGATGCAGCCTTGGTTCAGGTTGGCATGCTCGACAAGCGCACGCGGCGCAt                                                                                                                                          <  7:189282/142‑1 (MQ=255)
          gCCTGCCCATGACCGTGGACGATTTCATGGCCGCCATGTGCCAGCGCCGCCCGGCCTTCCTCGGCCTGTCGCGGCGCGTGAAGCCGGCCATTGATGCAGCCTTGGTTCAGGTTGGCATGCTCGACAAGCGCACGCGGCGCAt                                                                                                                                          >  8:189282/1‑142 (MQ=255)
             tGCCCATGACCGTGGACGATTTCATGGCCGCCATGTGCCAGCGCCGCCCGGCCTTCCTCGGCCTGTCGCGGCGCGTGAAGCCGGCCATTGATGCAGCCTTGGTTCAGGTTGGCATGCTCGACAAGCGCACGCGGCGCATGGGGGCACTg                                                                                                                                >  8:108705/1‑149 (MQ=255)
                  aTGACCGTGGACGATTTCATGGCCGCCATGTGCCAGCGCCGCCCGGCCTTCCTCGGCCTGTCGCGGCGCGTGAAGCCGGCCATTGATGCAGCCTTGGTTCAGGTTGGCATGCTCGACAAGCGCACGCGGCGCATGGGGGCACTGTCcgg                                                                                                                           >  1:299288/1‑149 (MQ=255)
                   tGACCGTGGACGATTTCATGGCCGCCATGTGCCAGCGCCGCCCGGCCTTCCTCGGCCTGTCGCGGCGCGTGAAGCCGGCCATTGATGCAGCCTTGGTTCAGGTTGGCATGCTCGACAAGCGCACGCGACGCATGGGGGCACTGTCcggc                                                                                                                          >  5:23426/1‑149 (MQ=255)
                   tGACCGTGGACGATTTCATGGCCGCCATGTGCCAGCGCCGCCCGGCCTTCCTCGGCCTGTCGCGGCGCGTGAAGCCGGCCATTGATGCAGCCTTGGTTCAGGTTGGCATGCTCGACAAGCGCACGCGACGCATGGGGGCACTGTCcggc                                                                                                                          >  5:27886/1‑149 (MQ=255)
                      ccGTGGACGATTTCATGGCCGCCATGTGCCAGCGCCGCCCGGCCTTCCTCGGCCTGTCGCGGCGCGTGAAGCCGGCCATTGATGCAGCCTTGGTTCAGGTTGGCATGCTCGACAAGCGCACGCGGCGCATGGGGGCACTGTCCGGCGGt                                                                                                                       >  2:270925/1‑149 (MQ=255)
                                     aTGGCCGCCATGTGCCAGCGCCGCCCGGCCTTCCTCGGCCTGTCGCGGCGCGTGAAGCCGGCCATTGATGCAGCCTTGGTTCAGGTTGGCATGCTCGACAAGCGCACGCGGCGCATGGGGGCACTGTCCGGCGGTGAGCGCCAGCGCGt                                                                                                         <  7:264129/149‑1 (MQ=255)
                                     aTGGCCGCCATGTGCCAGCGCCGCCCGGCCTTCCTCGGCCTGTCGCGGCGCGTGAAGCCGGCCATTGATGCAGCCTTGGTTCAGGTTGGCATGCTCGACAAGCGCACGCGACGCATGGGGGCACTGTCCGGCGGTGAGCGCCAGcgcg                                                                                                          <  6:23426/148‑1 (MQ=255)
                                     aTGGCCGCCATGTGCCAGCGCCGCCCGGCCTTCCTCGGCCTGTCGCGGCGCGTGAAGCCGGCCATTGATGCAGCCTTGGTTCAGGTTGGCATGCTCGACAAGCGCACGCGACGCATGGGGGCACTGTCCGGCGGTGAGCGCCAGcgcg                                                                                                          <  6:27886/148‑1 (MQ=255)
                                      tGGCCGCCATGTGCCAGCGCCGCCCGGCCTTCCTCGGCCTGTCGCGGCGCGTGAAGCCGGCCATTGATGCAGCCTTGGTTCAGGTTGGCATGCTCGACAAGCGCACGCGGCGCATGGGGGCACTGTCCGGCGGTGAGCGCCAGCGCGtg                                                                                                        <  6:122767/149‑1 (MQ=255)
                                      tGGCCGCCATGTGCCAGCGCCGCCCGGCCTTCCTCGGCCTGTCGCGGCGCGTGAAGCCGGCCATTGATGCAGCCTTGGTTCAGGTTGGCATGCTCGACAAGCGCACGCGGCGCATGGGGGCACTGTCCGGCGGTGAGCGCCAGCGCGtg                                                                                                        >  3:262956/1‑149 (MQ=255)
                                      tGGCCGCCATGTGCCAGCGCCGCCCGGCCTTCCTCGGCCTGTCGCGGCGCGTGAAGCCGGCCATTGATGCAGCCTTGGTTCAGGTTGGCATGCTCGACAAGCGCACGCGGCGCATGGGGGCACTGTCCGGCGGTGAGCGCCAGCGCGtg                                                                                                        >  8:142873/1‑149 (MQ=255)
                                           gccATGTGCCAGCGCCGCCCGGCCTTCCTCGGCCTGTCGCGGCGCGTGAAGCCGGCCATTGATGCAGCCTTGGTTCAGGTTGGCATGCTCGACAAGCGCACGCGGCGCATGGGGGCACTGTCCGGCGGTGAGCGCCAGCGCGTgctgct                                                                                                   >  6:12084/1‑149 (MQ=255)
                                             cATGTGCCCGCGCCGCCCGGCCGTCCGCGCCCGGGCGCGGCGCGTGAAGCCGGCCATTGACGCCGCCTCGGTTCAGGTTGGCATGCTCGACAAGCGCACGCGGCGCAGGGGGGCACTGTCCGGCGGTGAGCGCCAGCGCGTGCTGCTgg                                                                                                 <  8:125121/149‑1 (MQ=255)
                                             cATGTGCCAGCGCCGCCCGGCCTTCCTCGGCCTGTCGCGGCGCGTGAAGCCGGCCATTGATGCAGCCTTGGTTCAGGTTGGCATGCTCGACAAGCGCACGCGGCGCATGGGGGCAGTGTCCGGCGGTGAGCGCCAGCGCGTGCTGCTgg                                                                                                 <  2:41010/149‑1 (MQ=255)
                                             cATGTGCCAGCGCCGCCCGGCCTTCCTCGGCCTGTCGCGGCGCGTGAAGCCGGCCATTGATGCAGCCTTGGTTCAGGTTGGCATGCTCGACAAGCGCACGCGGCGCATGGGGGCACTGTCCGGCGGTGAGCGCCAGCGCGTGCTGCTgg                                                                                                 <  7:142873/149‑1 (MQ=255)
                                              aTGTGCCAGCGCCGCCCGGCCTTCCTCGGCCTGTCGCGGCGCGTGAAGCCGGCCATTGATGCAGCCTTGGTTCAGGTTGGCATGCTCGACAAGCGCACGCGGCGCATGGGGGCACTGTCCGGCGGTGAGCGCCAGCGCGTGCTGCTGGc                                                                                                <  8:74722/149‑1 (MQ=255)
                                               tgtgCCAGCGCCGCCCGGCCTTCCTCGGCCTGTCGCGGCGCGTGAAGCCGGCCATTGATGCAGCCTTGGTTCAGGTTGGCATGCTCGACAAGCGCACGCGGCGCATGGGGGCACTGTCCGGTGGTGAGCGCCAGCGCGTGCTGCTGGcc                                                                                               <  4:184599/149‑1 (MQ=255)
                                                   ccAGCGCCGCCCGGCCTTCCTCGGCCTGTCGCGGCGCGTGAAGCCGGCCATTGATGCAGCCTTGGTTCAGGTTGGCATGCTCGACAAGCGCACGCGGCGCATGGGGGCACTGTCCGGCGGTGAGCGCCAGCGCGTGCTGCTGGCCCAgg                                                                                           >  6:203882/1‑149 (MQ=255)
                                                    cAGCGCCGCCCGGCCTTCCTCGGCCTGTCGCGGCGCGTGAAGCCGGCCATTGATGCAGCCTTGGTTCAGGTTGGCATGCTCGACAAGCGCACGCGGCGCATGGGGGCACTGTCCGGCGGTGAGCGCCAGCGCGTGCTGCTGGCCCAggg                                                                                          >  7:157322/1‑149 (MQ=255)
                                                         ccgccCGGCCTTCCTCGGCCTGTCGCGGCGCGTGAAGCCGGCCATTGATGCAGCCTTGGTTCAGGTTGGCATGCTCGACAAGCGCACGCGGCGCATGGGGGCACTGTCCGGCGGTGAGCGCCAGCGCGTGCTGCTGGCCCAGGGCCTGa                                                                                     <  5:203882/149‑1 (MQ=255)
                                                         ccgccCGGCCTTCCTCGGCCTGTCGCGGCGCGTGAAGCCGGCCATTGATGCAGCCTTGGTTCAGGTTGGCATGCTCGACAAGCGCACGCGGCGCATGGGGGCACTGTCCGGCGGTGAGCGCCAGCGCGTGCTGCTGGCCCAGGGCCTGa                                                                                     <  8:233625/149‑1 (MQ=255)
                                                           gccCGGCCTTCCTCGGCCTGTCGCGGCGCGTGAAGCCGGCCATTGATGCAGCCTTGGTTCAGGTTGGCATGCTCGACAAGCGCACGCGGCGCATGGGGGCACTGTCCGGCGGTGAGCGCCAGCGCGTGCTGCTGGCCCAGGGCCTGAt                                                                                    >  5:148522/1‑148 (MQ=255)
                                                           gccCGGCCTTCCTCGGCCTGTCGCGGCGCGTGAAGCCGGCCATTGATGCAGCCTTGGTTCAGGTTGGCATGCTCGACAAGCGCACGCGGCGCATGGGGGCACTGTCCGGCGGTGAGCGCCAGCGCGTGCTGCTGGCCCAGGGCCTGAt                                                                                    <  6:148522/148‑1 (MQ=255)
                                                             ccGGCCTTCCTCGGCCTGTCGCGGCGCGTGAAGCCGGCCATTGATGCAGCCTTGGTTCAGGTTGGCATGCTCGACAAGCGCACGCGGCGCATGGGGGCACTGTCCGGCGGTGAGCGCCAGCGCGTGCTGCTGGCCCAGGGCCTGATccc                                                                                 <  2:3068/149‑1 (MQ=255)
                                                                gCCTTCCTCGGCCTGTCGCGGCGCGTGAAGCCGGCCATTGATGCAGCCTTGGTTCAGGTTGGCATGCTCGACAAGCGCACGCGGCGCATGGGGGCACTGTCCGGCGGTGAGCGCCAGCGCGTGCTGCTGGCCCAGGGCCTGATCCCCGa                                                                              <  1:11783/149‑1 (MQ=255)
                                                                 ccTTCCTCGGCCTGTCGCGGCGCGTGAAGCCGGCCATTGATGCAGCCTTGGTTCAGGTTGGCATGCTCGACAAGCGCACGCGGCGCATGGGGGCACTGTCCGGCGGTGAGCGCCAGCGCGTGCTGCTGGCCCAGGGCCTGATCCCCGAg                                                                             <  1:270925/149‑1 (MQ=255)
                                                                   ttCCTCGGCCTGTCGCGGCGCGTGAAGCCGGCCATTCATGCAGCCTTGGTTCAGGTTGGCATGCTCGACAAGCGCACGCGGCGAATGGTGGCACTGTCCGGCGGTGATCTCCAGCGCGTCCTGCTGGCCCAGGGCCTGATACCCGAGcc                                                                           >  6:248860/1‑149 (MQ=255)
                                                                              gTCGCGGCGCGTGAAGCCGGCCATTGATGCAGCCTTGGTTCAGGTTGGCATGCTCGACAAGCGCACGCGGCGCATGGGGGCACTGTCCGGCGGTGAGCGCCAGCGCGTGCTGCTGGCCCAGGGCCTGATCCCCGAGCCGCAGTtgctgg                                                                >  1:200953/1‑149 (MQ=255)
                                                                                     cgcgTGAAGCCGGCCATTGATGCAGCCTTGGTTCAGGGTGGCATGCTCGACAAGCGCACGCGGCGCATGGGGGCACTGTCCGGCGGTGAGCGCCAGCGCGTGCTGCTGGCCCAGGGCCTGATCCCCGAGCCGCAGTTGCTGGTGCTGGa                                                         <  8:230077/149‑1 (MQ=255)
                                                                                            aGCCGGCCATTGATGCAGCCTTGGTTCAGGTTGGCATGCTCGACAAGCGCACGCGGCGCATGGGGGCACTGTCCGGCGGTGAGCGCCAGCGCGTGCTGCTGGCCCAGGGCCTGATCCCCGAGCCGCAGTTGCTGGTGCTGGACGAACCg                                                  <  2:200953/149‑1 (MQ=255)
                                                                                                   cATTGATGCAGCCTTGGTTCAGGTTGGCATGCTCGACAAGCGCACGCGGCGCATGGGGGCACTGTCCGGCGGTGAGCGCCAGCGCGTGCTGCTGGCCCAGGGCCTGATCCCCGAGCCGCAGTTGCTGGTGCTGGACGAACCGATGTCgg                                           >  8:98021/1‑149 (MQ=255)
                                                                                                        aTGCAGCCTTGGTTCAGGTTGGCATGCTCGACAAGCGCACGCGGCGCATGGGGGCACTGTCCGGCGGTGAGCGCCAGCGCGTGCTGCTGGCCCAGGGCCTGATCCCCGAGCCGCAGTTGCTGGTGCTGGACGAACCGATGTCgg                                           >  8:222273/1‑144 (MQ=255)
                                                                                                        aTGCAGCCTTGGTTCAGGTTGGCATGCTCGACAAGCGCACGCGGCGCATGGGGGCACTGTCCGGCGGTGAGCGCCAGCGCGTGCTGCTGGCCCAGGGCCTGATCCCCGAGCCGCAGTTGCTGGTGCTGGACGAACCGATGTCgg                                           <  7:222273/144‑1 (MQ=255)
                                                                                                               cTTGGTTCAGGTTGGCATGCTCGACAAGCGCACGCGGCGCATGGGGGCACTGTCCGGCGGTGAGCGCCAGCGCGTGCTGCTGGCCCAGGGCCTGATCCCCGAGCCGCAGTTGCTGGTGCTGGACGAACCGATGTCGGCCCTCGATGAAg                               >  4:231337/1‑149 (MQ=255)
                                                                                                                   gTTCAGGTTGGCATGCTCGACAAGCGCACGCGGCGCATGGGGGCACTGTCCGGCGGTGAGCGCCAGCGCGTGCTGCTGGCCCAGGGCCTGATCCCCGAGCCGCAGTTGCTGGTGCTGGACGAACCGATGTCGGCCCTCGATGAAGCCg                            <  3:231337/148‑1 (MQ=255)
                                                                                                                         gTTGGCATGCTCGACAAGCGCACGCGGCGCATGGGGGCACTGTCCGGCGGTGAGCGCCAGCGCGTGCTGCTGGCCCAGGGCCTGATCCCCGAGCCGCAGTTGCTGGTGCTGGACGAACCGATGTCGGCCCTCGATGAAGCCGGAATc                       >  1:302006/1‑147 (MQ=255)
                                                                                                                          ttGGCATGCTCGACAAGCGCACGCGGCGCATGGGGGCACTGTCCGGCGGTGAGCGCCAGCGCGTGCTGCTGGCCCAGGGCCTGATCCCCGAGCCGCAGTTGCTGGTGCTGGACGAACCGATGTCGGCCCTCGATGAAGCCGGAATCCAg                    <  8:157322/149‑1 (MQ=255)
                                                                                                                          gtGGCATGCTCGACAAGCGCACGCGGCGCATGGGGGCACTGTCGGGCGGTGAGCGCCAGCGCGTGCTGCTGGCCCAGGGCCTGATCCCCGAGCCGCAGTTGCTGGTGCTGGACGAACCGATGTCGGCCCTCGATGAAGCCGGAATCCAg                    <  5:248860/148‑1 (MQ=255)
                                                                                                                               aTGCTCGACAAGCGCACGCGGCGCATGGGGGCACTGTCCGGCGGTGAGCGCCAGCGCGTGCTGCTgg                                                                                                 >  6:27710/1‑67 (MQ=255)
                                                                                                                               aTGCTCGACAAGCGCACGCGGCGCATGGGGGCACTGTCCGGCGGTGAGCGCCAGCGCGTGCTGCTGGCCCAGGGCCTGATCCCCGAGCCGCAGTTGCTGGTGCTGGACGAACCGATGTCGGCCCTCGATGAAGCCGGAATCCAGGTGtt               >  8:29430/1‑149 (MQ=255)
                                                                                                                                tGCTCGACAAGCGCACGCGGCGCATGGGGGCACTGTCCGGCGGTGAGCGCCAGCGCGTGCTGCTGGCCCAGGGCCTGATCCCCGAGCCGCAGTtgctg                                                                 >  5:194233/1‑98 (MQ=255)
                                                                                                                                tGCTCGACAAGCGCACGCGGCGCATGGGGGCACTGTCCGGCGGTGAGCGCCAGCGCGTGCTGCTGGCCCAGGGCCTGATCCCCGAGCCGCAGTtgctg                                                                 >  5:235467/1‑98 (MQ=255)
                                                                                                                                tGCTCGACAAGCGCACGCGGCGCATGGGGGCACTGTCCGGCGGTGAGCGCCAGCGCGTGCTGCTGGCCCAGGGCCTGATCCCCGAGCCGCAGTtgctg                                                                 <  6:235467/98‑1 (MQ=255)
                                                                                                                                tGCTCGACAAGCGCACGCGGCGCATGGGGGCACTGTCCGGCGGTGAGCGCCAGCGCGTGCTGCTGGCCCAGGGCCTGATCCCCGAGCCGCAGTtgctg                                                                 <  6:194233/98‑1 (MQ=255)
                                                                                                                                 gCTCGACAAGCGCACGCGGCGCATGGGGGCACTGTCCGGCGGTGAGCGCCAGCGCGTGCTGCTGGCCCAGGGCCTGATCCCCGAGCCGCAGTTGCTGGTGCTGGACGAACCGATGTCGGCCCTCGATGAAGCCGGAATCCAGGTGTTc              >  3:33468/1‑148 (MQ=255)
                                                                                                                                   tCGACAAGCGCACGCGGCGCATGGGGGCACTGTCCGGCGGTGAGCGCCAGCGCGTGCTGCTGGCCCAGGGCCTGATCCCCGAGCCGCAGTTGCTGGTGCTGGACGAACCGATGTCGGCCCTCGATGAAGCCGGAATc                       >  5:259820/1‑137 (MQ=255)
                                                                                                                                   tCGACAAGCGCACGCGGCGCATGGGGGCACTGTCCGGCGGTGAGCGCCAGCGCGTGCTGCTGGCCCAGGGCCTGATCCCCGAGCCGCAGTTGCTGGTGCTGGACGAACCGATGTCGGCCCTCGATGAAGCCGGAATc                       <  6:259820/137‑1 (MQ=255)
                                                                                                                                     gACAAGCGCACGCGGCGCATGGGGGCACTGTCCGGCGGTGAGCGCCAGCGCGTGCTGCTGGCCCAGGGCCTGATCCCCGAGCCGCAGTTGCTGGTGCTGGACGAACCGATGTCGGCCCTCGATGAAGCCGGAATCCAGGTGTTCGaa           <  4:262956/147‑1 (MQ=255)
                                                                                                                                          gcgcACGCGGCGCATGGGGGCACTGTCCGGCGGTGAGCGCCAGCGCGTGCTGCTGGCCCAGGGCCTGATCCCCGAGCCGCAGTTGCTGGTGCTGGACGAACCGATGTCGGCCCTCGATGAAGCCGGAATCCAGGTGTTCGAACAGTTAt    >  2:27125/1‑149 (MQ=255)
                                                                                                                                          gcgcACGCGGCGCATGGGGGCACTGTCCGGCGGTGAGCGCCAGCGCGTGCTGCTGGCCCAGGGCCTGATCCCCGAGCCGCAGTTGCTGGTGCTGGACGAACCCATGTCGGCCCTCGATGAAGCCGGAATCCAGGTGTTCGAACAGTTAt    >  3:102872/1‑149 (MQ=255)
                                                                                                                                           cgcACGCGGCGCATGGGGGCACGGTCCGGCGGTGAGCGCCAGCGCGTGCTGCTGGCCCAGGGCCTGATCCCCGAGCCGCAGTTGCTGGTGCTGGACGAACCGATGTCGGCCCTCGATGAAGCCGGAATCCAGGTGTTCGAACAGTTAtt   <  4:275031/149‑1 (MQ=255)
                                                                                                                                            gcACGCGGCGCATGGGGGCACTGTCCGGCGGTGAGCGCCAGCGCGTGCTGCTGGCCCAGGGCCTGATCCCCGAGCCGCAGTTGCTGGTGCTGGACGAACCGATGTCGGCCCTCGATGAAGCCGGAATCCAGGTGTTCGAACAGTTATTg  >  4:161794/1‑149 (MQ=255)
                                                                                                                                              aCGCGGCGCATGGGGGCACTGTCCGGCGGTGAGCGCCAGCGCGTGCTGCTGGCCCAGGGCCTGATCCCCGAGCCGCAGTTGCTGGTGCTGGACGAACCGATGTCGGCCCTCGATGAAGCCGGAATCCAGGTGTTCGAACAGTTATTg  >  8:150865/1‑147 (MQ=255)
                                                                                                                                                gcggcgCATGGGGGCACTGTCCGGCGGTGAGCGCCAGCGCGTGCTGCTGGCCCAGGGCCTGATCCCCGAGCCGCAGTTGCTGGTGCTGGACGAACCGATGTCGGCCCTCGAt                                   <  8:118683/112‑1 (MQ=255)
                                                                                                                                                gcggcgCATGGGGGCACTGTCCGGCGGTGAGCGCCAGCGCGTGCTGCTGGCCCAGGGCCTGATCCCCGAGCCGCAGTTGCTGGTGCTGGACGAACCGATGTCGGCCCTCGAt                                   >  7:118683/1‑112 (MQ=255)
                                                                                                                                                   |                                                                                                                                             
TTCGACCGCGGCCTGCCCATGACCGTGGACGATTTCATGGCCGCCATGTGCCAGCGCCGCCCGGCCTTCCTCGGCCTGTCGCGGCGCGTGAAGCCGGCCATTGATGCAGCCTTGGTTCAGGTTGGCATGCTCGACAAGCGCACGCGGCGCATGGGGGCACTGTCCGGCGGTGAGCGCCAGCGCGTGCTGCTGGCCCAGGGCCTGATCCCCGAGCCGCAGTTGCTGGTGCTGGACGAACCGATGTCGGCCCTCGATGAAGCCGGAATCCAGGTGTTCGAACAGTTATTG  >  NC_002947/4330822‑4331109

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