Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,406,056 C→G 7.2% A184A (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,0560CG 7.2% 208.7 / 5.7 69A184A (GCG→GCC) PP_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 = 8.59e-01

CACAAAGCCCTTCCACAACTCCGGCAACCATTGCATCCACATCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCA‑‑ACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTGGCTGAGGTGGGAATGGCGAT  >  NC_002947/5405910‑5406199
                                                                                                                                                    |                                                                                                                                               
cacaAAGCCCTTCCACAACTCCGGCAACCATTGCATCCACATCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCA‑‑ACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACcgccgccgc                                                                                                                                               >  2:155831/1‑149 (MQ=255)
  caAAGCCCTTCCACAACTCCGGCAACCATTGCATCCACATCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCA‑‑ACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCgg                                                                                                                                             <  6:213168/149‑1 (MQ=255)
   aaaGCCCTTCCACAACTCCGGCAACCATTGCATCCACATCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCA‑‑ACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGgt                                                                                                                                            >  6:364335/1‑149 (MQ=255)
   aaaGCCCTTCCACAACTCCGGCAACCATTGCATCCACATCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCA‑‑ACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGgt                                                                                                                                            >  7:24868/1‑149 (MQ=255)
   aaaGCCCTTCCACAACTCCGGCAACCATTGCATCCACATCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCA‑‑ACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGgt                                                                                                                                            >  5:38501/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                                                                                                                                      >  3:8170/1‑128 (MQ=255)
                         aaCCATTGCATCCACATCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCA‑‑ACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCGGCGGCGgtggtgta                                                                                                                                      <  4:8170/133‑6 (MQ=255)
                         aaCCATTGCATCCACATCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCA‑‑ACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAg                                                                                                                      <  5:364335/149‑1 (MQ=255)
                         aaCCATTGCATCCACATCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCA‑‑ACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAg                                                                                                                      >  7:155904/1‑149 (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                                                    >  8:303039/1‑148 (MQ=255)
                                                                                            ccA‑‑ACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGcc                                                    >  7:61281/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                                            >  6:45563/1‑111 (MQ=255)
                                                                                                                                           caccGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCtt                                            <  5:45563/111‑1 (MQ=255)
                                                                                                                                           caccGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTGGCTGAGGTGGGAATgg      <  2:146074/149‑1 (MQ=255)
                                                                                                                                           caccGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTGGCTGAGGTGGGAATgg      >  4:119567/1‑149 (MQ=255)
                                                                                                                                           caccGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTGGCTGAGGTGGGAATgg      <  4:192115/149‑1 (MQ=255)
                                                                                                                                              cgccgccgcGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTGGCTGAGGTGGGAATGGCGa   >  2:73506/1‑149 (MQ=255)
                                                                                                                                               gaatccgcGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTGGCTGAGGTGGGAATGGCGAt  <  1:222746/145‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
CACAAAGCCCTTCCACAACTCCGGCAACCATTGCATCCACATCGTCAGCAGTTCGCTCATCGGGCAATCCTCTGTTTCAACGTGGCCGACAGCCA‑‑ACGAGAGAACAAGGCCACTGGCACGCTGAGCACCAGGTACACCACCGCCGCCGCGGTGTAGACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTGGCTGAGGTGGGAATGGCGAT  >  NC_002947/5405910‑5406199

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