Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,973,308:1 +T 7.3% coding (609/1713 nt) PP_2600 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,973,3081.T 7.3% 374.9 / 22.5 96coding (609/1713 nt)PP_2600hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (41/48);  new base T (4/3);  total (45/51)
Fisher's exact test for biased strand distribution p-value = 7.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.60e-01

GCGACCCACGCCCCGCGCACCTGGCCACCCTCGATCAGCAGGCGTTTGGCCGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGG‑CGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCGGCACCGGCCATGATCGGCATGCCCATGAACGAGGTTTCGCGCATGGTCTT  >  NC_002947/2973159‑2973451
                                                                                                                                                      |                                                                                                                                               
gcgACCCACGCCCCGCGCACCTGGCCACCCTCGATCAGCAGGCGTTTGGCCGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATgg                                                                                                                                                  >  1:854530/1‑150 (MQ=255)
  gACCCACGCCCCGCGCACCTGGCCACCCTCGATCAGCAGGCGTTTGGCCGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGG‑Cg                                                                                                                                               >  1:334444/1‑150 (MQ=255)
  gACCCACGCCCCGCGCACCTGGCCACCCTCGATCAGCAGGCGTTTGGCCGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGG‑Cg                                                                                                                                               <  2:334444/150‑1 (MQ=255)
   aCCCACGCCCCGCGCACCTGGCCACCCTCGATCAGCAGGCGTTTGGCCGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGG‑c                                                                                                                                                >  1:1704852/1‑148 (MQ=255)
   aCCCACGCCCCGCGCACCTGGCCACCCTCGATCAGCAGGCGTTTGGCCGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGG‑c                                                                                                                                                <  2:1704852/148‑1 (MQ=255)
   aCCCACGCCCCGCGCACCTGGCCACCCTCGATCAGCAGGCGTTTGGCCGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGG‑Cgg                                                                                                                                              >  1:662425/1‑150 (MQ=255)
     ccACGCCCCGCGCACCTGGCCACCCTCGATCAGCAGGCGTTTGGCCGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGG‑CGGGc                                                                                                                                            >  1:142746/1‑150 (MQ=255)
     ccACGCCCCGCGCACCTGGCCACCCTCGATCAGCAGGCGTNTGGCCGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGG‑CGGGc                                                                                                                                            >  1:1358892/1‑150 (MQ=255)
     ccACGCCCCGCGCACCTGGCCACCCTCGACCAGCAGGCGTTTGGCCGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGG‑CGGGc                                                                                                                                            <  1:1095387/150‑1 (MQ=255)
        ggCCCCGCCCACCTCGCCACCCCCGATCAGCAGGCGTTTGGCCGACGCCGACGCGTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGGGCCACGAGGGCCACACCGTTGACCAGGTGCATGGCCCG‑CCATGGTCGGGCCat                                                                                                                                         <  1:244878/149‑2 (MQ=255)
         gCCCCGCGCACCTGGCCACCCTCGATCAGCAGGCGTTTGGCCGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGG‑CGGGCCAgg                                                                                                                                        <  2:854530/150‑1 (MQ=255)
          ccccGCGCACCTGGCCACCCTCGATCAGCAGGCGTTTGGCCGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCG‑CCATGGTCGGGCCatgc                                                                                                                                       <  2:228453/150‑4 (MQ=255)
           cccGCGCACCTGGCCACCCTCGATCAGCAGGCGTTTGGCCGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGG‑CGGGCCAggtg                                                                                                                                      <  2:972695/150‑1 (MQ=255)
               cgcACCTGGCCACCCTCGATCAGCAGGCGTTTGGCCGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATTG‑CGGGCCAGGTGGTAc                                                                                                                                  <  1:238546/150‑1 (MQ=255)
                         cACCCTCGATCAGCAGGCGTTTGGCCGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGG‑CGGGCCAggtg                                                                                                                                      <  1:861707/136‑1 (MQ=255)
                         cACCCTCGATCAGCAGGCGTTTGGCCGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGG‑CGGGCCAggtg                                                                                                                                      >  2:861707/1‑136 (MQ=255)
                          aCCCTCGATCAGCAGGCGTTTGGCCGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGG‑CGTGCCAGGTGGTACAGGTGGCGGcc                                                                                                                       >  2:933571/1‑150 (MQ=255)
                           cccTCGATCCGCAGGCGTTTGGCCGGCGCCGACGCTTGCAGGTGCACACCCAGTTCCTGTGCCGATCTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGG‑CGGGCCAGGTGGTACAGGTGGCGGCCa                                                                                                                      <  1:429142/150‑1 (MQ=255)
                            ccTCGATCAGCAGGCGTTTGGCCGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGG‑CGGGCCAGGTGGTACAGGTGGCGGCCaa                                                                                                                     >  1:1101524/1‑150 (MQ=255)
                             cTCGATCAGCAGGCGTTTGGCCGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGG‑CGGGCCAGGTGGTACAGGTGGCGGCCaaa                                                                                                                    <  2:1358892/150‑1 (MQ=255)
                             cTCGATCAGCAGGCGTTTGGCCGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGG‑CGGGCCAGGTGGTACAGGTGGCGGCCaaa                                                                                                                    >  1:1471856/1‑150 (MQ=255)
                             cTCGATCAGCAGGCGTTTGGCCGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGG‑CGGGCCAGGTGGTACAGGTGGCGGCCaaa                                                                                                                    <  2:1471856/150‑1 (MQ=255)
                             cTCGATCAGCAGGCGTTTGGCCGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGG‑CGGGCCAGGTGGTACAGGTGGCGGCCaaa                                                                                                                    <  1:534753/150‑1 (MQ=255)
                             cTCGATCAGCAGGCGTTTGGCCGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGG‑CGGGCCAGGTGGTACAGGTGGCGGCCaaa                                                                                                                    >  1:1637665/1‑150 (MQ=255)
                                 ggcagcagGCGTTTGGCCGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGG‑CGGGCCAGGTGGTACAGGTGGCGGCCAAAgcgc                                                                                                                >  1:1036972/3‑150 (MQ=255)
                                     gcagGCGTTTGGCCGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGG‑CGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGc                                                                                                            <  2:318259/150‑1 (MQ=255)
                                         gCGTTTGGCCGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGG‑CGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGTTTGCacac                                                                                                        <  2:662425/150‑1 (MQ=255)
                                           gTTTGGCCGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGG‑CGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGt                                                                                                      <  2:317152/150‑1 (MQ=255)
                                            tttGGCCGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGG‑CGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTg                                                                                                     <  2:1541272/150‑1 (MQ=255)
                                                 ccGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGG‑CGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTgcagca                                                                                                >  2:373725/1‑150 (MQ=255)
                                                  cGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGG‑CGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGcagcag                                                                                               <  2:1509488/150‑1 (MQ=255)
                                                  cGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCACGG‑CGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGcagcag                                                                                               >  1:1786800/1‑150 (MQ=255)
                                                   ggCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGG‑CGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAgg                                                                                              >  1:1924033/1‑150 (MQ=255)
                                                       ccGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGG‑CGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCtt                                                                                          <  2:1482740/150‑1 (MQ=255)
                                                       ccGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGG‑CGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCtt                                                                                          <  1:1355970/150‑1 (MQ=255)
                                                         gACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGG‑CGGGCCAggt                                                                                                                                       <  1:300324/103‑1 (MQ=255)
                                                         gACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGG‑CGGGCCAggt                                                                                                                                       >  2:300324/1‑103 (MQ=255)
                                                              ttGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGG‑CGGGCCAggtgg                                                                                                                                     >  2:1023809/1‑100 (MQ=255)
                                                              ttGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGG‑CGGGCCAggtgg                                                                                                                                     <  1:1023809/100‑1 (MQ=255)
                                                              ttGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGG‑CGGGCCAGGTGGTACAGGTggc                                                                                                                           <  1:1603840/110‑1 (MQ=255)
                                                              ttGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGG‑CGGGCCAGGTGGTACAGGTggc                                                                                                                           >  2:1603840/1‑110 (MQ=255)
                                                              ttGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGG‑CGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGa                                                                                   >  2:162905/1‑150 (MQ=255)
                                                                    gTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGG‑CGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGt                                                                             >  1:470971/1‑150 (MQ=255)
                                                                     tGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGG‑CGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTc                                                                            <  2:1924033/150‑1 (MQ=255)
                                                                      gCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCCGTGCCACCCAGTTGACCAGGTGCATGGCCCGACCATGG‑CGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCa                                                                           >  2:758575/1‑150 (MQ=255)
                                                                         cacCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGG‑CGGGCCAGGTGGTACAGGTGGCGGCCAAGGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGt                                                                        <  1:547189/150‑1 (MQ=255)
                                                                          acCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGG‑CGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCtt                                                                                                              <  2:267090/111‑1 (MQ=255)
                                                                          acCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGG‑CGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCtt                                                                                                              >  1:267090/1‑111 (MQ=255)
                                                                          acCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGG‑CGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGtt                                                                       <  1:680653/150‑1 (MQ=255)
                                                                                 tCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGG‑CGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAg                                                                <  1:373725/150‑1 (MQ=255)
                                                                                  ccTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGG‑CGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAgg                                                               <  1:933571/150‑1 (MQ=255)
                                                                                    tgtgCCGATTTAGCCAGGTGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGG‑CGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAggcg                                                             <  1:125101/150‑1 (MQ=255)
                                                                                    tgtgCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGG‑CGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAggcg                                                             <  2:1101524/150‑1 (MQ=255)
                                                                                      tgCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGG‑CGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAggcggc                                                           >  1:302326/1‑150 (MQ=255)
                                                                                       gCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGG‑CGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCa                                                                           <  1:1530344/133‑1 (MQ=255)
                                                                                       gCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGG‑CGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCa                                                                           >  2:1530344/1‑133 (MQ=255)
                                                                                        ccGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGG‑CGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCacac                                                                                                        >  2:1862569/1‑103 (MQ=255)
                                                                                        ccGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGG‑CGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCacac                                                                                                        <  1:1862569/103‑1 (MQ=255)
                                                                                         cGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGG‑CGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAg                                                        <  2:770750/150‑1 (MQ=255)
                                                                                                 ccAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGG‑CGGGCCAGGTGGTAAAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGGGGGGGTCATGTTGCACAAGGCGGCAAGGTCGGCAc                                                >  1:591213/1‑150 (MQ=255)
                                                                                                  cAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGG‑CGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCGGCAcc                                               >  1:1121070/1‑150 (MQ=255)
                                                                                                   aGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGG‑CGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCGGCACCg                                              >  2:1864937/1‑150 (MQ=255)
                                                                                                   aGGCGCGCCACCACTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGG‑CGGGCCCGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCGGCACCg                                              >  2:96620/1‑150 (MQ=255)
                                                                                                     gcgcgcCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGG‑CGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTTGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCGGCACCGGc                                            <  1:495757/150‑1 (MQ=255)
                                                                                                        cgcCACCAGTGCCCCACCGTTGACCAGGTGCATGGCCCGACCATGG‑CGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCGGCACCGGCCAt                                         <  1:1046878/150‑1 (MQ=255)
                                                                                                            accaGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGG‑CGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCGGCACCGGCCATGATc                                     >  2:162559/1‑150 (MQ=255)
                                                                                                                   ccACACCGTTGACCAGGTGCATGGCCCGACCATGG‑CGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCGGCACCGGCCATGATCGGCATGc                              <  2:142746/150‑1 (MQ=255)
                                                                                                                     acacCGTTGACCAGGTGCATGGCCCGACCATGG‑CGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATg                                                                         <  1:842382/105‑1 (MQ=255)
                                                                                                                     acacCGTTGACCAGGTGCATGGCCCGACCATGG‑CGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATg                                                                         >  2:842382/1‑105 (MQ=255)
                                                                                                                     acacCGTTGACCAGGTGCATGGCCCGACCATGG‑CGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAggcg                                                             <  2:151760/117‑1 (MQ=255)
                                                                                                                     acacCGTTGACCAGGTGCATGGCCCGACCATGG‑CGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAggcc                                                             >  1:151760/1‑116 (MQ=255)
                                                                                                                     acacCGTTGACCAGGTGCATGGCCCGACCATGG‑CGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCGGCACCGGCCATGATCGGCATGccc                            >  2:1117700/1‑150 (MQ=255)
                                                                                                                       acCGTTGACCAGGTGCATGGCCCGACCATGG‑CGGGCCAGGTGGTACAGGTGGAGGCCTAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCGGCACCGGCCATGATCGGCATGCCCAt                          <  2:1717759/150‑1 (MQ=255)
                                                                                                                         cGTTGACCAGGTGCATGGCCCGACCATGG‑CGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCGGCACCGGCCATGATCGGCATGCCCATGa                        <  1:162559/150‑1 (MQ=255)
                                                                                                                             gACCAGGTGCATGGCCCGACCATGG‑CGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCGGCACCGGCCATGATCGGCATGCCCATGAACGa                    <  1:872522/150‑1 (MQ=255)
                                                                                                                             gACCAGGTGCATGGCCCGACCATGG‑CGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCGGCACCGGCCATGATCGGCATGCCCATGAACGa                    <  2:302326/150‑1 (MQ=255)
                                                                                                                             gACCAGGTGCATGGCCCGACCATGG‑CGGGCCAGGTCGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGCGCGGGTCATGTTGAGGAAGGGGGCAAGGTCGGCACCGGCCATGATCGGCATGCCCATGAACGa                    <  2:1036972/150‑1 (MQ=255)
                                                                                                                              aCCAGGTGCATGGCCCGACCATGG‑CGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCGGCACCGGCCATGATCGGCATGCCCATGAACGAg                   >  1:107166/1‑150 (MQ=255)
                                                                                                                              aCCAGGTGCATGGCCCGACCATGG‑CGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCGGCACCGGCCATGATCGGCATGCCCATGAACGAg                   <  1:758575/150‑1 (MQ=255)
                                                                                                                              aCCAGGTGCATGGCCCGACCATGG‑CGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCGGCACCGGCCATGATCGGCATGCCCATGAACGAg                   <  2:107166/150‑1 (MQ=255)
                                                                                                                               ccAGGTGCATGGCCCGACCATGG‑CGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCGGCACCGGCCATGATCGGCATGCCCATGAACGAgg                  >  2:53334/1‑150 (MQ=255)
                                                                                                                               ccAGGTGCATGGCCCGACCATGG‑CGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCGGCACCGGCCATGATCGGCATGCCCATGAACGAgg                  >  1:539339/1‑150 (MQ=255)
                                                                                                                               ccAGGTGCATGGCCCGACCATGG‑CGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCGGCACCGGCCATGATCGGCATGCCCATGAACGAgg                  >  1:636271/1‑150 (MQ=255)
                                                                                                                               ccAGGTGCATGGCCCGACCATGG‑CGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCGGCACCGGCCATGATCGGCATGCCAATGAACGAgg                  >  1:148014/1‑150 (MQ=255)
                                                                                                                                   gTGCATGGCCCGACCATGG‑CGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAggcggc                                                           >  1:1346800/1‑105 (MQ=255)
                                                                                                                                   gTGCATGGCCCGACCATGG‑CGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAggcggc                                                           <  2:1346800/105‑1 (MQ=255)
                                                                                                                                    cacctGGCCCG‑CCATGGTCGGGCCATGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCGGCACCGGCCATGATc                                     <  1:174652/122‑1 (MQ=255)
                                                                                                                                    cacctGGCCCG‑CCATGGTCGGGCCATGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCGGCACCGGCCATGATc                                     >  2:174652/5‑126 (MQ=255)
                                                                                                                                    cacctGGCCCG‑CCATGGTCGGGCCATGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCGGCACCGGCCATGATCGGCATGCCCATGAACGAGGTTTcg             >  1:1222733/5‑150 (MQ=255)
                                                                                                                                      cctGGCCCG‑CCATGGTCGGGCCAGGTGGTACAGGTGGCGGCCAACGCGCTTGCCCGCGTGCAGCAGGGCGTGTGGCGAGCGTGTCATGTTGAGGAAGGCGGCAAGGTCGGCACCGGCCATGATCGGCATGCCCATGAACGAGGTTTcgcg           >  1:1206697/3‑150 (MQ=255)
                                                                                                                                       ctGGCCCG‑CCATGGTCGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCGGCACCGGCCATGATCGGCATGCCCATGAACGAGGTTTCgcgc          >  2:320757/2‑150 (MQ=255)
                                                                                                                                       aTGGCCCGACCATGG‑CGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCGGCACCGGCCATGATCGGc                                  >  1:1670706/1‑126 (MQ=255)
                                                                                                                                       aTGGCCCGACCATGG‑CGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCGGCACCGGCCATGATCGGc                                  <  2:1670706/126‑1 (MQ=255)
                                                                                                                                          gCCCTACCATGG‑CGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCGGCATCGGCCATGATCGGCATGCCCATGAACTAGGTTTCGCGCATg       <  1:859173/150‑1 (MQ=255)
                                                                                                                                           cccGACCATGG‑CGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGCGAACTATCAATGATAGCCGCCTTGACGTGTCCGGCCATGATCGGCATGCCCATGAACGAGGTTTCGCGCATgg      >  1:1684017/1‑150 (MQ=255)
                                                                                                                                            ccGACCATGG‑CGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCGGCACCGGCCATGATCGGCATGCCCATGAACGAGGTTTCGCGCATGGt     >  1:376274/1‑150 (MQ=255)
                                                                                                                                               aCCATGG‑CGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCGGCACCGGCCATGATCGGCATGCCCATGAACGAGGTTTCGCGCATGGTCtt  <  2:376274/150‑1 (MQ=255)
                                                                                                                                                   tGG‑CGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCt                                                                                                               >  2:904430/1‑37 (MQ=255)
                                                                                                                                                   tGG‑CGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCt                                                                                                               <  1:904430/37‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                               
GCGACCCACGCCCCGCGCACCTGGCCACCCTCGATCAGCAGGCGTTTGGCCGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGG‑CGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCGGCACCGGCCATGATCGGCATGCCCATGAACGAGGTTTCGCGCATGGTCTT  >  NC_002947/2973159‑2973451

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