Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 198,014 A→G 7.2% N1174D (AAT→GAT)  PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947198,0140AG 7.2% 196.0 / 5.2 69N1174D (AAT→GAT) PP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base A (31/33);  new base G (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 = 6.82e-01

GTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACT‑‑CCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGTCACCGTAAGTATTGTCAGCAATGGCAATGTGACCGAAGACCAGCAGCCTTCCTTTACTGTGAAAGTAAGCCAGGCGCTGGACCGTCCGTTGACGGTAACCCTGTCTAACGGCGACACCGTGACGATC  >  NC_002947/197868‑198160
                                                                                                                                                      |                                                                                                                                                   
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                                                                                                                            <  6:117818/148‑1 (MQ=255)
                        tCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACT‑‑CCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGTCACCGt                                                                                                                            <  8:236597/148‑1 (MQ=255)
                         ccAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACT‑‑CCGAATCCGCCCCCAGCCTCGACTGTCATCACTGCCAGCATCGATACCGTCAGCGt                                                                                                                            >  6:320889/1‑147 (MQ=255)
                               cccGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACT‑‑CCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGTCACCGTAAGTATTg                                                                                                                    <  7:128350/149‑1 (MQ=255)
                               cccGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACT‑‑CCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGTCACCGTAAGTATTg                                                                                                                    >  1:78375/1‑149 (MQ=255)
                               cccGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACT‑‑CCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGTCACCGTAAGTATTg                                                                                                                    >  5:26891/1‑149 (MQ=255)
                               cccGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACT‑‑CCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGTCACCGTAAGTATTg                                                                                                                    >  2:223888/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                                                           <  4:375802/149‑1 (MQ=255)
                                                                                       ccACGGGCGGCAACTTCGAACATCTGACT‑‑CCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGTCACCGTAAGTATTGTCAGCAATGGCAATGTGACCGAAGACCAGCAGCCTTCCTTTACTGTGAAAGTAAGc                                                           <  1:34438/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)
                                                                                                                                                      |                                                                                                                                                   
GTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACT‑‑CCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGTCACCGTAAGTATTGTCAGCAATGGCAATGTGACCGAAGACCAGCAGCCTTCCTTTACTGTGAAAGTAAGCCAGGCGCTGGACCGTCCGTTGACGGTAACCCTGTCTAACGGCGACACCGTGACGATC  >  NC_002947/197868‑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: