Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,346,270 C→T 13.8% intergenic (+35/‑121) 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,346,2700CT13.8% 239.5 / 25.7 94intergenic (+35/‑121)PP_3817/PP_3818polyamine ABC transporter ATP‑binding protein/phosphate‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (40/41);  new base T (2/11);  total (42/52)
Fisher's exact test for biased strand distribution p-value = 3.32e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

CGGGCATGCCGATGCAGCCGGGCGACAGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTTGCTGCTC  >  NC_002947/4346120‑4346417
                                                                                                                                                      |                                                                                                                                                   
cggGCATGCCGATGCAGCCGGGCGACAGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGc                                                                                                                                                     <  1:1415750/150‑1 (MQ=255)
 gggCATGCCGATGCAGCCGGGCGACAGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGc                                                                                                                                                     >  2:4858671/1‑150 (MQ=255)
 gggCATGCCGATGCAGCCGGGCGACAGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGc                                                                                                                                                     <  1:338381/150‑1 (MQ=255)
   gCATGCCGATGCAGCCGGGCGACAGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCAc                                                                                                                                                   <  1:2964569/150‑1 (MQ=255)
    cATGCCGATCCAGCCGGGCGACAGCGTGAAGGTGCGATGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCAcc                                                                                                                                                  <  2:441013/150‑1 (MQ=255)
     aTGCCGATGCAGCCGGGCGACAGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCAccg                                                                                                                                                 >  1:3907179/1‑150 (MQ=255)
          gATGCAGCCGGGCGACAGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCAccgcccgc                                                                                                                                            >  1:1266910/1‑150 (MQ=255)
           aTGCAGCCGGGCGACAGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCcgcg                                                                                                                                           >  1:3983605/1‑150 (MQ=255)
                gCCGGGCGACAGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGca                                                                                                                                      >  2:1883946/1‑150 (MQ=255)
                     gCGACAGCTTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGa                                                                                                                                 <  1:2806058/150‑1 (MQ=255)
                     gCGACAGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGTGTGCCGCgcgggcggtgctcga                                                                                                                                 <  1:1274445/150‑8 (MQ=38)
                     gCGACAGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGTGTGCCGCgcgggcggtgctcga                                                                                                                                 <  1:673983/150‑8 (MQ=38)
                     gCGACAGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGa                                                                                                                                 <  2:3983605/150‑1 (MQ=255)
                       gACAGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATc                                                                                                                               >  1:977357/1‑150 (MQ=255)
                        aCAGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCt                                                                                                                              <  2:2190143/150‑1 (MQ=255)
                         cAGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGTGTGCCGCGCGGGCGGTGCTCGATCTg                                                                                                                             <  1:55035/150‑1 (MQ=38)
                         cAGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGTGTGCCGCGCGGGCGGTGCTCGATCTg                                                                                                                             <  1:601581/150‑1 (MQ=38)
                         cAGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGTGTGCCGCGCGGGCGGTGCTCGATCTg                                                                                                                             <  2:3783944/150‑1 (MQ=38)
                          aGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGTGTGCCGCGCGGGCGGTGCTCGATCTgg                                                                                                                            <  1:4398709/150‑1 (MQ=37)
                          aGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCt                                                                                                                              <  1:2789203/148‑1 (MQ=255)
                          aGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCt                                                                                                                              >  2:2789203/1‑148 (MQ=255)
                          aGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCAGCCAGATCGAGCACCGCCCGCGCGGCCCACCATCTgg                                                                                                                            >  1:577384/1‑150 (MQ=255)
                           gCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGTGTGCCGCGCGGGCGGTGCTCGATCTGGc                                                                                                                           <  2:1635993/150‑1 (MQ=37)
                           gCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGTGTGCCGCGCGGGCGGTGCTCGATCTGGc                                                                                                                           <  2:1635994/150‑1 (MQ=37)
                           gCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGTGTGCCGCGCGGGCGGTGCTCGATCTGGc                                                                                                                           <  2:1636046/150‑1 (MQ=37)
                           gCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGTGTGCCGCGCGGGCGGTGCTCGATCTGGc                                                                                                                           <  2:3629765/150‑1 (MQ=37)
                           gCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGc                                                                                                                           <  2:2373822/150‑1 (MQ=255)
                            cGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCa                                                                                                                          >  1:2657987/1‑150 (MQ=255)
                             gTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCGTTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAg                                                                                                                         >  1:4565147/1‑150 (MQ=255)
                             gTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAg                                                                                                                         <  2:3041083/150‑1 (MQ=255)
                              tGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAgg                                                                                                                        <  1:1334377/150‑1 (MQ=255)
                               gAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGgc                                                                                                                       <  1:4858454/150‑1 (MQ=255)
                                aaGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGgcg                                                                                                                      <  2:3344918/150‑1 (MQ=255)
                                       gTTGGTCAGTGGCAGACGCTTGAATTTATACCGAATGGCCCGACAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGATATTTTACGTTGCAGCGCCTGCCAGATCGTGTGCCGCGCGGGCGGTGCTCTATCTGGCAGGCGCTGCaac                                                                                                               <  1:1601258/150‑3 (MQ=18)
                                           ggCTGTGTCTGACGCTTGCTTTTATACCTATTGGGCCGAGCGCGACCTTAGCATGTCGGCCGGGGCGATTTAACGATCTTTTACGTTGCAGGGCGTGCATGATGGAGCACCGCCCGGGGGGCACACTTTCTGGCAGGCGCTGCACTATTg                                                                                                           <  2:2261779/150‑1 (MQ=255)
                                             cAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCCCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGa                                                                                                         >  1:1987334/1‑150 (MQ=255)
                                                tGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGc                                                                                                      >  2:1855034/1‑150 (MQ=255)
                                                tGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGc                                                                                                      >  2:1854923/1‑150 (MQ=255)
                                                          ttGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCCGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCg                                                                                            >  2:69929/1‑150 (MQ=255)
                                                              aTTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTcac                                                                                        >  2:3758453/1‑150 (MQ=255)
                                                               tttATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTATTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTcaca                                                                                       <  2:2869410/150‑1 (MQ=255)
                                                                    aCCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTatc                                                                                  >  2:4005258/1‑150 (MQ=255)
                                                                      cGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTatcat                                                                                >  1:4669003/1‑150 (MQ=255)
                                                                               cGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGCGCACCGCCCGCGCGGCACAAGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACt                                                                       >  1:1195033/1‑150 (MQ=255)
                                                                               cGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCGCATTATCATCACGCCACt                                                                       >  2:2192261/1‑150 (MQ=255)
                                                                               cGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACt                                                                       >  2:226590/1‑150 (MQ=255)
                                                                                 agagCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGt                                                                     >  1:2538811/1‑150 (MQ=255)
                                                                                  gcgCGACCTGAGCACGTCGGCCGGGGCGCACTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGCGCACCGCCCGCGCGGCACACGTTCTGGCTGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTc                                                                    <  1:4693034/148‑1 (MQ=255)
                                                                                  gagCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTc                                                                    <  1:4692759/150‑1 (MQ=255)
                                                                                  gagCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTc                                                                    <  1:4692757/150‑1 (MQ=255)
                                                                                  gagCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTATTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTc                                                                    <  1:4695789/150‑1 (MQ=255)
                                                                                      gACCTGAGCATGTGGGCCGGGGCGAATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGTTCGCGCACCGCCCGCGGGGCACCCGATCTGGCTGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTaa                                                                <  1:1883946/150‑1 (MQ=255)
                                                                                      gACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTaa                                                                >  2:3376553/1‑150 (MQ=255)
                                                                                       tccTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCGGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAAt                                                               <  2:3821729/149‑1 (MQ=255)
                                                                                         cTGAGCACGTGGGCGGGGGCGAATTACCGCTTTTTTAGGTGGCCGCGCGTGCATGATCGAGCCCCGCCCGGGCGGAACACGATGTGGCTGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCt                                                             <  1:3498398/150‑1 (MQ=255)
                                                                                         cTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCt                                                             <  2:3769183/150‑1 (MQ=255)
                                                                                            aGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGcc                                                          <  1:4263358/150‑1 (MQ=255)
                                                                                            aGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGcc                                                          <  1:4263283/150‑1 (MQ=255)
                                                                                            aGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGAAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGcc                                                          <  1:4262916/150‑1 (MQ=255)
                                                                                               acgTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGc                                                       <  2:3099328/148‑1 (MQ=255)
                                                                                                       cGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTa                                               <  2:3310659/150‑1 (MQ=255)
                                                                                                       cGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCAGCCAGAAAGAGCACCGCCCCCGCGGCACACCATCTGGCAGGCGCGGCACTATTGGAGGCGTACCCCCCCTCACAATATCAAAAACCCACAGTCGTAATATTCCGGCCaaaacta                                               >  1:2211696/1‑145 (MQ=255)
                                                                                                                ttAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAgg                                      >  1:2532064/1‑150 (MQ=255)
                                                                                                                 tAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGt                                     >  2:4736984/1‑150 (MQ=255)
                                                                                                                  aaCGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTc                                    >  2:907004/1‑150 (MQ=255)
                                                                                                                    cGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCAc                                  <  1:3758453/150‑1 (MQ=255)
                                                                                                                     gCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCAcc                                 >  2:4589688/1‑150 (MQ=255)
                                                                                                                      ctctTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCAccc                                >  1:2455039/1‑150 (MQ=255)
                                                                                                                           gTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCAt                           <  2:3386075/150‑1 (MQ=255)
                                                                                                                                   cAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGcc                   <  1:3969580/150‑1 (MQ=255)
                                                                                                                                     gcgcCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGAAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTg                 <  1:829074/150‑1 (MQ=255)
                                                                                                                                      cgcCTGCCAGATCGAGCACCGCCCGCGCGGCACACGCTCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTAGGCCACAGGTCACCCGGCATGCTCCGCCTGc                <  2:2932154/150‑1 (MQ=255)
                                                                                                                                      cgcCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGc                <  2:374734/150‑1 (MQ=255)
                                                                                                                                      cgcCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGc                <  2:2932719/150‑1 (MQ=255)
                                                                                                                                         cTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTg             >  2:2125776/1‑150 (MQ=255)
                                                                                                                                          tGCCAGATCGTGTGCCGCGCGGGCGGTGCTCGATATGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCAACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGACACCCGGCATGCTCCGACTGCTTGt            >  1:4367160/1‑150 (MQ=25)
                                                                                                                                          tGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGt            >  1:3258506/1‑150 (MQ=255)
                                                                                                                                           gCCGGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCGCCCGGCATGCTCCGCCTGCTTGTc           >  2:4195634/1‑150 (MQ=255)
                                                                                                                                           gCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTc           >  1:4979219/1‑150 (MQ=255)
                                                                                                                                            ccAGATCGTGTGCCGCGCGGGCGGTGCTCGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCt          >  1:1991150/1‑150 (MQ=37)
                                                                                                                                            ccAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCt          >  1:141264/1‑150 (MQ=255)
                                                                                                                                            ccAGATCGAGCACCGCCCGCGCGGCACACGAGCTTGCAGCCGCTGCACTATAGGAGGAGTACGCCCCGTCACATTATCATCACGCCTCTGTCGTAATCTGCAGGCCAACTCTATGCCACAGGTCACGCGGCATGCTCCGCCTGCTGGTCt          >  2:4662020/1‑150 (MQ=255)
                                                                                                                                             cAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCtt         >  2:2052377/1‑150 (MQ=255)
                                                                                                                                              aGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTtg        <  2:4090951/150‑1 (MQ=255)
                                                                                                                                              aGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTtg        <  2:4565147/150‑1 (MQ=255)
                                                                                                                                               gATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTtgc       >  2:326614/1‑150 (MQ=255)
                                                                                                                                               gATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTtgc       >  1:20169/1‑150 (MQ=255)
                                                                                                                                               gATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTtgc       <  2:977357/150‑1 (MQ=255)
                                                                                                                                                 tCGAGCCCCGCCCGCGCGGCACACGATGTGGCTGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTtgctg     <  1:4858671/150‑1 (MQ=255)
                                                                                                                                                 tCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCTGGCATGCTCCGCCTGCTTGTCTtgctg     >  1:183572/1‑150 (MQ=255)
                                                                                                                                                 tCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTtgctg     <  2:3651938/150‑1 (MQ=255)
                                                                                                                                                  cGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTtgctgc    <  2:3741114/150‑1 (MQ=255)
                                                                                                                                                    aGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTTGCTGCTc  >  2:2546456/1‑150 (MQ=255)
                                                                                                                                                     gCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAActct                                                >  2:4534250/1‑103 (MQ=255)
                                                                                                                                                     gCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAActct                                                >  2:4534266/1‑103 (MQ=255)
                                                                                                                                                     gCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAActct                                                <  1:4534266/103‑1 (MQ=255)
                                                                                                                                                     gCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAActct                                                <  1:4534250/103‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                   
CGGGCATGCCGATGCAGCCGGGCGACAGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTTGCTGCTC  >  NC_002947/4346120‑4346417

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: