Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,303,020 G→C 5.1% G82A (GGT→GCT)  pvdQ → acyl‑homoserine lactone acylase PvdQ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,303,0200GC 5.1% 1060.1 / 33.1 312G82A (GGT→GCT) pvdQacyl‑homoserine lactone acylase PvdQ
Reads supporting (aligned to +/- strand):  ref base G (155/141);  new base C (5/11);  total (160/152)
Fisher's exact test for biased strand distribution p-value = 1.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.94e-01

Only 100 of 340 total aligned reads displayed.

GCCGACGCCAGTGCCCAGATCCGCTATACCCGTTACGGCGTACCGCATATCGTGGCCAAGGATGAGCGTGGCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGCGACCTGTTCTTCACCTGGCTCAATAGCCCAGCGGCGGTCGATGCCTTCCTGCAGGCGCAACCCGCGTCGGTGCAGGCATTGCTGG  >  NC_002947/3302872‑3303169
                                                                                                                                                    |                                                                                                                                                     
  cGACGCCAGTGCCCAGATCCGCTATACCCGTTACGGCGTACCGCATATCGTGGCCAAGGATGAGCGTGGCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCgtgagcggtga                                                                                                                                                    <  2:1258595/150‑1 (MQ=255)
      gCCAGTGCCCAGATCCGCTATACCCGTTACGGCGTACCGCATATCGTGGCCAAGGATGAGCGTGGCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAgcgc                                                                                                                                                >  1:6870988/1‑150 (MQ=255)
      gCCAGTGCCCAGATCCGCTATACCCGTTACGGCGTACCGCATATCGTGGCCAAGGATGAGCGTGGCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAgcgc                                                                                                                                                >  2:6068660/1‑150 (MQ=255)
       ccAGTGCCCAGATCCGCTATACCCGTTACGGCGTACCGCATATCGTGGCCAAGGATGAGCGTGGCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCt                                                                                                                                               <  2:6527090/150‑1 (MQ=255)
         aGTGCCCAGATCCGCTATCCCCGTTAGGGCGTACCGCATATCGTGGCCTAGGATGAGCGTGGCCTGGGGTATGGCGTCGGCTATGCCTACGCCCAGGACTACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGCTCACCGCTCa                                                                                                                                             <  2:256534/150‑1 (MQ=255)
         aGTGCCCAGATCCGCTATACCCGTTACGGCGTACCGCATATCGTGGCCAAGGATGAGCGTGGCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGCTCACCGCTCa                                                                                                                                             <  1:665863/150‑1 (MQ=255)
          gTGCCCAGATCCGCTATACCCGTTACGGCGTACCGCATATCGTGGCCAAGGATGAGCGTGGCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCAc                                                                                                                                            >  1:141125/1‑150 (MQ=255)
             cccAGATCCGCTATACCCGTTACGGCGTACCGCATATCGTGGCCAAGGATGAGCGTGGCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGCTCACCGCTCACggt                                                                                                                                         <  2:5319844/150‑3 (MQ=255)
              ccAGATCCGCTATACCCGTTACGGCGTACCGCATATCGTGGCCAAGGATGAGCGTGGCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTa                                                                                                                                        >  1:4459151/1‑150 (MQ=255)
              ccAGATCCGCTATACCCGTTACGGCGTACCGCATATCGTGGCCAAGGATGAGCGTGGCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTa                                                                                                                                        >  2:4865999/1‑150 (MQ=255)
              ccAGATCCGCTATACCCGTTACGGCGTACCGCATATCGTGGCCAAGGATGAGCGTGGCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTa                                                                                                                                        >  1:5890945/1‑150 (MQ=255)
              ccAGATCCGCTATACCCGTTACGGCGTACCGCATATCGTGGCCAAGGATGAGCGTGGCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTa                                                                                                                                        >  1:1599432/1‑150 (MQ=255)
              ccAGATCCGCTATACCCGTTACGGCGTACCGCATATCGTGGCCAAGGATGAGCGTGGCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTa                                                                                                                                        >  1:15024/1‑150 (MQ=255)
                   tCCGCTATACCCGTTACGGCGTACCGCATATCGTGGCCAAGGATGAGCGTGGCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCg                                                                                                                                   <  1:4773634/150‑1 (MQ=255)
                     cGCTATACCCGTTACGGCGTACCGCATATCGTGGCCAAGGATGAGCGTGGCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGgc                                                                                                                                 <  2:1062576/150‑1 (MQ=255)
                     cGCTATACCCGTTACGGCGTACCGCATATCGTGGCCAAGGATGAGCGTGGCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGgc                                                                                                                                 <  2:2272229/150‑1 (MQ=255)
                      gCTATACCCGTTACGGCGTACCGCATATCGTGGCCAAGGATGAGCGTGGCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGgcg                                                                                                                                >  2:2433163/1‑150 (MQ=255)
                         ataCCCGTTACGGCGTACCGCATATCGTGGCCAAGGATGAGCGTGGCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCgtgagcggt                                                                                                                                                      >  2:4743732/1‑125 (MQ=255)
                           aCCCGTTACGGCGTACCGCATATCGTGGCCAAGGATGAGCGTGGCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAg                                                                                                                           >  2:6847655/1‑150 (MQ=255)
                           aCCCGTTACGGCGTACCGCATATCGTGGCCAAGGATGAGCGTGGCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAg                                                                                                                           <  1:6188663/150‑1 (MQ=255)
                            cccGTTACGGCGTACCGCATATCGTGGCCAAGGATGAGCGTGGCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAgg                                                                                                                          <  2:3632099/150‑1 (MQ=255)
                             ccGTTACGGCGTACCGCATATCGTGGCCAAGGATGAGCGTGGCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAggg                                                                                                                         >  1:4935753/1‑150 (MQ=255)
                              cGTTACGGCGTACCGCATATCGTGGCCAAGGATGAGCGTGGCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGc                                                                                                                        <  2:3734116/150‑1 (MQ=255)
                              cGTTACGGCGTACCGCATATCGTGGCCAAGGATGAGCGTGGCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGc                                                                                                                        >  2:3892269/1‑150 (MQ=255)
                               gTTACGGCGTACCGCATATCGTGGCCAAGGATGAGCGTGGCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGTTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGcc                                                                                                                       <  2:1676442/150‑1 (MQ=255)
                                ttACGGCGTACCGCATATCGTGGCCAAGGATGAGCGTGGCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCa                                                                                                                      <  1:3710010/150‑1 (MQ=255)
                                   cGGCGTACCGCATATCGTGGCCAAGGATGAGCGTGGCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGAc                                                                                                                   >  1:1878672/1‑150 (MQ=255)
                                      cGTACCGCATATCGTGGCCAAGGATGAGCGTGGCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCt                                                                                                                >  1:5961785/1‑150 (MQ=255)
                                       gTACCGCATATCGTGGCCAAGGATGAGCGTGGCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTg                                                                                                               <  2:5551587/150‑1 (MQ=255)
                                         aCCGCATATCGTGGCCAAGGATGAGCGTGGCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGa                                                                                                             <  1:5346136/150‑1 (MQ=255)
                                         aCCGCATATCGTGGCCAAGGATGAGCGTGGCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGa                                                                                                             <  1:6466909/150‑1 (MQ=255)
                                         aCCGCATATCGTGGCCAAGGATGAGCGTGGCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGa                                                                                                             <  1:4395801/150‑1 (MQ=255)
                                         aCCGCATATCGTGGCCAAGGATGAGCGTGGCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGa                                                                                                             <  2:5261367/150‑1 (MQ=255)
                                         aCCGCATATCGTGGCCAAGGATGAGCGTGGCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGa                                                                                                             <  1:3916193/150‑1 (MQ=255)
                                           cGCATATCGTGGCCAAGGATGAGCGTGGCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGagc                                                                                                           >  2:2314302/1‑150 (MQ=255)
                                            gCATATCGTGGCCAAGGATGAGCGTGGCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGagca                                                                                                          >  1:6527329/1‑150 (MQ=255)
                                             cATATCGTGGCCAAGGATGAGCGTGGCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGagcag                                                                                                         >  2:5946248/1‑150 (MQ=255)
                                              atatCGTGGCCAAGGATGAGCGTGGCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGagcagc                                                                                                        <  1:5562418/150‑1 (MQ=255)
                                              atatCGTGGCCAAGGATGAGCGTGGCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGagcagc                                                                                                        <  1:2649382/150‑1 (MQ=255)
                                              atatCGTGGCCAAGGATGAGCGTGGCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGAGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGagcagc                                                                                                        <  1:876541/150‑1 (MQ=255)
                                                atCGTGGCCAAGGATGAGCGTGGCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGt                                                                                                      <  2:1811250/150‑1 (MQ=255)
                                                atCGTGGCCAAGGATGAGCGTGGCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGt                                                                                                      <  1:4582071/150‑1 (MQ=255)
                                                 tCGTGGCCAAGGATGAGCGTGGCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACTCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTg                                                                                                     <  1:5088415/150‑1 (MQ=255)
                                                 tCGTGGCCAAGGATGAGCGTGGCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTg                                                                                                     <  1:3280794/150‑1 (MQ=255)
                                                 tCGTGGCCAAGGATGAGCGTGGCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTg                                                                                                     >  1:6232723/1‑150 (MQ=255)
                                                   gTGGCCAAGGATGAGCGTGGCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGAc                                                                                                   >  2:5732858/1‑150 (MQ=255)
                                                             aTGAGCGTGGCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCa                                                                                         <  2:2231042/150‑1 (MQ=255)
                                                                aGCGTGGCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGCg                                                                                      <  2:4439203/150‑1 (MQ=255)
                                                                   gTGGCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGCGAcc                                                                                   <  2:1595326/150‑1 (MQ=255)
                                                                      gCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGGGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGCGACCTGt                                                                                >  2:5435534/1‑150 (MQ=255)
                                                                        cTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGCGACCTGttc                                                                              >  1:77017/1‑150 (MQ=255)
                                                                        cTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGCGACCTGttc                                                                              <  2:4940434/150‑1 (MQ=255)
                                                                         tGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGCGACCTGttct                                                                             >  2:6356611/1‑150 (MQ=255)
                                                                             cTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGCGACCTGTTCTTCAc                                                                         <  2:195923/150‑1 (MQ=255)
                                                                             cTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGCGACCTGTTCTTCAc                                                                         <  2:4041252/150‑1 (MQ=255)
                                                                             cTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGCGACCTGTTCTTCAc                                                                         <  2:77017/150‑1 (MQ=255)
                                                                              tATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGCGACCTGTTCTTCAcc                                                                        <  2:6125969/150‑1 (MQ=255)
                                                                                 ggCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGCGACCTGTTCTTCACCTgg                                                                     <  2:3371228/150‑1 (MQ=255)
                                                                                     tCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGCGACCTGTTCTTCACCTGGCTCa                                                                 >  2:1799893/1‑150 (MQ=255)
                                                                                       ggCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGCGACCTGTTCTTCACCTGGCTCAAt                                                               <  2:46072/150‑1 (MQ=255)
                                                                                             gCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGCGACCTGTTCTTCACCTGGCTCAATAGCCCa                                                         <  2:2083451/150‑1 (MQ=255)
                                                                                               cTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGCGACCTGTTCTTCACCTGGCTCAATAGCCCAgc                                                       >  2:873878/1‑150 (MQ=255)
                                                                                               cTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGCGACCTGTTCTTCACCTGGCTCAATAGCCCAgc                                                       >  2:2130218/1‑150 (MQ=255)
                                                                                                tACGCCCAGGACAACCAATGCCTGATCGCCAATGAAGTGCTGAACGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCTTGACAACCTGGCCAGCGACCTGTTCTTCAACTGGCTCAATAGCCCAgcg                                                      >  1:2260122/1‑150 (MQ=255)
                                                                                                 aCGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGCGACCTGTTCTTCACCTGGCTCAATAGCCCAgcgg                                                     >  2:4267892/1‑150 (MQ=255)
                                                                                                 aCGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGCGACCTGTTCTTCACCTGGCTCAATAGCCCAgcgg                                                     <  1:468567/150‑1 (MQ=255)
                                                                                                  cGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGCGACCTGTTCTTCACCTGGCTCAATAGCCCAgcggc                                                    <  1:2415725/150‑1 (MQ=255)
                                                                                                     ccAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGCGACCTGTTCTTCACCTGGCTCAATAGCCCAGCGGCGGt                                                 >  1:3020464/1‑150 (MQ=255)
                                                                                                     ccAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGCGACCTGTTCTTCACCTGGCTCAATAGCCCAGCGGCGGt                                                 >  2:6838504/1‑150 (MQ=255)
                                                                                                      cAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCCGCGTGACAACCTGGCCAGCGACCTGTTCTTCACCTGGCTCAATAGCCCAGCGGCGGTc                                                >  2:3415066/1‑150 (MQ=255)
                                                                                                      cAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGCGACCTGGTCTTCACCTGGCTCAATAGCCCAGCGGCGGTc                                                >  2:4630223/1‑150 (MQ=255)
                                                                                                          aCAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGCGACCTGTTCTTCACCTGGCTCAATAGCCCAGCGGCGGTCGATg                                            <  1:521966/150‑1 (MQ=255)
                                                                                                            aaCCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGCGACCTGTTCTTCACCTGGCACAATAGCCCAGCGGCGGTCGATGcc                                          >  1:4565487/1‑150 (MQ=255)
                                                                                                              ccTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTCTTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGCGACCTGTTCTTCACCTGGCTCAATAGCCCAGCGGCGGTCGATGCCtt                                        >  2:6181363/1‑150 (MQ=255)
                                                                                                               cTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACACTGGAGCAGCGTGACAACCTGGCCAGCGACCTGTTCTTCACCTGGCTCAATAGCCCAGCGGCGGTCGATGCCTTc                                       >  1:6741584/1‑150 (MQ=255)
                                                                                                                 aTGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGCGACCTGTTCTTCACCTGGCTCAATAGCCCAGCGGCGGTCGATGCCTTCCt                                     <  1:4820530/150‑1 (MQ=255)
                                                                                                                    ccTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGCGACCTGTTCTTCACCTGGCTCAATAGCCCAGCGGCGGTCGATGCCTTCCTGCa                                  <  1:965516/150‑1 (MQ=255)
                                                                                                                     cTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGCGACCTGTTCTTCACCTGGCTCAATAGCCCAGCGGCGGTCGATGCCTTCCTGCAg                                 <  2:4510098/150‑1 (MQ=255)
                                                                                                                     cTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGCGACCTGTTCTTCACCTGGCTCAATAGCCCAGCGGCGGTCGATGCCTTCCTGCAg                                 >  1:6240158/1‑150 (MQ=255)
                                                                                                                      tGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGCGACCTGTTCTTCACCTGGCTCAATAGCCCAGCGGCGGTCGATGCCTTCCTGCAgg                                <  1:499743/150‑1 (MQ=255)
                                                                                                                      tGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGCGACCTGTTCTTCACCTGGCTCAATAGCCCAGCGGCGGTCGATGCCTTCCTGCAgg                                >  1:1981192/1‑150 (MQ=255)
                                                                                                                      tGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGCGACCTGTTCTTCACCTGGCTCAATAGCCCAGCGGCGGTCGATGCCTTCCTGCAgg                                >  1:4058379/1‑150 (MQ=255)
                                                                                                                      tGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGCGACCTGTTCTTCACCTGGCTCAATAGCCCAGCGGCGGTCGATGCCTTCCTGCAgg                                >  1:3873162/1‑150 (MQ=255)
                                                                                                                           gCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGCGACCTGTTCTTCACCTGGCTCAATAGCCCAGCGGCGGTCGATGCCTTCCTGCAGGCGCaa                           >  1:3363105/1‑150 (MQ=255)
                                                                                                                             cAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGCGACCTGTTCTTCACCTGGCTCAATAGCCCAGCGGCGGTCGATGCCTTCCTGCAGGCGCAAcc                         >  1:3667689/1‑150 (MQ=255)
                                                                                                                             cAATGAAGGGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGCGACCTGTTCTTCACCTGGCTCAATAGCCCAGCGGCGGTCGATGCCTTCCTGCAGGCGCAAcc                         >  1:2838266/1‑150 (MQ=255)
                                                                                                                              aaTGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGCGACCTGTTCTTCACCTGGCTCAATAGCCCAGCGGCGGTCGATGCCTTCCTGCAGGCGCAAccc                        >  2:5034722/1‑150 (MQ=255)
                                                                                                                                tGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGCGACCTGTTCTTCACCTGGCTCAATAGCCCAGCGGCGGTCGATGCCTTCCTGCAGGCGCAACCcgc                      <  1:1930407/150‑1 (MQ=255)
                                                                                                                                  aaGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGCGACCTGTTCTTCACCTGGCTCAATAGCCCAGCGGCGGTCGATGCCTTCCTGCAGGCGCAACCCGCGt                    >  1:104526/1‑150 (MQ=255)
                                                                                                                                   aGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGCGACCTGTTCTTCACCTGGCTCGATAGCCCAGCGGCGGTCGATGCCTTCCTGCAGGCGCAACCCGCGTc                   <  2:2864736/150‑1 (MQ=255)
                                                                                                                                     tGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGCGACCTGTTCTTCACCTGGCTCAATAGCCCAGCGGCGGTCGATGCCTTCCTGCAGGCGCAACCCGCGTCgg                 <  1:2225096/150‑1 (MQ=255)
                                                                                                                                          aCCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGCGACCTGTTCTTCACCTGGCTCAATAGCCCAGCGGCGGTCGATGCCTTCCTGCAGGCGCAACCCGCGTCGGTGCAg            <  1:313163/150‑1 (MQ=255)
                                                                                                                                          aCCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGCGACCTGTTCTTCACCTGGCTCAATAGCCCAGCGGCGGTCGATGCCTTCCTGCAGGCGCAACCCGCGTCGGTGCAg            >  1:1339600/1‑150 (MQ=255)
                                                                                                                                          aCCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGCGACCTGTTCTTCACCTGGCTCAATAGCCCAGCGGCGGTCGATGCCTTCCTGCAGGCGCAACCCGCGTCGGTGCAg            <  2:141120/150‑1 (MQ=255)
                                                                                                                                           gcGTGAGCGCTCACCGCTCACGGTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGCGACCTGTTCTTCACCTGGCTCAATAGCCCAGCGGCGGTCGATGCCTTCCTGCAGGCGCAACCCGCGTCGGTGCAgg           >  2:6088581/2‑150 (MQ=255)
                                                                                                                                            cGTGAGCGCTCACCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGCGACCTGTTCTTCACCTGGCTCAATAGCCCAGCGGCGGTCGATGCCTTCCTGCAGGCGCAATCCGCGTCGGTGCAGGc          >  1:818806/1‑150 (MQ=255)
                                                                                                                                             gtgagcggtgagcgCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGCGACCTGTTCTTCACCTGGCTCAATAGCCCAGCGGCGGTCGATGCCTTCCTGCAGGCGCAACCCGCGTCGGTGCAGGCa         <  2:6854302/150‑1 (MQ=255)
                                                                                                                                              tgagcggtgagcgCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGCGACCTGTTCTTCACCTGGCTCAATAGCCCAGCGGCGGTCGATGCCTTCCTGCAGGCGCAACCCGCGTCGGTGCAggaag        <  2:3018488/150‑4 (MQ=255)
                                                                                                                                                agcggtgagcgCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGCGACCTGTTCTTCACCTGGCTCAATAGCCCAGCGGCGGTCGATGCCTTCCTGCAGGCGCAACCCGCGTCGGTGCAGGCATTg      >  1:4121363/1‑150 (MQ=255)
                                                                                                                                                 gcggtgagcgCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGCGACCTGTTCTTCACCTGGCTCAATAGCCCAGCGGCGGTCGATGCCTTCCTGCAGGCGCAACCCGCGTCGGTGCAGGCATTGc     >  2:4119811/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GCCGACGCCAGTGCCCAGATCCGCTATACCCGTTACGGCGTACCGCATATCGTGGCCAAGGATGAGCGTGGCCTGGGCTATGGCGTCGGCTATGCCTACGCCCAGGACAACCTATGCCTGCTCGCCAATGAAGTGCTGACCGTGAGCGGTGAGCGCTCACGCTATTTCGGCGCCAAGGGCCAGACCCTGGAGCAGCGTGACAACCTGGCCAGCGACCTGTTCTTCACCTGGCTCAATAGCCCAGCGGCGGTCGATGCCTTCCTGCAGGCGCAACCCGCGTCGGTGCAGGCATTGCTGG  >  NC_002947/3302872‑3303169

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: