Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,349,312 A→G 12.8% intergenic (+50/‑106) PP_3817 → / → PP_3818 polyamine ABC transporter ATP‑binding protein/phosphate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,349,3120AG12.8% 617.0 / 62.4 228intergenic (+50/‑106)PP_3817/PP_3818polyamine ABC transporter ATP‑binding protein/phosphate‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (100/97);  new base G (17/12);  total (118/110)
Fisher's exact test for biased strand distribution p-value = 5.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.72e-01

Only 100 of 239 total aligned reads displayed.

CCGGGCGACAGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTTGCTGCTCACCCTGCTACCACTCGC  >  NC_002947/4349164‑4349461
                                                                                                                                                    |                                                                                                                                                     
ccGGGCGACAGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGcac                                                                                                                                                      <  1:4673378/150‑1 (MQ=255)
    gCGACAGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGa                                                                                                                                                  <  2:4698067/150‑1 (MQ=255)
      gACAGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATc                                                                                                                                                <  1:2429380/150‑1 (MQ=255)
       aCAGCGTGTAGGTGCGTTGGGCAGTGGCTGACGCTTGCATTTATACCGAATGGGGCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTAGGTTGCTGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCt                                                                                                                                               <  2:5463568/150‑1 (MQ=255)
       aCAGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGACTGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGTGTGCCGCGCGGGCGGTGCTCGATCt                                                                                                                                               <  2:6051958/150‑1 (MQ=37)
       aCAGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGTGTGCCGCGCGGGCGGTGCTCGATCt                                                                                                                                               <  2:1013416/150‑1 (MQ=38)
       aCAGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCt                                                                                                                                               >  2:3117075/1‑150 (MQ=255)
        cAGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGTGTGCCGCGCGGGCGGTGCTCGATCTg                                                                                                                                              <  2:1060358/150‑1 (MQ=38)
           cGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCAGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCa                                                                                                                                           >  1:1120906/1‑150 (MQ=255)
           cGTGAAGGTGAGTTGTGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGTGTGCCGCGCGGGCGGTGCTCGATCTGGCa                                                                                                                                           <  2:2760671/150‑1 (MQ=37)
            gTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAg                                                                                                                                          <  1:3117075/150‑1 (MQ=255)
                  gTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTg                                                                                                                                    >  1:5432812/1‑150 (MQ=255)
                   tGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGc                                                                                                                                   >  1:492455/1‑150 (MQ=255)
                      gTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGAATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCaac                                                                                                                                <  2:2675169/150‑3 (MQ=255)
                        tGGGCTGTGGCTGACGCTTTCTTTTATACCGAAGGGGCCGAGCGCGACCTGAGCAAGTCGGCGGGGGCGAATTAAGGCTTTGGTACGTGGCTGCGCCTGCCTTTGATTGCCCGTCCCGCGCGGAACACGATCGTGCAGGCGCTGCACTAt                                                                                                                              <  2:4578720/150‑1 (MQ=255)
                        tGGGCAGTGGCAGACGCTTGCATTTATACCGATTGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCATATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTAt                                                                                                                              <  2:3450817/150‑1 (MQ=255)
                        tGGGCAGTGGCAGACGCTTGCATTTATACCGACTGGGCCGAGAGCGACCTGAGCATGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTAt                                                                                                                              <  2:3688805/150‑1 (MQ=255)
                         gtgCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCGGCACTAtt                                                                                                                             >  2:6992932/3‑150 (MQ=255)
                         gggCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGATGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTAtt                                                                                                                             <  2:6824478/150‑1 (MQ=255)
                                 gCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGt                                                                                                                     >  2:5566984/1‑150 (MQ=255)
                                 gCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGt                                                                                                                     >  2:3281177/1‑150 (MQ=255)
                                 gCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGt                                                                                                                     >  2:3281333/1‑150 (MQ=255)
                                 gCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCAGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCGGCCCTATTGGAGGCgg                                                                                                                     >  1:748411/1‑149 (MQ=255)
                                  cAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCCGGCGCGGCACTAATGGAGGCGTa                                                                                                                    >  1:6903806/1‑150 (MQ=255)
                                     aCGCTTGCATTTATACCGACGGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGc                                                                                                                 <  2:5673446/150‑1 (MQ=255)
                                         ttGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCg                                                                                                             >  1:6044592/1‑150 (MQ=255)
                                            cATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGAATTAACGCTCTTGTACGTTGCAGCGCCTGCAAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTca                                                                                                          <  2:1351547/150‑1 (MQ=255)
                                             aTTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTcac                                                                                                         >  2:7336313/1‑150 (MQ=255)
                                              tttATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGgcg                                                                                                                                       >  2:4238340/1‑119 (MQ=255)
                                                   aCCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCCCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTatc                                                                                                   >  1:2791903/1‑150 (MQ=255)
                                                    ccGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTatca                                                                                                  >  1:6127914/1‑150 (MQ=255)
                                                      gAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTGGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTatcatc                                                                                                <  2:3361035/150‑1 (MQ=255)
                                                          gggCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGc                                                                                            >  2:4792191/1‑150 (MQ=255)
                                                            gCCGAGAGCGCCCTGAGCATGTCGGCGGGGGCGCATTAACGCTCTTGTCGGTTGCAGCGCCTGCCTGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCa                                                                                          <  1:5384492/150‑1 (MQ=255)
                                                            gCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGGCAGCCAGATAGAGCACCGCCCGCGCCGCACACGATCAGGCAGGCGCTGCCCGATTGGGGGCGTACGCCCCCGCACAATATCATCACGccc                                                                                          >  1:4982712/1‑149 (MQ=255)
                                                                 gagCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTc                                                                                     <  2:7186254/150‑1 (MQ=255)
                                                                      aCCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAAt                                                                                >  2:300410/1‑150 (MQ=255)
                                                                          gagaAAGTCGGCCGGGGCGAATTAACGCTCTTGTACGTTGCAGCGCCTGCCATATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGc                                                                            <  1:6584466/150‑1 (MQ=255)
                                                                          gAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGca                                                                                                                                                       <  1:4750169/75‑1 (MQ=255)
                                                                          gAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGca                                                                                                                                                       >  2:4750169/1‑75 (MQ=255)
                                                                          gAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGc                                                                            >  2:5279531/1‑150 (MQ=255)
                                                                                   ggCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGTCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAAct                                                                   >  2:7044631/1‑150 (MQ=255)
                                                                                      cGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTa                                                                >  2:6788758/1‑150 (MQ=255)
                                                                                       ggggggCTTTACGGCTGTTGTGGGTGGCTGCGCGTGCCAGATGGCGCACGGCCGGGGGGGCACTCTTTGTGGCTGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTAt                                                               <  1:4517488/150‑1 (MQ=255)
                                                                                       ggggCGATTTAACGATTTTGTACTTTGAAGCGCCTGCAATATCGAGCACCGCCCGGGCGGAAAACGTTGTGGCTGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTAt                                                               <  1:1408631/150‑1 (MQ=255)
                                                                                          gcgcATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCAACACGCCACTGTCGTAATCTGCCGGCCAACTCTATGcc                                                            >  2:5237325/1‑150 (MQ=255)
                                                                                           cgcATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACTATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCca                                                           <  2:701452/150‑1 (MQ=255)
                                                                                            gcATTAACGCTCTTGTACGTTACAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCcac                                                          >  2:2509995/1‑150 (MQ=255)
                                                                                             cATTACCGCTCTTGTAGGTTGCCGCGCCTGCCTGATCGAGCCCCGCCCGCGCGGAACACTATGTGGCTGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCcaca                                                         <  1:3875938/150‑1 (MQ=255)
                                                                                              aTTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAg                                                        <  2:6339328/150‑1 (MQ=255)
                                                                                               ttAACGCTGTTTTACGTTGCAGCGCCTGCCTGATCGAGCACCGCCCGGGCGGCACACTATCTGGCTGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAgg                                                       <  1:5982844/150‑1 (MQ=255)
                                                                                                  aCGCTCTTGTACGTTGCAGCGCCTGCCCGATCCAGCACCCCCCCCGCGGCACACCATCTGGCAGGCGCTGCACTAATGGAGGCGTACGCCCCGGCACAATATCATCACGCCACAGTCGTAATATGCCGGCCCACTCTATGCCACAGGTCa                                                    >  1:5841547/1‑150 (MQ=255)
                                                                                                     ctctTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCAccc                                                 <  2:1204642/150‑1 (MQ=255)
                                                                                                                tGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCg                                      >  1:4280091/1‑150 (MQ=255)
                                                                                                                 gCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGc                                     <  1:5188992/150‑1 (MQ=255)
                                                                                                                   aGCGCCTGCCAGATCGAGCACCGCCCGCGGGGCACACGATCTGGCTGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCt                                   <  1:865223/150‑1 (MQ=255)
                                                                                                                     cgcCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGc                                 <  2:6300202/150‑1 (MQ=255)
                                                                                                                      gcCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCt                                >  2:7306375/1‑150 (MQ=255)
                                                                                                                         tGCCAGATCGTGTGCCGCGCGGGCGGTGCTCGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCa                                                                                          >  1:5222869/1‑89 (MQ=17)
                                                                                                                         tGCCAGATCGTGTGCCGCGCGGGCGGTGCTCGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCa                                                                                          <  2:5222869/89‑1 (MQ=17)
                                                                                                                          gCCAGATCGTGTGCCGCGCGGGCGGTGCTCGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCtt                               >  2:4968521/1‑147 (MQ=37)
                                                                                                                          gCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTc                            >  1:4869532/1‑150 (MQ=255)
                                                                                                                           ccAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCg                                                                          >  2:2084960/1‑103 (MQ=255)
                                                                                                                           ccAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCt                           >  2:4471364/1‑150 (MQ=255)
                                                                                                                            cAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTCTGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCtt                          >  2:5749804/1‑150 (MQ=255)
                                                                                                                            cAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCtt                          <  1:5506697/150‑1 (MQ=255)
                                                                                                                            cAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCtt                          >  1:5580838/1‑150 (MQ=255)
                                                                                                                            cAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCtt                          >  1:3247448/1‑150 (MQ=255)
                                                                                                                             cgATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTtg                         <  1:5979660/149‑1 (MQ=255)
                                                                                                                             aGATCGTGTGCCGCGCGGGCGGTGCTCGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTa                                                                <  2:1859534/111‑1 (MQ=25)
                                                                                                                             aGATCGTGTGCCGCGCGGGCGGTGCTCGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTtg                         >  1:7289832/1‑150 (MQ=38)
                                                                                                                             aGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTCTCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTtg                         >  2:5991645/1‑150 (MQ=255)
                                                                                                                             aGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCg                                                >  1:204415/1‑127 (MQ=255)
                                                                                                                             aGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCtt                               <  1:2854117/144‑1 (MQ=255)
                                                                                                                             aGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTtg                         >  1:90689/1‑150 (MQ=255)
                                                                                                                              gATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCTTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCtt                          <  1:4226946/148‑1 (MQ=255)
                                                                                                                              gATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCtt                          >  2:4226946/1‑148 (MQ=255)
                                                                                                                              gATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTtgc                        >  2:1602709/1‑150 (MQ=255)
                                                                                                                               aTCGCGCCCCGCCCGCGCGGCACCCGCTCGGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTCTGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTtgct                       <  2:3245420/150‑1 (MQ=255)
                                                                                                                               aTCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTtgct                       <  1:6815184/150‑1 (MQ=255)
                                                                                                                               aTCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTtgct                       >  1:5788155/1‑150 (MQ=255)
                                                                                                                                tCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACGAATGGAGGCGTACGCCCCCTCACATTATCATCACCCCACTGTCGTAATCAGCCGGCCAACTCGATGCCACAGGTCACCCGGCATGCTCCGCCGGCGGGTCTtgctg                      >  2:6868095/1‑150 (MQ=255)
                                                                                                                                 cgtgtgccgcgcggGCGGTGCTCGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTtgctgc                     >  2:2522053/7‑150 (MQ=38)
                                                                                                                                  gtgtgccgcgcggGCGGTGCTCGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTTgctgct                    >  1:6792473/6‑150 (MQ=38)
                                                                                                                                   tgtgccgcgcggGCGGTGCTCGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTTGCTGCTc                   >  1:3524905/5‑150 (MQ=39)
                                                                                                                                   aGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTTGCTGCTc                   >  2:5884523/1‑150 (MQ=255)
                                                                                                                                    gCACCGCCCTAGCTGCACACGATCTGGAAGTCTCTGCAATATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCAGGCCAACTCTATGCCACAGGTCACACGTCATGCTCCGCCTGCTTGTCTTGCTGCTCa                  >  2:2609631/1‑150 (MQ=255)
                                                                                                                                      aCCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTa                                                                                  >  1:2244701/1‑84 (MQ=255)
                                                                                                                                      aCCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTa                                                                                  <  2:2244701/84‑1 (MQ=255)
                                                                                                                                      aCCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAAct                                                                   <  2:5946273/99‑1 (MQ=255)
                                                                                                                                      aCCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTTGCTGCTCAcc                >  1:5274149/1‑150 (MQ=255)
                                                                                                                                       ccgcccgcGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTTGCTGCTCAccc               >  1:5945003/1‑150 (MQ=255)
                                                                                                                                        cgcccgcGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTTGCTGCTCACCCt              >  2:6199331/1‑150 (MQ=255)
                                                                                                                                          cccgcGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTTGCTGCTCACCCTGc            <  2:4337954/150‑1 (MQ=255)
                                                                                                                                            cgcgcgGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTTGCTGCTCACCCTGCTa          <  1:1519596/150‑1 (MQ=255)
                                                                                                                                             gcgcgGCACACTATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACGGGTCACCCGGCATGCTCCGCCTGCTTGTCTTGCTGCTCACCCTGCTAc         <  1:7090481/150‑1 (MQ=255)
                                                                                                                                             gcgcgGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTTGCTGCTCACCCTGCTAc         >  2:929067/1‑150 (MQ=255)
                                                                                                                                                cgGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTTGCTGCTCACCCTGCTACCAc      <  2:3492563/150‑1 (MQ=255)
                                                                                                                                                   cacaCGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCaa                                                                     >  1:3982571/1‑84 (MQ=255)
                                                                                                                                                   cacaCGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCaa                                                                     <  2:3982571/84‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CCGGGCGACAGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTTGCTGCTCACCCTGCTACCACTCGC  >  NC_002947/4349164‑4349461

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