Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,349,295 A→T 12.9% intergenic (+33/‑123) 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,2950AT12.9% 557.5 / 69.0 224intergenic (+33/‑123)PP_3817/PP_3818polyamine ABC transporter ATP‑binding protein/phosphate‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (95/98);  new base T (13/16);  total (108/116)
Fisher's exact test for biased strand distribution p-value = 6.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.10e-01

Only 100 of 228 total aligned reads displayed.

GCGGGCATGCCGATGCAGCCGGGCGACAGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTTGCTGCTC  >  NC_002947/4349146‑4349444
                                                                                                                                                     |                                                                                                                                                     
gcggGCATGCCGATGCAGCCGGGCGACAGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGa                                                                                                                                                       >  1:3492563/1‑150 (MQ=255)
gcggGCATGCCGATGCAGCCGGGCGACAGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGa                                                                                                                                                       <  2:288968/150‑1 (MQ=255)
  gggCATGCCGATGCAGCCGGGCGACAGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGc                                                                                                                                                     >  1:1817689/1‑150 (MQ=255)
   ggCATGCCGATGCAGCCGGGCGACAGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCa                                                                                                                                                    >  2:5390177/1‑150 (MQ=255)
   ggCATGCCGATGCAGCCGGGCGACAGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCa                                                                                                                                                    <  1:5007132/150‑1 (MQ=255)
    gCATGCCGATGCAGCCGGGCGACAGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCAc                                                                                                                                                   >  2:7090481/1‑150 (MQ=255)
       tGCCGATGCAGCCGGGCGACAGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCAccgc                                                                                                                                                <  2:2955471/150‑1 (MQ=255)
            aTGCAGCCGGGCGACAGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCcgcg                                                                                                                                           <  1:6442430/150‑1 (MQ=255)
               cAGCCGGGCGACAGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCgg                                                                                                                                        >  1:5042994/1‑150 (MQ=255)
                aGCCGGGCGACAGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAg                                                                                                                                                      >  1:2637837/1‑135 (MQ=255)
                aGCCGGGCGACAGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAg                                                                                                                                                      <  2:2637837/135‑1 (MQ=255)
                     ggCGACAGCGTGAAGGTGCGTTGGGCTGTGGCAGACGCTTGCATTTATACCGATTGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTTCGTTGCAGCGCCTGCCAGATCGTGTGCCGCgcgggcggtgctcg                                                                                                                                  <  1:1107360/150‑7 (MQ=37)
                      gCGACAGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGa                                                                                                                                 <  2:4698067/150‑1 (MQ=255)
                        gACAGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATc                                                                                                                               <  1:5390177/150‑1 (MQ=255)
                        gACAGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATc                                                                                                                               <  1:2429380/150‑1 (MQ=255)
                         aCAGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGACTGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGTGTGCCGCGCGGGCGGTGCTCGATCt                                                                                                                              <  2:6051958/150‑1 (MQ=37)
                         aCAGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGTGTGCCGCGCGGGCGGTGCTCGATCt                                                                                                                              <  2:1013551/150‑1 (MQ=38)
                             cGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCa                                                                                                                          >  2:3211120/1‑150 (MQ=255)
                             cGTGAAGGTGAGTTGTGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGTGTGCCGCGCGGGCGGTGCTCGATCTGGCa                                                                                                                          <  2:2760671/150‑1 (MQ=37)
                                     tGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGc                                                                                                                  >  1:492455/1‑150 (MQ=255)
                                        gTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGAATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCaac                                                                                                               <  2:2675169/150‑3 (MQ=255)
                                          tGGGCTGTGGCAGACGCTTGCTTTTATACCGATTGGGCCGAGAGCGACCTGAGAATGTGGGCGGGGGCGCATTAACGCTCTTGTACGTGGCAGCGCGTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGATGGCGCTGCACTAt                                                                                                             <  2:3450944/150‑1 (MQ=255)
                                           gtgCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCGGCACTAtt                                                                                                            >  2:6992932/3‑150 (MQ=255)
                                           gggCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGATGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTAtt                                                                                                            <  2:6824478/150‑1 (MQ=255)
                                              cAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCACCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGa                                                                                                         <  2:4600574/150‑1 (MQ=255)
                                                   gCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGt                                                                                                    >  2:3281194/1‑150 (MQ=255)
                                                   gCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGt                                                                                                    >  2:5566984/1‑150 (MQ=255)
                                                   gCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGt                                                                                                    >  2:23382/1‑150 (MQ=255)
                                                    cAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCCGGCGCGGCACTAATGGAGGCGTa                                                                                                   >  1:6903806/1‑150 (MQ=255)
                                                         gCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGccc                                                                                              >  1:4337954/1‑150 (MQ=255)
                                                          cTTGCATTTAGCCCGACGGGGCCGAGCGCGACCTGAGCACGGCGGCCGGGGCGCATTACCGCTCTTGTCCGGTGCAGCGCCGGCCAGATCGAGCACCGCCCGCGCGGCCCCCGATCGGGCCGGCGCTGCACTATTGGAGGCGTACGcccc                                                                                             <  2:4110787/150‑1 (MQ=255)
                                                           ttGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCg                                                                                            >  1:6044592/1‑150 (MQ=255)
                                                            tGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGt                                                                                           >  1:1190579/1‑150 (MQ=255)
                                                            tGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGt                                                                                           >  1:6099824/1‑150 (MQ=255)
                                                              cATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGAATTAACGCTCTTGTACGTTGCAGCGCCTGCAAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTca                                                                                         <  2:1351547/150‑1 (MQ=255)
                                                                tttATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGgcg                                                                                                                      >  2:4238340/1‑119 (MQ=255)
                                                                tttATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGgcg                                                                                                                      <  1:4238340/119‑1 (MQ=255)
                                                                 ttATACCGAATGGGCCGAGAGCGACCTGTGCACGTCGGCGGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCGGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACAt                                                                                      <  1:2976459/150‑1 (MQ=255)
                                                                     aCCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCCCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTatc                                                                                  >  1:2791903/1‑150 (MQ=255)
                                                                            gggCCGAGAGCGACCTGAGCAAGTCGGCGGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGc                                                                           <  1:7237670/150‑1 (MQ=255)
                                                                            gggCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGc                                                                           >  2:4792191/1‑150 (MQ=255)
                                                                                   gcgCGACCTGAGCAAGTCGGCCGGGGCGCATTCACGCTCTTTTCCGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTc                                                                    <  2:3432594/148‑1 (MQ=255)
                                                                                   gagCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTc                                                                    <  2:7186254/150‑1 (MQ=255)
                                                                                           tGAGCAAGTCGGCCGGGGCGCATTACCGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGCCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTg                                                            <  2:2527483/150‑1 (MQ=255)
                                                                                            gAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGc                                                           >  2:5279531/1‑150 (MQ=255)
                                                                                             aGCAAGTCGGCGGGGGCGCATTAACGCTCTTGTACGTGGCAGCGCCTGCCAGATCGAGCCCCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGcc                                                          <  2:3298512/150‑1 (MQ=255)
                                                                                              gCATGTCGGCGGGGGCGCATTACCGCTCTTGTAGGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATGTGGCTGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCg                                                         <  1:3813560/150‑1 (MQ=255)
                                                                                              gCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCt                                                                                                                              <  2:808888/81‑1 (MQ=255)
                                                                                              gCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCt                                                                                                                              >  1:808888/1‑81 (MQ=255)
                                                                                               cAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGACACTGTCGTAATCTGCCgg                                                        >  1:7337294/1‑150 (MQ=255)
                                                                                                  gTCGGCCGGGGCGCATTAACGCTCTGGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCGGCACGAGTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCa                                                     >  2:3664822/1‑150 (MQ=255)
                                                                                                     ggCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCCCCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCCCAGTATCATCACGCCACTGGCGTAATCTGCCGGCCAAcg                                                  >  2:5929860/1‑149 (MQ=255)
                                                                                                     ggCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAAct                                                  >  2:5465347/1‑150 (MQ=255)
                                                                                                     ggCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAAct                                                  <  1:6622276/150‑1 (MQ=255)
                                                                                                     ggCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAAct                                                  <  1:4838162/150‑1 (MQ=255)
                                                                                                        cGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTa                                               >  2:6788758/1‑150 (MQ=255)
                                                                                                         ggggCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCCActctct                                              >  2:7316972/1‑148 (MQ=255)
                                                                                                         ggggCGATTTAACGATTTTGTACTTTGAAGCGCCTGCAATATCGAGCACCGCCCGGGCGGAAAACGTTGTGGCTGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTAt                                              <  1:1408631/150‑1 (MQ=255)
                                                                                                               ctttAACGCTCTTGTACGTTGCAGCGCCTGCCTTATGGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCcaca                                        <  1:5249847/148‑1 (MQ=255)
                                                                                                               cATTACCGCTCTTGTAGGTTGCCGCGCCTGCCTGATCGAGCCCCGCCCGCGCGGAACACTATGTGGCTGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCcaca                                        <  1:3875938/150‑1 (MQ=255)
                                                                                                                aTTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAg                                       <  2:6339328/150‑1 (MQ=255)
                                                                                                                    aCGCTCTTGTACGTTGCAGCGCCTGCCCGATCCAGCACCCCCCCCGCGGCACACCATCTGGCAGGCGCTGCACTAATGGAGGCGTACGCCCCGGCACAATATCATCACGCCACAGTCGTAATATGCCGGCCCACTCTATGCCACAGGTCa                                   >  1:5841547/1‑150 (MQ=255)
                                                                                                                        tctTGTACGTTGCAGCGCCTGCCAGATAGATCACCGCCCGCGCGGCCCACGATCTGGCAGGCGCTGCCCTATTTGAGGCGTAAGCCCCCGCACAATATAATAACGCCCCAGTCGTAATCTGCCGGCCAACACTAAGCCACAGGTCACCCg                               >  1:4602721/1‑150 (MQ=255)
                                                                                                                           tGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCa                            <  2:3067239/150‑1 (MQ=255)
                                                                                                                             aatagtGCAGCGCCTGCCAGATCGTGTGCCGCGCGGGCGGTGCTCGATCTGGCAGGCGCTGCACTATTGGAGGCGGACGCCCCGTCACAATATCATCACGCCACTGTCGTAATCAGCCGGCCAACGCAATGCCACAGGTCACCCGGCATg                          >  2:3263006/6‑150 (MQ=18)
                                                                                                                                gTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTc                       <  2:5227147/150‑1 (MQ=255)
                                                                                                                                   gCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTTCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGc                    <  2:4872320/150‑1 (MQ=255)
                                                                                                                                   gCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGc                    <  1:5188992/150‑1 (MQ=255)
                                                                                                                                   gCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGc                    <  1:5491802/150‑1 (MQ=255)
                                                                                                                                    cAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGcc                   <  2:6242915/150‑1 (MQ=255)
                                                                                                                                     aGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCt                  <  2:7020731/150‑1 (MQ=255)
                                                                                                                                       cgcCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGc                <  2:6300202/150‑1 (MQ=255)
                                                                                                                                        gcCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCt               >  2:7306375/1‑150 (MQ=255)
                                                                                                                                        gcCTGCCAGACCGAGCACCGCCCGCGCGGCACCCGATCTGGCCGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCt               <  1:921252/150‑1 (MQ=255)
                                                                                                                                           tGCCAGATCGTGTGCCGCGCGGGCGGTGCTCGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGt            >  1:3216835/1‑150 (MQ=37)
                                                                                                                                           tGCCAGATCGTGTGCCGCGCGGGCGGTGCTCGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACAGTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGt            >  1:890330/1‑150 (MQ=37)
                                                                                                                                           tGCCAGATCGAGCACCGCCCGCGCGGCACACTATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGt            <  1:5369105/150‑1 (MQ=255)
                                                                                                                                           tGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGt            >  2:1063457/1‑150 (MQ=255)
                                                                                                                                            gCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTc           >  1:4869532/1‑150 (MQ=255)
                                                                                                                                             ccAGATCGTGTGCCGCGCGGGCGGTGCTCGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCt          >  1:6660706/1‑150 (MQ=37)
                                                                                                                                             ccAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCg                                                         <  1:2084960/103‑1 (MQ=255)
                                                                                                                                             ccAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCg                                                         >  2:2084960/1‑103 (MQ=255)
                                                                                                                                             ccAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCt          >  2:1242203/1‑150 (MQ=255)
                                                                                                                                              cAGATCGTGTGCCGCGCGGGCGGTGCTCGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGt                                                                  <  2:5219723/93‑1 (MQ=18)
                                                                                                                                              cAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTCTGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCtt         >  2:5749804/1‑150 (MQ=255)
                                                                                                                                              cAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCtt         >  1:5580222/1‑150 (MQ=255)
                                                                                                                                               cgATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTtg        <  1:5979660/149‑1 (MQ=255)
                                                                                                                                               aGATCGTGTGCCGCGCGGGCGGTGCTCGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTa                                               <  2:1859534/111‑1 (MQ=25)
                                                                                                                                               aGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTtg        >  1:3112495/1‑150 (MQ=255)
                                                                                                                                               aGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTtg        >  1:90689/1‑150 (MQ=255)
                                                                                                                                                gATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCtt         >  2:4226946/1‑148 (MQ=255)
                                                                                                                                                gATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTtgc       >  2:1602709/1‑150 (MQ=255)
                                                                                                                                                 aTCGAGCACCGCCCGCGCGGCACCCGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTtgct      <  1:3988286/150‑1 (MQ=255)
                                                                                                                                                 aTCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTtgct      <  1:6815184/150‑1 (MQ=255)
                                                                                                                                                 aTCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTtgct      >  1:5788155/1‑150 (MQ=255)
                                                                                                                                                 aTCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTtgct      <  2:4472852/150‑1 (MQ=255)
                                                                                                                                                  tCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACGAATGGAGGCGTACGCCCCCTCACATTATCATCACCCCACTGTCGTAATCAGCCGGCCAACTCGATGCCACAGGTCACCCGGCATGCTCCGCCGGCGGGTCTtgctg     >  2:6868095/1‑150 (MQ=255)
                                                                                                                                                   cGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTtgctgc    >  2:434291/1‑150 (MQ=255)
                                                                                                                                                    gAGCACCGCCCGCGCGGCACACGATCGGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTTgctgct   <  1:4702343/150‑1 (MQ=255)
                                                                                                                                                     aGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTTGCTGCTc  >  2:5884435/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCGGGCATGCCGATGCAGCCGGGCGACAGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTTGCTGCTC  >  NC_002947/4349146‑4349444

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: