Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,857,389 T→A 8.2% Q312L (CAG→CTG)  PP_2505 ← GAF domain/GGDEF domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,857,3890TA 8.2% 1110.3 / 64.4 352Q312L (CAG→CTG) PP_2505GAF domain/GGDEF domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base T (175/147);  new base A (16/13);  total (191/161)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.29e-01

Only 100 of 353 total aligned reads displayed.

CCGAATGTGATGCGTACCCTCCCGGATCCGGAGATCTTAACTTCCTCAGGATAACGCTGGCACATAGACAACTTGGTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATTGACTGACGCAAGCGCTCCGCAATCTGCGCCAATTCTTGGCTATCTTTAACCCGTAGCAAGGCTACGAATTCC  >  NC_002947/2857241‑2857538
                                                                                                                                                    |                                                                                                                                                     
     tgtgATGCGTACCCTCCCGGATCCGGAGATCTTAACTTCCTCAGGATAACGCTGGCACATAGACAACTTGGTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGa                                                                                                                                                 <  1:270840/150‑1 (MQ=255)
        gATGCGTACCCTCCCGGATCCGGAGATCTTAACTTCCTCAGGATAACGCTGGCACATAGACAACTTGGTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCa                                                                                                                                              >  2:3271404/1‑150 (MQ=255)
          tGCGTACCCTCCCGGATCCGGAGATCTTAACTTCCTCAGGATAACGCTGGCACATAGACAACTTGGTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAgg                                                                                                                                            <  2:102198/150‑1 (MQ=255)
          tGCGTACCCTCCCGGATCCGGAGATCTTAACTTCCTCAGGATAACGCTGGCACATAGACAACTTGGTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAgg                                                                                                                                            <  1:5016324/150‑1 (MQ=255)
          tGCGTACCCTCCCGGATCCGGAGATCTTAACTTCCTCAGGATAACGCTGGCACATAGACAACTTGGTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCAGGGTGAGCAga                                                                                                                                            <  1:1466612/150‑2 (MQ=255)
          tGCGTACCCTCCCGGATCCGGAGATCTTAACTTCCTCAGGATAACGCTGGCACATAGACAACTTGGTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAACCTGCTCACCCAGGGTGAGCAga                                                                                                                                            <  1:3301144/150‑2 (MQ=255)
          tGCGTACCCTCCCGGATCCGGAGATCTTAACTTCCTCAGGATAACGCTGGCACATAGACAACTTGGTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAACCTGCTCACCCAGGGTGAGCAga                                                                                                                                            <  2:7130772/150‑2 (MQ=255)
          tGCGTACCCTCCCGGATCCGGAGATCTTAACGTCCTCAGGATAACGCTGGCACATAGACAACTTGGTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAgg                                                                                                                                            <  1:5065732/150‑1 (MQ=255)
          tGCGTACCCTCCCGGATCCGGAGATCTTAACGTCCTCAGGATAACGCTGGCACATAGACAACTTGGTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAgg                                                                                                                                            <  1:3767365/150‑1 (MQ=255)
           gCGTACCCTCCCGGATCCGGAGATCTTAACTTCCTCAGGATAACGCTGGCACATAGACAACTTGTTCAAGTCCACCATGTTTGGCCGCATATAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAggg                                                                                                                                           <  1:5157806/150‑1 (MQ=255)
           gCGTACCCTCCCGGATCCGGAGATCTTAACTTCCTCAGGATAACGCTGGCACATAGACAACTTGGTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAggg                                                                                                                                           <  1:5634233/150‑1 (MQ=255)
               aCCCTCCCGGATCCGGAGATCTTAACTTCCTCAGGATAACGCTGGCACATAGACAACTTGGTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCa                                                                                                                                       <  2:1230962/150‑1 (MQ=255)
                   tCCCGGATCCGGAGATCTTAACTTCCTCAGGATAACGCTGGCACATAGACAACTTGGTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCa                                                                                                                                   >  1:2535110/1‑150 (MQ=255)
                   tCCCGGATCCGGAGATCTTAACTTCCTCAGGATAACGCTGGCACATAGACAACTTGGTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCa                                                                                                                                   >  1:78859/1‑150 (MQ=255)
                   tCCCGGATCCGGAGATCTTAACTTCCTCAGGATAACGCTGGCACATAGACAACTTGGTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCa                                                                                                                                   >  1:5341902/1‑150 (MQ=255)
                      cGGATCCGGAGATCTTAACTTCCTCAGGATAACGCTGGCACATAGACAACTTGGTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATg                                                                                                                                >  2:4253858/1‑150 (MQ=255)
                      cGGATCCGGAGATCTTAACTTCCTCAGGATAACGCTGGCACATAGACAACTTGGTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATg                                                                                                                                >  2:4254193/1‑150 (MQ=255)
                      cGGATCCGGAGATCTTAACTTCCTCAGGATAACGCTGGCACATAGACAACTTGGTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATg                                                                                                                                <  2:5893372/150‑1 (MQ=255)
                       ggATCCGGAGATCTTAACTTCCTCAGGATAACGCTGGCACATAGACAACTTGGTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGCGCCCTGCTCACCCagggt                                                                                                                               <  1:5957992/150‑5 (MQ=255)
                        gATCCGGAGATCTTAACTTCCTCAGGATAACGCTGGCACAGAGACAACTTGGTCACGTCCTCCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCt                                                                                                                              <  2:5746709/150‑1 (MQ=255)
                              gagaTCTTAACTTCCTCAGGATAACGCTGGCACATAGACAACTTGGTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGt                                                                                                                        >  2:4504683/1‑150 (MQ=255)
                                   cTTAACTTCCTCAGGATAACGCTGGCACATAGACAACTTGGTCACGTCCACCATGTTTGGCCGCAGAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCt                                                                                                                   >  1:4109003/1‑150 (MQ=255)
                                    ttAACTTCCTCAGGATAACGCTGGCACATAGACAACTTGGTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCtt                                                                                                                  >  2:835428/1‑150 (MQ=255)
                                     tAACTTCCTCAGGATAACGCTGGCACATAGACAACTTGGTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCAGGGTGAGCAgattga                                                                                                                                        >  2:1163080/1‑122 (MQ=255)
                                     tAACTTCCTCAGGATAACGCTGGCACATAGACAACTTGGTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCAGGGTGAGCAgattga                                                                                                                                        <  1:1163080/127‑6 (MQ=255)
                                         ttCCTCAGGATAACGCTGGCACATAGACAACTTGGTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAg                                                                                                             >  2:4363593/1‑150 (MQ=255)
                                         ttCCTCAGGATAACGCTGGCACATAGACAACTTGGTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAg                                                                                                             >  1:6363653/1‑150 (MQ=255)
                                           ccTCAGGATAACGCTGGCACATAGACAACTTGGTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAggg                                                                                                           <  1:112934/150‑1 (MQ=255)
                                           ccTCAGGATAACGCTGGCACATAGACAACTTGGTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAggg                                                                                                           <  1:4759073/150‑1 (MQ=255)
                                            cTCAGGATATCGCTGGCACATAGACTACTTGGTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGgc                                                                                                          <  2:5453095/150‑1 (MQ=255)
                                            cTCAGGATAACGCTGGCACATAGACAACTTGGTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGgc                                                                                                          >  1:3902700/1‑150 (MQ=255)
                                            cTCAGGATAACGCTGGCACATAGACAACTTGGTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGgc                                                                                                          >  2:6839080/1‑150 (MQ=255)
                                            cTCAGGATAACGCTGGCACATAGACAACTTGGTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGgc                                                                                                          >  1:1814715/1‑150 (MQ=255)
                                             tCAGGATAACGCTGGCACATAGACAACTTGGTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGgcc                                                                                                         >  1:5547152/1‑150 (MQ=255)
                                             tCAGGATAACGCTGGCACATAGACAACTTGGTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGgcc                                                                                                         >  2:2904880/1‑150 (MQ=255)
                                                 gATAACGCTGGCACATAGACAACTTGGTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGccgccg                                                                                                     >  2:3379266/1‑150 (MQ=255)
                                                 gATAACGCTGGCACATAGACAACTTGGTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGccgccg                                                                                                     >  1:1958528/1‑150 (MQ=255)
                                                        cTGGCACATAGACAACTTGGTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACa                                                                                              <  1:5797008/150‑1 (MQ=255)
                                                          ggCACATAGACAACTTGGTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGa                                                                                            >  1:1376328/1‑150 (MQ=255)
                                                           gCACATAGACAACTTGGTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGaa                                                                                           <  2:4466435/150‑1 (MQ=255)
                                                               aTAGACAACTTGGTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCGGAGGGCCGCCGGGTAACAGAAACgg                                                                                       >  2:7010867/1‑150 (MQ=255)
                                                                tAGACAACTTGGTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCTCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGgc                                                                                      <  2:1508000/150‑1 (MQ=255)
                                                                   aCAACTTGGTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGgcgcg                                                                                   >  1:4815235/1‑150 (MQ=255)
                                                                       cTTGGTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGc                                                                               <  2:6463198/150‑1 (MQ=255)
                                                                       cTTGGTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGc                                                                               <  2:2402208/150‑1 (MQ=255)
                                                                          ggTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGa                                                                            <  2:3651340/150‑1 (MQ=255)
                                                                           gTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGAt                                                                           <  2:730894/150‑1 (MQ=255)
                                                                               cGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTCTTGCATCAGCCATCTGCACACCCTGGGTGCGCAGGGCGCACCCACTGCGCATGGTCATCTGCCGAGGGCCGCCGGGTACCAGAAACGGCGCGCTGCCGATtgac                                                                       >  2:6683964/1‑150 (MQ=255)
                                                                                gTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATtgact                                                                      <  2:5018644/150‑1 (MQ=255)
                                                                                    accaTGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGc                                                                                                                                               >  2:2586798/1‑73 (MQ=255)
                                                                                    accaTGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATTGACTGACg                                                                  >  1:893997/1‑150 (MQ=255)
                                                                                    accaTGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATTGACTGACg                                                                  >  2:4794181/1‑150 (MQ=255)
                                                                                    accaGGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATTGACTGACg                                                                  >  2:5594362/1‑150 (MQ=255)
                                                                                      caTGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATTGACTGACGCa                                                                >  2:1782055/1‑150 (MQ=255)
                                                                                      caTGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATTGACTGACGCa                                                                >  1:6188965/1‑150 (MQ=255)
                                                                                      caTGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCAGGGTGAGCAga                                                                                                                                            <  2:1603718/74‑2 (MQ=255)
                                                                                      caTGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCAGGGTGAGCAga                                                                                                                                            >  1:1604428/1‑73 (MQ=255)
                                                                                       aTGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGgcc                                                                                                         <  1:7050591/108‑1 (MQ=255)
                                                                                       aTGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGgcc                                                                                                         >  2:7050591/1‑108 (MQ=255)
                                                                                           ttGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATTGACTGACGCAAgcgc                                                           >  2:3453842/1‑150 (MQ=255)
                                                                                            tGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATTGACTGACGCAAGCGCt                                                          <  2:1646722/150‑1 (MQ=255)
                                                                                                cGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATTGACTGACGCAAGCGCTCCGc                                                      >  1:6875640/1‑150 (MQ=255)
                                                                                                 gCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATTGACTGACGCAAGCGCTCCGCa                                                     >  1:6963091/1‑150 (MQ=255)
                                                                                                   aTAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATTGACTGACGCAAGCGCTCCGCAAt                                                   >  1:1561586/1‑150 (MQ=255)
                                                                                                   aTAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATTGACTGACGCAAGCGCTCCGCAAt                                                   <  1:3166344/150‑1 (MQ=255)
                                                                                                   aTAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATTGACTGACGCAAGCGCTCCGCAAt                                                   >  1:3669167/1‑150 (MQ=255)
                                                                                                      gagaCCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATTGACTGACGCAAGCGCTCCGCAATCTg                                                >  1:4492163/1‑150 (MQ=255)
                                                                                                       agaCCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATTGACTGACGCAAGCGCTCCGCAATCTgc                                               <  1:3505327/150‑1 (MQ=255)
                                                                                                         aCCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATTGACTGACGCAAGCGCTCCGCAATCTgcgc                                             <  1:1340232/150‑1 (MQ=255)
                                                                                                         aCCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATTGACTGACGCAAGCGCTCCGCAATCTgcgc                                             <  2:4725687/150‑1 (MQ=255)
                                                                                                                tCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATTGACTGACGCAAGCGCTCCGCAATCTGCGCCAATTCt                                      >  2:3578746/1‑150 (MQ=255)
                                                                                                                   gcgcgcTTTATTGCATCATCACTCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATTGACTGACGCAAGCGCTCCGCAATCTGCGCCAATTCTTgg                                   <  2:6212588/150‑1 (MQ=255)
                                                                                                                   gcgcgcTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCTTGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATTGACTGACGCAAGCGCTCCGCAATCTGCGCCAATTCTTgg                                   <  1:7196298/150‑1 (MQ=255)
                                                                                                                    cgcgcTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATTGACTGACGCAAGCGCTCCGCAATCTGCGCCAATTCTTGGc                                  <  1:2019265/150‑1 (MQ=255)
                                                                                                                       gcTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATTGACTGACGCAAGCGCTCCGCAATCTGCGCCAATTCTTGGCTAt                               <  1:1782055/150‑1 (MQ=255)
                                                                                                                       gcTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATTGACTGACGCAAGCGCTCCGCAATCTGCGCCAATTCTTGGCTAt                               >  1:6049749/1‑150 (MQ=255)
                                                                                                                         tttATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATTGACTGACGCAAGCGCTCCGCAATCTGCGCCAATTCTTGGCTATCt                             >  1:6806723/1‑150 (MQ=255)
                                                                                                                         tttATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATTGACTGACGCAAGCGCTCCGCAATCTGCGCCAATTCTTGGCTATCt                             <  2:134644/150‑1 (MQ=255)
                                                                                                                              tGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATTGACTGACGCAAGCGCTCCGCAATCTGCGCCAATTCTTGGCTATCTTTAAc                        >  2:2141738/1‑150 (MQ=255)
                                                                                                                               gCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACACAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATTGACTGACGCAAGCGCTCCGCAATCTGCGCAAATTCTTGGCTATCTTTAAcc                       >  2:4144794/1‑150 (MQ=255)
                                                                                                                                catcatCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATTGACTGAAGCAAGCGCTCCGCAATCTGCGCCAATTCtt                                     <  2:3875807/135‑1 (MQ=255)
                                                                                                                                 atcatcAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATTGACTGACGCAAGCGCTCCGCAATCTGCGCCAATTCTTGGCTATCTTTAACCCg                     >  1:243299/1‑150 (MQ=255)
                                                                                                                                 atcatcAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATTGACTGACGCAAGCGCTCCGCAATCTGCGCCAATTCTTGGCTATCTTTAACCCg                     >  1:6471950/1‑150 (MQ=255)
                                                                                                                                       gcccTGCTCACCCAGGGTGAGCAGAGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATTGACTGACGCAAGCGCTCCGCACTCTGCGCCAATTCTTGGCTATCTTTAACCCGTAGCaa               >  1:3203081/4‑150 (MQ=255)
                                                                                                                                        aTCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATTGACTGACGCAAGCGCTCCGCAATCTGCGCCAATTCTTGGCTATCTTTAACCCGTAGCAAg              <  1:2622691/150‑1 (MQ=255)
                                                                                                                                        aTCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATTGACTGACGCAAGCGCTCCGCAATCTGCGCCAATTCTTGGCTATCTTTAACCCGTAGCAAg              <  1:4071299/150‑1 (MQ=255)
                                                                                                                                         tCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATTGACTGACGCAAGCGCTCCGCAATCTGCGCCAATTCTTGGCTATCTTTAACCCGTAGCAAgg             >  1:1248232/1‑150 (MQ=255)
                                                                                                                                         ccTGCTCACCCAGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATTGACTGACGCAAGCGCTCCGCAATCTGCGCCAATTCTTGGCTATCTTTAACCCGTAGCAAgg             >  1:4405606/2‑150 (MQ=255)
                                                                                                                                         ccTGCTCACCCAGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATTGACTGACGCAAGCGCTCCGCAATCTGCGCCAATTCTTGGCTATCTTTAACCCGTAGCAAgg             >  2:6920230/2‑150 (MQ=255)
                                                                                                                                           tGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATTGACTGACGCAAGCGCTCCGCAATCTGCGCCAATTCTTGGCTATCTTTAACCCGTAGCAAGGCt           <  2:4667663/150‑1 (MQ=255)
                                                                                                                                           tGCTCACCCAGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATTGACTGACGCAAGCGCTCCGCACTCTGCGCCAATTCTTGGCTATCTTTAACCCGTAGCAAGGCt           >  2:6019129/1‑150 (MQ=255)
                                                                                                                                           tGCTCACCCAGGGTGAGCAGAGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATTGACTGACGCAAGCGCTCCGCAATCTGCGCCAATTCTTGGCTATCTTTAACCCGTAGCAAGGCt           >  1:672435/1‑150 (MQ=255)
                                                                                                                                              tCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATTGACTGACGCAAGCGCTCCGCAATCTGCGCCAATTCTTGGCTATCTTTAACCCGTAGCAAGGCTACg        >  2:2649075/1‑150 (MQ=255)
                                                                                                                                               cACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATTGACTGACGCAAGCGCTCCGCAATCTGCGCCAATTCTTGGCTATCTTTAACCCGTAGCAAGGCTACGa       >  2:581555/1‑150 (MQ=255)
                                                                                                                                                aCCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATTGACTGACGCAAGCGCTCCGCAATCTGCGCCAATTCTTGGCTATCTTTAACCCGTAGCAAGGCTACGaa      >  1:3608298/1‑150 (MQ=255)
                                                                                                                                                 cccTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATTGACTGACGCAAGCGCTCCGCAATCTGCGCCAATTCTTGGCTATCTTTAAc                        >  2:508649/1‑131 (MQ=255)
                                                                                                                                                  ccTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATTGACTGACGCAAGCGCTCCGCAATCTGCGCCAATTCTTGGCTATCTTTAACCCGTAGCAAGGCTACGAAtt    >  2:3607877/1‑150 (MQ=255)
                                                                                                                                                  ccTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATTGACTGACGCAAGCGCTCCGCAATCTGCGCCAATTCTTGGCTATCTTTAACCCGTAGCAAGGCTACGAAtt    >  2:2714158/1‑150 (MQ=255)
                                                                                                                                                   cTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATTGACTGACGCAAGCGCTCCGCAATCTGCGCCAATTCTTGGCTATCTTTAACCCGTAGCAAGGCTACGAATtc   >  2:3686809/1‑150 (MQ=255)
                                                                                                                                                   cTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATTGACTGACGCAAGCGCTCCGCAATCTGCGCCAATTCTTGGCTATCTTTAACCCGTAGCAAGGCTACGAATtc   >  1:301388/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CCGAATGTGATGCGTACCCTCCCGGATCCGGAGATCTTAACTTCCTCAGGATAACGCTGGCACATAGACAACTTGGTCACGTCCACCATGTTTGGCCGCATAGAGACCGCCGTCAGCGCGCTTTATTGCATCATCAATCTGCTCACCCTGGGTGAGCAGGGCGCACCCAATGCTCATGGTCATCTTCCGAGGGCCGCCGGGTAACAGAAACGGCGCGCTGCCGATTGACTGACGCAAGCGCTCCGCAATCTGCGCCAATTCTTGGCTATCTTTAACCCGTAGCAAGGCTACGAATTCC  >  NC_002947/2857241‑2857538

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: