Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,349,314 A→T 13.0% intergenic (+52/‑104) 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,3140AT13.0% 589.0 / 69.6 228intergenic (+52/‑104)PP_3817/PP_3818polyamine ABC transporter ATP‑binding protein/phosphate‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (101/92);  new base T (17/13);  total (118/110)
Fisher's exact test for biased strand distribution p-value = 6.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.14e-01

Only 100 of 240 total aligned reads displayed.

GCGACAGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTTGCTGCTCACCCTGCTACCACTCGCGA  >  NC_002947/4349168‑4349463
                                                                                                                                                  |                                                                                                                                                     
gCGACAGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGa                                                                                                                                                    <  2:4698067/150‑1 (MQ=255)
  gACAGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATc                                                                                                                                                  <  1:2429380/150‑1 (MQ=255)
   aCAGCGTGTAGGTGCGTTGGGCAGTGGCTGACGCTTGCATTTATACCGAATGGGGCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTAGGTTGCTGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCt                                                                                                                                                 <  2:5463568/150‑1 (MQ=255)
   aCAGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGTGTGCCGCGCGGGCGGTGCTCGATCt                                                                                                                                                 <  2:1013416/150‑1 (MQ=38)
   aCAGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCt                                                                                                                                                 >  2:3117075/1‑150 (MQ=255)
       cGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCa                                                                                                                                             >  2:3211120/1‑150 (MQ=255)
       cGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCa                                                                                                                                             <  1:4211461/150‑1 (MQ=255)
       cGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCAGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCa                                                                                                                                             >  1:1120906/1‑150 (MQ=255)
       cGTGAAGGTGAGTTGTGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGTGTGCCGCGCGGGCGGTGCTCGATCTGGCa                                                                                                                                             <  2:2760671/150‑1 (MQ=37)
        gTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGTGTGCCGCGCGGGCGGTGCTCGATCTGGCAg                                                                                                                                            <  2:6555288/150‑1 (MQ=37)
               tGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGc                                                                                                                                     >  1:492455/1‑150 (MQ=255)
                  gTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGAATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCaac                                                                                                                                  <  2:2675169/150‑3 (MQ=255)
                    tGGGCAGTGGCAGACGCTTGCATTTATACCGATTGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCATATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTAt                                                                                                                                <  2:3450817/150‑1 (MQ=255)
                    tGGGCAGTGGCAGACGCTTGCATTTATACCGACTGGGCCGAGAGCGACCTGAGCATGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTAt                                                                                                                                <  2:3688805/150‑1 (MQ=255)
                     gtgCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCGGCACTAtt                                                                                                                               >  2:6992932/3‑150 (MQ=255)
                        cAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCACCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGa                                                                                                                            <  2:4600574/150‑1 (MQ=255)
                             gCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGt                                                                                                                       >  2:3281177/1‑150 (MQ=255)
                             gCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGt                                                                                                                       >  2:3281333/1‑150 (MQ=255)
                              cAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCGGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTa                                                                                                                      <  2:6799192/150‑1 (MQ=255)
                                   gCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGccc                                                                                                                 >  1:4337954/1‑150 (MQ=255)
                                        cATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGAATTAACGCTCTTGTACGTTGCAGCGCCTGCAAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTca                                                                                                            <  2:1351547/150‑1 (MQ=255)
                                         aTTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTcac                                                                                                           >  2:7336313/1‑150 (MQ=255)
                                          tttATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGgcg                                                                                                                                         <  1:4238340/119‑1 (MQ=255)
                                          tttATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGgcg                                                                                                                                         >  2:4238340/1‑119 (MQ=255)
                                          tttATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTcaca                                                                                                          >  2:5506697/1‑150 (MQ=255)
                                                ccGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTatca                                                                                                    >  1:6127914/1‑150 (MQ=255)
                                                      gggCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGc                                                                                              >  2:4792191/1‑150 (MQ=255)
                                                        gCCGAGAGCGCCCTGAGCATGTCGGCGGGGGCGCATTAACGCTCTTGTCGGTTGCAGCGCCTGCCTGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCa                                                                                            <  1:5384492/150‑1 (MQ=255)
                                                         ccGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGcc                                                                                                                  >  2:2384334/1‑127 (MQ=255)
                                                             gagCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTc                                                                                       <  2:7186254/150‑1 (MQ=255)
                                                                    cTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCt                                                                                >  2:6079175/1‑150 (MQ=255)
                                                                      gagaatgTCGGCCGGGGCGCTTTAACGCTCTTGTACGTTGCCGCGCCTGCCATATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGc                                                                              <  1:2836548/144‑1 (MQ=255)
                                                                        gCATGTCGGCGGGGGCGCATTACCGCTCTTGTAGGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATGTGGCTGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCg                                                                            <  1:3813560/150‑1 (MQ=255)
                                                                        gCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCt                                                                                                                                                 >  1:808888/1‑81 (MQ=255)
                                                                         cAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGACACTGTCGTAATCTGCCgg                                                                           >  1:7337294/1‑150 (MQ=255)
                                                                               ggCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGTCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAAct                                                                     >  2:7044631/1‑150 (MQ=255)
                                                                               ggCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCCCCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCCCAGTATCATCACGCCACTGGCGTAATCTGCCGGCCAAcg                                                                     >  2:5929860/1‑149 (MQ=255)
                                                                               ggCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAAct                                                                     <  1:6622276/150‑1 (MQ=255)
                                                                                   ggggCGATTTAACGATTTTGTACTTTGAAGCGCCTGCAATATCGAGCACCGCCCGGGCGGAAAACGTTGTGGCTGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTAt                                                                 <  1:1408631/150‑1 (MQ=255)
                                                                                       cgcATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACTATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCca                                                             <  2:701452/150‑1 (MQ=255)
                                                                                          tttAACGCTGTTGTATTTTGATGCGCGTGCCATTTCGAGCACGGCCGGGGGGGCAAACTTTGTGTCTGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCTCTGTCGTAATCTGCCGGCCAACTCTATGCCACAg                                                          <  1:3944921/149‑1 (MQ=255)
                                                                                           ttAACGCTGTTTTACGTTGCAGCGCCTGCCTGATCGAGCACCGCCCGGGCGGCACACTATCTGGCTGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAgg                                                         <  1:5982844/150‑1 (MQ=255)
                                                                                           ttAACGCTCTTGTACGTTGCAGCGCCAGCCAGAACGAGCACCCCCCCCGCGGCACACCATCTGGCAGGCGCTGCACTATTTGAGGCGTAAGCCCCCTCACATTATCAACACCCCACAGTCGTAATAATCCGGCCAACACTAATCCCCAgg                                                         >  1:7020731/1‑150 (MQ=255)
                                                                                             aaCGCTCTTGTACGTTGCAGCGCCTGCCAGATAGAGCACCGCCCGCGCGGCACACGATATGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGGCACATTATCATCACGCCCCTGTCGTAATATTCCGGGCAACTCTATGCCACAgggc                                                       >  1:909186/1‑148 (MQ=255)
                                                                                              aCGCTCTTGTACGTTGCAGCGCCTGCCCGATCCAGCACCCCCCCCGCGGCACACCATCTGGCAGGCGCTGCACTAATGGAGGCGTACGCCCCGGCACAATATCATCACGCCACAGTCGTAATATGCCGGCCCACTCTATGCCACAGGTCa                                                      >  1:5841547/1‑150 (MQ=255)
                                                                                                gCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCAcc                                                    >  2:6832379/1‑150 (MQ=255)
                                                                                                  tctTGTACGTTGCAGCGCCTGCCAGATCGAGCACCCCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCg                                                  >  2:2560973/1‑150 (MQ=255)
                                                                                                  tctTGTACGTTGCAGCGCCTGCCAGATAGATCACCGCCCGCGCGGCCCACGATCTGGCAGGCGCTGCCCTATTTGAGGCGTAAGCCCCCGCACAATATAATAACGCCCCAGTCGTAATCTGCCGGCCAACACTAAGCCACAGGTCACCCg                                                  >  1:4602721/1‑150 (MQ=255)
                                                                                                          gTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTc                                          <  2:5227208/150‑1 (MQ=255)
                                                                                                           ttGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTcc                                         <  1:6079175/150‑1 (MQ=255)
                                                                                                             gCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTTCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGc                                       <  2:4872320/150‑1 (MQ=255)
                                                                                                             gCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGc                                       <  1:5491802/150‑1 (MQ=255)
                                                                                                             gCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGc                                       <  1:5188992/150‑1 (MQ=255)
                                                                                                              cAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGcc                                      <  2:6242915/150‑1 (MQ=255)
                                                                                                               aGCGCCTGCCAGATCGAGCACCGCCCGCGGGGCACACGATCTGGCTGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCt                                     <  1:865223/150‑1 (MQ=255)
                                                                                                                  gcCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCt                                  <  2:6323700/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)
                                                                                                                     tGCCAGATCGTGTGCCGCGCGGGCGGTGCTCGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGt                               >  1:3216835/1‑150 (MQ=37)
                                                                                                                     tGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGt                               <  1:5369019/150‑1 (MQ=255)
                                                                                                                     tGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGt                               >  2:1063457/1‑150 (MQ=255)
                                                                                                                      gCCAGATCGTGTGCCGCGCGGGCGGTGCTCGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTACTCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTc                              >  2:5702748/1‑150 (MQ=37)
                                                                                                                      gCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTc                              >  1:4869532/1‑150 (MQ=255)
                                                                                                                       ccAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCg                                                                            >  2:2084960/1‑103 (MQ=255)
                                                                                                                       ccAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCg                                                                            <  1:2084960/103‑1 (MQ=255)
                                                                                                                       ccAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCt                             >  2:1242203/1‑150 (MQ=255)
                                                                                                                        cAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCtt                            <  1:5506697/150‑1 (MQ=255)
                                                                                                                        cAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCtt                            >  1:5580838/1‑150 (MQ=255)
                                                                                                                        cAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCtt                            >  2:5749815/1‑150 (MQ=255)
                                                                                                                         cgATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTtg                           <  1:5979660/149‑1 (MQ=255)
                                                                                                                         aGATCGTGTGCCGCGCGGGCGGTGCTCGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTa                                                                  >  1:1859534/1‑111 (MQ=25)
                                                                                                                         aGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTtg                           >  1:3112495/1‑150 (MQ=255)
                                                                                                                         aGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTtg                           >  2:6873015/1‑150 (MQ=255)
                                                                                                                          gATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTtgc                          >  2:1602709/1‑150 (MQ=255)
                                                                                                                          gATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTtgc                          <  2:4116383/150‑1 (MQ=255)
                                                                                                                           aTCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTtgct                         <  1:6815184/150‑1 (MQ=255)
                                                                                                                           aTCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTtgct                         <  2:4472852/150‑1 (MQ=255)
                                                                                                                            tCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTtgctg                        >  2:3333770/1‑150 (MQ=255)
                                                                                                                             cgtgtgccgcgcggGCGGTGCTCGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTtgctgc                       >  2:2522053/7‑150 (MQ=38)
                                                                                                                             cGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTtgctgc                       >  2:2908490/1‑150 (MQ=255)
                                                                                                                              gtgtgccgcgcggGCGGTGCTCGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTTgctgct                      >  1:6792473/6‑150 (MQ=38)
                                                                                                                              gAGCACCGCCCGCGCGGCACACGATCGGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTTgctgct                      <  1:4702343/150‑1 (MQ=255)
                                                                                                                               tgtgccgcgcggGCGGTGCTCGATCTGGCGGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTTGCTGCTc                     >  1:5804699/5‑150 (MQ=39)
                                                                                                                               tgtgccgcgcggGCGGTGCTCGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTTGCTGCTc                     >  1:3524905/5‑150 (MQ=39)
                                                                                                                                  aCCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTa                                                                                    <  2:2244701/84‑1 (MQ=255)
                                                                                                                                  aCCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTa                                                                                    <  2:5994771/84‑1 (MQ=255)
                                                                                                                                  aCCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAAct                                                                     <  2:5946405/99‑1 (MQ=255)
                                                                                                                                  aCCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAAct                                                                     >  1:5946273/1‑99 (MQ=255)
                                                                                                                                  aCCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTTGCTGCTCAcc                  >  1:2438218/1‑150 (MQ=255)
                                                                                                                                      cccgcGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTTGCTGCTCACCCTGc              <  2:4337954/150‑1 (MQ=255)
                                                                                                                                        cgcgcgGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTTGCTGCTCACCCTGCTa            <  2:1817689/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                                                                       <  2:3982571/84‑1 (MQ=255)
                                                                                                                                               cacaCGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTTGCTGCTCACCCTGCTACCACTcg     >  2:5273311/1‑150 (MQ=255)
                                                                                                                                                acacGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACATGTCACCCGGCATGCTCAGCCTGCTTGTCTTGCTGCTCACCCTGCTACCACTcgc    >  1:5121669/1‑150 (MQ=255)
                                                                                                                                                 cacGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTTGCTGCTCACCCTGCTACCACTcgcg   >  2:5597729/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
GCGACAGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTTGCTGCTCACCCTGCTACCACTCGCGA  >  NC_002947/4349168‑4349463

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: