Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 197,572 G→A 13.9% K1026K (AAG→AAA PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947197,5720GA13.9% 181.0 / 22.2 86K1026K (AAG→AAAPP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base G (41/33);  new base A (4/8);  total (45/41)
Fisher's exact test for biased strand distribution p-value = 2.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.34e-01

TGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACC  >  NC_002947/197428‑197716
                                                                                                                                                  |                                                                                                                                                
tGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGAccc                                                                                                                                              <  3:272365/149‑1 (MQ=37)
  ccaccaCCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCa                                                                                                                                             <  1:195886/148‑1 (MQ=37)
  ccaccaCCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAg                                                                                                                                            <  5:90214/149‑1 (MQ=37)
   caccacCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAgg                                                                                                                                           >  4:352742/1‑149 (MQ=37)
   caccacCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAgg                                                                                                                                           <  7:297338/149‑1 (MQ=37)
   caccacCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAgg                                                                                                                                           >  2:281363/1‑149 (MQ=37)
    accaccGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAggg                                                                                                                                          >  2:235125/1‑149 (MQ=37)
      caccGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAggg                                                                                                                                          >  1:118783/1‑147 (MQ=37)
      caccGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAggg                                                                                                                                          >  5:79159/1‑147 (MQ=37)
         cgcgACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCg                                                                                                                                     >  4:80652/1‑149 (MQ=37)
          gcgACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGt                                                                                                                                    <  1:328851/149‑1 (MQ=37)
           cgACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCAGCAAGACCCAGGGCAGCGTc                                                                                                                                   >  1:90755/1‑149 (MQ=37)
              cccTGACGGCGCGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCCTGATCCTGAGCAATCCTGCCCAGACCCCGGTAACCGAGACCCTGACCAACGGCAAACCCATTACCTGTGAAGCCGGCACGACCCAGGCCAGCGTc                                                                                                                                   >  8:166593/1‑146 (MQ=11)
               cgtGAGGGCGAGCCCGTCGGTCAGCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAAGCGTGACGCTGTCCAACGGCCATACCTTTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGAtt                                                                                                                               <  7:166593/147‑1 (MQ=11)
                cTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTCACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGAtt                                                                                                                               >  6:78032/1‑148 (MQ=37)
                 tgccGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGCTCCTGAGCAATCCTGCCCAGCCGCCGGTCACCGTGACCCTGTCCCCCGGCCACACCATGACCGTTGCAGCCGGCAAGCCCCCGGGCAGCGGCGAttgc                                                                                                                             >  8:134821/4‑147 (MQ=11)
                     ggCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTcc                                                                                                                            <  2:238199/146‑1 (MQ=37)
                     ggCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTcc                                                                                                                            >  1:238199/1‑146 (MQ=37)
                      gCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGAGCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGCCAGCGTCGATTTCCAGAc                                                                                                                        <  5:199601/149‑1 (MQ=21)
                            ccGGCGGTCACCGAAGGCGGCGGGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGCCGCCGGTAACCGGGACCCTGTCCCACGGCCAACCCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGc                                                                                                                  <  7:134821/149‑1 (MQ=11)
                                 ggTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAAt                                                                                                              <  7:162000/148‑1 (MQ=37)
                                 ggTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAAt                                                                                                              <  3:80652/148‑1 (MQ=37)
                                 ggTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATg                                                                                                             <  1:315768/149‑1 (MQ=37)
                                    cACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACg                                                                                                          >  6:315896/1‑149 (MQ=37)
                                    cACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACAGTGACCCTGGCCCACGGCCAAACAATTACCGTTGAAGCCGGCAAGACCCAGGGGAGCGGCGAGTTGCAGCCCCCGGCGATTGACg                                                                                                          >  6:131162/1‑149 (MQ=11)
                                      ccGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTc                                                                                                        >  6:128665/1‑149 (MQ=37)
                                       cGAAGGCGGCTTGATCACCTTCACCGTGATCATGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCGTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCt                                                                                                       <  6:79159/149‑1 (MQ=255)
                                         aaGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTa                                                                                                      >  5:282065/1‑148 (MQ=37)
                                                 gataTCACCTAACACGGTATCCTGACCAATCCTGCCCAGACGCCGGTAACCGTTACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACg                                                                                              <  7:51693/145‑1 (MQ=11)
                                                    aTCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAATCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTc                                                                                          >  5:151769/1‑149 (MQ=37)
                                                    aTCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAATCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTc                                                                                          <  6:151769/149‑1 (MQ=37)
                                                       aCCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGGCGATTTGCCGACCCCGGCGAATGAAGGGGACAACAACggagaacg                                                                                       >  8:51693/1‑143 (MQ=11)
                                                        ccTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGAcc                                                                                      <  1:281363/149‑2 (MQ=25)
                                                         cTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGAccg                                                                                     >  6:129616/1‑147 (MQ=25)
                                                             cacCGTTATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCa                                                                                  <  3:280145/149‑1 (MQ=255)
                                                               ccGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCAGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGc                                                                                <  2:90755/149‑1 (MQ=255)
                                                                  tGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTc                                                                             >  4:297112/1‑149 (MQ=21)
                                                                    aTCCTGAGCAATCCTGCCCAGACGCCGGTAACCCTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTc                                                                             <  2:111331/147‑1 (MQ=18)
                                                                       cTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGCCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCAt                                                                        >  2:108118/1‑149 (MQ=17)
                                                                        tGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATc                                                                       >  7:280644/1‑149 (MQ=37)
                                                                        tGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATc                                                                       >  2:81843/1‑149 (MQ=17)
                                                                        tGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATc                                                                       >  6:88057/1‑149 (MQ=17)
                                                                        tGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACAATc                                                                       >  2:155729/1‑149 (MQ=17)
                                                                        tGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACACAGGCGAATCACGTCTACAACAACGTTAGTACCGTAAGCTTCACCAt                                                                        >  8:167971/1‑148 (MQ=11)
                                                                        tGAGCAATCCTGCCAAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGCCGTCTACAACAACGGTACGACTGTCAGCGTCACCATc                                                                       >  7:134172/1‑149 (MQ=255)
                                                                          aGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATc                                                                       >  2:67538/1‑147 (MQ=37)
                                                                              aTCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGt                                                                                                         <  7:70377/109‑1 (MQ=35)
                                                                              aTCCTGCCCAGACGCCGGTAACCGTGACCCTGTCAAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGt                                                                                                         >  8:70377/1‑109 (MQ=25)
                                                                                  tGCCCAGACGCCGGTAACCGTGCCCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGCACGCca                                                             >  7:55731/1‑149 (MQ=255)
                                                                                     ccAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACaacaac                                                                                               >  1:60121/1‑112 (MQ=25)
                                                                                      cAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAgg                                                         <  6:322483/149‑1 (MQ=255)
                                                                                      cAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGGCACCATCGAGAACGCCACAgg                                                         >  5:353608/1‑149 (MQ=11)
                                                                                      cAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAgg                                                         <  5:128665/149‑1 (MQ=17)
                                                                                       aGACGCCGGTAACCGTGACCCTGTCCAACCGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCCATATCAAGACCCCCGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATAGATAACGCCACAgtc                                                        >  4:12763/1‑147 (MQ=11)
                                                                                        gACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCg                                                       <  5:315896/149‑1 (MQ=37)
                                                                                        gACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCg                                                       <  6:323199/149‑1 (MQ=37)
                                                                                         aCGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCgg                                                      <  3:29466/149‑1 (MQ=17)
                                                                                              ggTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAAtt                                                 <  5:108045/149‑1 (MQ=255)
                                                                                              ggTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAgg                                                         <  3:28993/141‑1 (MQ=17)
                                                                                              ggTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAgg                                                         >  4:28993/1‑141 (MQ=17)
                                                                                               gTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAAttt                                                <  1:67538/149‑1 (MQ=37)
                                                                                                tAACCGTGACCCTGTCCAACGGCCAAACCATTACAGTTGTAGCCGTCAAGACCCAGGTCAGCGTCGATTTCCATACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCCTCACCATCGAGAACGCCACAGGCGGTAAtt                                                 <  5:78032/147‑1 (MQ=11)
                                                                                                tAACCGTGACCCTGCCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTc                                               <  1:108118/149‑1 (MQ=14)
                                                                                                   ccGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAg                                            <  5:223314/149‑1 (MQ=16)
                                                                                                    cGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCATGACCCAGGCCAGCGTCGATTTCCTGACCGGGGCGTATGAGGTCTACATCAAGGGTTCGACTGTCAGCGTCACCATCGAGAACGGCTCTGGCGGTAATTTCGTGc                                           <  6:236340/149‑1 (MQ=11)
                                                                                                     gTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACACCGGGTCGACTGTCCGCGTCACGATCGAGAACGCCACACGCGGTAATTTGGAGc                                           >  4:15908/1‑148 (MQ=11)
                                                                                                         cccTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACt                                       >  1:329796/1‑148 (MQ=12)
                                                                                                         cccTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTg                                      >  3:255285/1‑149 (MQ=32)
                                                                                                           cTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGAc                                    <  3:352742/149‑1 (MQ=32)
                                                                                                                cAACGGCCAAACGATTACCGTTGAAGACGGCAAAGAACCAGTGCAGATTCGACTTCCAGACACAGGCGAATTAAGTCTACAACAACGGTTCGAATGTCAGAGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTTAcc                                   >  1:235330/1‑146 (MQ=11)
                                                                                                                  aCGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAAc                             <  2:302650/149‑1 (MQ=12)
                                                                                                                      ccAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGa                         >  6:126000/1‑149 (MQ=34)
                                                                                                                         aaCGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGa                                             <  7:167759/126‑1 (MQ=12)
                                                                                                                         aaCGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTAGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGa                                             >  8:167759/1‑126 (MQ=255)
                                                                                                                           cgaTTACCGTTGAAGCCGGCAAAACCCAGGGCAGCTTCTACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCACCGTCACCATCGAGAACGCCACAGGCTGTAATTTCGAGCAACTGACCCCGAACCCGACCCCg                    <  6:263684/147‑1 (MQ=11)
                                                                                                                             aTTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCAGGc                  >  3:250657/1‑149 (MQ=37)
                                                                                                                             aTTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACAGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGc                  >  1:57758/1‑149 (MQ=21)
                                                                                                                                 ccGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCa               <  5:126000/148‑1 (MQ=37)
                                                                                                                                 ccGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCa               <  8:280644/148‑1 (MQ=37)
                                                                                                                                 ccGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGGCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACCCCGGCTCAg              <  4:210745/149‑1 (MQ=255)
                                                                                                                                   gTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGAc            >  2:174242/1‑149 (MQ=37)
                                                                                                                                    ttGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGa                                             >  7:325723/1‑115 (MQ=32)
                                                                                                                                    ttGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGa                                             <  8:325723/115‑1 (MQ=32)
                                                                                                                                          ccGGCAAAACCCAGGCCAGCGTCTACTTCCATACCCCCTCGAATGACTTCTACAACAACGTTTCGACTTTCAGCGTCACCATCGAGAACGCCACATGCTGTAATATCGAGCAAAAGACCCCGAACACGACCCCGGCTCAGACCACGAta     <  2:235330/149‑2 (MQ=11)
                                                                                                                                              cAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAcc  >  6:69902/1‑148 (MQ=38)
                                                                                                                                              cAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAAAGGTTCGACTGTCAGCGTCACCATCGAGAAAGCCACACGCGGGAATTTCGAGCAACTGAAACCGAGCCAGAACCCGGCTCAGACCAcagatcgg  >  8:163566/1‑141 (MQ=11)
                                                                                                                                                  |                                                                                                                                                
TGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACC  >  NC_002947/197428‑197716

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