Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,160,4320TG36.1% 176.0 / 39.3 161L159R (CTG→CGG) PP_3661membrane protein
Reads supporting (aligned to +/- strand):  ref base T (52/49);  new base G (42/15);  total (97/64)
Fisher's exact test for biased strand distribution p-value = 7.10e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 164 total aligned reads displayed.

GTGTCGGTGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGT  >  NC_002947/4160284‑4160579
                                                                                                                                                    |                                                                                                                                                    
gtgtCGGTGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTg                                                                                                                                                     >  2:5100167/1‑150 (MQ=255)
 tgtCGGTGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAACGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCAGTCGAATCGGCGCGgc                                                                                                                                                    >  2:3269843/1‑150 (MQ=255)
 tgtCGGTGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTgc                                                                                                                                                    >  1:2828701/1‑150 (MQ=255)
  gtCGGTGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATcggcgcggcg                                                                                                                                                   >  1:2820952/1‑150 (MQ=255)
      gTGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAgc                                                                                                                                               >  2:2795669/1‑150 (MQ=255)
      gTGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGcgccgc                                                                                                                                               >  1:2673493/1‑150 (MQ=255)
       tGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAgcg                                                                                                                                              >  2:2447634/1‑150 (MQ=255)
        gCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAgcg‑c                                                                                                                                            <  1:501249/150‑1 (MQ=255)
         cTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑cc                                                                                                                                           <  2:5446033/150‑1 (MQ=255)
         cTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑cc                                                                                                                                           <  1:4938631/150‑1 (MQ=255)
         cTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑cc                                                                                                                                           <  1:3121985/150‑1 (MQ=255)
          tGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCggcgcggcgcccccccc                                                                                                                                           >  1:863437/1‑143 (MQ=255)
          tGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGCGGCGCCGCG‑ccc                                                                                                                                          >  1:1101983/1‑150 (MQ=255)
          tGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑ccc                                                                                                                                          <  1:948343/150‑1 (MQ=255)
           ggTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGCGGCgccccg‑cccc                                                                                                                                         >  1:4117472/1‑150 (MQ=255)
           ggTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑ccga                                                                                                                                         <  2:3488179/150‑3 (MQ=255)
           ggTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑ccga                                                                                                                                         <  2:925009/150‑3 (MQ=255)
            gTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGCGGCGCCGCG‑ccccc                                                                                                                                        >  1:3505539/1‑150 (MQ=255)
            gTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAGCG‑ccccc                                                                                                                                        >  2:3040636/1‑150 (MQ=255)
                tCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGCGGCGCCGCG‑cccccccaa                                                                                                                                    >  1:1603278/1‑147 (MQ=255)
                tCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGCGGCGCCCCG‑cccccccaa                                                                                                                                    >  1:4426569/1‑147 (MQ=255)
                tCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGGGGCGCCGCG‑cccccccaa                                                                                                                                    >  1:3141281/1‑147 (MQ=255)
                  aCCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAGCG‑CCCCCCAGaaa                                                                                                                                  >  2:1348361/1‑150 (MQ=255)
                   ccAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACGGCCTTGCTTACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGAATTCGAATCGGCGCTGGGCAGCG‑CCCCCCATaaaa                                                                                                                                 <  2:3771596/150‑1 (MQ=255)
                    cAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCAAAAAAg                                                                                                                                >  2:1335892/1‑150 (MQ=255)
                     aGTACCCTGCTTATCGGCGTGTGGTGTGGCATGAATGTTTTCAAGCTGGACCGCGATACCGCCTTGCTGACCGCCGCGGGAAGCGCCATCTGCGGCGCCGCCGCGGTGCTGGCATTGTAATGGGCGCTGCGCAGCG‑CCCCCCATAAAAgc                                                                                                                               <  2:86349/150‑1 (MQ=255)
                         cccTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCa                                                                                                                           >  2:4237645/1‑150 (MQ=255)
                          ccTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCCGCG‑CCCCCCAAAAAAGCGCCaa                                                                                                                          >  1:395262/1‑149 (MQ=255)
                          ccTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCaa                                                                                                                          >  2:2084567/1‑149 (MQ=255)
                          ccTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCAAAAAAGCGCCaa                                                                                                                          >  2:2084669/1‑149 (MQ=255)
                           ctgctgATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCCCCG‑CCCCCCCGAAAAGCGcccgg                                                                                                                         >  1:2365365/1‑147 (MQ=255)
                            tgctgATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAGCG‑CCCCCCCAAAAATCGCCAAgg                                                                                                                        >  1:4800851/1‑150 (MQ=255)
                             gctgATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCCGCG‑CCCCCCATAAAAGCGCCaatgc                                                                                                                       >  2:1288974/1‑146 (MQ=255)
                             gctgATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCC‑CCCCCCCAAAAAACGcccaggc                                                                                                                       >  1:1717211/1‑145 (MQ=255)
                             gctgATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCCGCG‑CCCCCCCAAAAAGCCCCAAGGc                                                                                                                       >  1:1718132/1‑150 (MQ=255)
                                 aTCGGCCTGTGGTGTGGCATGAAAGGGTTCAAGCTCGACCGCGATACCGCCTTGCTGCCCGCCGCCGGCTGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTc                                                                                                                   <  2:4805887/150‑1 (MQ=255)
                                  tCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCg                                                                                                                  >  2:2266048/1‑150 (MQ=255)
                                    ggCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCAGCGCAGCG‑CCCCCCAAAAAAGCGCCAAGGCAGTCGGa                                                                                                                >  2:1391473/1‑150 (MQ=255)
                                     gCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGCGGGGCCGCG‑CCCCCCCCAAAAGCCCCCTGGCcggcgggg                                                                                                               >  1:197318/1‑148 (MQ=255)
                                      ccTGTGGTGTGGCATTAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGATTTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAg                                                                                                              <  2:1650999/150‑1 (MQ=255)
                                        tgtggtgtggCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCCTAAAAGCCCCAAGGCAGTCGGAATCg                                                                                                            >  2:3095768/1‑150 (MQ=255)
                                             tgtggCTTTAATGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTg                                                                                                       <  2:2609735/150‑1 (MQ=255)
                                             tgtggCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCCCCG‑CCCCCCCAAAAACCGCCCGGGCAGGAGGAAGCGTAGTg                                                                                                       >  1:4293458/1‑150 (MQ=255)
                                                 gCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAGCG‑CCCCCCCAAAAATCGCCAAGGCAGTCGGAAGCGTAGTGCGGt                                                                                                   >  1:5664917/1‑150 (MQ=255)
                                                 gCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAGCG‑CCCCCCCAAAAAAAGCCCGGGCCGTCGGAAGCGGAGttgggt                                                                                                   >  1:2623857/1‑145 (MQ=255)
                                                  cttGAATGTTTTCAAGCTCGACCGCTATACCGGCTTGCTGACGGCCGGGGGGTGGGCCTTGTGGGGCGCCGCCGCGGTGCTGGAATTCGAATGGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGtt                                                                                                  <  2:60540/148‑1 (MQ=255)
                                                   aTGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCACCG‑CCCCCCCGAAAAGCGCCCACGCAGGCGGAAGCGTAGTGCTGTTc                                                                                                 >  1:2313341/1‑150 (MQ=255)
                                                    ttaatgtgtTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGAATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCg                                                                                                <  2:695886/145‑1 (MQ=255)
                                                    tGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGCAGCGCCGCG‑CCCCCCCCAAAAAAGCCAAGGCAGGCGGAAGCCGAGTGGTGttag                                                                                                >  1:605045/1‑148 (MQ=255)
                                                        aGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTAc                                                                                            <  2:5009192/150‑1 (MQ=255)
                                                            ttCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCAGCGCCCCG‑CCCCCCAGAAAAACCCCAACGGAGGCCGAAACCTACTGCTGTTggggaacatg                                                                                        >  1:3618807/1‑143 (MQ=255)
                                                             tCAAGCTCGACGGCTTTACCGCGTTGCTGCCCGCCGCGGGCTGCGCATTGTGGGGCGCCGCCGCCGTGGTGGCTTTCGAAGCGGGGGTGCGCCGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGt                                                                                       <  1:4618868/150‑1 (MQ=255)
                                                              cAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAGCG‑CCCCCCAGAAAAGCGCCATGGCAGTCGGAAGCGTAGAGCTGTTCGGTACCTTGTc                                                                                      >  2:4981641/1‑150 (MQ=255)
                                                               aaGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCCGCC‑CCCCCCAAAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTCCCGTGTCg                                                                                     >  2:844447/1‑150 (MQ=255)
                                                               aaGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAGCG‑CCCCCCCTAAAAGCCCCCAGGCAGTCGGCAGCGGAGTGGGGTTCGGTGCCgtgtgg                                                                                     >  1:5092794/1‑148 (MQ=255)
                                                                 gCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAAGCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGAt                                                                                   <  2:302569/150‑1 (MQ=255)
                                                                   tCGACCGGGATACCGCCTTGCTGACCGCCGCGGGCTGCGCCATCGGCGGGGGCGCCGCCTGGCTGGCATTCGAATCGGGGCTGCCCAGCG‑CCCCCCATAAAAGCGCCATGTCAGTCGGAAGCGTAGTGCTGTTCGGTAACTTGTCGATGt                                                                                 <  2:5071690/150‑1 (MQ=255)
                                                                     gACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCGTGTGGATGTTc                                                                               >  2:3663775/1‑150 (MQ=255)
                                                                     gACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCGCCC‑CCCCCCCGAAAAGACCCATGGCAGTCGGAAGCGTAGTTGTGTTGGGGACCTTGTCGATGttg                                                                               >  1:3053906/1‑149 (MQ=255)
                                                                      aCCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGATTCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTcc                                                                              <  2:4373650/150‑1 (MQ=255)
                                                                        cgcgATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGACGCGTAGTGCTGTTCGGTACCGTGTCGATGTTCCtt                                                                            >  2:1463475/1‑150 (MQ=255)
                                                                         gcgATACGGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGGTGGCATTCGAATCGGGGCTGCGCCGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCttt                                                                           <  1:2329861/150‑1 (MQ=255)
                                                                            aTACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAGCG‑CCCCCCATAAAAGAGCCAAGGGCGTCGGAAGCGGAGTGCTGTTCGGTACCTTGTGGATGTTACTTTaac                                                                        >  1:1792417/1‑148 (MQ=255)
                                                                                cGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCCGCG‑CCCCCCCAAAAAGCGCCCAGGCCGTCGGAAGCCTAGTGCTGTTCGGTACCATGTAGATGTTGCGTTAACCAca                                                                    >  1:5072342/1‑149 (MQ=255)
                                                                                 gCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCCGAAAAGCGCCCTGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCGTGTCGATGTTCCGTTACCCACTg                                                                   >  2:5453327/1‑150 (MQ=255)
                                                                                 gCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAGCG‑CCCCCCCGAAAAGCGCCCGGGCAGTCGGAAGCGTAGTGCTGTTCGGTTCCGGGTCGCTGTTGCTTTACCCaagg                                                                   >  1:1505316/1‑147 (MQ=255)
                                                                                      gCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGAt                                                              <  2:2341600/150‑1 (MQ=255)
                                                                                       cTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATc                                                             <  2:5219821/150‑1 (MQ=255)
                                                                                       cTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGAATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATc                                                             <  2:1761069/150‑1 (MQ=255)
                                                                                        tGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGCTGCGCCGCGCCCCCCCCAAAAAGCCCCA‑GGCAGGCGGTAGCGGAGTGGTGTTGGGTACCGTGTGGGTGTTGCTTTTACCAAAGGCGGTCa                                                            >  1:2600139/1‑150 (MQ=255)
                                                                                         gACCGCCGCGGGCAGCGCCATCGGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCaa                                                           <  2:2318634/150‑1 (MQ=255)
                                                                                           ccgccgccgGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATg                                                         <  2:3688121/150‑1 (MQ=255)
                                                                                            cgccgccgGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGCCTCGGCGCGGCGCAGCG‑CCCCCCCGAAAAGCGCCCGGGCAGTCGGACGCGGCGGGCGGGTCGGTACCTTGTCGCTGTGCCTGTACCCACTGGCGAGCCATGc                                                        >  2:1951895/1‑150 (MQ=255)
                                                                                            cgccgccgGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGc                                                        >  2:71349/1‑150 (MQ=255)
                                                                                                gccgGCTGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCggc                                                    <  2:4903269/150‑1 (MQ=255)
                                                                                                gccgGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGCGGCGCCGCG‑CCCCCCCAAAAAGCGCCATGGCAGTCGGAAGCGTAGTGGGGTTCGGTACCGTGTCGATGTTCCATTAACCACAGGCGCTCAAAGCCcgc                                                    >  1:380432/1‑150 (MQ=255)
                                                                                                gccgGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCGGTTCGGGACCGTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCggc                                                    >  2:1331974/1‑150 (MQ=255)
                                                                                                 gcgGCCGCGCCCTGGGCGGCGCCGCCGCCGTGGTGGAATGCGACGCGGCGCGGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCggct                                                   <  1:5486478/149‑1 (MQ=255)
                                                                                                         gcctcgggcggcgCCGCCGCGGGGCGGGCATCCGAAGCGGGGGTGCGCCGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCAt                                           <  1:1011469/143‑1 (MQ=255)
                                                                                                         cccATCGGGGGCGCCGCCGCCGTGCTGGCATTCTACGCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCAt                                           <  1:4145342/149‑1 (MQ=255)
                                                                                                           cctcgggggcgCCGCCGCGGGGCTGGCATGCGCATGGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCGTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCt                                         <  1:1815907/143‑1 (MQ=255)
                                                                                                           cATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAGCG‑CCCCCCAGAAAAGCGCCAAGGCAGTCGGAAGCGGAGTGCTGTTGGGTACCAGGTCGTTGTTGCTTTACCCACAGGCGATCAATTCCGGCTGGGTGCATCt                                         >  1:1747441/1‑150 (MQ=255)
                                                                                                            aTCCGCGGGGCCGCCGCCGTGGGGGCATTCGCATCGGGGCGGCGCCGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTg                                        <  2:2126816/150‑1 (MQ=255)
                                                                                                             tCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTgg                                       >  2:1384080/1‑150 (MQ=255)
                                                                                                             tCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCCGCG‑CCCCCCCTACAAGCGCCCAGGCAGTCGGAAGCGTAGTGCTGGTCGGTACCGTGTCGATGGTGCGTTACCCACTGGCGATCAATGCCGGCTGGCTGCagcgag                                       >  2:1491354/1‑145 (MQ=255)
                                                                                                                        gccgccgTGCTGGCATTCGAATCGGCGCTGCGCCGCG‑CCCCCCCTAAAAGCGCCCTGGCAGTCGGACGCGTAGTGCTGTTCGGTACCTTGTCGATGTTGCATTACCCACTGGCGATCAAATCCGGCTGGCTGCATCTGGATAACAAGGgc                            >  1:4556947/1‑150 (MQ=255)
                                                                                                                            ccgTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCGTGTCGATGTTGCTTTAACCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGccc                        >  2:4528247/1‑149 (MQ=255)
                                                                                                                                 ctggatttCGACTCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGgctg                   <  2:197659/144‑1 (MQ=255)
                                                                                                                                    gCATTCGAATCGGCGCGGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTAACCACAGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGgctgctg                >  2:2162397/1‑150 (MQ=255)
                                                                                                                                       ttCGACGCGGCGCGGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGctgctgctg             <  1:3955284/150‑1 (MQ=255)
                                                                                                                                       ttCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAATCGCCAAGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTAACCACAGGCGATCAATGCCGGCTGGCTGCATCTGGATACCAAGGGCGCCGGGctgctgctg             >  1:5706831/1‑150 (MQ=255)
                                                                                                                                         cgcctCGGCGCGGCGCCCCC‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTggg           <  1:4816845/146‑1 (MQ=255)
                                                                                                                                          tggggggCGCTGCGCAGCG‑CCCCCCATAAAATCGCCAAGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTTTCGATGTTACTTTAACCACAGGCGATCAAATCCGGCTGGCTGCAACTGGATACCAAGGGCGCCGGGCTGCTGCTGggc          >  2:5727585/6‑150 (MQ=255)
                                                                                                                                          gactCGGCGCGGCGCAGCG‑CCCCCCATAAAAGCGCCAGGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCGGGTCGATGTTCCTTTACCCACTGGCGATCACTGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGggc          <  2:3994779/147‑1 (MQ=255)
                                                                                                                                          gAAGCGGCGCGGCGCCGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGggc          <  1:3558396/150‑1 (MQ=255)
                                                                                                                                           ggggggCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGAt                                                              <  2:2318538/93‑1 (MQ=255)
                                                                                                                                            aTCGGCGCTGCGCAGCG‑CCCCCCCTAACAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGggcgg        <  2:4882618/150‑1 (MQ=255)
                                                                                                                                             ggggCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGggcggc       >  2:5111412/3‑150 (MQ=255)
                                                                                                                                              gggcgcggcgcAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCTGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCa      <  1:5100167/149‑1 (MQ=255)
                                                                                                                                              cGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCa      <  2:2828701/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GTGTCGGTGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGT  >  NC_002947/4160284‑4160579

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