Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 197,572 G→A 17.0% 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,5720GA17.0% 233.6 / 37.0 107K1026K (AAG→AAAPP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base G (43/45);  new base A (9/9);  total (53/54)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.23e-01

Only 100 of 107 total aligned reads displayed.

GCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCG‑‑TTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACC  >  NC_002947/197429‑197716
                                                                                                                                                  |                                                                                                                                                
gCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCG‑‑TTGAAGCCGGCAAGAcc                                                                                                                                               >  8:35516/1‑147 (MQ=37)
 ccaccaCCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCG‑‑TTGAAGCCGGCAAGACCCa                                                                                                                                             <  6:231365/148‑1 (MQ=37)
  caccacCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCG‑‑TTGAAGCCGGCAAGACCCAgg                                                                                                                                           >  3:138349/1‑149 (MQ=37)
   accaccGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCATGTCCAACGGCCAAACCATTACCG‑‑TTGAAGCCGGCAAGACCCAggg                                                                                                                                          <  1:233061/149‑1 (MQ=37)
     caccGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGAGCCTGTCCAACGGCCAAACCATTACCG‑‑TTGAAGCCGGCAAGACCCAggg                                                                                                                                          >  7:170187/1‑147 (MQ=37)
        cgcgACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCG‑‑TTGAAGCCGGAAAGCACCAGGGCAGCg                                                                                                                                     >  6:263612/1‑149 (MQ=255)
          cgACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCG‑‑TTGAAGCCGGCAAGACCCAGGGCAGCGTc                                                                                                                                   <  3:98470/149‑1 (MQ=37)
                tGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCG‑‑TTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTc                                                                                                                             >  7:262446/1‑149 (MQ=37)
                 gACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCG‑‑TTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTc                                                                                                                             >  2:96446/1‑148 (MQ=37)
                  aCGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCG‑‑TTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTcc                                                                                                                            <  6:143677/148‑1 (MQ=37)
                  aCGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCG‑‑TTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTcc                                                                                                                            <  2:288094/148‑1 (MQ=37)
                  aCGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCG‑‑TTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCa                                                                                                                           <  4:133431/149‑1 (MQ=37)
                     gCGAGCCCTTCGGTCACCGAAGGCGGCGTGATCCCCTACACCGTGATCCTGAGCAATGCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCG‑‑TTGAAGCCGTCAAGACCCAGGGCAGCGTCGATTTCCAGAc                                                                                                                        <  3:191126/149‑1 (MQ=255)
                       gAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCG‑‑TTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGAccc                                                                                                                      >  6:313208/1‑149 (MQ=37)
                           ccGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCG‑‑TTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTc                                                                                                                             >  6:285230/1‑138 (MQ=37)
                           ccGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCG‑‑TTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTc                                                                                                                             <  5:285230/138‑1 (MQ=37)
                                ggTCACCTAAGGCGGCGTGATCACCTACACCGTTATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCG‑‑TTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATg                                                                                                             <  5:263612/149‑1 (MQ=255)
                                ggTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCG‑‑TTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAAt                                                                                                              <  6:24330/148‑1 (MQ=37)
                                  tCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCG‑‑TTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGAc                                                                                                           <  5:304417/149‑1 (MQ=37)
                                    aCCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGTCGCCGGTAACCGTGACCCTGTCCAACGGCCAAAGCATTACCG‑‑TTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACg                                                                                                          <  7:261282/148‑1 (MQ=255)
                                      cGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACAATGCCCAACGGAAAAAAAATAACCT‑‑TTGAAGCAGTCAAGAACAAGGGCAGCGTCGATTTCAAGAAAAAGGAGAATGACGTc                                                                                                        >  4:252330/1‑148 (MQ=2)
                                      cGAAGGCGGCGTAATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCG‑‑TTGAAGCCGGCAAGACCCa                                                                                                                                             <  2:162476/111‑1 (MQ=12)
                                      cGAAGGCGGCGTAATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACAG‑‑TTGAAGCCGGCAAGACCCa                                                                                                                                             >  1:162476/1‑111 (MQ=255)
                                         aGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCG‑‑TTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTac                                                                                                     >  5:301524/1‑148 (MQ=37)
                                           gcggcgTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCG‑‑TTGAAGCCGGCAAGACCCAGGGCAGCGTCGAtt                                                                                                                               <  7:203374/120‑1 (MQ=39)
                                           gcggcgTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCAATACCG‑‑TTGAAGCCGGCAAGACCCAGGGCAGCGTCGAtt                                                                                                                               >  8:203374/1‑120 (MQ=255)
                                               cgTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGCCATGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTAcaacaa                                                                                                <  1:96446/149‑1 (MQ=255)
                                               cgTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCG‑‑TTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACg                                                                                              <  5:82951/149‑1 (MQ=37)
                                               cgTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCG‑‑TTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACg                                                                                              <  6:223122/149‑1 (MQ=37)
                                                 tGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCG‑‑TTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGt                                                                                            <  2:96417/149‑1 (MQ=37)
                                                   aTCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCG‑‑TTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTc                                                                                          <  7:69762/149‑1 (MQ=37)
                                                     cACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGGTGACCCTGTCCAACGGCCAAACCATTACCG‑‑TTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGtt                                                                                           >  1:71202/1‑147 (MQ=255)
                                                     cACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCG‑‑TTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGtt                                                                                           <  2:71202/146‑1 (MQ=37)
                                                      aCCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCG‑‑TTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGa                                                                                        <  8:170187/148‑1 (MQ=37)
                                                            acCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCG‑‑TTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAg                                                                                 >  7:76663/1‑149 (MQ=21)
                                                              cGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCG‑‑TTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGtt                                                                                           <  8:113270/137‑1 (MQ=37)
                                                              cGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCG‑‑TTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGtt                                                                                           >  7:113270/1‑137 (MQ=37)
                                                                     cTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCG‑‑TTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGAc                                                                                                           <  6:110236/114‑1 (MQ=35)
                                                                     cTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCG‑‑TTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGAc                                                                                                           >  5:110236/1‑114 (MQ=35)
                                                                      tgcgCAACCCTGCCCAGACGCCGGTGCCCGTGACCCTGTCCAACGGCCAAACGATTACCG‑‑TTGAAGCCGGCAACCCCCAGGGCAGCGTCGACTTCCCGCCCCCGGCGCATGACGGCTACAACAACGGTTCGCCTGTCAGCGGCACCATc                                                                       >  1:193594/4‑149 (MQ=11)
                                                                        aGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCG‑‑TTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCg                                                                      <  1:225331/148‑1 (MQ=17)
                                                                        aGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCG‑‑TTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTACAGCCCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCCCCCTCga                                                                     >  6:24501/1‑149 (MQ=11)
                                                                         gCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCG‑‑TTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCgag                                                                    <  5:313145/149‑1 (MQ=37)
                                                                             tCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCG‑‑TTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACg                                                                <  6:240000/149‑1 (MQ=37)
                                                                              ccggCCCAGACGCCGGTACCCGTGACCCTGTCCAACGGCCAAACCATTACCG‑‑TTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGc                                                               <  4:81570/146‑1 (MQ=255)
                                                                               cTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCG‑‑TTGAAGCAGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCAGGCGAATGACGTCTACAACAACGGTTCTACTGTCAGCGTCACCAACGAGAACGc                                                               >  3:191321/1‑148 (MQ=255)
                                                                                  cccAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCG‑‑TTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCcac                                                            <  7:54430/148‑1 (MQ=37)
                                                                                  cccAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCG‑‑TTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCca                                                             >  6:282011/1‑147 (MQ=17)
                                                                                  cccAGACGCCGGTAACCGTGACCCTATCCAACGGCCAAACCATTACCG‑‑TTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCca                                                             >  2:142101/1‑147 (MQ=17)
                                                                                   ccAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCG‑‑TTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAg                                                          >  2:11405/1‑149 (MQ=37)
                                                                                    cAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCG‑‑TTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAgg                                                         <  8:262446/149‑1 (MQ=37)
                                                                                     aGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCG‑‑TTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGc                                                                                >  4:181546/1‑125 (MQ=17)
                                                                                     aGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCG‑‑TTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGc                                                                                <  3:181546/125‑1 (MQ=17)
                                                                                      gACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCG‑‑TTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCg                                                       <  4:138349/149‑1 (MQ=37)
                                                                                       aCGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCG‑‑TTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCgg                                                      >  2:183455/1‑149 (MQ=37)
                                                                                        cGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCG‑‑TTGAAGCCGGCAAGACCCa                                                                                                                                             >  7:235504/1‑61 (MQ=40)
                                                                                        cGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCG‑‑TTGAAGCCGGCAAGACCCa                                                                                                                                             <  8:235504/61‑1 (MQ=40)
                                                                                         gCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCG‑‑TTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGt                                                     >  5:311867/1‑148 (MQ=255)
                                                                                             gTAACCGTGACCCTGTCCAACGGCCAAACCATTACCG‑‑TTGAAGCCGGCAAGACCCATGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCcaca                                                           <  5:272950/138‑1 (MQ=37)
                                                                                             gTAACCGTGACCCTGTCCAACGGCCAAACCATTACCG‑‑TTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTAcaacaa                                                                                                >  5:178793/1‑101 (MQ=35)
                                                                                             gTAACCGTGACCCTGTCCAACGGCCAAACCATTACCG‑‑TTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTAcaacaa                                                                                                <  6:178793/101‑1 (MQ=35)
                                                                                             gTAACCGTGACCCTGTCCAACGGCCAAACCATTACCG‑‑TTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCcaca                                                           >  6:272950/1‑138 (MQ=37)
                                                                                              tgaCCGTGACCCTGTCCAACGGCCAAACGATTACCG‑‑TTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTc                                               <  8:3833/147‑1 (MQ=255)
                                                                                              tgaCCGTGACCCTGTCCAACGGCCAAACGATTACCG‑‑TTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTc                                               <  6:222263/147‑1 (MQ=255)
                                                                                              tAACCGTGACCCTGTCCAACGGCCAAACCATTACCG‑‑TTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTc                                               >  8:136955/1‑149 (MQ=35)
                                                                                                aCCGTGACCCTGTCCAACGGCCAAACCATTACCG‑‑TTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCg                                              <  1:11405/148‑1 (MQ=34)
                                                                                                 ccGTGACCCTGTCCAACGGCCAAACCATTACCG‑‑TTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAg                                            >  2:76852/1‑149 (MQ=33)
                                                                                                 ccGTGACCCTGTCCAACGGCCAAACCATTACCG‑‑TTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAg                                            >  4:277961/1‑149 (MQ=33)
                                                                                                  cGTGACCCTGTCCAACGGCCAAACGATTACCG‑‑TTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGc                                           <  6:161494/149‑1 (MQ=12)
                                                                                                  cGTGACCCTGTCCAACGGCCAAACCATTACCG‑‑TTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTaa                                                   <  5:108699/141‑1 (MQ=32)
                                                                                                  cGTGACCCTGTCCAACGGCCAAACCATTACCG‑‑TTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTaa                                                   >  6:108699/1‑141 (MQ=32)
                                                                                                       aCCCTGTCCAACGGCCAAACGATTACCG‑‑TTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAAc                                        <  5:201017/148‑1 (MQ=12)
                                                                                                         ccTGTCCAACGGCCAAACGATTATCG‑‑TTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCg                                                       <  8:172302/131‑1 (MQ=255)
                                                                                                         ccTGTCCAACGGCCAAACGATTATCG‑‑TTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCg                                                       >  7:172302/1‑131 (MQ=255)
                                                                                                           tGTCCAACGGCCAAACGATTACCG‑‑TTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGAcc                                   >  6:90396/1‑149 (MQ=12)
                                                                                                                aaCGGCCAAACCATTACCG‑‑TTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAAc                                                                 >  1:284259/1‑114 (MQ=32)
                                                                                                                   gtccAAACCAATACCG‑‑TTGAAGCCTGCAAGACCCAGGGCATCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGTCGGTAATTTCGAGCAACTGAcc                                   <  5:68667/139‑1 (MQ=255)
                                                                                                                   ggCCAAACGATTACCG‑‑TTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGGTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAAccc                           <  8:293950/149‑1 (MQ=255)
                                                                                                                   ggCCAAACCATTACAG‑‑TTGAAGACGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCTACTGTCAGCGTCACCATCGAGAACGCAACAGGCGGTAATTTCGAGCACCTGAcc                                   >  6:68667/1‑141 (MQ=11)
                                                                                                                    gCCAAACCATTACCG‑‑TTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACt                                                                                      <  6:34202/89‑1 (MQ=31)
                                                                                                                    gCCAAACCATTACCG‑‑TTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACt                                                                                      >  5:34202/1‑89 (MQ=31)
                                                                                                                     ccAAACGATTACCG‑‑TTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGaa                              <  5:304691/144‑1 (MQ=12)
                                                                                                                       aaaCGATTACCG‑‑TTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGAcc                       >  2:55759/1‑148 (MQ=255)
                                                                                                                       aaaCCATTACCG‑‑TTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGATTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCACCTGACCCCGAAAccgccg                       >  4:69192/1‑149 (MQ=11)
                                                                                                                           cATTACCG‑‑TTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCgg                   <  5:313208/149‑1 (MQ=37)
                                                                                                                            aTTACCG‑‑TTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACCCCGGc                  <  5:90396/149‑1 (MQ=12)
                                                                                                                             ttACCG‑‑TTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCTAATGACGTCTACAACAACTGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCt                 <  3:252330/149‑1 (MQ=255)
                                                                                                                             ttACCG‑‑TTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCt                 >  2:235435/1‑149 (MQ=37)
                                                                                                                              tACCG‑‑TTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACCCCGGCTc                >  6:138188/1‑149 (MQ=12)
                                                                                                                                ccG‑‑TTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAg              <  7:136955/149‑1 (MQ=37)
                                                                                                                                  g‑‑TTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGAc            >  7:205260/1‑149 (MQ=37)
                                                                                                                                      tGAAGCCGTCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGTCGAAGGACGTCTCCAACACCGGTTCGACTGGCAGCGTCACCATCTAGAACGCCACAGGCGGTAATTTCGAGCAACGGACCCCGAACCCGACCCCGGCTCAGACCa          <  2:193594/149‑1 (MQ=11)
                                                                                                                                        aaGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCAc         >  4:275266/1‑148 (MQ=37)
                                                                                                                                          gCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGAt      >  1:56761/1‑149 (MQ=37)
                                                                                                                                           ccGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATc     <  8:205260/149‑1 (MQ=37)
                                                                                                                                            cGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCCACAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGCACCCGCCGCCGCCTCAGACCACGATCa    >  3:24740/1‑149 (MQ=11)
                                                                                                                                            cGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGGAATTTCGAGCAACTGACCCCGAACCCGACACCGGCTCAGACCACGATc     >  4:241251/1‑148 (MQ=255)
                                                                                                                                             ggCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAc   >  2:292717/1‑149 (MQ=38)
                                                                                                                                              gCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAcc  <  3:277961/149‑1 (MQ=38)
                                                                                                                                              gCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACCCCGGCTCAGACCACGATCaac  >  7:264608/1‑147 (MQ=2)
                                                                                                                                                  |                                                                                                                                                
GCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCG‑‑TTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACC  >  NC_002947/197429‑197716

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: