Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,954,932 A→G 5.0% V24A (GTG→GCG)  PP_5219 ← membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,954,9320AG 5.0% 234.4 / 2.3 79V24A (GTG→GCG) PP_5219membrane protein
Reads supporting (aligned to +/- strand):  ref base A (37/38);  new base G (2/2);  total (39/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.74e-01

CAATTGCAGCAGCGCCAGGTGCTCGAAGGCCAGAGCATAGGTGACCAGCCCTACGCCTTTGGTCAGCAGGCCGGTACTGGTCGACAGTGCGCCGGTACTCAGGCTGTTGATCGAGTCGGCCAGGGTGAAATTGCGCTGGCCGCGCACGCGGTCGGCGACCAGCTCTACGGCTATCAGCACAAAAAAGAATGGTACGGCCAGCAGTATCAGGTCCATGGGCGCGACCTCAGGCGTTATCCACCAAGATTAGGCCGCACCCACAGGAATCCCTATGGCGACATCTGCCAAACTAGA  >  NC_002947/5954787‑5955080
                                                                                                                                                 |                                                                                                                                                    
cAATTGCAGCAGCGCCAGGTGCTCGAAGGCCAGAGCATAGGTGACCAGCCCTACGCCTTTGGTCAGCAGGCCGGTACTGGTCGACAGTGCGCCGGTACTCAGGCTGTTGATCGAGTCGGCCAGGGTGAAATTGCGCTGGCCGCGCAcgc                                                                                                                                                   <  5:135530/149‑1 (MQ=255)
   ttGCAGCAGCGCCAGGTGCTCGAAGGCCTGAGCATAGGTGACCAGCCCTACGCCTTTGGTCAGCAGGCCGGTACTGGTCGACAGTGCGCCGGTACTCAGGCTGTTGATCGAGTCGGCCAGGGTGAAATTGCGCTGGCCGCGCACGCGGt                                                                                                                                                >  1:243653/1‑149 (MQ=255)
     gcagcaGCGCCAGGTGCTCGAAGGCCAGAGCATAGGTGACCAGCCCTACGCCTTTGGTCAGCAGGCCGGTACTGGTCGACAGTGCGCCGGTACTCAGGCTGTTGATCGAGTCGGCCAGGGTGAAATTGCGCTGGCCGCGCACGCGGTc                                                                                                                                               <  7:106918/148‑1 (MQ=255)
      cagcagTGCCAGGTGCTCGAAGGCCAGAGCATAGGTGACCAGCCCTACGCCTTTGGTCAGCAGGCCGGTACTGGTCGACAGTGCGCCGGTACTCAGGCTGTTGATCGAGTCGGCCAGGGTGAAATTGCGCTGGCCGCGCACGCGGTCgg                                                                                                                                             >  2:157823/1‑149 (MQ=255)
       agcagcGCCAGGTGCTCGAAGGCCAGAGCATAGGTGACCAGCCCTACGCCTTTGGTCAGCAGGCCGGTACTGGTCGACAGTGCGCCGGTACTCAGGCTGTTGATCGAGTCGGCCAGGGTGAAATTGCGCTGGCCGCGCACGCGGTCGGc                                                                                                                                            >  3:138777/1‑149 (MQ=255)
           gcgcCAGGTGCTCGAAGGCCAGAGCATAGGTGACCAGCCCTACGCCTTTGGTCAGCAGGCCGGTACTGGTCGACAGTGCGCCGGTACTCAGGCTGTTGATCGAGTCGGCCAGGGTGAAATTGCGCTGGCCGCGCACGCGGTCGGCGAcc                                                                                                                                        <  1:275157/149‑1 (MQ=255)
            cgcCAGGTGCTCGAAGGCCAGAGCATAGGTGACCAGCCCTACGCCTTTGGTCAGCAGGCCGGTACTGGTCGACAGTGCGCCGGTACTCAGGCTGTTGATCGAGTCGGCCAGGGTGAAATTGCGCTGGCCGCGCACGCGGTCGGCGAc                                                                                                                                         >  7:250674/1‑147 (MQ=255)
               cAGGTGCTCGAAGGCCAGAGCATAGGTGACCAGCCCTACGCCTTTGGTCAGCAGGCCGGTACTGGTCGACAGTGCGCCGGTACTCAGGCTGTTGATCGAGTCGGCCAGGGTGAAATTGCGCTGGCCGCGCAcgc                                                                                                                                                   >  6:208819/1‑134 (MQ=255)
               cAGGTGCTCGAAGGCCAGAGCATAGGTGACCAGCCCTACGCCTTTGGTCAGCAGGCCGGTACTGGTCGACAGTGCGCCGGTACTCAGGCTGTTGATCGAGTCGGCCAGGGTGAAATTGCGCTGGCCGCGCAcgc                                                                                                                                                   <  5:208819/134‑1 (MQ=255)
                 ggTGCTCGAAGGCCAGAGCATAGGTGACCAGCCCTACGCCTTTGGTCAGCAGGCCGGTACTGGTCGACAGTGCGCCGGTACTCAGGCTGTTGATCGAGTCGGCCAGGGTGAAATTGCGCTGGCCGCGCACGCGGTCGGCGACCAGctct                                                                                                                                  <  4:302849/149‑1 (MQ=255)
                     cTCGAAGGCCAGAGCATAGGTGACCAGCCCTACGCCTTTGGTCAGCAGGCCGGTACTGGTCGACAGTGCGCCGGTACTCAGGCTGTTGATCGAGTCGGCCAGGGTGAAATTGCGCTGGCCGCGCACGCGGTCGGCGACCAGCTCTACgg                                                                                                                              <  4:317076/149‑1 (MQ=255)
                     cTCGAAGGCCAGAGCATAGGTGACCAGCCCTACGCCTTTGGTCAGCAGGCCGGTACTGGTCGACAGTGCGCCGGTACTCAGGCTGTTGATCGAGTCGGCCAGGGTGAAATTGCGCTGGCCGCGCACGCGGTCGGCGACCAGCTCTACg                                                                                                                               <  8:219438/148‑1 (MQ=255)
                      tCGAAGGCCAGAGCATAGGTGACCAGCCCTACGCCTTTGGTCAGCAGGCCGGTACTGGTCGACAGTGCGCCGGTACTCAGGCTGTTGATCGAGTCGGCCAGGGTGAAATTGCGCTGGCCGCGCACGCGGTCGGCGACCAGCTCTACGGc                                                                                                                             >  1:344555/1‑149 (MQ=255)
                         aaGGCCAGAGCATAGGTGACCAGCCCTACGCCTTTGGTCAGCAGGCCGGTACT‑GTCGACAGTGCGCCGGTACTCAGGCTGTTGATCGAGTCGTCCAGTGTGAAATTGCGCTGGCCGCGCACGCGGTCGGCGACCAGCTCTACGGCAATc                                                                                                                         <  5:303702/149‑1 (MQ=255)
                           ggCCAGAGCATAGGTGACCAGCCCTACGCCTTTGGTCAGCAGGCCGGTACTGGTCGACAGTGCGCCGGTACTCAGGCTGTTGATCGAGTCGGCCAGGGTGAAATTGCGCTGGCCGCGCACGCGGTCGGCGACCAGCTCTACGGCTATCa                                                                                                                        <  3:220713/149‑1 (MQ=255)
                                       ggTGACCAGCCCTACGCCTTTGGTCAGCAGGCCGGTACTGGTCGACAGTGCGCCGGTACTCAGGCTGTTGATCGAGTCGGCCAGGGTGAAATTGCGCTGGCCGCGCACGCGGTCGGCGACCAGCTCTACGGCTATCAGCACAAAAAAg                                                                                                             <  3:400769/148‑1 (MQ=255)
                                       ggTGACCAGCCCTACGCCTTTGGTCAGCAGGCCGGTACTGGTCGACAGTGCGCCGGTACTCAGGCTGTTGATCGAGTCGGCCAGGGTGAAATTGCGCTGGCCGCGCACGCGGTCGGCGACCAGCTCTACGGCTATCAGCACAAAAAAGa                                                                                                            <  7:31407/149‑1 (MQ=255)
                                        gTGACCAGCCCTACGCCTTTGGTCAGCAGGCCGGTACTGGTCGACAGTGCGCCGGTACTCAGGCTGTTGATCGAGTCGGCCAGGGTGAAATTGCGCTGGCCGCGCACGCGGTCGGCGACCAGCTCTACGGCTATCAGCACAAAAAAGaa                                                                                                           >  1:19213/1‑149 (MQ=255)
                                        atctccAGGCCTACGCCTTTGGTCAGCAGGCCGGTACTGGTCGACAGTGCGCCGGTACTCAGGCTGTTGATCGAGTCGGCCAGGGTGAAATTGCGCTGGCCGCGCACGCGGTCGGCGACCAGCTCTACGGCTATCAGCACAAAAAAGaa                                                                                                           <  3:265681/145‑1 (MQ=255)
                                           aCCAGCCCTACGCCTTTGGTCAGCAGGCCGGTACTGGTCGACAGTGCGCCGGTACTCAGGCTGTTGATCGAGTCGGCCAGGGTGAAATTGCGCTGGCCGCGCACGCGGTCGGCGACCAGCTCTACGGCTATCAGCACAAAAAAGAATgg                                                                                                        <  5:4521/149‑1 (MQ=255)
                                           aCCAGCCCTACGCCTTTGGTCAGCAGGCCGGTACTGGTCGACAGTGCGCCGGTACTCAGGCTGTTGATCGAGTCGGCCAGGGTGAAATTGCGCTGGCCGCGCACGCGGTCGGCGACCAGCTCTACGGCTATCAGCACAAAAAAGAATg                                                                                                         <  1:157823/148‑1 (MQ=255)
                                             cAGCCCTACGCCTTTGGTCAGCAGGCCGGTACTGGTCGACAGTGCGCCGGTACTCAGGCTGTTGATCGAGTCGGCCAGGGTGAAATTGCGCTGGCCGCGCACGCGGTCGGCGACCAGCTCTACGGCTATCAGCACAAAAAAGaa                                                                                                           >  4:265681/1‑144 (MQ=255)
                                              aGCCCTACGCCTTTGGTCAGCAGGCCGGTACTGGTCGACAGTGCGCCGGTACTCAGGCTGTTGATCGAGTCGGCCAGGGTGAAATTGCGCTGCCCGCGCACGCGGTCGGCGACCAGCTCTACGGCTATCAGCACAAAAAAGAATGGTAc                                                                                                     >  2:344180/1‑149 (MQ=255)
                                               gCCCTACGCCTTTGGTCAGCAGGCCGGTACTGGTCGACAGTGCGCCGGTACTCAGGCTGTTGATCGAGTCGGCCAGGGTGAAATTGCGCTGGCCGCGCACGCGGTCGGCGACCAGCTCTACGGCTATCAGCACAAAAAAGAATGGTACg                                                                                                    >  3:365966/1‑149 (MQ=255)
                                               gCCCTACGCCTTTGGTCAGCAGGCCGGTACTGGTCGACAGTGCGCCGGTACTCAGGCTGTTGATCGAGTCGGCCAGGGTGAAATTGCGCTGGCCGCGCACGCGGTCGGCGACCAGCTCTACGGCTATCAGCACAAAAAAGAATGGTACg                                                                                                    >  8:59021/1‑149 (MQ=255)
                                                cccTACGCCTTTGGTCAGCAGGCCGGTACTGGTCGACAGTGCGCCGGTACTCAGGCTGTTGATCGAGTCGGCCAGGGTGAAATTGCGCTGGCCGCGCACGCGGTCGGCGACCAGCTCTACGGCTATCAGCACAAAAAAGAATGGTACgg                                                                                                   >  1:274679/1‑149 (MQ=255)
                                                  cTACGCCTTTGGTCAGCAGGCCGGTACTGGTCGACAGTGCGCCGGTACTCAGGCTGTTGATCGAGTCGGCCAGGGTGAAATTGCGCTGGCCGCGCACGCGGTCGGCGACCAGCTCTACGGCTATCAGCACAAAAAAGAATGGTACGGc                                                                                                  >  2:37907/1‑148 (MQ=255)
                                                      gCCTTTGGTCAGCAGGCCGGTACTGGTCGACAGTGCGCCGGTACTCAGGCTGTTGATCGAGTCGGCCAGGGTGAAATTGCGCTGGCCGCGCACGCGGTCGGCGACCAGCTCTACGGCTATCAGCAGAAAAAAGAATGGTACGGCcag                                                                                               >  3:391929/1‑147 (MQ=255)
                                                      gCCTTTGGTCAGCAGGCCGGTACTGGTCGACAGTGCGCCGGTACTCAGGCTGTTGATCGAGTCGGCCAGGGTGAAATTGCGCTGGCCGCGCACGCGGTCGGCGACCAGCTCTACGGCTATCAGCACAAAAAAGAATGGTACGGCcag                                                                                               >  7:226842/1‑147 (MQ=255)
                                                        cTTTGGTCAGCAGGCCGGTACTGGTCGACAGTGCGCCGGTACTCAGGCTGTTGATCGAGTCGGCCAGGGTGAAATTGCGCTGGCCGCGCACGCGGTCGGCGACCAGCTCTACGGCTATCAGCACAAAAAAGAATGGTACGGCcagcag                                                                                            >  2:253068/1‑148 (MQ=255)
                                                         tttGGTCAGCAGGCCGGTACTGGTCGACAGTGCGCCGGTACTCAGGCTGTTGATCGAGTCGGCCAGGGTGAAATTGCGCTGGCCGCGCACGCGGTCGGCGACCAGCTCTACGGCTATCAGCACaa                                                                                                                  <  1:150694/125‑1 (MQ=255)
                                                         tttGGTCAGCAGGCCGGTACTGGTCGACAGTGCGCCGGTACTCAGGCTGTTGATCGAGTCGGCCAGGGTGAAATTGCGCTGGCCGCGCACGCGGTCGGCGACCAGCTCTACGGCTATCAGCACaa                                                                                                                  >  2:150694/1‑125 (MQ=255)
                                                            ggTCAGCAGGCCGGTACTGGTCGACAGTGCGCCGGTACTCAGGCTGTTGATCGAGTCGGCCAGGGTGAAATTGCGCTGGCCGCGCACGCGGTCGGCGACCAGCTCTACGGCTATCAGCACAAAAAAGAATGGTACGGCCAGCAGTATc                                                                                        <  6:1421/148‑1 (MQ=255)
                                                             gTCAGCAGGCCGGTACTGGTCGACAGTGCGCCGGTACTCAGGCTGTTGATCGAGTCGGCCAGGGTGAAATTGCGCCGACCGCGTGCGCGGCCAGCgcaat                                                                                                                                       >  2:343687/1‑96 (MQ=37)
                                                             gTCAGCAGGCCGGTACTGGTCGACAGTGCGCCGGTACTCAGGCTGTTGATCGAGTCGGCCAGGGTGAAATTGCGCCGACCGCGTGCGCGGCCAGCgcaat                                                                                                                                       >  2:349158/1‑96 (MQ=37)
                                                             gTCAGCAGGCCGGTACTGGTCGACAGTGCGCCGGTACTCAGGCTGTTGATCGAGTCGGCCAGGGTGAAATTGCGCCGACCGCGTGCGCGGCCAGCgcaat                                                                                                                                       <  1:343687/100‑5 (MQ=37)
                                                             gTCAGCAGGCCGGTACTGGTCGACAGTGCGCCGGTACTCAGGCTGTTGATCGAGTCGGCCAGGGTGAAATTGCGCCGACCGCGTGCGCGGCCAGCgcaat                                                                                                                                       <  1:349158/100‑5 (MQ=37)
                                                              tCAGCAGGCCGGTACTGGTCGACAGTGCGCCGGTACTCAGGCTGTTGATCGAGTCGGCCAGGGTGAAATTGCGCTGGCCGCGCACGCGGTCGGCGACCAGCTCTACGGCTATCAGcac                                                                                                                    >  6:361404/1‑118 (MQ=255)
                                                              tCAGCAGGCCGGTACTGGTCGACAGTGCGCCGGTACTCAGGCTGTTGATCGAGTCGGCCAGGGTGAAATTGCGCTGGCCGCGCACGCGGTCGGCGACCAGCTCTACGGCTATCAGcac                                                                                                                    <  5:361404/118‑1 (MQ=255)
                                                              tCAGCAGGCCGGTACTGGTCGACAGTGCGCCGGTACTCAGGCTGTTGATCGAGTCGGCCAGGGTGAAATTGCGCTGGCCGCGCACGCGGTCGGCGACCAGCTCTACGGCTATCAGCACAAAAAAGAATGGTACg                                                                                                    <  1:314900/134‑1 (MQ=255)
                                                              tCAGCAGGCCGGTACTGGTCGACAGTGCGCCGGTACTCAGGCTGTTGATCGAGTCGGCCAGGGTGAAATTGCGCTGGCCGCGCACGCGGTCGGCGACCAGCTCTACGGCTATCAGCACAAAAAAGAATGGTACg                                                                                                    >  2:314900/1‑134 (MQ=255)
                                                              tCAGCAGGCCGGTACTGGTCGACAGTGCGCCGGTACTCAGGCTGTTGATCGAGTCGGCCAGGGTGAAATTGCGCTGGCCGCGCACGCGGTCGGCGACCAGCTCTACGGCTATCAGCACAAAAAAGAATGGTACGGCCAGCAGTATCAgg                                                                                     >  8:82375/1‑149 (MQ=255)
                                                               cagcagGCCGGTACTGGTCGACAGTGCGCCGGTACTCAGGCTGTTGATCGAGTCGGCCAGGGTGAAATTGCGCTGGCCGCGCACGCGGTCGGCGACCAGCTCTACGGCTATCAGCACAAAAAAGAATGGTACGGCCAGCAGTATCAGGt                                                                                    >  7:44484/1‑149 (MQ=255)
                                                                agcagGCCGGTACTGGTCGACAGTGCGCCGGTACTCAGGCTGTTGATCGAGTCGGCCAGGGTGAAATTGCGCTGGCCGCGCACGCGGTCGGCGACCAGCTCTACGGCTATCAGCACAAAAAAGAATGGTACGGCCAGCAGTATCAGGt                                                                                    <  1:93690/148‑1 (MQ=255)
                                                                   agtccGGTACTTGTCGACAGTGCGCCGGTACTCAGGCAGTTGATCGAGTCGGCCATGGTGAAATTGCGCTGGCCGCGCACGCGGTCGGCGACCAGCTCTACGGCTATCAGCACAAAAAAGAATGGTACGGCCAGCAGTATCAGGTCCAt                                                                                <  5:364959/146‑1 (MQ=255)
                                                                   agGCCGGTACTGGTCGACAGTGCGCCGGTACTCAGGCTGTTGATCGAGTCGGCCAGGGTGAAATTGCGCTGGCCGCGCACGCGGTCGGCGACCAGCTCTACGGCTATCAGCACAAAAAAGAATGGTACGGCCAGCAGTATCAGGTCCAt                                                                                <  4:41215/149‑1 (MQ=255)
                                                                   agGCCGGTACTGGTCGACAGTGCGCCGGTACTCAGGCTGTTGATCGAGTCGGCCAGGGTGAAATTGCGCTGGCCGCGCACGCGGTCGGCGACCAGCTCTACGGCTATCAGCACAAAAAAGAATGGTACGGCCAGCAGTATCAGGTCCAt                                                                                >  8:109844/1‑149 (MQ=255)
                                                                       cGGTACTGGTCGACAGTGCGCCGGTACTCAGGCTGTTGATCGAGTCGGCCAGGGTGAAATTGCGCTGGCCGCGCACGCGGTCGGCGACCAGCTCTTCGGCTATCAGCACAAAAAAGAATGGTACGGCCAGCAGTATCAGGTCCATGGgc                                                                            <  8:250674/149‑1 (MQ=255)
                                                                        ggTACTGGTCGACAGTGCGCCGGTACTCAGGCTGTTGATCGAGTCGGCCAGGGTGAAATTGCGCTGGCCGCGCACGCGGTCGGCGACCAGCTCTACGGCTATCATCACAAAAAAGAATGGTACGGCCAGCAGTATCAGGTCCATGGgcg                                                                           <  1:90694/149‑1 (MQ=255)
                                                                            cTGGTCGACAGTGCGCCGGTACTCAGGCTGTTGATCGAGTCGGCCAGGGTGAAATTGCGCTGGCCGCGCACGCGGTCGGCGACCAGCTCTACGGCTATCAGCACAAAAAAGAATGGTACGGCCAGCAGTATCAGGTCCATGGGCGCGa                                                                        <  4:365966/148‑1 (MQ=255)
                                                                            cTGGTCGACAGTGCGCCGGTACTCAGGCTGTTGATCGAGTCGGCCAGGGTGAAATTGCGCTGGCCGCGCACGCGGTCGGCGACCAGCTCTACGGCTATCAGCACAAAAAAGAATGGTACGGCCAGCAGTATCAGGTCCATGGGCGCGAc                                                                       <  2:344555/149‑1 (MQ=255)
                                                                             tGGTCGACAGTGCGCCGGTACTCAGGCTGTTGATCGAGTCGGCCAGGGTGAAATTGCGCTGGCCGCGCACGCGGTCGGCGACCAGCTCTACGGCTATCAGCACAAAAAAGAATGGTACGGCCAGCAGTATCAGGTCCATGGGCGCGAcc                                                                      >  5:353761/1‑149 (MQ=255)
                                                                               gTCGACAGTGCGCCGGTACTCAGGCTGTTGATCGAGTCGGCCAGGGTGAAATTGCGCTGGCCGCGCACGCGGTCGGCGACCAGCTCTACGGCTATCAGCACAAAAAAGAATGGTACGGCCAGCAGTATCAGGTCCATGGGCGCGACCTc                                                                    <  2:243653/149‑1 (MQ=255)
                                                                                 cGACAGTGCGCCGGTACTCAGGCTGTTGATCGAGTCGGCCAGGGTGAAATTGCGCTGGCCGCGCACGCGGTCGGCGACCAGCTCTACGGCTATCAGCACAAAAAAGAATGGTACGGCCAGCAGTATCAGGTCCATGGGCGCGACCTCAg                                                                  <  8:226842/149‑1 (MQ=255)
                                                                                  gACAGTGCGCCGGTACTCAGGCTGTTGATCGAGTCGGCCAGGGTGAAATTGCGCTGGCCGCGCACGCGGTCGGCGACCAGCTCTACGGCTATCAGCACAAAAAAGAATGGTACGGCCAGCAGTATCAGGTCCATGGGCGCGACCTCAgg                                                                 <  2:19213/149‑1 (MQ=255)
                                                                                    cAGTGCGCCGGTACTCAGGCTGTTGATCGAGTCGGCCAGGGTGAAATTGCGCTGGCCGCGCACGCGGTCGGCGACCAGCTCTACGGCTATCAGCACAAAAAAGAATGGTACGGCCAGCAGTATCAGGTCCATGGGCGCGACCTCAGGCg                                                               >  3:418691/1‑149 (MQ=255)
                                                                                     aGTGCGCCGGTACTCAGGCTGTTGATCGAGTCGGCCAGGGTGATATTGCGCTGGCCGCGCACGCGGTCGGCGACCAGCTCTACGGCTATCAGCACAAAAAAGAATGGTACGGCCAGCAGTATCAGGTCCATGGGCGCGACCTCAGTCGt                                                              >  6:402872/1‑149 (MQ=255)
                                                                                      gTGCGCCGGTACTCAGGCTGTTGATCGACTCGGCCAGGGTGAAATTGCGCTGGCCGCGCACGCGGTCGGCGACCAGCTCTACGGCTATCAGCACAAAAAAGAATGGTACGGCCAGCAGTATCAGGTCCATGGGCGCGACCTCAGGCGt                                                              >  2:417627/1‑148 (MQ=255)
                                                                                        gcgcCGGTACTCAGGCTGTTGATCGAGTCGGCCAGGGTGAAATTGCGCTGGCCGCGCACGCGGTCGGCGACCAGCTCTACGGCTATCAGCACAAAAAAGAATGGTACGGCCAGCAGTATCAGGTCCATGGGCGCGACCTCAGGCGTTAt                                                           >  5:160453/1‑149 (MQ=255)
                                                                                         cgcCGGTACTCAGGCTGTTGATCGAGTCGGCCAGGGTGAAATTGCGCTGGCCGCGCACGCGGTCGGCGACCAGCTCTACGGCTATCAGCACAAAAAAGAATGGTACGGCCAGCAGTATCAGGTCCATGGGCGCGACCTCAGGCGTTATc                                                          >  3:30474/1‑149 (MQ=255)
                                                                                         cgcCGGTACTCAGGCTGTTGATCGAGTCGGCCAGGGTGAAATTGCGCTGGCCGCGCACGCGGTCGGCGACCAGCTCTACGGCTATCAGCACAAAAAAGAATGGTACGGCCAGCAGTATCAGGTCCATGGGCGCGACCTCAGGCGTTATc                                                          <  2:211159/149‑1 (MQ=255)
                                                                                          gcCGGTACTCAGGCTGTTGATCGAGTCGGCCAGGGTGAAATTGCGCTGGCCGCGCACGCGGTCGGCGACCAGCTCTACGGCTATCAGCACAAAAAAGAATGGTACGGCCAGCAGTATCAGGTCCATGGGCGCGACCTCAGGCGTTATcc                                                         >  5:107945/1‑149 (MQ=255)
                                                                                          gcCGGTACTCAGGCTGTTGATCGAGTCGGCCAGGGTGAAATTGCGCAGGCCGCGCACGCGGTCGGCGACCAGCTCTACGGCTATCAGCACAAAAAAGAATGGTACGGCCAGCAGTATCAGGTCCATGGGCGCGACCTCAGGCGTTATc                                                          >  1:315604/1‑148 (MQ=255)
                                                                                           ccGGTACTCAGGCTGTTGATCGAGTCGGCCAGGGTGAAATTGCGCTGGCCGCGCACGCGGTCGGCGACCAGCTCTACGGCTATCAGCACAAAAAAGAATGGTACGGCCAGCAGTATCAGGTCCATGGGCGCGACCTCAGGCGTTATcca                                                        >  5:233504/1‑149 (MQ=255)
                                                                                              gTACTCAGGCTGTTGATCGAGTCGGCCAGGGTGAAATTGCGCTGGCCGCGGACGCGGTCGGCGACCAGCTCTACGGCTATCAGCACAAAAAAGAATGGTACGGCCAGCAGTATCAGGTCCATGGGCGCGACCTCAGGCGTTATccacca                                                     <  6:125377/149‑1 (MQ=255)
                                                                                              gTACTCAGGCTGTTGATCGAGTCGGCCAGGGTGAAATTGCGCTGGCCGCGCACGCGGTCGGCGACCAGCTCTACGGCTATCAGCACAAAAAAGAATGGTACGGCCAGCAGTATCAGGTCCATGGGCGCGACCTCAGGCGTTATccacca                                                     <  8:44484/149‑1 (MQ=255)
                                                                                                 cTCAGGCTGTTGATCGAGTCGGCCAGGGTGAAATTGCGCTGGCCGCGCACGCGGTCGGCGACCAGCTCTACGGCTATCAGCACAAAAAAGAATGGTACGGCCAGCAGTATCAGGTCCATGGGCGCGACCTCAGGCGTTATCCACCAAg                                                   >  6:239530/1‑148 (MQ=255)
                                                                                                 cTCAGGCTGTTGATCGAGTCGGCCAGGGTGAAATTGCGCTGGCCGCGCACGCGGTCGGCGACCAGCTCTACGGCTATCAGCACAAAAAAGAATGGTACGGCCAGCAGTATCAGGTCCATGGGCGCGACCTCAGGCGTTATCCACCAAGa                                                  <  7:82375/149‑1 (MQ=255)
                                                                                                    aGGCTGTTGATCGAGTCGGCCAGGGTGAAATTGCGCTGGCCGCGCACGCGGTCGGCGACCAGCGCTACGGCTATCAGCACAAAAAAGAATGGTACGGCCAGCAGTATCAGGTCCATGGGCGCGACCTCAGGCGTTATCCACCAAGAt                                                 >  7:7631/1‑147 (MQ=255)
                                                                                                      gCTGTTGATCGAGTCGGCCAGGGTGAAATTGCGCTGGCCGCGCACGCGGTCGGCGACCAGCTCTACGGCTATCAGCACAAAAAAGAATGGTACGGCCAGCAGTATCAGGTCCATGGGCGCGACCTCAGGCGTTATCCACCAAGATTAgg                                             >  1:322368/1‑149 (MQ=255)
                                                                                                         gTTGATCGAGTCGGCCAGGGTGAAATTGCGCTGGCCGCGCACGCGGTCGGCGACCAGCTCTACGGCTATCAGCACAAAAAAGAATGGTACGGCCAGCAGTATCAGGTCCATGGGCGCGACCTCAGGCGTTATCCACCAAGATTAGGCCg                                          >  5:8306/1‑149 (MQ=255)
                                                                                                          ttGATCGAGTCGGCCAGGGTGAAATTGCGCTGGCCGCGCACGCGGTCGGCGACCAGCTCTACGGCTATCAGCACAAAAAAGAATGGTACGGCCAGCAGTATCAGGTCCATGGGCGCCACCTCAGGCGTTATCCACCAAGATTAGGcc                                           <  1:37907/147‑1 (MQ=255)
                                                                                                                         aGGGTGAAATTGCGCTGGCCGCGCACGCGGTCGGCGACCAGCTCTACGGCTATCAGCACAAAAAAGAATGGTACGGCCAGCAGTATCAGGTCCATGGGCGCGACCTCAGGCGTTATCCACCAAGATTAGGCCGCACCCACAGGAATccc                          <  2:315604/149‑1 (MQ=255)
                                                                                                                            gTGAAATTGCGCTGGCCGCGCACGCGGTCGGCGACCAGCTCTACGGCTATCAGCACAAAAAAGAATGGTACGGCCAGCAGTATCAGGTCCATGGGCGCGACCTCAGGCGTTATCCACCAAGATTAGGCCGCACCCACAGGAATCCCTAt                       <  3:35251/149‑1 (MQ=255)
                                                                                                                             tGAAATTGCGCTGGCCGCGCAGGCGGTCGGCGACCAGCTCTACGGCTATCAGCACAAAAAAGAATGGTACGGCCAGCAGTATCAGGTCCATGGGCGCGACCTCAGGCGTTATCCACCAAGATTAGGCCGCACCCACAGGAATCCCTAt                       <  1:417627/148‑1 (MQ=255)
                                                                                                                                       cTGGCCGCGCACGCGGTCGGCGACCAGCTCTACGGCTATCAGCACAAAAAAGAATGGTACGGCCAGCAGTATCAGGTCCATGGGCGCGACCTCAGGCGTTATCCACCAAGATTAGGCCGCACCCACAGGAATCCCTATGGCGACATCTg            >  7:11255/1‑149 (MQ=255)
                                                                                                                                        tGGCCGCGCACGCGGTCGGCGACCAGCTCTACGGCTATCAGCACAAAAAAGAATGGTACGGCCAGCAGTATCAGGTCCATGGGCGCGACCTCAGGCGTTATCCACCAAGATTAGGCCGCACCCACAGGAATCCCTATGGCGACATCTGc           >  5:243452/1‑149 (MQ=255)
                                                                                                                                          gCCGCGCACGCGGTCGGCGACCAGCTCTACGGCTATCAGCACAAAAAAGAATGGTACGGCCAGCAGTATCAGGTCCATGGGCGCGACCTCAGGCGTTATCCACCAAGATTAGGCCGCACCCACAGGAATCCCTATGGCTACATCTGcc          <  5:239530/148‑1 (MQ=255)
                                                                                                                                                cACGCGGTCGGCGACCAGCTCTACGGCTATCAGCACAAAAAAGAATGGTACGGCCAGCAGTATCAGGTCCATGGGCGCGACCTCAGGCGTTATCCACCAAGATTAGGCCGCACCCACAGGAATCCCTATGGCGACATCTGCCAAACTag   <  2:319670/149‑1 (MQ=255)
                                                                                                                                                 aCGCGGTCGGCGACCAGCTCTACGGCTATCAGCACAAAAAAGAATGGTACGGCCAGCAGTATCAGGTCCATGGGCGCGACCTCAGGCGTTATCCACCAAGATTAGGCCGCACCCACAGGAATCCCTATGGCGACATCTGCCAAACTaga  <  2:322368/149‑1 (MQ=255)
                                                                                                                                                 aCGCGGTCGGCGACCAGCTCTACGGCTATCAGCACAAAAAAGAATGGTACGGCCAGCAGTATCAGGTCCATGGGCGCGACCTCAGGCGTTATCCACCAAGATTAGGCCGCACCCACAGGAATCCCTATGGCGACATCTGCCAAACTaga  <  1:344180/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
CAATTGCAGCAGCGCCAGGTGCTCGAAGGCCAGAGCATAGGTGACCAGCCCTACGCCTTTGGTCAGCAGGCCGGTACTGGTCGACAGTGCGCCGGTACTCAGGCTGTTGATCGAGTCGGCCAGGGTGAAATTGCGCTGGCCGCGCACGCGGTCGGCGACCAGCTCTACGGCTATCAGCACAAAAAAGAATGGTACGGCCAGCAGTATCAGGTCCATGGGCGCGACCTCAGGCGTTATCCACCAAGATTAGGCCGCACCCACAGGAATCCCTATGGCGACATCTGCCAAACTAGA  >  NC_002947/5954787‑5955080

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