Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,908,169 C→T 13.5% intergenic (‑100/‑182) 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,1690CT13.5% 180.6 / 20.7 88intergenic (‑100/‑182)PP_4317/PP_4318hypothetical protein/transposase
Reads supporting (aligned to +/- strand):  ref base C (36/39);  new base T (4/8);  total (40/48)
Fisher's exact test for biased strand distribution p-value = 5.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.20e-01

GCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGAGCGATGTGGGCTGATGTAC  >  NC_002947/4908020‑4908315
                                                                                                                                                     |                                                                                                                                                  
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)
         ggaggCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATACCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGatgatg                                                                                                                                           >  1:935734/1‑150 (MQ=17)
            ggCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGa                                                                                                                                        <  2:9575/150‑1 (MQ=32)
             gCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAg                                                                                                                                       <  1:475176/150‑1 (MQ=32)
             gCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAg                                                                                                                                       <  2:1095629/150‑1 (MQ=32)
               cGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAgcg                                                                                                                                     >  1:973949/1‑150 (MQ=32)
               cGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAgcg                                                                                                                                     >  2:536947/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                                                                                                                                        >  1:397713‑M2/6‑136 (MQ=255)
                          cttgcGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGa                                                                                                                                        <  2:397713‑M2/131‑1 (MQ=255)
                          cttgcGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGtttttt                                                                                                                          <  2:1152546‑M2/145‑1 (MQ=255)
                          cttgcGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGtttttt                                                                                                                          <  1:536947‑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:712540/1‑150 (MQ=32)
                                                                                     ccTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAg                                                               <  1:624217/150‑1 (MQ=32)
                                                                                      cTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCCa                                                                          >  2:996172/1‑138 (MQ=33)
                                                                                      cTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCCa                                                                          <  1:996172/138‑1 (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                                    <  1:442402/139‑1 (MQ=40)
                                                                                                                           ggggACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCa                                    >  2:442402/1‑139 (MQ=40)
                                                                                                                           ggggACCCCGGAGCGCAGCGCAGGGGCCGGATGACGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGAc                         >  2:564457/1‑150 (MQ=255)
                                                                                                                            gggACCCCGGAGCGTAGCGTAGGGGTcccctcactccgggcttgctcccgggaggaccgcgctgcaggccccatcctggggcccagcgcttgacgggcatccatgcc                                                                   >  2:181519‑M2/1‑26 (MQ=255)
                                                                                                                            gggACCCCGGAGCGTAGCGTAGGGGTcccctcactccgggcttgctcccgggaggaccgcgctgcaggccccatcctggggcccagcgcttgacgggcatccatgcc                                                                   <  1:181519‑M2/107‑82 (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)
                                                                                                                                             cgcagGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGAGCGATGTGGGCTGa       <  2:1245927/150‑1 (MQ=35)
                                                                                                                                                agGGGCcccctcactccgggcttgctcccgggaggaccgcgctgcaggccccatcctggggcccagcgcttgacgggcatccatgcccgtcacctccctccgcaagccctgcgctcggcctcctgaagtcgcaatctgcatcgcctgaac    >  2:228133‑M2/1‑6 (MQ=255)
                                                                                                                                                agGGGCCccctcactccggcctcgctcccggcaggaccgcgccgaaggccccatcctggggcctcagcgcttgacgggcatccatgcccgtcacctgcctccgcaaggcctgcgttcggcctcctgaagtcgcgaagatcaagatcaaga    >  1:107019‑M2/1‑7 (MQ=255)
                                                                                                                                                  gggCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGAGCGATGTGGGCTGATGTAc  <  2:712540/150‑1 (MQ=35)
                                                                                                                                                     |                                                                                                                                                  
GCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCGAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGAGCGATGTGGGCTGATGTAC  >  NC_002947/4908020‑4908315

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: