Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,349,298 A→G 12.8% intergenic (+36/‑120) 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,2980AG12.8% 581.4 / 65.0 242intergenic (+36/‑120)PP_3817/PP_3818polyamine ABC transporter ATP‑binding protein/phosphate‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (101/96);  new base G (8/22);  total (116/126)
Fisher's exact test for biased strand distribution p-value = 1.75e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.74e-03

Only 100 of 250 total aligned reads displayed.

GGCATGCCGATGCAGCCGGGCGACAGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTTGCTGCTCACC  >  NC_002947/4349149‑4349447
                                                                                                                                                     |                                                                                                                                                     
ggCATGCCGATGCAGCCGGGCGACAGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCa                                                                                                                                                       <  2:3470291/150‑1 (MQ=255)
 gCATGCCGATGCAGCCGGGCGACAGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCAc                                                                                                                                                      <  2:468199/150‑1 (MQ=255)
  cATGCCGATGCAGCCGGGCGACAGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCAcc                                                                                                                                                     <  2:8659346/150‑1 (MQ=255)
  cATGCCGATGCAGCCGGGCGACAGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCAcc                                                                                                                                                     <  2:359177/150‑1 (MQ=255)
   aTGCCGATGCAGCCGGGCGACAGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCAccg                                                                                                                                                    >  1:5647732/1‑150 (MQ=255)
           gCAGCCGGGCGACAGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCcgcgcg                                                                                                                                            <  2:810469/150‑1 (MQ=255)
             aGCCGGGCGACAGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGc                                                                                                                                          <  1:4015559/150‑1 (MQ=255)
             aGCCGGGCGACAGCGTGAAGGTGCGTTGGCCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGc                                                                                                                                          <  2:6900913/150‑1 (MQ=255)
              gCCGGGCGACAGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGca                                                                                                                                         <  1:1003510/150‑1 (MQ=255)
               ccGGGCGACAGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGcac                                                                                                                                        <  2:3843193/150‑1 (MQ=255)
                      aCAGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGTGTGCCGCGCGGGCGGTGCTCGATCt                                                                                                                                 <  1:5762320/150‑1 (MQ=38)
                      aCAGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCt                                                                                                                                 <  1:316327/150‑1 (MQ=255)
                      aCAGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGGTGCAGCGCCTGCCAGATCGTGTGCCGCGCGGGCGGTGCTCGATCt                                                                                                                                 <  1:6831476/150‑1 (MQ=37)
                       cAGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGGGTGCCGCGCGGGCGGTGCTCGATCTg                                                                                                                                <  2:1720143/150‑1 (MQ=38)
                       cAGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTg                                                                                                                                <  1:8235288/150‑1 (MQ=255)
                        aGCGTTAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTgg                                                                                                                               <  2:491517/150‑1 (MQ=255)
                        aGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTgg                                                                                                                               <  1:8704456/150‑1 (MQ=255)
                        aGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTgg                                                                                                                               <  1:2282152/150‑1 (MQ=255)
                        aGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCGTGCATTTATACCGACGGGGCCGAGAGCGACCTGAGCACGTCGGCCGGGGCGCATTACCGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTgg                                                                                                                               <  2:7498402/150‑1 (MQ=255)
                         gCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGTGTGCCGCGCGGGCGGTGCTCGATCTGGc                                                                                                                              <  2:929704/150‑1 (MQ=37)
                         gCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGTGTGCCGCGCGGGCGGTGCTCGATCTGGc                                                                                                                              <  2:6607449/150‑1 (MQ=37)
                           gTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGTGTGCCGCGCGGGCGGTGCTCGATCTGGCAg                                                                                                                            <  1:7517447/150‑1 (MQ=37)
                           gTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAg                                                                                                                            <  2:8478411/150‑1 (MQ=255)
                            tGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCCCGCGGCACACGATCTGGCAgg                                                                                                                           >  1:7728590/1‑150 (MQ=255)
                                 gTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTACCGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTg                                                                                                                      >  2:4472748/1‑150 (MQ=255)
                                 gTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTg                                                                                                                      >  2:182041/1‑150 (MQ=255)
                                  tGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGc                                                                                                                     >  2:376005/1‑150 (MQ=255)
                                  tGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGc                                                                                                                     >  2:217420/1‑150 (MQ=255)
                                   gCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGACTGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGGTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCa                                                                                                                    <  2:4973330/150‑1 (MQ=255)
                                   gCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCa                                                                                                                    >  2:2044407/1‑150 (MQ=255)
                                      ttGGGCTGTGGCTGACGGTTTCTTTTTTACCGATTGGGCCGAGAGGTACCTGAGCATGTCGGCGGGGGCGCATTAACGCTCTTGTAGGTTGCAGCGCGTGCCAGATCGAGCCCCGCCCGGGCGGCACACGATCTGGCAGGCGCTGCACTa                                                                                                                 <  2:6079166/150‑1 (MQ=255)
                                        gggCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTAtt                                                                                                               <  2:2071307/150‑1 (MQ=255)
                                           cAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGa                                                                                                            <  2:1342178/150‑1 (MQ=255)
                                           cAGGGGCCGACGCTTGCATTTAGACCGACTGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTGGCCGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGa                                                                                                            <  2:7973316/150‑1 (MQ=255)
                                              tGGCAGACGCTTGAATTTATACCGAATGGGCCGAGAGCGCCCTGAGCAAGGCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGc                                                                                                         <  1:7930349/150‑1 (MQ=255)
                                                gCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGt                                                                                                       >  1:6465194/1‑150 (MQ=255)
                                                       cTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTGGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGcccc                                                                                                <  1:8797102/150‑1 (MQ=255)
                                                        ttGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCg                                                                                               >  2:8440733/1‑150 (MQ=255)
                                                        ttGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCg                                                                                               >  2:5830726/1‑150 (MQ=255)
                                                              ttATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCGGCCAGATCGAGCCCCCCCCCCGCGGCCAACGATCTGGCAGGCGCTGCCCGAGTGGAGGCGTACCCCCCGGCACAt                                                                                         >  1:391159/1‑150 (MQ=255)
                                                                 taCCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTgg                                                                                                                               >  1:8189994/1‑109 (MQ=255)
                                                                 taCCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCCCGCGGCACACGATCTGGCAGGCGCTGCACTATTTGAGGCGTACGCCCCGTCACATTat                                                                                      >  1:5185968/1‑150 (MQ=255)
                                                                  aCCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCCCGCGGCACACGATCTTGCAGGCGCTGCACTATTGGAGGCGTACGCCCCCGCACATTata                                                                                     >  1:143228/1‑149 (MQ=255)
                                                                   ccGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTatca                                                                                    >  2:2666219/1‑150 (MQ=255)
                                                                     gAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCGGCCAGATCGAGCCCCCCCCGCGCGGCACACGCGCGGGCAGGCGCGGCACGATTGGCGGCGTACGCCCCGGCCCAGTaacatc                                                                                  >  1:94004/1‑150 (MQ=255)
                                                                     catggggCCGAGTGCGACCTGAGCAAGGCGGCGGGGGCGAATTAACGCTCTTGTAGTTTGCAGCGCCTGCCAGAGCGCGCACCGCCCGCGCGGCACACTATGTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTatcatc                                                                                  <  1:2343875/146‑1 (MQ=255)
                                                                       aTGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCAc                                                                                >  2:1592204/1‑150 (MQ=255)
                                                                         gggCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGc                                                                              >  1:3432516/1‑150 (MQ=255)
                                                                          ggCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCGGCCCGATAGAGCCCCGCCCCCGCGGCACACCATCTGGCAGGCGCGGCACTCTTTGAGGCGTACGCCCCGGCACAATATCATCACGcc                                                                             >  1:2912986/1‑150 (MQ=255)
                                                                           gCCGAGAGCGACCTGAGCAAGTCGGCCGGAGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCa                                                                            >  2:1319182/1‑150 (MQ=255)
                                                                             cgtgagCGACCTGAGCATGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCTGCGCCTGCCATATCGAGCACCGCCCGCGGGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACt                                                                          <  1:1389814/147‑1 (MQ=255)
                                                                                         gAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCCCCCCCGCGGCACACCATCTGGCAGGCGCGGCCCGATTTGAGGCGTACGCCCCGGCACAATATAATAAAGCCCCAGTCGTAATCTGc                                                              >  1:4836289/1‑150 (MQ=255)
                                                                                          aGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGcc                                                             >  2:6396963/1‑150 (MQ=255)
                                                                                           gCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACCCGATCTGGCAGGCGCGGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCg                                                            >  2:7719753/1‑150 (MQ=255)
                                                                                           gCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGa                                                                                                            <  2:6452254/102‑1 (MQ=255)
                                                                                                tgggCGGGGGCGATTTAACGCTCTTGTAGGTTGAAGCGCCTGCATTATCGCGCACCGCCCGCGCGGCACACTATGTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCaa                                                       <  1:1890256/148‑1 (MQ=255)
                                                                                                tCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCaa                                                       >  2:8702887/1‑150 (MQ=255)
                                                                                                   gCCGGGGCGCATTAACGCTCTTGTCCGTTGCAGCGCCTGCCAGATCGAGCCCCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAActc                                                    <  1:733992/150‑1 (MQ=255)
                                                                                                    ccGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAActct                                                   >  2:4007314/1‑150 (MQ=255)
                                                                                                     cGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTa                                                  <  2:767867/150‑1 (MQ=255)
                                                                                                       ggggGCTTTAACGATGTTTTACGTGGCAGGGCGTGCCATATGGTGCACCGCCCGCGGGGCACACTTTCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATg                                                <  1:1543282/150‑1 (MQ=255)
                                                                                                           gctttAACGGTCTTTTAGGTTGCAGCGCCTGCCAGATCGAGCCCCGCCCGCGCGGAACACGATCGGGCTGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCcac                                            <  1:4589921/147‑1 (MQ=255)
                                                                                                           gcATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCTGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCcac                                            <  2:856776/150‑1 (MQ=255)
                                                                                                            cATTAACGCTGTTGTAGGTTGCCGCGCCTGCCAGTTCGCGCACCGCCCGCGGGGCACACGTTCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCcaca                                           <  1:2504984/150‑1 (MQ=255)
                                                                                                            cATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCcaca                                           <  2:5677725/150‑1 (MQ=255)
                                                                                                             aTTATCGCTCTTGTAGTTTGCAGCGCCTGCCATATGTAGCACGGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAg                                          <  1:4565667/150‑1 (MQ=255)
                                                                                                             aTTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAg                                          >  2:3135155/1‑150 (MQ=255)
                                                                                                              ttAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCCCCGCCCGCGCGGCACACTATCTGGCTGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAgg                                         <  1:1188187/150‑1 (MQ=255)
                                                                                                              ttAACGCTCTTGTACGTTGCAGCGCCAGCCAGATCGATCACCGCCCCCCCGGCACACGGTCTGGCAGGCGCTGCACTAATTGAGGCGTACGCCCCCTCACAATAAAAACACGCCACTGTCGTAATATTCCGGCCAACTCAATTCCACAgg                                         >  1:1250798/1‑150 (MQ=255)
                                                                                                                   gCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCAcc                                    >  1:7191277/1‑150 (MQ=255)
                                                                                                                    ctctTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCAccc                                   <  2:7866081/150‑1 (MQ=255)
                                                                                                                     tctTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCg                                  >  1:6567516/1‑150 (MQ=255)
                                                                                                                       ttGTCCTTTGCCGCGCCTGCCAGATCGCGCACCGCCCGGGGGGCACACGATCTGGCTGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGc                                <  1:2357367/150‑1 (MQ=255)
                                                                                                                          tACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATg                             >  1:5003455/1‑150 (MQ=255)
                                                                                                                          aatagtGCAGCGCCAGCCAGATCCTGTTCCCCGCCGGCGGTGCACGATCTGGCAGGCGCGGCACAATTGGAGGCCTACGCCCCCTCACCATATAATCACCCCCCAGTCGTAAACTGCCGGCCAACTCTATGCCACAGGTCACCCGGCaag                             >  2:4546273/6‑148 (MQ=255)
                                                                                                                                 cAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGcc                      <  1:4249783/150‑1 (MQ=255)
                                                                                                                                   gcgcCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTg                    >  1:1685025/1‑150 (MQ=255)
                                                                                                                                    cgcCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGc                   <  2:6229919/150‑1 (MQ=255)
                                                                                                                                      ccTGCCAGTTCGAGCACCGCCCGCGCGCAACACTTTCTGGCAGGCGCTGCACTATGGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCtt                 <  1:6776169/150‑1 (MQ=255)
                                                                                                                                       cagacAGATCGAGCAAAGCCAGCGAGGCACACGAAATGGCAGGCGCTGCACTATTCGAGGAGTACGCCCCGTCACATTATAATCACGCCACTGTCGTAATCTGACGCCCAACTCTATGCCACAGGTCACCCTGCATGCTCCGCCAGCTTg                >  2:5128324/5‑150 (MQ=255)
                                                                                                                                       cTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTg                >  2:2251813/1‑150 (MQ=255)
                                                                                                                                        cgCCAGTCCGAGCCCCGCCCGCGCGCCCCACGCCGGGGCAGCCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGt               <  1:3685599/149‑1 (MQ=255)
                                                                                                                                         gCCAGATCGTGTGCCGCGCGGGCGGTGCTCGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTc              >  1:7048218/1‑150 (MQ=37)
                                                                                                                                          ccAGATCGTGTGCCGCGCGGGCGGTGCTCGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCt             >  2:4203222/1‑150 (MQ=37)
                                                                                                                                          ccAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCt             >  1:5112480/1‑150 (MQ=255)
                                                                                                                                           cAGATCGTGTGCCGCGCGGGCGGTGCTCGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCtt            >  2:1925485/1‑150 (MQ=38)
                                                                                                                                           cAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTATCtt            >  1:6215451/1‑150 (MQ=255)
                                                                                                                                            aGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTtg           >  2:2594592/1‑150 (MQ=255)
                                                                                                                                            aGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTtg           >  1:6251485/1‑150 (MQ=255)
                                                                                                                                             gATCGCGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTtgc          <  2:3070094/150‑1 (MQ=255)
                                                                                                                                              aTCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTtgct         >  1:8666336/1‑150 (MQ=255)
                                                                                                                                              aTCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTAGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTtgct         >  1:6030757/1‑150 (MQ=255)
                                                                                                                                               tCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTtgctg        <  1:8014994/150‑1 (MQ=255)
                                                                                                                                                 gAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTCTCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTTgctgct      <  2:7183867/150‑1 (MQ=255)
                                                                                                                                                   gCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTTGCTGCTCa    >  2:5301399/1‑150 (MQ=255)
                                                                                                                                                    cACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCcaca                                           <  1:4357368/110‑1 (MQ=255)
                                                                                                                                                    cACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCcaca                                           >  2:4357368/1‑110 (MQ=255)
                                                                                                                                                    cACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTTGCTGCTCAc   >  1:6016230/1‑150 (MQ=255)
                                                                                                                                                     aCCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTTGCTGCTCAcc  >  1:6676333/1‑150 (MQ=255)
                                                                                                                                                     aCCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTTGCTGCTCAcc  >  2:8043242/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGCATGCCGATGCAGCCGGGCGACAGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTTGCTGCTCACC  >  NC_002947/4349149‑4349447

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: