Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,887,817 C→T 5.1% intergenic (+67/‑158) oprN → / → PP_3428 multidrug RND transporter outer membrane protein OprN/TPR repeat‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,887,8170CT 5.1% 951.8 / 21.9 294intergenic (+67/‑158)oprN/PP_3428multidrug RND transporter outer membrane protein OprN/TPR repeat‑containing protein
Reads supporting (aligned to +/- strand):  ref base C (145/134);  new base T (7/8);  total (152/142)
Fisher's exact test for biased strand distribution p-value = 7.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 309 total aligned reads displayed.

TCGACCTTTATCGCGGCATCGTCTCGATCTACAAGGCCCTGGGTGGCGGCTGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCCGGTTCGGTTGCATTTACCCCCCGGCTGGCCCAAGCTCTGCCCTTCCCGCCGCCACGG  >  NC_002947/3887668‑3887966
                                                                                                                                                     |                                                                                                                                                     
    ccTTTATCGCGGCATCGTCTCGATCTACAAGGCCCTGGGTGGCGGCTGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCGGGCCATGGCCTGCgg                                                                                                                                                   <  1:99679/150‑1 (MQ=255)
    ccTTTATCGCGGCATCGTCTCGATCTACAAGGCCCTGGGTGGCGGCTGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCgg                                                                                                                                                   <  1:2619231/150‑1 (MQ=255)
     cTTTCTCGCGGCATCGTCTCGATCTACAAGGCCCTGGGTGGCGGCTGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACGGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCGGGCCATGGCCTGCggg                                                                                                                                                  <  1:637615/150‑1 (MQ=255)
     cTTTATCGCGGCATCGTCTCGATCTACAAGGCCCTGGGTGGCGGCTGGCAACCGGTGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTTGCTCCTCCCCCAACCGTTTGCCCCGCGGGCCATGGCCTGCggg                                                                                                                                                  <  2:3861816/150‑1 (MQ=255)
     cTTTATCGCGGCATCGTCTCGATCTACAAGGCCCTGGGTGGCGGCTGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCGGGCCATGGCCTGCggg                                                                                                                                                  <  1:4289744/150‑1 (MQ=255)
     cTTTATCGCGGCATCGTCTCGATCTACAAGGCCCTGGGTGGCGGCTGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCGGGCCATGGCCTGCggg                                                                                                                                                  <  1:5245493/150‑1 (MQ=255)
      tttATCGCGGCATCGTCTCGATCTACAAGGCCCTGGGTGGCGGCTGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCgggg                                                                                                                                                 <  2:3629080/150‑1 (MQ=255)
          tCGCGGCATCGTCTCGATCTACAAGGCCCTGGGTGGCGGCTGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCGGGCCATGGCCTGCGGGGcaaa                                                                                                                                             <  2:6668333/150‑4 (MQ=255)
           cgcgGCATCGTCTCGATCTACAAGGCCCTGGGTGGCGGCTGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCtttt                                                                                                                                            <  2:2423499/150‑1 (MQ=255)
            gcgGCATCGTCTCGATCTACAAGGCCCTGGGTGGCGGCTGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCttttt                                                                                                                                           <  1:2802989/150‑1 (MQ=255)
            gcgGCATCGTCTCGATCTACAAGGCCCTGGGTGGCGGCTGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCttttt                                                                                                                                           <  1:3266304/150‑1 (MQ=255)
             cgGCATCGTCTCGATCTACAAGGCCCTGGGTGGCGGCTGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCtttttt                                                                                                                                          <  2:4454356/150‑1 (MQ=255)
             cgGCATCGTCTCGATCTACAAGGCCCTGGGTGGCGGCTGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCtttttt                                                                                                                                          <  2:2298248/150‑1 (MQ=255)
               gcctCGTCTCGATCTACAAGGCCCTGGGTGGCGGCTGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTAt                                                                                                                                        <  2:3476804/147‑1 (MQ=255)
               gCATCGTCTCGATCTACAAGGCCCTGGGTGGCGGCTGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTAt                                                                                                                                        <  2:6134117/150‑1 (MQ=255)
                  tCGTCTCGATCTACAAGGCCCTGGGTGGCGGCTGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGtt                                                                                                                                     >  1:2366899/1‑150 (MQ=255)
                         gATCTACAAGGCCCTGGGTGGCGGCTGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATc                                                                                                                              <  1:6233227/150‑1 (MQ=255)
                         gATCTACAAGGCCCTGGGTGGCGGCTGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATc                                                                                                                              <  2:2915783/150‑1 (MQ=255)
                         gATCTACAAGGCCCTGGGTGGCGGCTGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATc                                                                                                                              >  1:5786948/1‑150 (MQ=255)
                         gATCTACAAGGCCCTGGGTGGCGGCTGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATc                                                                                                                              >  1:504277/1‑150 (MQ=255)
                                aaGGCCCTGGGTGGCGGCTGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCa                                                                                                                       >  2:6045637/1‑150 (MQ=255)
                                  ggCCCTGGGTGGCGGCTGGCAACCGGAGACGGTGGCCAGCGCGCACTGAGTCCACTGCCTCAACGTCTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATTGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCAcc                                                                                                                     <  1:4694714/150‑1 (MQ=255)
                                   gCCCTGGGTGGCGGCTGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCAccc                                                                                                                    <  1:5167193/150‑1 (MQ=255)
                                   gCCCTGGGTGGCGGCTGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCAccc                                                                                                                    >  1:1175097/1‑150 (MQ=255)
                                    cccTGGGTGGCGGCTGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCAcccc                                                                                                                   >  2:508378/1‑150 (MQ=255)
                                     ccTGGGTGGCGGCTGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCa                                                                                                                  >  2:4626898/1‑150 (MQ=255)
                                        gggTGGCGGCTGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATa                                                                                                               <  1:3198164/150‑1 (MQ=255)
                                            ggcggcTGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCa                                                                                                           <  1:1797316/150‑1 (MQ=255)
                                            ggcggcTGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCa                                                                                                           <  1:6010545/150‑1 (MQ=255)
                                             gcggcTGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAg                                                                                                          >  1:1664421/1‑150 (MQ=255)
                                             gcggcTGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAg                                                                                                          <  1:1104986/150‑1 (MQ=255)
                                               ggcTGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTg                                                                                                        >  1:5097243/1‑150 (MQ=255)
                                                  tGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCAt                                                                                                     <  1:5410752/150‑1 (MQ=255)
                                                  tGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCAt                                                                                                     >  1:1478439/1‑150 (MQ=255)
                                                  tGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCATATTCACCCCAATACGCAGTGCAt                                                                                                     >  2:1875764/1‑150 (MQ=255)
                                                    gCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATgg                                                                                                   >  2:6758705/1‑150 (MQ=255)
                                                     cAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGa                                                                                                  <  1:3625132/150‑1 (MQ=255)
                                                         cGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCCTGGAATGc                                                                                              >  2:281961/1‑150 (MQ=255)
                                                          ggAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGcc                                                                                             >  1:6627382/1‑150 (MQ=255)
                                                            agaCGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCAc                                                                                           <  1:2536272/150‑1 (MQ=255)
                                                                ggTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCgaa                                                                                       <  2:921784/150‑3 (MQ=255)
                                                                ggTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCa                                                                                       <  2:2604760/150‑1 (MQ=255)
                                                                      cAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGgcctg                                                                                 >  2:2374770/1‑150 (MQ=255)
                                                                        gcgcgcACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGgcctggc                                                                               <  2:3123554/150‑1 (MQ=255)
                                                                              aCTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGgcctggcgcctgn                                                                         <  1:4752496/150‑2 (MQ=255)
                                                                              aCTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTGTGTTATGTTTGCCATCGTATTCACCCCACTACGCAGTGCATGGAATGCCACCGCAGgcctggcgcctgg                                                                         >  2:5141902/1‑150 (MQ=255)
                                                                                tGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGcc                                                                       <  2:224245/150‑1 (MQ=255)
                                                                                    tCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACg                                                                   >  2:5262582/1‑150 (MQ=255)
                                                                                    tCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACg                                                                   >  2:207302/1‑150 (MQ=255)
                                                                                        cTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGt                                                               >  1:2702597/1‑150 (MQ=255)
                                                                                             tCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCa                                                                                                           >  1:5897554/1‑101 (MQ=255)
                                                                                             tCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCcg                                                          <  1:3821846/150‑1 (MQ=255)
                                                                                                aCGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCcggtt                                                       >  1:2972247/1‑150 (MQ=255)
                                                                                                aCGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCcggtt                                                       <  2:1710068/150‑1 (MQ=255)
                                                                                                 cGACTCCTTTGGTTGGCTCCTCCCCCCCCCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCcggttc                                                      <  2:193893/150‑1 (MQ=255)
                                                                                                  gACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCcggttcg                                                     >  2:2719761/1‑150 (MQ=255)
                                                                                                     tCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCtttttt                                                                                                                                          >  1:2090590/1‑62 (MQ=255)
                                                                                                     tCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCcggttcggtt                                                  >  2:6276122/1‑150 (MQ=255)
                                                                                                     tCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCcggttcggtt                                                  <  1:3314829/150‑1 (MQ=255)
                                                                                                     tCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCcggttcggtt                                                  >  1:1429408/1‑150 (MQ=255)
                                                                                                       cTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCCGGTTCGGTTGc                                                >  2:2959376/1‑150 (MQ=255)
                                                                                                        tttGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCCGGTTCGGTTGCa                                               <  2:734880/150‑1 (MQ=255)
                                                                                                        tttGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCCGGTTCGGTTGCa                                               >  2:3803573/1‑150 (MQ=255)
                                                                                                        tttGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCCGGTTCGGTTGCa                                               <  1:2892325/150‑1 (MQ=255)
                                                                                                           ggttggCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCCGGTTCGGTTGCAttt                                            >  1:4500114/1‑150 (MQ=255)
                                                                                                               ggCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCCGGTTCGGTTGCATTTAccc                                        <  1:798257/150‑1 (MQ=255)
                                                                                                                 ctcctcCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCCGGTTCGGTTGCATTTAccccc                                      >  2:2190236/1‑150 (MQ=255)
                                                                                                                   cctccCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCCGGTTCGGTTGCATTTACCCCCCg                                    >  2:1363001/1‑150 (MQ=255)
                                                                                                                    ctccCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCCGGTTCGGTTGCATTTACCCCCCgg                                   >  1:4544258/1‑150 (MQ=255)
                                                                                                                     tccCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCCACGGCGTGCCCGGTTCGGTTGCATTTACCCCCCTGc                                  >  1:1002011/1‑150 (MQ=255)
                                                                                                                     tccCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCCGGTTCGGTTGCATTTACCCCCCGGc                                  >  1:3911397/1‑150 (MQ=255)
                                                                                                                     tccCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCCGGTTCGGTTGCATTTACCCCCCGGc                                  >  2:2499192/1‑150 (MQ=255)
                                                                                                                     ccccccAACCGTGTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCCGGTTCGGTTGCATTTACCCCCCGGc                                  <  1:1203997/149‑1 (MQ=255)
                                                                                                                        cccAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCCGGTTCGGTTGCATTTACCCCCCGGCTgg                               >  1:3584568/1‑150 (MQ=255)
                                                                                                                        cccAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCCGGTTCGGTTGCATTTACCCCCCGGCTgg                               >  2:5487992/1‑150 (MQ=255)
                                                                                                                         ccAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCCGGTTCGGTTGCATTTACCCCCCGGCTGGc                              >  1:1946204/1‑150 (MQ=255)
                                                                                                                           aaCCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCCGGTTCGGTTGCATTTACCCCCCGGCTGGccc                            <  1:1043456/150‑1 (MQ=255)
                                                                                                                             ccGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCCGGTTCGGTTGCATTTACCCCCCGGCTGGCCCaa                          <  1:3081507/150‑1 (MQ=255)
                                                                                                                               gTTTGCCCCGCAGGCCATGGCCCGCGGGGCTATTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCCGGTTCGGTTGCATTTACCCCCCGGCTGGCCCAAGc                        <  1:2252456/150‑1 (MQ=255)
                                                                                                                               aaaagCCCCGCGGGCCATGGCCTGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCCGGTTCGGTTGCATTTACCCCCCGGCTGGCCCAAGc                        >  1:500883/5‑150 (MQ=255)
                                                                                                                                  tGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCCGGTTCGGTTGCATTTACCCCCCGGCTGGCCCAAGctct                     <  2:5857562/150‑1 (MQ=255)
                                                                                                                                   gCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCCGGTTCGGTTGCATTTACCCCCCGGCTGGCCCAAGCTCTg                    <  2:5841296/150‑1 (MQ=255)
                                                                                                                                   gCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCCGGTTCGGTTGCATTTACCCCCCGGCTGGCCCAAGCTCTg                    <  2:4534293/150‑1 (MQ=255)
                                                                                                                                     cccGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCCGGTTCGGTTGCATTTACCCCCCGGCTGGCCCAAGCTCTGcc                  >  2:4943576/1‑150 (MQ=255)
                                                                                                                                     cccGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCCGGTTCGGTTGCATTTACCCCCCGGCTGGCCCAAGCTCTGcc                  <  1:6810062/150‑1 (MQ=255)
                                                                                                                                     cccGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCCGGTTCGGTTGCATTTACCCCCCGGCTGGCCCAAGCTCTGcc                  <  1:3612444/150‑1 (MQ=255)
                                                                                                                                     cccGCAGGCCAGGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCGGGCCAACGGCGTGCCCGGTTCGGTTGCATTTACCCCCCGGCTGGCCCAAGCTCTGcc                  <  1:629656/150‑1 (MQ=255)
                                                                                                                                      ccGCGGGCCATGGCCTGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCCGGTTCGGTTGCATTTACCCCCCGGCTGGCCCAAGCTCTGccc                 >  1:553529/1‑150 (MQ=255)
                                                                                                                                      ccGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCCGGTTCGGTTGCATTTACCCCCCGGCTGGCCCAAGCTCTGccc                 >  1:1935972/1‑150 (MQ=255)
                                                                                                                                       cgcgGGCCATGGCCTGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCCGGTTCGGTTGCATTTACCCCCCGGCTGGCCCAAGCTCTGCCCt                >  1:4582075/1‑150 (MQ=255)
                                                                                                                                       cGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCCGGTTCGGTTGCATTTACCCCCCGGCTGGCCCAAGCTCTGCCCt                >  2:5955924/1‑150 (MQ=255)
                                                                                                                                       cGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCCGGTTCGGTTGCATTTACACCCCGGCTGGCCCAAGCTCTGCCCt                >  1:244432/1‑150 (MQ=255)
                                                                                                                                        gCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCCGGTTCGGTTGCATTTACCCCCCGGCTGGCCCAAGCTCTGCCCtt               >  2:1105162/1‑150 (MQ=255)
                                                                                                                                          aGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCCGGTTCGGTTGCATTTACCCCCCGGCTGGCCCAAGCTCTGCCCTTcc             >  1:1244293/1‑150 (MQ=255)
                                                                                                                                              cATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCCGGTTCGGTTGCATTTACCCCCCGGCTGGCCCAAGCTCTGCCCTTCccgcc         >  1:6687826/1‑150 (MQ=255)
                                                                                                                                              cATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCCGGTTCGGTTGCATTTACCCCCCGGCTGGCCCAAGCTCTGCCCTTCccgcc         >  2:2987546/1‑150 (MQ=255)
                                                                                                                                              cATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATAAGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCCGGTTCGGTTGCATTTACCCCCCGGCTGGCCCAAGCTCTGCCCTTCccgcc         >  1:3040160/1‑150 (MQ=255)
                                                                                                                                                    ccGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCCGGTTCGGTTGCATTTACCCCCCGGCTGGCCCAAGCTCTGCCCTTCCCGCCGCCACg   >  2:4074385/1‑150 (MQ=255)
                                                                                                                                                    ccGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCCGGTTCGGTTGCATTTACCCCCCGGCTGGCCCAAGCTCTGCCCTTCCCGCCGCCACg   >  2:1631830/1‑150 (MQ=255)
                                                                                                                                                     cgcgGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCCGGTTCGGTTGCATTTACCCCCCGGCTGGCCCAAGCTCTGCCCTTCCCGCCGCCACgg  >  2:3431562/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TCGACCTTTATCGCGGCATCGTCTCGATCTACAAGGCCCTGGGTGGCGGCTGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCCGGTTCGGTTGCATTTACCCCCCGGCTGGCCCAAGCTCTGCCCTTCCCGCCGCCACGG  >  NC_002947/3887668‑3887966

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