Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,421,080:1 +G 5.2% intergenic (+29/+90) PP_3901 → / ← PP_3902 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,421,0801.G 5.2% 1889.7 / 87.6 423intergenic (+29/+90)PP_3901/PP_3902hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (150/251);  new base G (7/15);  total (157/266)
Fisher's exact test for biased strand distribution p-value = 6.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

Only 100 of 433 total aligned reads displayed.

GGAATCCCAGGAAGACTGGGGCGCGCTCGTCGAGGCCATGGAAGGCACCGAGGGCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGCATACTCATACTCGCTCCATCACCGGATAGACACCGGCCAGCAACAAAGGCAGCGAAAGTCGCCA  >  NC_002947/4420931‑4421229
                                                                                                                                                      |                                                                                                                                                     
ggAATCCCAGGAAGACTGGGGCGCGCTCGTCGAGGCCATGGAAGGCACCGAGGGCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGa                                                                                                                                                        <  2:1116957/150‑1 (MQ=255)
ggAATCCCAGGAAGACTGGGGCGCGCTCGTCGAGGCCATGGAAGGCACCGAGGGCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGa                                                                                                                                                        >  1:2491518/1‑150 (MQ=255)
 gAATCCCAGGAAGACTGGGGCGCGCTCGTCGAGGCCATGGAAGGCACCGAGGGCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAg                                                                                                                                                      >  2:6050520/1‑150 (MQ=255)
      ccAGGAAGACTGGGGCGCGCTCGTCGAGGCCATGGAAGGCACCGAGGGCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGc                                                                                                                                                 >  1:5690897/1‑150 (MQ=255)
       cAGGAAGACTGGGGCGCGCTCGTCGAGGCCATGGAAGGCACCGAGGGCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCt                                                                                                                                                >  1:6293623/1‑150 (MQ=255)
           aaGACTGGGGCGCGCTCGTCGAGGCCATGGAAGGCACCGAGGGCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCTCGAGGCGGGcaaggt                                                                                                                                            <  1:3920083/150‑6 (MQ=255)
           aaGACTGGGGCGCGCTCGTCGAGGCCATGGAAGGCACCGAGGGCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCTCGAGGCGGGcaaggt                                                                                                                                            <  2:2350057/150‑6 (MQ=255)
           aaGACTGGGGCGCGCTCGTCGAGGCCATGGAAGGCACCGAGGGCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCTCGAGGCGGGcaaggt                                                                                                                                            <  1:2388759/150‑6 (MQ=255)
            aGACTGGGGCGCGCTCGTCGAGGCCATGGAAGGCACCGAGGGCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCtttttt                                                                                                                                           <  2:5433398/150‑1 (MQ=255)
            aGACTGGGGCGCGCTCGTCGAGGCCATGGAAGGCACCGAGGGCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCTCGAGGCGGGcaaggtg                                                                                                                                           <  1:96130/150‑7 (MQ=255)
             gACTGGGGCGCGCTCGTCGAGGCCATGGAAGGCACCGAGGGCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCttttttt                                                                                                                                          <  2:3367231/150‑1 (MQ=255)
             gACTGGGGCGCGATCGTCGAGGCCATGGAAGGCACCGAGGGCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCttttttt                                                                                                                                          <  1:3210468/150‑1 (MQ=255)
                      gcgcTCGTCGAGGCCATGGAAGGCACCGAGGGCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAg                                                                                                                                 >  2:3375977/1‑150 (MQ=255)
                      gcgcTCGTCGAGGCCATGGAAGGCACCGAGGGCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAg                                                                                                                                 >  1:4393812/1‑150 (MQ=255)
                       cgcTCGTCGAGGCCATGGAAGGCACCGAGGGCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGa                                                                                                                                >  2:5349516/1‑150 (MQ=255)
                          tcgtcgAGGCCATGGAAGGCACCGAGGGCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGaaaa                                                                                                                             <  1:3685600/150‑1 (MQ=255)
                          tcgtcgAGGCCATGGAAGGCACCGAGGGCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGaaaa                                                                                                                             <  1:3567330/150‑1 (MQ=255)
                                  gCCATGGAAGGCACCGAGGGCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTa                                                                                                                     >  2:3611142/1‑150 (MQ=255)
                                   ccATGGAAGGCACCGAGGGCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTAt                                                                                                                    <  1:3782055/150‑1 (MQ=255)
                                    cATGGAAGGCACCGAGGGCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATg                                                                                                                   <  2:205000/150‑1 (MQ=255)
                                    cATGGAAGGCACCGAGGGCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATg                                                                                                                   <  1:1011903/150‑1 (MQ=255)
                                     aTGGAAGGCACCGAGGGCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATgg                                                                                                                  <  1:5828806/150‑1 (MQ=255)
                                     aTGGAAGGCACCGAGGGCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATgg                                                                                                                  <  1:5245632/150‑1 (MQ=255)
                                            gCACCGAGGGCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTaccac                                                                                                           >  2:4087536/1‑150 (MQ=255)
                                             cACCGAGGGCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTaccacc                                                                                                          <  2:4371519/150‑1 (MQ=255)
                                              aCCGAGGGCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTAccacca                                                                                                         >  2:4807284/1‑150 (MQ=255)
                                               ccGAGGGCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCAt                                                                                                        <  2:5519446/150‑1 (MQ=255)
                                                cGAGGGCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATa                                                                                                       >  1:4288904/1‑150 (MQ=255)
                                                 gAGGGCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATaa                                                                                                      >  2:4195787/1‑150 (MQ=255)
                                                   gggCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACa                                                                                                    <  1:6094294/150‑1 (MQ=255)
                                                      cGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACAttt                                                                                                 <  1:222803/150‑1 (MQ=255)
                                                        tAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTAt                                                                                               <  1:3962130/150‑1 (MQ=255)
                                                          aGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGa                                                                                             <  2:6753322/150‑1 (MQ=255)
                                                          aGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGa                                                                                             <  1:1511396/150‑1 (MQ=255)
                                                          aGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGa                                                                                             <  2:3489622/150‑1 (MQ=255)
                                                           gCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGATCGTGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAg                                                                                            <  1:5590270/150‑1 (MQ=255)
                                                            ccTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGc                                                                                           >  1:1600242/1‑150 (MQ=255)
                                                            ccTGGCCATCGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCGAGAAGACTGATCTCGACAAACCCCCCGTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGc                                                                                           <  2:1526448/150‑1 (MQ=255)
                                                              tGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCaa                                                                                         >  1:444411/1‑150 (MQ=255)
                                                               ggCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCaaa                                                                                        <  2:4416199/150‑1 (MQ=255)
                                                                gCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAAc                                                                                       <  1:2584759/150‑1 (MQ=255)
                                                                     tGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCAtt                                                                                  <  1:656565/150‑1 (MQ=255)
                                                                     tGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCAtt                                                                                  <  2:6556937/150‑1 (MQ=255)
                                                                     tGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCAtt                                                                                  <  2:4298280/150‑1 (MQ=255)
                                                                        aGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGac                                                                               <  2:1730864/150‑1 (MQ=255)
                                                                          atgatgGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGacac                                                                             >  1:692382/1‑150 (MQ=255)
                                                                             atgGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACaaa                                                                          <  2:1630083/150‑1 (MQ=255)
                                                                              tgGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACaaaa                                                                         <  1:6409507/150‑1 (MQ=255)
                                                                              tgGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACaaaa                                                                         >  1:2083986/1‑150 (MQ=255)
                                                                               ggcggcGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAat                                                                        >  2:5196005/1‑150 (MQ=255)
                                                                                 cggcgTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAatat                                                                      <  2:1865553/150‑1 (MQ=255)
                                                                                   gcgTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGa                                                                    <  2:712810/150‑1 (MQ=255)
                                                                                   gcgTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGa                                                                    <  2:1259178/150‑1 (MQ=255)
                                                                                   gcgTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGa                                                                    <  2:216874/150‑1 (MQ=255)
                                                                                   gcgTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGa                                                                    <  2:2944612/150‑1 (MQ=255)
                                                                                   gcgTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGa                                                                    <  2:3732232/150‑1 (MQ=255)
                                                                                   gcgTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGa                                                                    <  1:829673/150‑1 (MQ=255)
                                                                                   gcgTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGa                                                                    <  1:4313702/150‑1 (MQ=255)
                                                                                   gcgTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGa                                                                    <  1:545828/150‑1 (MQ=255)
                                                                                         gtgtGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACATTGCCCGCATCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGTGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATTATcatac                                                              >  2:6364808/1‑150 (MQ=255)
                                                                                            tgACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGcatactca                                                           >  1:2940014/1‑150 (MQ=255)
                                                                                             gACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGcatactcat                                                          >  2:5892689/1‑150 (MQ=255)
                                                                                                 tGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCttttttt                                                                                                                                          >  1:5332135/1‑66 (MQ=255)
                                                                                                 tGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGCatactcatactc                                                      >  1:374414/1‑150 (MQ=255)
                                                                                                   ggAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGCATACTCATACTCGc                                                    <  2:2768260/150‑1 (MQ=255)
                                                                                                   ggAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGCATACTCATACTCGc                                                    <  2:5987525/150‑1 (MQ=255)
                                                                                                      aaaGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGCATACTCATACTCGCTcc                                                 <  1:1337768/150‑1 (MQ=255)
                                                                                                      aaaGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGCATACTCATACTCGCTcc                                                 >  1:6010588/1‑150 (MQ=255)
                                                                                                            gCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCATTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCAGTGACACAAAATATGAGCATACTCATACTCGCTCCATCAcc                                           <  1:1976421/150‑1 (MQ=255)
                                                                                                                gaagaCTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGCATACTCATACTCGCTCCATCACCGGAt                                       <  1:2803924/150‑1 (MQ=255)
                                                                                                                 aagaCTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGCATACTCATACTCGCTCCATCACCGGATa                                      <  1:1192670/150‑1 (MQ=255)
                                                                                                                 aagaCTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGCATACTCATACTCGCTCCATCACCGGATa                                      <  2:464763/150‑1 (MQ=255)
                                                                                                                 aagaCTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGCATACTCATACTCGCTCCATCACCGGATa                                      <  1:412681/150‑1 (MQ=255)
                                                                                                                 aagaCTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGCATACTCATACTCGCTCCATCACCGGATa                                      <  1:463305/150‑1 (MQ=255)
                                                                                                                  agaCTGATTTCGACAAACCCACCTTGCCAGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTTAATTATGGGTACCACCATAACATTTATGAGCAAAACCATTGAAACAAAATATTAGCATACTCATACTCTCTCCATCACCGGATAg                                     >  1:4734674/1‑150 (MQ=255)
                                                                                                                   gaCTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGa                                                                                             >  1:3308565/1‑93 (MQ=255)
                                                                                                                      tGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGCATACTCATACTCGCTCCATCACCGGATAGacac                                 <  1:6774464/150‑1 (MQ=255)
                                                                                                                        aTTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGCATACTCATACTCGCTCCATCACCGGATAGACACCg                               >  2:6084675/1‑150 (MQ=255)
                                                                                                                         tttCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGCATACTCATACTCGCTCCATCACCGGATAGACACCgg                              >  2:2860143/1‑150 (MQ=255)
                                                                                                                          ttCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGCATACTCATACTCGCTCCATCACCGGATAGACACCGGc                             <  1:4796106/150‑1 (MQ=255)
                                                                                                                          ttCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGCATACTCATACTCGCTCCATCACCGGATAGACACCGGc                             <  1:3435983/150‑1 (MQ=255)
                                                                                                                            ctaCAAACCCACCTTGCCCGCCTCGAGCGGGATTTTTTTGCGCCATAGAAAAGTGAATCATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGCATACTCATACTCGCTCCATCACCGGATAGACACCGGCCa                           <  1:5644564/148‑1 (MQ=255)
                                                                                                                             gACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGCATACTCATACTCGCTCCATCACCGGATAGACACCGGCCAg                          >  1:1133800/1‑150 (MQ=255)
                                                                                                                             gACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGCATACTCATACTCGCTCCATCACCGGATAGACACCGGCCAg                          >  2:2440907/1‑150 (MQ=255)
                                                                                                                              aCAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGCATACTCATACTCGCTCCATCACCGGATAGACACCGGCCAGc                         >  1:5387089/1‑150 (MQ=255)
                                                                                                                                    ccACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGCATACTCATACTCGCTCCATCACCGGATAGACACCGGCCAGcaacaca                   >  2:1354858/1‑148 (MQ=255)
                                                                                                                                    ccACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGCATACTCATACTCGCTCCATCACCGGATAGACACCGGCCAGCAACaaa                   <  1:3097763/150‑1 (MQ=255)
                                                                                                                                    ccACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGCATACTCATACTCGCTCCATCACCGGATAGACACCGGCCAGCAACaaa                   >  2:2588980/1‑150 (MQ=255)
                                                                                                                                    ccACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGCATACTCATACTCGCTCCATCACCGGATAGACACCGGCCAGCAACaaa                   >  1:4963048/1‑150 (MQ=255)
                                                                                                                                    ccACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGCATACTCATACTCGCTCCATCACCGGATAGACACCGGCCAGCAACaaa                   <  2:1936448/150‑1 (MQ=255)
                                                                                                                                     cACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGCATACTCATACTCGCTCCATCACCGGATAGACACCGGCCAGCAACACAg                  <  1:2912457/150‑1 (MQ=255)
                                                                                                                                      aCCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGCATACTCATACTCGCTCCATCACCGGATAGACACCGGCCAGCAACAAAgg                 <  1:5102126/150‑1 (MQ=255)
                                                                                                                                        cTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGCATACTCATACTCGCTCCATCACCGGATAGACACCGGCCAGCAACAAAGGCa               <  2:3157807/150‑1 (MQ=255)
                                                                                                                                        cTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGCATACTCATACTCGATCCATCACCGGATAGACACCGGCCAGCAACAAAGGCa               >  2:6231687/1‑150 (MQ=255)
                                                                                                                                        aaagCCCGCTCGAGGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGCATACTCATACTCGCTCCATCACCGGATAGACACCGGCCAGCAACAAAGGCa               >  2:5165583/4‑150 (MQ=255)
                                                                                                                                         aagCCCGCTCGAGGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGCATACTCATACTCGCTCCATCACCGGATAGACACCGGCCAGCAACAAAGGCAg              >  1:3367033/3‑150 (MQ=255)
                                                                                                                                             ccGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGCATACTCATACTCGCTCCATCACCGGATAGACACCGGCCAGCAACaaa                   >  1:5218972/1‑141 (MQ=255)
                                                                                                                                                 cTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGCATACTCATACTCGCTCCATCACCGGATAGACACCGGCCAGCAACAAAGGCAGCGAAAGTc      <  2:5006131/150‑1 (MQ=255)
                                                                                                                                                     aGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGCATACTCATACTCGCTCCATCACCGGATAGACACCGGCCAGCGACAAAGGCAGCGAAAGTCGCCa  <  1:4195787/150‑1 (MQ=255)
                                                                                                                                                     aGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGCATACTCATACTCGCTCCATCACCGGATAGACACCGGCCAGCAACAAAGGCAGCGAAAGTCGCCa  <  1:2037289/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
GGAATCCCAGGAAGACTGGGGCGCGCTCGTCGAGGCCATGGAAGGCACCGAGGGCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGCATACTCATACTCGCTCCATCACCGGATAGACACCGGCCAGCAACAAAGGCAGCGAAAGTCGCCA  >  NC_002947/4420931‑4421229

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