Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,406,053 C→G 7.2% A185A (GCG→GCC PP_4750 ← amino acid ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,406,0530CG 7.2% 214.7 / 4.8 69A185A (GCG→GCCPP_4750amino acid ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base C (30/34);  new base G (2/3);  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 = 9.17e-01

CACAAAGCCCTTCCACAACTCCGGCAACCATTGCATCCACATCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCA‑‑ACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTGGCTGAGGTGGGAATGGCGA  >  NC_002947/5405910‑5406198
                                                                                                                                                 |                                                                                                                                                 
cacaAAGCCCTTCCACAACTCCGGCAACCATTGCATCCACATCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCA‑‑ACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACcgccgccgc                                                                                                                                              >  2:155831/1‑149 (MQ=255)
  caAAGCCCTTCCACAACTCCGGCAACCATTGCATCCACATCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCA‑‑ACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCgg                                                                                                                                            <  6:213168/149‑1 (MQ=255)
   aaaGCCCTTCCACAACTCCGGCAACCATTGCATCCACATCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCA‑‑ACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGgt                                                                                                                                           >  5:38501/1‑149 (MQ=255)
   aaaGCCCTTCCACAACTCCGGCAACCATTGCATCCACATCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCA‑‑ACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGgt                                                                                                                                           >  7:24868/1‑149 (MQ=255)
   aaaGCCCTTCCACAACTCCGGCAACCATTGCATCCACATCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCA‑‑ACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGgt                                                                                                                                           >  6:364335/1‑149 (MQ=255)
       cccTTCCACAACTCCGGCAACCATTGCATCCACATCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCA‑‑ACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCGGCGGCGgtggtg                                                                                                                                       <  1:162905/149‑4 (MQ=255)
       cccTTCCACAACTCCGGCAACCATTGCATCCACATCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCA‑‑ACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAg                                                                                                                                       <  1:64662/149‑1 (MQ=255)
         cTTCCACAACTCCGGCAACCATTGCATCCACATCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCA‑‑ACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGAc                                                                                                                                     >  6:232707/1‑149 (MQ=255)
             cacaACTCCGGCAACCATTGCATCCACATCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCA‑‑ACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCt                                                                                                                                  >  7:51802/1‑148 (MQ=255)
                 aCTCCGGCAACCATTGCATCCACATCGTCAGCACTTCCCTCATCGGGAAATCCTCTCTTTCAACGTGGCCGACAGCCA‑‑ACGAGAGAACAAGGCCACTTGCACGCAGAGCAACACGTACACCAACGCCGCCGCGGTGTAGAAGCTCAGGc                                                                                                                             >  7:84794/1‑149 (MQ=255)
                   tCCGGCAACCATTGCATCCACATCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCA‑‑ACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGccc                                                                                                                           >  1:15710/1‑149 (MQ=255)
                     cGGCAACCATTGCATCCACATCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCA‑‑ACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGccc                                                                                                                           <  1:30706/147‑1 (MQ=255)
                     cGGCAACCATTGCATCCACATCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCA‑‑ACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCt                                                                                                                          <  8:51802/148‑1 (MQ=255)
                         aaCCATTGCATCCACATCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCA‑‑ACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCGGCGGCGgtggtgta                                                                                                                                     <  4:8170/133‑6 (MQ=255)
                         aaCCATTGCATCCACATCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCA‑‑ACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCGGCGGCGgtggtgta                                                                                                                                     >  3:8170/1‑128 (MQ=255)
                         aaCCATTGCATCCACATCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCA‑‑ACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAg                                                                                                                     >  7:155904/1‑149 (MQ=255)
                         aaCCATTGCATCCACATCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCA‑‑ACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAg                                                                                                                     <  5:364335/149‑1 (MQ=255)
                            cATTGCATCCACATCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCA‑‑ACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGtt                                                                                                                  <  2:145078/149‑1 (MQ=255)
                            cATTGCATCCACATCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCA‑‑ACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGtt                                                                                                                  >  1:145078/1‑149 (MQ=255)
                                 cATCCACATCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCCACGTGTCCGACAGCCA‑‑ACGAGAGTACAAGGGCACTGGCACTGTGAGCAGCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCt                                                                                                              <  5:266037/148‑1 (MQ=255)
                                  aTCCACATCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCA‑‑ACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCtgt                                                                                                            <  2:15710/149‑1 (MQ=255)
                                   tCCACATCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCAAGACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCtg                                                                                                             >  1:299104/1‑149 (MQ=255)
                                    ccACATCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCA‑‑ACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCtgtg                                                                                                           >  6:58088/1‑148 (MQ=255)
                                     cacaTCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCA‑‑ACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGa                                                                                                          >  1:146074/1‑148 (MQ=255)
                                                cagTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCA‑‑ACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGc                                                                                              <  5:330180/149‑1 (MQ=255)
                                                   ttCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCA‑‑ACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCACCGCCGCGGTGTAGACGCTCAGGCCCt                                                                                                                          >  4:174955/1‑118 (MQ=255)
                                                   ttCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCA‑‑ACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCACCGCCGCGGTGTAGACGCTCAGGCCCt                                                                                                                          <  3:174955/118‑1 (MQ=255)
                                                      gCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCA‑‑ACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTg                                                                                                    <  8:257285/137‑1 (MQ=255)
                                                      gCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCA‑‑ACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTg                                                                                                    >  7:257285/1‑137 (MQ=255)
                                                      gCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCA‑‑ACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGcc                                                                                             >  4:201202/1‑144 (MQ=255)
                                                      gCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCA‑‑ACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGcc                                                                                             <  3:201202/144‑1 (MQ=255)
                                                      gCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCA‑‑ACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGaaa                                                                                        <  5:363939/149‑1 (MQ=255)
                                                       cTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCA‑‑ACGAGAGAACAAGGCCACTGTCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAg                                                                                       <  5:232707/149‑1 (MQ=255)
                                                          aTCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCA‑‑ACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCt                                                                                     <  8:24868/148‑1 (MQ=255)
                                                                cAATCCTCTGTTTCAACGTGGCCGACAGCCA‑‑ACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCa                                                                              <  5:58088/149‑1 (MQ=255)
                                                                     ctctGTTTCAACGTGGCCGACAGCCA‑‑ACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCaa                                                                          <  1:155831/148‑1 (MQ=255)
                                                                     ctctGTTTCAACGTGGCCGACAGCCA‑‑ACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAAc                                                                         <  6:38501/149‑1 (MQ=255)
                                                                         gTTTCAACGTGGCCGACAGCCA‑‑ACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGcc                                                                      >  3:192115/1‑148 (MQ=255)
                                                                          tttCAACGTGGCCGACAGCCA‑‑ACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGc                                                                       <  3:9492/146‑1 (MQ=255)
                                                                                  tcgCCGACAGCCA‑‑ACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCTCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAg                                                            <  8:123202/147‑1 (MQ=255)
                                                                                            ccA‑‑ACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGcc                                                   >  7:61281/1‑148 (MQ=255)
                                                                                            ccA‑‑ACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGcc                                                   >  8:303039/1‑148 (MQ=255)
                                                                                             cA‑‑ACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACggac                                                   >  7:43453/1‑145 (MQ=255)
                                                                                              a‑‑aCGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAg                                                <  6:77906/149‑1 (MQ=255)
                                                                                                 aCGAGAGAAGAAGGCGACTTGCACGCTTAGCACGAGGTACACCACCGCCGCCGCGGTGTAGACGCTATGGCCCTGTAAGGTTTGCTGTGAAACCTTGTAAGCCTGATTGCTTATTTCACCAACTCCAATCGTCGAGTGCACGGCCGatt                                               <  8:84794/149‑3 (MQ=255)
                                                                                                       gaACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGa                                         <  8:61281/149‑1 (MQ=255)
                                                                                                       gaACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCCGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGa                                         <  2:130699/149‑1 (MQ=255)
                                                                                                        aaCAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGcc                                                         <  8:178956/132‑1 (MQ=255)
                                                                                                        aaCAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGcc                                                         >  7:178956/1‑132 (MQ=255)
                                                                                                                cACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCaa                                >  7:214558/1‑149 (MQ=255)
                                                                                                                       aCGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTa                         <  7:303039/149‑1 (MQ=255)
                                                                                                                         gCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAgg                       >  5:220966/1‑149 (MQ=255)
                                                                                                                           tGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGAAGGTg                     >  6:152390/1‑149 (MQ=255)
                                                                                                                             aGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTgg                    <  8:155904/148‑1 (MQ=255)
                                                                                                                               cACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTGGCt                  <  6:220966/148‑1 (MQ=255)
                                                                                                                                  cAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGa                                                                                   <  2:181206/80‑1 (MQ=255)
                                                                                                                                  cAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGa                                                                                   >  1:181206/1‑80 (MQ=255)
                                                                                                                                  cAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTGGCTGAgg              >  5:65847/1‑149 (MQ=255)
                                                                                                                                     gTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTGGCTGAGGTgg           >  7:354596/1‑149 (MQ=255)
                                                                                                                                      tACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTGGCTGAGGTggg          <  8:354596/149‑1 (MQ=255)
                                                                                                                                       tctacaccGCGGCGGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTGGCTGAGGTggg          >  6:158199/5‑148 (MQ=255)
                                                                                                                                         accaccGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTgg                    <  5:179819/136‑1 (MQ=255)
                                                                                                                                         accaccGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTgg                    >  6:179819/1‑136 (MQ=255)
                                                                                                                                         accaccGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTAAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTGGCTGAGGTGGGAAt       >  3:229390/1‑149 (MQ=255)
                                                                                                                                           caccGCGGCGGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTGGCTGAGGTGGGAATgg     <  5:158199/149‑1 (MQ=255)
                                                                                                                                           caccGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCtt                                           <  5:45563/111‑1 (MQ=255)
                                                                                                                                           caccGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCtt                                           >  6:45563/1‑111 (MQ=255)
                                                                                                                                           caccGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTGGCTGAGGTGGGAATgg     >  4:119567/1‑149 (MQ=255)
                                                                                                                                           caccGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTGGCTGAGGTGGGAATgg     <  2:146074/149‑1 (MQ=255)
                                                                                                                                           caccGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTGGCTGAGGTGGGAATgg     <  4:192115/149‑1 (MQ=255)
                                                                                                                                              cgccgccgcGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTGGCTGAGGTGGGAATGGCGa  >  2:73506/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
CACAAAGCCCTTCCACAACTCCGGCAACCATTGCATCCACATCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCA‑‑ACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTGGCTGAGGTGGGAATGGCGA  >  NC_002947/5405910‑5406198

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: