Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 338,065 G→T 12.4% 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,0650GT12.4% 185.5 / 13.4 69A168A (GCC→GCAPP_0280amino acid ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base G (25/35);  new base T (4/5);  total (29/40)
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.97e-01

GCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCACCAGGTAGAACACGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCACGCGACATGCCCAGCGCCCGTGCCGCTTCAATCTCACCTTTAGGGATCGC  >  NC_002947/337920‑338211
                                                                                                                                                  |                                                                                                                                                  
gcagcTTGAAGCCTTGCACCAGCATGAACGAGATCACCAGGTAGAACACGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCa                                                                                                                                                 >  2:164330/1‑149 (MQ=255)
   gcTTGAAGCCTTGCACCAGCATGAACGAGATCACCAGGTAGAACACGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTg                                                                                                                                              <  4:177615/149‑1 (MQ=255)
       gAAGCCTTGCACCAGCATGAACGAGATCACCAGGTAGAACACGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGcc                                                                                                                                           <  4:311620/148‑1 (MQ=255)
       gAAGCCTTGCACCAGCATGAACGAGATCACCAGGTAGAACACGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGc                                                                                                                                            <  6:294564/147‑1 (MQ=255)
        aaGCCTTGCACCAGCATGAACGAGATCACCAGGTAGAACACGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGTGCCCTGGcca                                                                                                                                          <  4:207084/148‑2 (MQ=255)
        aaGCCTTGCACCAGCATGAACGAGATCACCAGGTAGAACACGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCtt                                                                                                                                         <  7:264191/149‑1 (MQ=255)
          gCCTTGCACCAGCATGAACGAGATCACCTGGTAGAACACGCCTGCGGGGAAGAAGATTTCCACGGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTGCCAGCAGGGTCACGGTACTGGCCAGTGCCCTGGccagta                                                                                                                                       <  5:267465/149‑5 (MQ=255)
          gCCTTGCACCAGCATGAACGAGATCACCAGGTAGAACACGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGa                                                                                                                                       <  7:252757/149‑1 (MQ=255)
          gCCTTGCACCAGCATGAACGAGATCACCAGGTAGAACACGCCGGCGGCGAAGAAGATTTCCACCGGCATGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGTGCCCTGGccagta                                                                                                                                       <  5:261434/149‑5 (MQ=255)
               gCACCAGCATGAACGAGATCACCAGGTAGAACACGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGTGCCCTGGccagtaccg                                                                                                                                    <  6:108250/147‑8 (MQ=255)
                cACCAGCATGAACGAGATCACCAGGTAGAACACGCCTGCGGCGAAGAAGATTTTCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCa                                                                                                                                 <  5:435316/149‑1 (MQ=255)
                 aCCAGCATGAACGAGATCACCAGGTAGAACACGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAg                                                                                                                                <  6:122657/149‑1 (MQ=255)
                    agcctcaaCGAGATCACCAGGTTCAACACGCCTGCGTCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGAt                                                                                                                             <  4:324450/143‑1 (MQ=255)
                    aGCATGAACGAGATCACCAGGTAGAACACGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGAt                                                                                                                             <  8:243416/149‑1 (MQ=255)
                         gAACGAGATCACCAGGTAGAACACGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGGAGCAGGATCAcgt                                                                                                                        >  8:161134/1‑147 (MQ=255)
                          aaCGAGATCACCAGGTAGAACACGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTc                                                                                                                       <  2:112729/149‑1 (MQ=255)
                           aCGAGATCACCAGTTAGAACACGCCTGCGGCGAAGAACATTTCCACCGGCAGGTAGGTGCTGGCGATGATTGTACGGGCCATGCCGGTCAGTTCCAGCAGTGTCACGGTACTGGCGAGGGCACTGGCCTTGATCATCAGGATCACTTc                                                                                                                       <  8:157063/148‑1 (MQ=255)
                            cGAGATCACCAGGTAGAACACGCCTGCGGCGAGAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCg                                                                                                                      >  1:267048/1‑149 (MQ=255)
                                 tCACCAGGTAGAACACGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTg                                                                                                                <  2:192452/149‑1 (MQ=255)
                                  cACCAGGTAGAACACGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCa                                                                                                                                    >  2:422922/1‑128 (MQ=255)
                                  cACCAGGTAGAACACGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCa                                                                                                                                    <  1:422922/128‑1 (MQ=255)
                                        gTAGAACACGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCg                                                                                                         <  6:271661/149‑1 (MQ=255)
                                        gTAGAACACGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCg                                                                                                         <  7:85452/149‑1 (MQ=255)
                                        gTAGAACACGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCg                                                                                                         >  3:300943/1‑149 (MQ=255)
                                         tAGAACACGCCTGCGGCGAAGAAGATTTCCCCCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGGCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCgg                                                                                                        <  7:384834/149‑1 (MQ=255)
                                              cacGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGc                                                                                                   <  7:317653/149‑1 (MQ=255)
                                                cGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCg                                                                                                 <  5:195599/149‑1 (MQ=255)
                                                      cggcgAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATgcgc                                                                                           <  1:210516/149‑1 (MQ=255)
                                                         cgAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATgcgcg                                                                                          >  6:142082/1‑147 (MQ=255)
                                                             agaagaTTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAgcg                                                                                     >  8:248918/1‑149 (MQ=255)
                                                             agaagaTTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAgcg                                                                                     <  4:74567/149‑1 (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                               <  2:180407/133‑1 (MQ=255)
                                                                                                                                   cACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCACGCGACATGCCCAGCGCCCGt                               >  1:180407/1‑133 (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:321116/149‑1 (MQ=255)
                                                                                                                                                aGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCACGCGACATGCCCAGCGCCCGTGCCGCTTCAATCTCACCTTTAGGGATCGc  <  6:219949/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
GCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCACCAGGTAGAACACGCCTGCGGCGAAGAAGATTTCCACCGGCAGGTAGGTGCGGGCGATGATGGTACGGGCCATGCCGGTCAGTTCCAGCAGGGTCACGGTACTGGCCAGGGCACTGGCCTTGAGCATCAGGATCACTTCGTTGCTGTACGCCGGCAGGCCGATGCGCGCAGCGCGTGGCAGCATGATGTAGAACAGCGTTTTGCCACGCGACATGCCCAGCGCCCGTGCCGCTTCAATCTCACCTTTAGGGATCGC  >  NC_002947/337920‑338211

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: