Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,908,158 C→T 6.9% intergenic (‑89/‑193) PP_4317 ← / → PP_4318 hypothetical protein/transposase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,908,1580CT 6.9% 228.1 / 6.1 87intergenic (‑89/‑193)PP_4317/PP_4318hypothetical protein/transposase
Reads supporting (aligned to +/- strand):  ref base C (40/41);  new base T (2/4);  total (42/45)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01

CTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGAGCGATGTGGG  >  NC_002947/4908010‑4908306
                                                                                                                                                    |                                                                                                                                                    
cttgatcttCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATAGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGagcgca                                                                                                                                                     >  2:1171996/1‑150 (MQ=17)
       ttCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCGGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAgg                                                                                                                                              <  2:447460/150‑1 (MQ=17)
       ttCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAgg                                                                                                                                              >  2:475176/1‑150 (MQ=32)
          gcgACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGc                                                                                                                                           >  1:1082499/1‑150 (MQ=32)
          gcgACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGc                                                                                                                                           >  1:9575/1‑150 (MQ=32)
             aCTTCAGGAGGCCGAACGCAGGCCTTGCCGAGGCAGCTGTCGAGAGCTACTTGCTGCAAAGCTGTCAGGCGCGATTATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCgg                                                                                                                                        >  1:348515/1‑150 (MQ=2)
                 cAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGNGAGGCCNGAGTGAGGGGACNCCGGAGCGCAGCGCAGGGGCCGgatga                                                                                                                                    <  1:1171996/150‑1 (MQ=17)
                  aggaggCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAgg                                                                                                                                              >  1:493106/1‑139 (MQ=33)
                  aggaggCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAgg                                                                                                                                              <  2:493106/139‑1 (MQ=33)
                   ggaggCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATACCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGatgatg                                                                                                                                  >  1:935734/1‑150 (MQ=17)
                      ggCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGa                                                                                                                               <  2:9575/150‑1 (MQ=32)
                       gCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAg                                                                                                                              <  2:1095629/150‑1 (MQ=32)
                       gCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAg                                                                                                                              <  1:475176/150‑1 (MQ=32)
                         cGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAgcg                                                                                                                            >  2:536947/1‑150 (MQ=32)
                         cGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAgcg                                                                                                                            >  1:973949/1‑150 (MQ=32)
                           aaCGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCa                                                                                                                          <  2:1016825/150‑1 (MQ=32)
                              gCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGGCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGCTGGGGGCGCAGCg                                                                                                                       >  1:429397/1‑150 (MQ=12)
                                aGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGCCACAAGAGCCCAGCCCATGGGCCGGATGATGGGAGCGCAGCGtt                                                                                                                     >  2:782687/1‑150 (MQ=11)
                                    cttgcGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGa                                                                                                                               <  2:397713‑M2/131‑1 (MQ=255)
                                    cttgcGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGa                                                                                                                               >  1:397713‑M2/6‑136 (MQ=255)
                                    cttgcGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGtttttt                                                                                                                 <  1:536947‑M2/145‑1 (MQ=255)
                                    cttgcGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGtttttt                                                                                                                 <  2:1152546‑M2/145‑1 (MQ=255)
                                     ttgcGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGACCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTg                                                                                                                >  2:268284‑M2/5‑150 (MQ=255)
                                       gcGGAGGCAGGTGACGGGCATGTATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGt                                                                                                              <  2:348515‑M2/148‑1 (MQ=255)
                                         ggAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTa                                                                                                            >  2:1043838/1‑150 (MQ=32)
                                                 ttGACGGGCATCCATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGt                                                                                                    >  2:614572/2‑150 (MQ=12)
                                                   gACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTcg                                                                                                  >  2:88661/1‑150 (MQ=32)
                                                   gACGGGCATCCATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTcg                                                                                                  >  1:1129944/1‑150 (MQ=12)
                                                      gggCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGt                                                                                               >  1:709057/1‑150 (MQ=32)
                                                       ggcatccaTGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGtt                                                                                              <  2:1129944/143‑1 (MQ=17)
                                                        gCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAgcgcc                                                                                                                          >  1:643455/1‑120 (MQ=18)
                                                        gCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCa                                                                                                                          <  2:643455/121‑1 (MQ=33)
                                                          aTGGATGCCCGTCAAGCGCTGAGGCCCGAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGa                                                                                           <  2:206468/150‑1 (MQ=17)
                                                          aTGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGa                                                                                           <  2:973949/150‑1 (MQ=32)
                                                            ggATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACa                                                                                         >  2:624217/1‑150 (MQ=32)
                                                                  ccGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGt                                                                                   <  2:1082499/150‑1 (MQ=32)
                                                                  ccGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAACGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGt                                                                                   >  2:568791/1‑150 (MQ=17)
                                                                   cGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTg                                                                                  >  2:126746/1‑150 (MQ=32)
                                                                   cGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTATTGCGTTTGACCAAAAGTg                                                                                  >  2:713255/1‑150 (MQ=12)
                                                                       aaGCGCTGAGGCCCGAGGATGGGACCTTCGGCGCGGTCCTGCCGGGAGCGGGGCCGGAGAGAGGGCACCCCGGTGCGCAGCGCGGGGGCGGGATGATTGGAGCGGCGCGTTTGTTGGGTACTTTTTGTCGCGTTTGACAATAAGTAACtt                                                                              >  2:486228/1‑149 (MQ=11)
                                                                              gAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTaa                                                                       <  2:206569/150‑1 (MQ=32)
                                                                                 gCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAggg                                                                    >  2:1269279/1‑150 (MQ=32)
                                                                                   cccAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGttt                                                                                             <  1:54234/123‑1 (MQ=33)
                                                                                   cccAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGttt                                                                                             >  2:54234/1‑123 (MQ=33)
                                                                                     cAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCCAt                                                                >  2:254282/1‑150 (MQ=32)
                                                                                        gATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCCATccc                                                             >  1:864834/1‑150 (MQ=32)
                                                                                               ccTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGTAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCAGGTCTATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAg                                                      <  1:1043838/150‑1 (MQ=11)
                                                                                               ccTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAg                                                      <  1:624217/150‑1 (MQ=32)
                                                                                               ccTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAg                                                      >  1:712540/1‑150 (MQ=32)
                                                                                                cTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCCa                                                                 <  1:996172/138‑1 (MQ=33)
                                                                                                cTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCCa                                                                 >  2:996172/1‑138 (MQ=33)
                                                                                                     gcgcGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGt                                                                                                    >  2:1148568/1‑98 (MQ=34)
                                                                                                     gcgcGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGt                                                                                                    <  1:1148568/98‑1 (MQ=34)
                                                                                                         ggTCCTGCCGGGAGCGAGGCCGGAGTGAGGTGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTa                                            >  2:519861/1‑150 (MQ=37)
                                                                                                               gCCGGGAGCGAGGCCGGAGTGAGGTGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTAGAATTc                                      <  1:519861/150‑1 (MQ=255)
                                                                                                                      gCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGaaa                               <  1:268284/150‑1 (MQ=255)
                                                                                                                      gCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGaaa                               <  1:713255/150‑1 (MQ=255)
                                                                                                                        gAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATc                             >  1:277689/1‑150 (MQ=255)
                                                                                                                         aGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCccctcactccggccttgctcccggcaggaccgcgccgaaggccccatcctggggcctcagcgcttgacgggcatccatgcccgtcacctgcctccgcaaggcctgcgttc                            <  1:842914‑M2/150‑111 (MQ=255)
                                                                                                                          ggCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCa                           <  1:1216301/150‑1 (MQ=255)
                                                                                                                          ggCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCa                           <  1:568791/150‑1 (MQ=255)
                                                                                                                             cGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCCACCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTc                        >  1:272430/1‑150 (MQ=255)
                                                                                                                               gAGTGAGGGGACCCCGGAGCGCAGCGGAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAg                      <  2:277689/150‑1 (MQ=255)
                                                                                                                               gAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCccctcactccggccttgctcccggcaggaccgcgccgaaggccccatcctggggcctcagcgctt                                                                         >  2:946258‑M2/1‑34 (MQ=255)
                                                                                                                               gAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCccctcactccggccttgctcccggcaggaccgcgccgaaggccccatcctggggcctcagcgctt                                                                         <  1:946258‑M2/99‑66 (MQ=255)
                                                                                                                                 gTGAGGGGACCCCGGAGCGTAGCGTAGGGGTcccctcactccgggcttgctcccgggaggaccgcgctgcaggccccatcctggggcccagcgcttgacgggcatccatgcccgtcacctccctccgcaagccctgcgctcggcctcctg                    <  2:934640‑M2/150‑120 (MQ=255)
                                                                                                                                 gTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGa                    <  1:614572/150‑1 (MQ=255)
                                                                                                                                  tGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGaa                   >  1:1243807/1‑150 (MQ=255)
                                                                                                                                    aGGGGACCCCGGAGCGTAGCGTAGGGGCcccctcactccgggcttgctcccgggaggaccgcgctgcaggccccatcctggggcccagcgcttgacgggcatccatgcccgtcacctccctccgcaagccctgcgctcggcctcctgaag                 <  1:228133‑M2/150‑123 (MQ=255)
                                                                                                                                     ggggACCCCGGAGCGCAGCGCAGGGGTCccctcactccggccttgctcccggcaggaccgcgccgaaggccccatcctggggcctcagcgcttgacgggcatccatgcccgtcacctgcctccgcaaggcctgcgttcggcctcctgaag                <  2:504682‑M2/150‑123 (MQ=255)
                                                                                                                                     ggggACCCCGGAGCGCAGCGCAGGGGCCccctcactccggcctcgctcccggcaggaccgcgccgaaggccccatcctggggcctcagcgcttgacgggcatccatgcccgtcacctgcctccgcaaggcctgcgttcggcctcctgaag                <  2:107019‑M2/150‑123 (MQ=255)
                                                                                                                                     ggggACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCa                           >  2:442402/1‑139 (MQ=40)
                                                                                                                                     ggggACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCa                           <  1:442402/139‑1 (MQ=40)
                                                                                                                                     ggggACCCCGGAGCGCAGCGCAGGGGCCGGATGACGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGAc                >  2:564457/1‑150 (MQ=255)
                                                                                                                                      gggACCCCGGAGCGTAGCGTAGGGGTcccctcactccgggcttgctcccgggaggaccgcgctgcaggccccatcctggggcccagcgcttgacgggcatccatgcc                                                          <  1:181519‑M2/107‑82 (MQ=255)
                                                                                                                                      gggACCCCGGAGCGTAGCGTAGGGGTcccctcactccgggcttgctcccgggaggaccgcgctgcaggccccatcctggggcccagcgcttgacgggcatccatgcc                                                          >  2:181519‑M2/1‑26 (MQ=255)
                                                                                                                                      gggACCCCGGAGCGCAGCGCAGGGGCCccctcactccggcctcgctcccggcaggaccgcgccgaaggccccatcctggggcctcagcgcttgacgggcatccatgcccgtcacctgcctccgcaaggcctgcgttcggcctcctgaagt               <  2:9860‑M2/150‑124 (MQ=255)
                                                                                                                                      gggACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCCACCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGAcc               <  2:272430/150‑1 (MQ=255)
                                                                                                                                      ggaACCCCGGAGCGCAGCGCAGGGGTTccctcactccggccttgctcccggcaggaccgcgccgaaggccccatcctggggcctcagcgcttgacgggcatccatgcctgtcacctgcctccgcaaggcctgcgttcggcctcctgaagt               <  2:165964‑M2/150‑124 (MQ=255)
                                                                                                                                      ggaACCCCGGAGCGCAGCGCAGGGGTTccctcactccggccttgctcccggcaggaccgcgccgaaggccccatcctggggcctcagcgcttgacgg                                                                    >  1:1250978‑M2/1‑27 (MQ=255)
                                                                                                                                      ggaACCCCGGAGCGCAGCGCAGGGGTTccctcactccggccttgctcccggcaggaccgcgccgaaggccccatcctggggcctcagcgcttgacgg                                                                    <  2:1250978‑M2/97‑71 (MQ=255)
                                                                                                                                       ggACCCCGGAGCGCAGCGCAGGGGCCccctcactccggccttgctcccggtaggacggggccgaaggccccatcctggggcctcagcgcttgacgggcatccatgcccgtcacctgcctccgcaaggcctgcgttcggccgcctgaagtc              <  1:837831‑M2/150‑125 (MQ=255)
                                                                                                                                       ggACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGAccc              >  1:1223485/1‑150 (MQ=40)
                                                                                                                                           cccGGAGCGCAGCGCAGGGGTTccctcactccggccttgctcccgacaggaccgcgccgaaggccccatcctggggcctcatcgcttgacgggcatccatgcccgtcacctgcctccgcaaggcctgcgttcggccacctgaagtcgtga          <  1:110293‑M2/150‑129 (MQ=255)
                                                                                                                                             cGGAGCGTAGCGCAGGGGTcccctcgctccgggcttgctcccgggaggaccgcgctgcaggccccatcctggggcccagcgcttgacgggcatccatgcccgtcacctccctccgcaagccctgcgctcggcctcctgaagtcgcaatct        <  2:679801‑M2/150‑132 (MQ=255)
                                                                                                                                             cGGAGCGCAGCGCAGGGGTCccctcactccggccttgctcccggcaggaccgcgccgaaggccccatcctggggcctcagcgcttgacgggcatccatgcccgtcacctgccttcgcaaggcct                                  >  2:821093‑M2/1‑20 (MQ=255)
                                                                                                                                             cGGAGCGCAGCGCAGGGGTCccctcactccggccttgctcccggcaggaccgcgccgaaggccccatcctggggcctcagcgcttgacgggcatccatgcccgtcacctgccttcgcaaggcct                                  <  1:821093‑M2/124‑105 (MQ=255)
                                                                                                                                               gAGCGTAGCGCAGGGGTcccctcgctccgggcttgctcccgggaggaccgcgctgcaggccccatcctggggcccagcgcttgacgggcatccatgcccgtcacctccctccgcaagccctgcgctcggcctcctgaagtcgcaatctgc      >  1:679801‑M2/1‑17 (MQ=255)
                                                                                                                                                agcgcagcgcAGGGGACGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGAGCGAtgt     <  2:935734/150‑1 (MQ=38)
                                                                                                                                                 gcgcagcgcaGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGAGCGAtgtg    <  2:709057/150‑1 (MQ=39)
                                                                                                                                                   gcagcgcagGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGAGCGATGTggg  <  1:1269279/150‑1 (MQ=38)
                                                                                                                                                    |                                                                                                                                                    
CTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGAGCGATGTGGG  >  NC_002947/4908010‑4908306

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