Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 338,068 A→C 12.5% S167R (AGT→AGG)  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,0680AC12.5% 191.3 / 17.7 71S167R (AGT→AGG) PP_0280amino acid ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base A (26/36);  new base C (4/5);  total (30/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCACCAGGTAGAACACGCCTGCGGCGA‑AGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCACGCGACATGCCCAGCGCCCGTGCCGCTTCAATCTCACCTTTAGGGATCGCCTGCA  >  NC_002947/337920‑338216
                                                                                                                                                     |                                                                                                                                                    
gcagcTTGAAGCCTTGCACCAGCATGAACGAGATCACCAGGTAGAACACGCCTGCGGCGA‑AGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCa                                                                                                                                                      >  2:164330/1‑149 (MQ=255)
   gcTTGAAGCCTTGCACCAGCATGAACGAGATCACCAGGTAGAACACGCCTGCGGCGA‑AGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTg                                                                                                                                                   <  4:177615/149‑1 (MQ=255)
       gAAGCCTTGCACCAGCATGAACGAGATCACCAGGTAGAACACGCCTGCGGCGA‑AGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGcc                                                                                                                                                <  4:311620/148‑1 (MQ=255)
       gAAGCCTTGCACCAGCATGAACGAGATCACCAGGTAGAACACGCCTGCGGCGA‑AGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGc                                                                                                                                                 <  6:294564/147‑1 (MQ=255)
        aaGCCTTGCACCAGCATGAACGAGATCACCAGGTAGAACACGCCTGCGGCGA‑AGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGTGCCCTGGcca                                                                                                                                               <  4:207084/148‑2 (MQ=255)
        aaGCCTTGCACCAGCATGAACGAGATCACCAGGTAGAACACGCCTGCGGCGA‑AGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCtt                                                                                                                                              <  7:264191/149‑1 (MQ=255)
          gCCTTGCACCAGCATGAACGAGATCACCTGGTAGAACACGCCTGCGGGGA‑AGAAGATTTCCACGGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTGCCAGCAGGGTCACGGTACTGGCCAGTGCCCTGGccagta                                                                                                                                            <  5:267465/149‑5 (MQ=255)
          gCCTTGCACCAGCATGAACGAGATCACCAGGTAGAACACGCCTGCGGCGA‑AGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGa                                                                                                                                            <  7:252757/149‑1 (MQ=255)
          gCCTTGCACCAGCATGAACGAGATCACCAGGTAGAACACGCCGGCGGCGA‑AGAAGATTTCCACCGGCATGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGTGCCCTGGccagta                                                                                                                                            <  5:261434/149‑5 (MQ=255)
               gCACCAGCATGAACGAGATCACCAGGTAGAACACGCCTGCGGCGA‑AGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGTGCCCTGGccagtaccg                                                                                                                                         <  6:108250/147‑8 (MQ=255)
                cACCAGCATGAACGAGATCACCAGGTAGAACACGCCTGCGGCGA‑AGAAGATTTTCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCa                                                                                                                                      <  5:435316/149‑1 (MQ=255)
                 aCCAGCATGAACGAGATCACCAGGTAGAACACGCCTGCGGCGA‑AGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAg                                                                                                                                     <  6:122657/149‑1 (MQ=255)
                    agcctcaaCGAGATCACCAGGTTCAACACGCCTGCGTCGA‑AGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGAt                                                                                                                                  <  4:324450/143‑1 (MQ=255)
                    aGCATGAACGAGATCACCAGGTAGAACACGCCTGCGGCGA‑AGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGAt                                                                                                                                  <  8:243416/149‑1 (MQ=255)
                         gAACGAGATCACCAGGTAGAACACGCCTGCGGCGA‑AGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGGAGCAGGATCAcgt                                                                                                                             >  8:161134/1‑147 (MQ=255)
                          aaCGAGATCACCAGGTAGAACACGCCTGCGGCGA‑AGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTc                                                                                                                            <  2:112729/149‑1 (MQ=255)
                           aCGAGATCACCAGTTAGAACACGCCTGCGGCGA‑AGAACATTTCCACCGGCAGGTAGGTGCTGGCGATGATTGTACGGGCCATGCCGGTCAGTTCCAGCAGTGTCACGGTACTGGCGAGGGCACTGGCCTTGATCATCAGGATCACTTc                                                                                                                            <  8:157063/148‑1 (MQ=255)
                            cGAGATCACCAGGTAGAACACGCCTGCGGCGAGAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCg                                                                                                                           >  1:267048/1‑149 (MQ=255)
                                 tCACCAGGTAGAACACGCCTGCGGCGA‑AGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTg                                                                                                                     <  2:192452/149‑1 (MQ=255)
                                  cACCAGGTAGAACACGCCTGCGGCGA‑AGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCa                                                                                                                                         >  2:422922/1‑128 (MQ=255)
                                  cACCAGGTAGAACACGCCTGCGGCGA‑AGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCa                                                                                                                                         <  1:422922/128‑1 (MQ=255)
                                        gTAGAACACGCCTGCGGCGA‑AGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCg                                                                                                              <  7:85452/149‑1 (MQ=255)
                                        gTAGAACACGCCTGCGGCGA‑AGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCg                                                                                                              >  3:300943/1‑149 (MQ=255)
                                        gTAGAACACGCCTGCGGCGA‑AGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCg                                                                                                              <  6:271661/149‑1 (MQ=255)
                                         tAGAACACGCCTGCGGCGA‑AGAAGATTTCCCCCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGGCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCgg                                                                                                             <  7:384834/149‑1 (MQ=255)
                                              cacGCCTGCGGCGA‑AGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGc                                                                                                        <  7:317653/149‑1 (MQ=255)
                                                cGCCTGCGGCGA‑AGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCg                                                                                                      <  5:195599/149‑1 (MQ=255)
                                                      cggcgA‑AGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATgcgc                                                                                                <  1:210516/149‑1 (MQ=255)
                                                         cgA‑AGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATgcgcg                                                                                               >  6:142082/1‑147 (MQ=255)
                                                             agaagaTTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAgcg                                                                                          <  4:74567/149‑1 (MQ=255)
                                                             agaagaTTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAgcg                                                                                          >  8:248918/1‑149 (MQ=255)
                                                                                 taggtGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGaa                                                                      <  3:62034/149‑1 (MQ=255)
                                                                                  aggtGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAAc                                                                     >  4:285840/1‑149 (MQ=255)
                                                                                   ggtGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAAc                                                                     <  2:267048/148‑1 (MQ=255)
                                                                                      gCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGCACAGCg                                                                 >  4:80831/1‑149 (MQ=255)
                                                                                           cGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTTTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTg                                                            <  4:416012/149‑1 (MQ=255)
                                                                                             atgatgGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGcc                                                          <  4:349601/149‑1 (MQ=255)
                                                                                                  ggTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCAcgc                                                      >  7:58563/1‑148 (MQ=255)
                                                                                                         gCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCACGCGACATGcc                                              >  8:366723/1‑149 (MQ=255)
                                                                                                         gCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCACGCGACATGcc                                              >  1:207449/1‑149 (MQ=255)
                                                                                                          ccATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCACGCGACATGcc                                              >  6:278951/1‑148 (MQ=255)
                                                                                                              gCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGtttt                                                             >  5:361121/1‑129 (MQ=255)
                                                                                                              gCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGtttt                                                             <  6:361121/129‑1 (MQ=255)
                                                                                                                cGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCACGCGACATGCCCAGCGcc                                       <  4:404893/149‑1 (MQ=255)
                                                                                                                cGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTAGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCATCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCACGCGACATCCCCAGCGcc                                       >  1:305517/1‑149 (MQ=255)
                                                                                                                 ggTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCACGCGACATGCCCAGCGccc                                      >  6:291007/1‑149 (MQ=255)
                                                                                                                 ggTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCACGCGACATGCCCAGCGccc                                      >  6:297193/1‑149 (MQ=255)
                                                                                                                         ccAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCACGCGACATGCCCAGCGCCCGTGCCGc                               <  7:248918/148‑1 (MQ=255)
                                                                                                                         ccAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCACGCGACATGCCCAGCGCCCGTGCCGCt                              >  1:73071/1‑149 (MQ=255)
                                                                                                                          cagcagGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCACGCGACATGCCCAGCGCCCGTGCCGCtt                             <  2:207449/149‑1 (MQ=255)
                                                                                                                          cagcagGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCACGCGACATGCCCAGCGCCCGTGCCGCtt                             >  8:278361/1‑149 (MQ=255)
                                                                                                                              agGGGCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGCTCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCACGCGACATGCCCAGCGCCCGTGCCGCTTCAAt                         <  7:161134/149‑1 (MQ=255)
                                                                                                                                 gTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCACGCGACATGCCCAGCGCCCGTGCCGCTTCAAtctc                      <  1:164330/149‑1 (MQ=255)
                                                                                                                                  tCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCACGCGACATGCCCAGCGCCCGTGCCGCTTCAATCTCa                     >  5:321116/1‑149 (MQ=255)
                                                                                                                                  tCACGGTACTGGCCAGGGCACAGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCAGATGCGCCCAGCGCGTGGCAGCATGATGTAGAACAGCGATATGCCACGCGACATGCCCAGCCACCGTGCCGCTTCAATCTCa                     >  5:184505/1‑149 (MQ=255)
                                                                                                                                   cACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCACGCGACATGCCCAGCGCCCGt                                    >  1:180407/1‑133 (MQ=255)
                                                                                                                                   cACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCACGCGACATGCCCAGCGCCCGt                                    <  2:180407/133‑1 (MQ=255)
                                                                                                                                   cACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCACGCGACATGCCCAGCGCCCGTGCCGCTTCAATCTCAc                    >  8:280026/1‑149 (MQ=255)
                                                                                                                                       ctcaaggCCAGTGCCCTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCACGCGACATGCCCAGCGCCCGTGCCTCTTCAATCTCACCtt                 >  7:185368/6‑148 (MQ=255)
                                                                                                                                        tcaaggCCAGTGCCCTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCACGCGACATGCCCAGCGCCCGTGCCGCTTCAATCTCACCtt                 >  7:205711/5‑147 (MQ=255)
                                                                                                                                        tcaaggCCAGTGCCCTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCACGCGACATGCCCAGCGCCCGTGCCGCTTCAATCTCACCtt                 >  1:37283/5‑147 (MQ=255)
                                                                                                                                         caaggCCAGTGCCCTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCACGCGACATGCCCAGCGCCCGTGCCGCTTCAATCTCACCTTTAg              >  3:428511/4‑149 (MQ=255)
                                                                                                                                         aCTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCACGCGACATGCCCAGCGCCCGTGCCGCTTCAATCTCACCTTTAg              >  3:155702/1‑149 (MQ=255)
                                                                                                                                          cTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCACGCGACATGCCCAGCGCCCGTGCCGCTTCAATCTCACCTTTAg              <  5:142082/148‑1 (MQ=255)
                                                                                                                                          aaggCCAGTGCCCTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCACGCGACATGCCCAGCGCCCGTGCCGCTTCAATCTCACCTTTa               <  2:37283/145‑1 (MQ=255)
                                                                                                                                             gCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCatg                                                                              <  2:273234/81‑1 (MQ=255)
                                                                                                                                             gCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCatg                                                                              >  1:273234/1‑81 (MQ=255)
                                                                                                                                                aGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCACGCGACATGCCCAGCGCCCGTGCCGCTTCAATCTCACCTTTAGGGATCGc       <  6:219949/149‑1 (MQ=255)
                                                                                                                                                aGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCACGCGACATGCCCAGCGCCCGTGCCGCTTCAATCTCACCTTTAGGGATCGc       <  6:321116/149‑1 (MQ=255)
                                                                                                                                                   gCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCACGCGACATGCCCAGCGCCCGTGCCGCTTCAATCTCACCTTTAGGGATCGCCTg    >  4:86103/1‑149 (MQ=255)
                                                                                                                                                   gCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCACGCGACATGCCCAGCGCCCGTGCCGCTTCAATCTCACCTTTAGGGATCGCCTg    <  7:366723/149‑1 (MQ=255)
                                                                                                                                                   gCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCACGCGACATGCCCAGCGCCCGTGCCGCTTCAATCTCACCTTTAGGGATCGCCTg    >  1:73186/1‑149 (MQ=255)
                                                                                                                                                     aCTGGCCTTGAGCATCAGGATGACTTCGTTGCTTTACGGCGGCATGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCACGCGACATGCCCAGCGCCCGTGCCGCTTCAATCTCACCTTTAGGGATCGCCTGCa  <  2:305517/149‑1 (MQ=255)
                                                                                                                                                     aCTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCACGCGACATGCCCAGCGCCCGTGCCGCTTCAATCTCACCTTTAGGGATCGCCTg    >  7:190189/1‑147 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCACCAGGTAGAACACGCCTGCGGCGA‑AGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCACGCGACATGCCCAGCGCCCGTGCCGCTTCAATCTCACCTTTAGGGATCGCCTGCA  >  NC_002947/337920‑338216

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