Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 198,013 C→T 7.3% I1173I (ATC→ATT PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947198,0130CT 7.3% 196.9 / 3.7 69I1173I (ATC→ATTPP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base C (30/33);  new base T (2/3);  total (33/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.06e-01

CGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACT‑‑CCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGTCACCGTAAGTATTGTCAGCAATGGCAATGTGACCGAAGACCAGCAGCCTTCCTTTACTGTGAAAGTAAGCCAGGCGCTGGACCGTCCGTTGACGGTAACCCTGTCTAACGGCGACACCGTGACGATC  >  NC_002947/197865‑198160
                                                                                                                                                        |                                                                                                                                                    
cGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACT‑‑CCGAATCCGACCCCAGCCTCGACTGtcatc                                                                                                                                                      <  1:5785/149‑1 (MQ=255)
   gTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCTAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACT‑‑CCGAATCCGACCCCAGCCTCGACTGTCATCAAt                                                                                                                                                   <  8:230900/149‑1 (MQ=255)
    tAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACT‑‑CCGAATCCGACCCCAGCCTCGACTGTCATCAATg                                                                                                                                                  <  6:255087/149‑1 (MQ=255)
     aaaaCCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACT‑‑CCGAATCCGACCCCAGCCTCGACTGTCATCAATGa                                                                                                                                                 >  1:161347/1‑149 (MQ=255)
      aaaCCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACT‑‑CCGAATCCGACCCCAGCCTCGACTGTCATCAATGAc                                                                                                                                                <  5:263794/149‑1 (MQ=255)
       aaCCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACT‑‑CCGAATCCGACCCCAGCCTCGACTGTCATCAATGACa                                                                                                                                               <  6:122583/149‑1 (MQ=255)
            aGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACT‑‑CCGAATCCGACCCCAGCCTCGACTGTCATTGATGACAgtcga                                                                                                                                          <  1:94945/149‑5 (MQ=255)
            aGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACT‑‑CCGAATCCGACCCCAGCCTCGACTGTCATTGATGACAgtcga                                                                                                                                          <  1:88878/149‑5 (MQ=255)
               gCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACT‑‑CCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGAt                                                                                                                                       <  8:158278/149‑1 (MQ=255)
                       gATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACT‑‑CCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGTCAc                                                                                                                               <  2:171461/149‑1 (MQ=255)
                        aTTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACT‑‑CCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGTCAcc                                                                                                                              >  5:117818/1‑149 (MQ=255)
                         tttCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACT‑‑CCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGTCACCg                                                                                                                             >  2:295296/1‑149 (MQ=255)
                           tCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACT‑‑CCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGTCACCGt                                                                                                                            <  8:236597/148‑1 (MQ=255)
                           tCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACT‑‑CCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGTCACCGt                                                                                                                            <  6:117818/148‑1 (MQ=255)
                            ccAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACT‑‑CCGAATCCGCCCCCAGCCTCGACTGTCATCACTGCCAGCATCGATACCGTCAGCGt                                                                                                                            >  6:320889/1‑147 (MQ=255)
                                  cccGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACT‑‑CCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGTCACCGTAAGTATTg                                                                                                                    <  7:128350/149‑1 (MQ=255)
                                  cccGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACT‑‑CCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGTCACCGTAAGTATTg                                                                                                                    >  5:26891/1‑149 (MQ=255)
                                  cccGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACT‑‑CCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGTCACCGTAAGTATTg                                                                                                                    >  2:223888/1‑149 (MQ=255)
                                  cccGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACT‑‑CCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGTCACCGTAAGTATTg                                                                                                                    >  1:78375/1‑149 (MQ=255)
                                   ccGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACT‑‑CCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGTCACCGTAAGTATTGt                                                                                                                   >  8:349615/1‑149 (MQ=255)
                                          aTGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACT‑‑CCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGTCACCGTAAGTATTGTCAgcaat                                                                                                           <  2:161347/149‑1 (MQ=255)
                                             aCGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACT‑‑CCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGTCACCGTAAGTATTGTCAgcaatgg                                                                                                         <  2:78375/148‑1 (MQ=255)
                                             aCGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACT‑‑CCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGTCACCGTAAGTATTGTCAgcaatgg                                                                                                         >  3:137309/1‑148 (MQ=255)
                                               gTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACT‑‑CCGAATCCGACCCCAGCCTCGACTGTCATTGATGACAgtcgag                                                                                                                                         <  7:291381/115‑6 (MQ=255)
                                               gTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACT‑‑CCGAATCCGACCCCAGCCTCGACTGTCATTGATGACAgtcgag                                                                                                                                         >  8:291381/1‑110 (MQ=255)
                                                tCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACT‑‑CCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGTCACCGTAAGTATTGTCAgcaatggcaat                                                                                                     <  5:320889/149‑1 (MQ=255)
                                                tCTACAACAACGGCTCGACCGGCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACT‑‑CCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGTCACCGTAAGTATTGTCAgcaatggcaat                                                                                                     <  4:198956/149‑1 (MQ=255)
                                                   acaacaACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACT‑‑CCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGTCACCGTAAGTATTGTCAGCAATGGCAAtgt                                                                                                   >  4:257730/1‑148 (MQ=255)
                                                    ctacaacGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACT‑‑CCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGTCACCGTAAGTATTGTCAGCAATGGCAATGTGa                                                                                                 <  4:312131/147‑1 (MQ=255)
                                                    caacaaCGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACT‑‑CCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGTCACCGTAAGTATTGTCAGCAATGGCAATGTGa                                                                                                 <  7:107354/149‑1 (MQ=255)
                                                     aacaacGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACT‑‑CCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGTCACCGTAAGTATTGTCAGCAATGGCAATGTGAc                                                                                                >  1:341313/1‑149 (MQ=255)
                                                      acaacGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAAGCTTCGAACATCTGACT‑‑CCGAATCCGACCCCACCCTCGACTGTCATCAATGACAGCATCGATACCGTCCCCGTAAGTATTGTCAGCAATGGCAATGTGAc                                                                                                >  1:289221/1‑149 (MQ=255)
                                                           ggCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTAGGCGAATCCGCCCCCAGCCTCGACTGTCATAAATGACAGCAACGATACCGTCACCGTAAGTATTGTAAGCAATGGCAATGTGACCGaa                                                                                            >  1:112462/1‑149 (MQ=37)
                                                            gCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACT‑‑CCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGTCACCGTAAGTATTGTCAGCAATGGCAATGTGACCGAAGAc                                                                                         <  1:109084/149‑1 (MQ=255)
                                                            gCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACT‑‑CCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGTCACCGTAAGTATTGTCAGCAATGGCAATGTGACCGAAGAc                                                                                         <  7:101082/149‑1 (MQ=255)
                                                            gCTCCACCGTCAGCGCCACCAACGCGACCGCCCCGCGCGGCAACACCGCACCTCTGACC‑‑CCGACTCCGACCCCAGCCTCGCCCGTCACCAATCACAGCAGCGATACCGTCACCGTCAGCACTGTCAGCAATGGCAATGTGACCGAacac                                                                                         <  4:137309/149‑1 (MQ=255)
                                                             cTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACT‑‑CCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGTCACCGTAAGTATTGTCAGCAATGGCAATGTGACCGAAGAc                                                                                         >  7:45479/1‑148 (MQ=255)
                                                              tCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACT‑‑CCGAATCCGACCCCAGCCTCGGCTGTCATCAATGACAGCATCGATACCGTCACCGTAAGTATTGTCAGCAATGGCAATGTGACCGAAGAc                                                                                         >  4:16959/1‑147 (MQ=255)
                                                                gACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACT‑‑CCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGTCACCATAAGTATTGTCAGCAATGGCAATGTGACCGAAGACcag                                                                                      >  7:235443/1‑148 (MQ=255)
                                                                 aCCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACT‑‑CCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGTCACCGTAAGTATTGTCAGCAATGGCAATGTGACCGAAGACcagca                                                                                    <  7:349615/149‑1 (MQ=255)
                                                                       aGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACT‑‑CCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGTCACCGTAAGTATTGTCAGCAATGGCAATGTGACCGAAGACCAGCAGccttc                                                                              >  6:253800/1‑149 (MQ=255)
                                                                           tCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACT‑‑CCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGTCACCGTAAGTATTGTCAGCAATGGCAATGTGACCGAAGACCAGCAGCCTTCCttt                                                                          <  6:26891/149‑1 (MQ=255)
                                                                               cATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACT‑‑CCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGTCACCGTAAGTATTGTCAGCAATGGCAATGTGACCGAAGACCAGCAGCCTTCCTTTACtg                                                                      >  4:258038/1‑149 (MQ=255)
                                                                               cATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACT‑‑CCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGTCACCGTAAGTATTGTCAGCAATGGCAATGTGACCGAAGACCAGCAGCCTTCCTTTACtg                                                                      >  6:45526/1‑149 (MQ=255)
                                                                                    agagCGCCACGGGCGGCAACTTCGAACATCTGACT‑‑CCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGTCACCGTAAGTATTGTCAGCAATGGCAATGTGACCGAAGACCAGCAGCCTTCCTTTACTGTGaaa                                                                 <  1:261689/149‑1 (MQ=255)
                                                                                    agagCGCCACGGGCGGCAACTTCGAACATCTGACT‑‑CCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGTCACCGTAAGTATTGTCAGCAATGGCAATGTGACCGAAGACCAGCAGCCTTCCTTTACTGTGa                                                                   >  8:225769/1‑147 (MQ=255)
                                                                                    agagCGCCACGGGCGGCAACTTCGAACATCTGACT‑‑CCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGGCACCGTAAGTATTGTCAGCAATGGCAATGTGCCCGAAGACCAGCAGCCTTCCTTTAGTGTGaaa                                                                 >  6:359653/1‑149 (MQ=255)
                                                                                          ccACGGGCGGCAACTTCGAACATCTGGCT‑‑CCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGTCACCGTAAGTATTGTCAGCAATGGCAATGTGACCGAAGACCAGCAGCCTTCCTTTACTGTGAAAGTAAGc                                                           >  4:187486/1‑149 (MQ=255)
                                                                                          ccACGGGCGGCAACTTCGAACATCTGACT‑‑CCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGTCACCGTAAGTATTGTCAGCAATGGCAATGTGACCGAAGACCAGCAGCCTTCCTTTACTGTGAAAGTAAGc                                                           <  1:34438/149‑1 (MQ=255)
                                                                                          ccACGGGCGGCAACTTCGAACATCTGACT‑‑CCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGTCACCGTAAGTATTGTCAGCAATGGCAATGTGACCGAAGACCAGCAGCCTTCCTTTACTGTGAAAGTAAGc                                                           <  4:375802/149‑1 (MQ=255)
                                                                                           cACGGGCGGCAACTTCGAACATCTGACT‑‑CCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGTCACCGTAAGTATTGTCAGCAATGGCAATGTGACCGAAGACCAGCAGCCTTCCTTTACTGTGAAAGTAAGc                                                           >  5:96341/1‑148 (MQ=255)
                                                                                           cACCGGCGGCAAATTCGAACATCTGACT‑‑CCGAAACCGACCCCAGCCTCGACTTTCATCAATGACATCATCGATACCGTCACCGTAATTATTGTCACCATATGGCAATGTGACCGAAGACCAGCAGCCTTCCTTTACTGTGAAAGTAAGc                                                           <  1:295296/149‑1 (MQ=255)
                                                                                             cGGGCGGCAACTTCGAACATCTGACT‑‑CCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGTCACCATAAGTATTGTCAGCAATGGCAATGTGACCGAAGACCAGCAGCCTTCCTTTACTTTCACAGTaaaggaa                                                        <  8:235443/149‑6 (MQ=255)
                                                                                                gcggcAACTTCGAACATCTGACT‑‑CCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGTCACCGTAAGTATTGTCAGCAATGGCAATGTGACCGAAGACCAGCAGCCTTCCTTTACTGTGAAAGTAAGCCAGgcg                                                     <  6:134806/149‑1 (MQ=255)
                                                                                                   gcAACTTCGAACATCTGACT‑‑CCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGTCACCGTAAGTATTGTCAGCAATGGCAATGTGACCGAAGACCAGCAGCCTTCCTTTACTGTGAAAGTAAGCCAGGCGCTg                                                  <  8:45479/149‑1 (MQ=255)
                                                                                                     aaCTACGAACATCTGACT‑‑CCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGTCACCGTAAGTATTGTCAGCAATGGCAATGTGACCGAAGACCAGCAGCCTTCCTTTACTGTGAAAGTAAGCCAGGCGCTGGa                                                <  8:348005/149‑1 (MQ=255)
                                                                                                            gAACATCTGACT‑‑CCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGTCACCGTAAGTATTGTCAGCAATGGCAATGTGACCGAAGACCAGCAGCCTTCCTTTACTGTGAAAGTAAGCCAGGCGCTGGAccgtc                                           <  5:45526/148‑1 (MQ=255)
                                                                                                              aCATCTGACT‑‑CCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGTCACCCAAAGTATTGTCAGCAAAGGCAATGTGACCGAAGACCATCACCCTTCCTTTACTGTGAAAGTAATCCAGCCGCTGCAccgtacgt                                        >  2:342128/1‑149 (MQ=255)
                                                                                                                     aCT‑‑CCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGTCACCGTAAGTATTGTCAGCAATGGCAATGTGACCGAAGACCAGCAGCCTTCCTTTACTGTGAAAGTAAGCCAGGCGCTGGACCGTCCGTTGa                                     <  3:258038/145‑1 (MQ=255)
                                                                                                                      cT‑‑CCGAATCCGACCCCAGCCTCGAACTGTCATCAATGACAGCATCGATACCGTCACCGTAAGTATTGTCAGCAATGGCAATGTGACCGAAGACCAGCAGCCTTCCTTTACTGTGAAAGTAAGCCAGGCGCTGGACCGTCCGTTGACgg                                  >  7:57920/1‑148 (MQ=255)
                                                                                                                       t‑‑CCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGTCACCGTAAGTATTGTCAGCAATGGCAATGTGACCGAAGACCAGCAGCCTTCCTTTACTGTGAAAGTAAGACAGGCGCTGGACCGTCCGTTGACGGTaa                               >  8:91924/1‑149 (MQ=255)
                                                                                                                           cGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGTCACCGTAAGTATTGTCAGCAATGGCAATGTGACCGAAGACCAGCAGCCTTCCTTTACTGTGAAAGTAAGCCAGGCGCTGGACCGTCCGTTGACGGTAAcc                             >  4:122088/1‑149 (MQ=255)
                                                                                                                                  gACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGTCACCGTAAGTATTGTCAGCAATGGCAATGTGACCGAAGACCAGCAGCCTTCCTTTACTGTGAAAGTAAGCCAGGCGCTGGACCGTCCGTTGACGGTAACCCTGTc                        <  3:131528/147‑1 (MQ=255)
                                                                                                                                  gACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGTCACCGTAAGTATTGTCAGCAATGGCAATGTGACCGAAGACCAGCAGCCTTCCTTTACTGTGAAAGTAAGCCAGGCGCTGGACCGTCCGTTGACGGTAACCCTGTc                        >  4:131528/1‑147 (MQ=255)
                                                                                                                                       cAGCCTCGACTGTCATCAATGACAGCATCGATACCGTCACCGTAAGTATTGTCAGCAATGGCAATGTGACCGAAGACCAGCAGCCTTCCTTTACTGTGAAAGTAAGCCAGGCGCTGGACCGTCCGTTGACGGTAACCCTGTCTAACGGc                 <  1:366568/149‑1 (MQ=255)
                                                                                                                                       cAGCCTCGACTGTCATCAATGACAGCATCGATACCGTCACCGTAAGTATTGTCAGCAATGGCAATGTGACCGAAGACCAGCAGCCTTCCTTTACTGTGAAAGTAAGCCAGGCGCTGGACCGTCCGTTGACGGTAACCCTGTCTAACGGc                 <  2:341313/149‑1 (MQ=255)
                                                                                                                                            tCGACTGTCATCAATGACAGCATCGATACCGTCACCGTAAGTATTGTCAGCAATGGCAATGTGACCGAAGACAAGCAGCCTTCCTTTACTGTGAAAGTAAGCCAGGCGCTGGACCGTCCGTTGACGGTAACCCTGTCTAACGGCGacac            >  1:338508/1‑149 (MQ=255)
                                                                                                                                              atgcTGTCATTGATGACAGCATCGATACCGTCACCGTAAGTATTGTCAGCAATGGCAATGTGACCGAAGACCAGCAGCCTTCCTTTACTGTGAAAGTAAGCCAGGCGCTGGACCGTCCGTTGACGGTAACCCTGTCTAACGGCGACAcc           >  1:188446/4‑149 (MQ=255)
                                                                                                                                                 cTGTCATCAATGACAGCATCGATACCGTCACCGTAAGTATTGTCAGCAATGGCAATGTGACCGAAGACCAGCAGCCTTCCTTTACTGTGAAAGTAAGCCAGGCGCTGGACCGTCCGTTGACGGTAACCCTGTCTAACGGCGACACCGt         >  4:1087/1‑148 (MQ=255)
                                                                                                                                                  tGTCATCAATGACAGCATCGATACCGTCACCGTAAGTATTGTCAGCAATGGCAATGTGACCGAAGACCAGCAGCCTTCCTTTACTGTGAAAGTAAGCCAGGCGCTGGACCGTCCGTTGACGGTAACCCTGTCTAACGGCGACACCGTGa       >  2:233649/1‑149 (MQ=255)
                                                                                                                                                    tcatcaATGACAGCATCGATACCGTCACCGTAAGTATTGTCAGCAATGGCAATGTGACCGAAGACCAGCAGCCTTCCTTTACTGTGAAAGTAAGCCAGGCGCTGGACCGTCCGTTGACGGTAACCCTGTCTAACGGCGACACCGTGAc      >  5:264286/1‑148 (MQ=255)
                                                                                                                                                      atcaATGACAGCATCGATACCGTCACCGTAAGTATTGTCAGCAATGGCAATGTGACCGAAGACCAGCAGCCTTCCTTTACTGTGAAAGTAAGCCAGGCGCTGGACCGTCCGTTGACGGTAACCCTGTCTAACGGCGACACCGTGACGAt   >  2:82733/1‑149 (MQ=255)
                                                                                                                                                       tcaATGACAGCATCGATACCGTCACCGTAAGTATTGTCAGCAATGGCAATGTGACCGAAGACCTGCAGCCTTCCTTTACTGTGAAAGTAAGCCAGGCGCTGGACCGTCCGTTGACGGTAACCCTGTCTAACGGCGACACCGTGACGATc  <  2:289221/149‑1 (MQ=255)
                                                                                                                                                        |                                                                                                                                                    
CGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACT‑‑CCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGTCACCGTAAGTATTGTCAGCAATGGCAATGTGACCGAAGACCAGCAGCCTTCCTTTACTGTGAAAGTAAGCCAGGCGCTGGACCGTCCGTTGACGGTAACCCTGTCTAACGGCGACACCGTGACGATC  >  NC_002947/197865‑198160

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: