Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 338,065 G→T 7.5% A168A (GCC→GCA)  PP_0280 ← amino acid ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947338,0650GT 7.5% 209.0 / 4.9 66A168A (GCC→GCA) PP_0280amino acid ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base G (24/37);  new base T (3/2);  total (27/39)
Fisher's exact test for biased strand distribution p-value = 3.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCACCAGGTAGAACACGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCACGCGACATGCCCAGCGCCCGTGCCGCTTCAATCTCACCTTTAGGGAT  >  NC_002947/337919‑338208
                                                                                                                                                  |                                                                                                                                               
agcagcTTGAAGCCTTGCACCAGCATGAACGAGATCACCAGGTAGAACACGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGTGc                                                                                                                                               >  6:266268/1‑149 (MQ=255)
         aaGCCTTGCACCAGCATGAACGAGATCACCAGGTAGAACACGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCtt                                                                                                                                      <  8:50673/149‑1 (MQ=255)
         aaGCCTTGCACCAGCATGAACGAGATCACCAGGTAGAACACGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCtt                                                                                                                                      <  6:90019/149‑1 (MQ=255)
         aaGCCTTGCACCAGCATGAACGAGATCACCAGGTAGAACACGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGACCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCtt                                                                                                                                      <  7:275548/149‑1 (MQ=255)
           gCCTTGCACCAGCATGAACGAGATCACCAGGTAGAACACGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGTGCCCTGGccagta                                                                                                                                    <  5:266268/149‑5 (MQ=255)
             cTTGCACCAGCATGAACGAGATCACCAGGTAGAACACGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGc                                                                                                                                  <  4:41682/149‑1 (MQ=255)
                 caccacCATGATCGAGATCACCAGGTAGAACACGCCTGCGGCGAAGAAGATTTCGACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATc                                                                                                                               <  6:4319/148‑1 (MQ=255)
                 cACCAGCATGAACGAGATCACCAGGTAGAACACGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATc                                                                                                                               <  3:30075/148‑1 (MQ=255)
                   ccAGCATGAACGAGATCACCAGGTAGAACACGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAgg                                                                                                                            <  7:227146/149‑1 (MQ=255)
                   ccAGCATGAACGAGATCACCAGGTAGAACACGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAgg                                                                                                                            <  3:22846/149‑1 (MQ=255)
                       cATGAACGAGATCACCAGGTAGAACACGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCa                                                                                                                        >  4:50288/1‑149 (MQ=255)
                            aCGAGATCACCAGGTAGAACACGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGTGCGATGATGGTACGGGCCATGCCGGGCAGTTCCATCAGGGTCACGGTGCTGGCCATGGCACTGGCCTTGAGCATCAGGATCACTTCg                                                                                                                   <  2:45208/149‑1 (MQ=255)
                             cGAGATCACCAGGTAGAACACGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGt                                                                                                                  <  4:61369/149‑1 (MQ=255)
                              gagaTCACCAGGTAGAACACGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGCCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGtt                                                                                                                 >  8:265046/1‑149 (MQ=255)
                               agaTCACCAGGTAGAACACGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTg                                                                                                                <  4:269499/149‑1 (MQ=255)
                                gaTCACCAGGTAGAACACGCCTGCGGCGAAGAAGATTTCCACCGGCTGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGATCCAGCAGGCTCACGGTACAGGCCAGGGCACTGGCCTAGAGCATCAGGATCACTTCGTTGc                                                                                                               >  7:73117/1‑149 (MQ=255)
                                gaTCACCAGGTAGAACACGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGACACTGGCCTTGAGCATCAGGATCACTTCGTTGc                                                                                                               >  6:228101/1‑149 (MQ=255)
                                  tCACCAGGTAGAACACGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTg                                                                                                             <  2:137737/149‑1 (MQ=255)
                                    aCCAGGTAGAACACGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCa                                                                                                                        <  8:20067/136‑1 (MQ=255)
                                    aCCAGGTAGAACACGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCa                                                                                                                        >  7:20067/1‑136 (MQ=255)
                                    aCCAGGTAGAACACGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTa                                                                                                           >  6:299714/1‑149 (MQ=255)
                                      cAGGTAGAACACGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACg                                                                                                         >  6:20019/1‑149 (MQ=255)
                                        ggTAGAACACGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGcc                                                                                                       <  3:242999/149‑1 (MQ=255)
                                         gTAGAACACGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCg                                                                                                      >  7:129234/1‑149 (MQ=255)
                                            gAACACGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCgg                                                                                                     >  3:82369/1‑147 (MQ=255)
                                               gacGCCTGCGGCGAAGTAGATTTCCCCCGGCAGGTAGGTGGGGGCGATGATGGTACGGCCCTTGCCGGTCAGTTCCAGGAGGGTCACGGTACTGGCCTGGGCTCTGGCCTTGAGCTTCAGGATCACTTCGTTGCTGTTCGCCGGCAGGc                                                                                                <  8:73117/148‑1 (MQ=255)
                                               cacGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCTGCAGGc                                                                                                <  6:147862/149‑1 (MQ=255)
                                               cacGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTAGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGc                                                                                                >  6:96516/1‑149 (MQ=255)
                                                    cTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATg                                                                                           <  1:5820/149‑1 (MQ=255)
                                                       cggcgAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATgcgc                                                                                        >  8:135995/1‑149 (MQ=255)
                                                             agaagaTTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAgcg                                                                                  >  8:225791/1‑149 (MQ=255)
                                                                 gaTTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTg                                                                              >  4:196412/1‑149 (MQ=255)
                                                                        aCCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGgcagc                                                                         >  7:206874/1‑147 (MQ=255)
                                                                        aCCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGgcagc                                                                         <  8:206874/147‑1 (MQ=255)
                                                                                gtaggtGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAg                                                                >  5:142911/1‑148 (MQ=255)
                                                                                 taggtGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTTGCCTTGAGCTTCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGaa                                                              <  5:96516/149‑1 (MQ=255)
                                                                                   ggtGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAAc                                                             <  2:51666/148‑1 (MQ=255)
                                                                                         ggCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGttt                                                      <  3:196412/149‑1 (MQ=255)
                                                                                           cGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTg                                                    <  5:17286/149‑1 (MQ=255)
                                                                                            gatgatGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGCCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGc                                                   <  7:265046/149‑1 (MQ=255)
                                                                                             atgatgGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGTAGCATGATGTAGAACAGCGTTTTGcc                                                  <  8:129234/149‑1 (MQ=255)
                                                                                               gatgGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCAc                                                >  4:189927/1‑149 (MQ=255)
                                                                                                    tACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCAGTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCACGCGa                                            >  7:9161/1‑148 (MQ=255)
                                                                                                      cGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCACGCGACAt                                         <  3:189927/149‑1 (MQ=255)
                                                                                                         gCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGtt                                                       >  5:268319/1‑132 (MQ=255)
                                                                                                         gCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGtt                                                       <  6:268319/132‑1 (MQ=255)
                                                                                                         gCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCACGCGACATGcc                                      <  1:242140/149‑1 (MQ=255)
                                                                                                                cGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCAcg                                               >  7:173300/1‑133 (MQ=255)
                                                                                                                cGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCAcg                                               <  8:173300/133‑1 (MQ=255)
                                                                                                                     aGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCACGCGACATGCCCAGCGCCCGTGc                          <  3:50288/149‑1 (MQ=255)
                                                                                                                     aGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCACGCGACATGCCCAGCGCCCGTGc                          <  7:225791/149‑1 (MQ=255)
                                                                                                                            gcagGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCACGCGACATGCCCAGCGCCCGTGCCGCTTCa                   <  5:299714/149‑1 (MQ=255)
                                                                                                                             cagGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCACGCGACATGCCCAGCGCCCGTGCCGCTTCaa                  <  7:85745/149‑1 (MQ=255)
                                                                                                                             cagGGGCACGGTACTGGCCAGGGCACTGGCCGTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCACGCGACATGCCCAGCGCCCGTGCCGCTTCaa                  <  7:146363/149‑1 (MQ=255)
                                                                                                                                 gTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCACGCGACATGCCCAGCGCCCGTGCCGCTTCa                   <  7:229122/144‑1 (MQ=255)
                                                                                                                                 gTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCACGCGACATGCCCAGCGCCCGTGCCGCTTCa                   >  8:229122/1‑144 (MQ=255)
                                                                                                                                  tCACGGTTCTGGCCAGGGCACTGGCCTTGATCATCAGGATCACTTCTTTGCTGTACGCCGGGAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTATAACAGCGTTTTGCCACGCGACATGCCCAGCGCCCGTGCCGCTTCAATCTCa             <  8:299611/149‑1 (MQ=255)
                                                                                                                                  tCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGTCACGCGACATGCCCAGCGCCCGTGCCGCTTCAATCTCa             <  4:82369/149‑1 (MQ=255)
                                                                                                                                   cACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGt                                                                  <  7:279340/95‑1 (MQ=255)
                                                                                                                                   cACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGt                                                                  >  8:279340/1‑95 (MQ=255)
                                                                                                                                    aCGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCACGCGACATGCCCAGCGCCCGTGCCGCTTCAATCTCAcc           >  4:1347/1‑149 (MQ=255)
                                                                                                                                     cGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCACGCGACATGCCCAGCGCCCGTGCCGCTTCAATCTCACCt          >  5:92221/1‑149 (MQ=255)
                                                                                                                                       ctcaaggCCAGTGCCCTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCACGCGACATGCCCAGCGCCCGTGCCGCTTCAATCTCACCtt         >  2:120586/6‑148 (MQ=255)
                                                                                                                                          aaggCCAGTGCCCTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCACGCGACATGCCCAGCGCCCGTGCCGCTTCAATCTCACCTTTAgg     >  4:163708/3‑149 (MQ=255)
                                                                                                                                           aggCCAGTGCCCTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCACGCGACATGCCCAGCGCCCGTGCCGCTTCAATCTCACCTTTAggg    <  1:120586/148‑1 (MQ=255)
                                                                                                                                             gCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGACCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCACGCGACATGCCCAGCGCCCGTGCCGCTTCAATCTCACCTTTAGGGAt  >  1:84188/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
AGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCACCAGGTAGAACACGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCACGCGACATGCCCAGCGCCCGTGCCGCTTCAATCTCACCTTTAGGGAT  >  NC_002947/337919‑338208

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