Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,300,198:1 +C 5.5% intergenic (+53/+29) pvdQ → / ← lpxH acyl‑homoserine lactone acylase PvdQ/UDP‑2,3‑diacylglucosamine hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,300,1981.C 5.5% 375.0 / 14.7 91intergenic (+53/+29)pvdQ/lpxHacyl‑homoserine lactone acylase PvdQ/UDP‑2,3‑diacylglucosamine hydrolase
Reads supporting (aligned to +/- strand):  ref base . (44/42);  new base C (2/3);  total (46/45)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.50e-01

TTCCGTTTACCGAGGCGCAGATCAAGGCAGATCCTGAGTACCGGGAAGTGGTGATCAGTGAGCGGGACAAGGGGGCGGTAGTCAGCCAGCCCTGAGGCGTCAGCCTGCGCCGGCCCATTCGCGGGTAAACCCGCTCCTACAGGGGCCTGCGGGCTTGCCCGCGAGGCCCTGTAATGTCAGGAAAACTCGAACGGCGCCAACTGAAAGCCCTGTTCGTCAACCTGCAAGGCCCAACCGCGGCGGTCCCAGTCGCCCAGCACGATACGCCGTGCCGGCTGGCCATCGACCATCAGCTTG  >  NC_002947/3300051‑3300347
                                                                                                                                                    |                                                                                                                                                     
ttCCGTTTACCGAGGCGCAGATCAAGGCAGATCCTGAGTACCGGGAAGTGGTGATCAGTGAGCGGGACAAGGGGGCGGTAGTCAGCCAGCCCTGAGGCGTCAGCCTGCGCCGGCCCATTCGCGGGTAAACCCGCTCCTACAGGGGCCTGc                                                                                                                                                     <  1:408702/150‑1 (MQ=255)
 tCCGTTTACCGAGGCGCAGATCAAGGCAGATCCTGAGTACCGGGAAGTGGTGATCAGTGAGCGGGACAAGGGGGCGGTAGTCAGCCAGCCCTGAGGCGTCAGCCTGCGCCGGCCCATTCGCGGGTAAACCCGCTCCTACAGGGGCCTGCg                                                                                                                                                    <  2:86032/150‑1 (MQ=255)
 tCCGTTTACCGAGGCGCAGATCAAGGCAGATCCTGAGTACCGGGAAGTGGTGATCAGTGAGCGGGACAAGGGGGCGGTAGTCAGCCAGCCCTGAGGCGTCAGCCTGCGCCGGCCCATTCGCGGGTAAACCCGCTCCTACAGGGGCCTGCg                                                                                                                                                    <  2:1991051/150‑1 (MQ=255)
    gTTTACCGAGGCGCAGATCAAGGCAGATCCTGAGTACCGGGAAGTGGTGATCAGTGAGCGGGACAAGGGGGCGGTAGTCAGCCAGCCCTGAGGCGTCAGCCTGCGCCGGCCCATTCGCGGGTAAACCCGCTCCTACAGGGGCCTGCGGGc                                                                                                                                                 >  1:504320/1‑150 (MQ=255)
    gTTTACCGAGGCGCAGATCAAGGCAGATCCTGAGTACCGGGAAGTGGTGATCAGTGAGCGGGACAAGGGGGCGGTAGTCAGCCAGCCCTGAGGCGTCAGCCTGCGCCGGCCCATTCGCGGGTAAACCCGCTCCTACAGGGGCCTGCGGGc                                                                                                                                                 >  2:487440/1‑150 (MQ=255)
             ggCGCAGATCAAGGCAGATCCTGAGTACCGGGAAGTGGTGATCAGTGAGCGGGACAAGGGGGCGGTAGTCAGCCAGCCCTGAGGCGTCAGCCTGCGCCGGCCCATTCGCGGGTAAACCCGCTCCTACAGGGGCCTCGCGGGCAAGCCcgc                                                                                                                                         <  2:518363/150‑1 (MQ=255)
                gcAGATCAAGGCAGATCCTGAGTACCGGGAAGTGGTGATCAGTGAGCGGGGCAAGGGGGCGGTGGTCAGCCAGCCCTGAGGCGTCAGCCTGCGCCGGCCCATTCGCGGGTAAACCCGCTCCTACAGGGGCCTGCAGGATTGCCCGCGAgg                                                                                                                                     >  2:1325539/1‑150 (MQ=255)
                gcAGATCAAGGCAGATCCTGAGTACCGGGAAGTGGTGATCAGTGAGCGGGACAAGGGGGCGGTAGTCAGCCAGCCCTGAGGCGTCAGCCTGCGCCGGCCCATTCGCGGGTAAACCCGCTCCTACAGGGGCCTGCGGGCTTGCCCGCGAgg                                                                                                                                     >  2:1302608/1‑150 (MQ=255)
                 cAGATCAAGGCAGATCCTGAGTACCGGGAAGTGGTGATCAGTGAGCGGGACAAGGGGGCGGTAGTCAGCCAGCCCTGAGGCGTCAGCCTGCGCCGGCCCATTCGCGGGTAAACCCGCTCCTACAGGGGCCTGCGGGCTTGCCCGCGAGGc                                                                                                                                    <  1:102373/150‑1 (MQ=255)
                 cAGATCAAGGCAGATCCTGAGTACCGGGAAGTGGTGATCAGTGAGCGGGACAAGGGGGCGGTAGTCAGCCAGCCCTGAGGCGTCAGCCTGCGCCGGCCCATTCGCGGGTAAACCCGCTCCTACAGGGGCCTGCGGGCTTGCCCGCGAGGc                                                                                                                                    >  1:1710515/1‑150 (MQ=255)
                 cAGATCAAGGCAGATCCTGAGTACCGGGAAGTGGTGATCAGTGAGCGGGACAAGGGGGCGGTAGTCAGCCAGCCCTGAGGCGTCAGCCTGCGCCGACCCATTCGCGGGTAAACCCGCTCTTACCGGGGCCTGCGGGCTTGCCCGCGAGGc                                                                                                                                    <  1:1243669/150‑1 (MQ=255)
                    aTCAAGGCAGATCCTGAGTACCGGGAAGTGGTGATCAGTGAGCGGGACAAGGGGGCGGTAGTCAGCCAGCCCTGAGGCGTCAGCCTGCGCCGGCCCATTCGCGGGTAAACCCGCTCCTACAGGGGCCTGCGGGCTTGCCCGCGAGGCCCt                                                                                                                                 <  1:929514/150‑1 (MQ=255)
                    aTCAAGGCAGATCCTGAGTACCGGGAAGTGGTGATCAGTGAGCGGGACAAGGGGGCGGTAGTCAGCCAGCCCTAAGGCGTCAGCCTGCGCCGGCCCATTCGCGGGTAAACCCGCTCCTACAGGGGCCTGCGGGCTTGCCCGCGAGGCCCt                                                                                                                                 <  1:1578754/150‑1 (MQ=255)
                     tCAAGGCAGATCCTGAGTACCGGGAAGTGGTGATCAGTGAGCGGGACAAGGGGGCGGTAGTCAGCCAGCCCTGAGGCGTCAGCCTGCGCCGGCCCATTCGCGGGTAAACCCGCTCCTACAGGGGCCTGCGGGCTTGCCCGCGAGGCCCTg                                                                                                                                <  1:269103/150‑1 (MQ=255)
                     tCAAGGCAGATCCTGAGTACCGGGAAGTGGTGATCAGTGAGCGGGACAAGGGGGCGGTAGTCAGCCAGCCCTGAGGCGTCAGCCTGCGCCGGCCCATTCGCGGGTAAACCCGCTCCTACAGGGGCCTGCGGGCTTGCCCGCGAGGCCCTg                                                                                                                                <  2:756662/150‑1 (MQ=255)
                       aaGGCAGATCCTGAGTACCGGGAAGTGGTGATCAGTGAGCGGGACAAGGGGGCGGTAGTCAGCCAGCCCTGAGGCGTCAGCCTGCGCCGGCCCATTCGCGGGTAAACCCGCTCCTACAGGGCCTCGCGGGCAAGCCCGCAGGcccctgta                                                                                                                             <  2:1993149/150‑6 (MQ=38)
                         ggCAGATCCTGAGTACCGGGAAGTGGTGATCAGTGAGCGGGACAAGGGGGCGGTAGTCAGCCAGCCCTGAGGCGTCAGCCTGCGCCGGCCCATTCGCGGGTAAACCCGCTCCTACAGGGGCCTGCGGGCTTGCCCGCGAGGCCCTGTAAt                                                                                                                            <  1:1042506/150‑1 (MQ=255)
                          gCAGATCCTGAGTCACGGGAAGTGGTGATCAGTGAGCGGGACAAGGGGGCGGTAGTCAGCCAGCCCTGAGGCGTCAGCCTGCGCCGGCCCATTCGCGGGTAAACCCGCTCCTACAGGGGCCTGCGGGCTTGCCCGCGAGGCCCTGTAATg                                                                                                                           >  2:2434608/1‑150 (MQ=255)
                             gATCCTGAGTACCGGGAAGTGGTGATCAGTGAGCGGGACAAGGGGGCGGTAGTCAGCCAGCCCTGAGGCGTCAGCCTGCGCCGGCCCATTCGCGGGTAAACCCGCTCCTACAGGGGCCTGCGGGCTTGCCCGCGAGGCCCTGTAATGTCa                                                                                                                        >  2:139997/1‑150 (MQ=255)
                               tCCTGAGTACCGGGAAGTGGTGATCAGTGAGCGGGACAAGGGGGCGGTAGTCAGCCAGCCCTGAGGCGTCAGCCTGCGCCGGCCCATTCGCGGGTAAACCCGCTCCTACAGGGGCCTGCGGGCTTGCCCGCGAGGCCCTGTAATGTCAgg                                                                                                                      >  2:1070691/1‑150 (MQ=255)
                                  tGAGTACCGGGAAGTGGTGATCAGTGAGCGGGACAAGGGGGCGGTAGTCAGCCAGCCCTGAGGCGTCAGCCTGCGCCGGCCCATTCGCGGGTAAACCCGCTCCTACAGGGGCCTGc                                                                                                                                                     >  1:1957616/1‑116 (MQ=255)
                                  tGAGTACCGGGAAGTGGTGATCAGTGAGCGGGACAAGGGGGCGGTAGTCAGCCAGCCCTGAGGCGTCAGCCTGCGCCGGCCCATTCGCGGGTAAACCCGCTCCTACAGGGGCCTGc                                                                                                                                                     <  2:1957616/116‑1 (MQ=255)
                                  tGAGTACCGGGAAGTGGTGATCAGTGAGCGGGACAAGGGGGCGGTAGTCAGCCAGCCCTGAGGCGTCAGCCTGCGCCGGCCCATTCGCGGGTAAACCCGCTCCTACAGGGGCCTGCGGGCTTGCCCGCGAGGCCCTGTAATGTCAGGaaa                                                                                                                   >  2:761107/1‑150 (MQ=255)
                                  tGAGTACCGGGAAGTGGTGATCAGTGAGCGGGACAAGGGGGCGGTAGTCAGCCAGCCCTGAGGCGTCAGCCTGCGCCGGCCCATTCGCGGGTAAACCCGCTCCTACAGGGGCCTGCGGGCTTGCCCGCGAGGCCCTGTAATGTCAGGaaa                                                                                                                   >  2:611263/1‑150 (MQ=255)
                                   gAGTACCGGGAAGTGGTGATCAGTGAGCGGGACAAGGGGGCGGTAGTCAGCCAGCCCTGAGGCGTCAGCCTGCGCCGGCCCATTCGCGGGTAAACCCGCTCCTACAGGGGCCTGCGGGCTTGCCCGCGAGGCCCTGTAATGTCAGGaaaa                                                                                                                  <  2:1703535/150‑1 (MQ=255)
                                        ccGGGAAGTGGTGATCAGTGAGCGGGACAAGGGGGCGGTAGTCAGCCAGCCCTGAGGCGTCAGCCTGCGCCGGCCCATTCGCGGGTAAACCCGCTCCTACAGGGGCCTGCGGGCTTGCCCGCGAGGCCCTGTAATGTCAGGAAAACTCGa                                                                                                             <  1:761107/150‑1 (MQ=255)
                                         cGGGAAGTGGTGATCAGTGAGCGGGACAAGGGGGCGGTAGTCAGCCAGCCCTGAGGCGTCAGCCTGCGCCGGCCCATTCGCGGGTAAACCCGCTCCTACAGGGGCCTGCGGGCTTGCCCGCGAGGCCCTGTAATGTCAGGAAAACTCGaa                                                                                                            <  1:2119904/150‑1 (MQ=255)
                                         cGGGAAGTGGTGATCAGTGAGCGGGACAAGGGGGCGGTAGTCAGCCAGCCCTGAGGCGTCAGCCTGCGCCGGCCCATTCGCGGGTAAACCCGCTCCTACAGGGGCCTGCGGGCTTGCCCGCGAGGCCCTGTAATGTCAGGAAAACTCGaa                                                                                                            >  1:794214/1‑150 (MQ=255)
                                         cGGGAAGTGGTGATCAGTGAGCGGGACAAGGGGGCGGTAGTCAGCCAGCCCTGAGGCGTCAGCCTGCGCCGGCCCATTCGCGGGTAAACCCGCTCCTACAGGGGCCTGCGGGCTTGCCCGCGAGGCCCTGTAATGTCAGGAAAACTCGaa                                                                                                            >  1:2315705/1‑150 (MQ=255)
                                             aaGTGGTGATCAGTGAGCGGGACAAGGGGGCGGTAGTCAGCTAGCCCTGAGGCGTCAGCCTGCGCCGGCCCATTCGCGGGTAAACCCGCTCCTACAGGGGCCTGCGGGCTTGCCCGCGAGGCCCTGTAATGTCAGGAAAACTCGAACGgc                                                                                                        >  2:1646395/1‑150 (MQ=255)
                                             aaGTGGTGATCAGTGAGCGGGACAAGGGGGCGGTAGTCAGCCAGCCCTGAGGCGTCAGCCTGCGCCGGCCCATTCGCGGGTAAACCCGCTCCTACAGGGGCCTGCGGGCTTGCCCGCGAGGCCCTGTAATGTCAGGAAAACTCGAACGgc                                                                                                        <  1:1951105/150‑1 (MQ=255)
                                              aGTGGTGATCAGTGAGCGGGACAAGGGGGCGGTAGTCAGCCAGCCCTGAGGCGTCAGCCTGCGCCGGCCCATTCGCGGGTAAGCCCGCTCCTACAGGGGCCTGCGGGCTTGCCCGCGAGGCCCTGTAATGTCAGGAAAACTCGAACGgcg                                                                                                       <  1:86375/150‑1 (MQ=255)
                                                        aGTGAGCGGGACAAGGGGGCGGTAGTCAGCCAGCCCTGAGGCGTCAGCCTGCGCCGGCCCATTCTCGGGTAAACCCGCTCCTACAGGGGCCTGCGGGCTTGCCCGCGAGGCCCTGTAATGTCAGGAAACCTCGAACGGGGCCAACTGaaa                                                                                             <  2:611683/150‑1 (MQ=255)
                                                         gTGAGCGGGACAAGGGGGCGGTAGTCAGCCAGCCCTGAGGCGTCAGCCTGCGCCGGCCCATTCGCGGGTAAACCCGCTCCTACAGGGGCCTGCGGGCTTGCCCGCGAGGCCCTGGAATGTCAGGAAAACTCGAACGGCGCCAACTGAAAg                                                                                            >  2:1686082/1‑150 (MQ=255)
                                                             gCGGGACAAGGGGGCGGTAGTCAGCCAGCCCTGAGGCGTCAGCCTGCGCCGGCCCATTCGCGGGTAAACCCGCTCCTACAGGGGCCTGCGGGCTTGCCCGCGAGGCCCTGTAATGTCAGGAAAACTCGAACGGCGCCAACTGAAAGCCCt                                                                                        >  1:1498411/1‑150 (MQ=255)
                                                               gggACAAGGGGGCGGTAGTCAGCCAGCCCTGAGGCGTCAGCCTGCGCCGGCCCATTCGCGGGTAAACCCGCTCCTACAGGGGCCTCGCGGGCAAGCCcgcaggc                                                                                                                                     <  2:124022/104‑5 (MQ=37)
                                                               gggACAAGGGGGCGGTAGTCAGCCAGCCCTGAGGCGTCAGCCTGCGCCGGCCCATTCGCGGGTAAACCCGCTCCTACAGGGGCCTCGCGGGCAAGCCcgcaggc                                                                                                                                     >  1:124022/1‑100 (MQ=37)
                                                                 gACAAGGGGGCGGTAGTCAGCCAGCCCTGAGGCGTCAGCCTGCGCCGGCCCATTCGCGGGTAAACCCGCTCCTACAGGGGCCTGCgg                                                                                                                                                   <  1:1937962/87‑1 (MQ=255)
                                                                 gACAAGGGGGCGGTAGTCAGCCAGCCCTGAGGCGTCAGCCTGCGCCGGCCCATTCGCGGGTAAACCCGCTCCTACAGGGGCCTGCgg                                                                                                                                                   >  2:1937962/1‑87 (MQ=255)
                                                                    aaGGGGGCGGTAGTCAGCCAGCCCTGAGGCGTCAGCCTGCGCCGGCCCATTCGCGGGTAAACCCGCTCCTACAGGGGCCTGCGGGCTTGCCCGCGAGGCCCTGTAATGTCAGGaa                                                                                                                    <  1:2083007/115‑1 (MQ=255)
                                                                    aaGGGGGCGGTAGTCAGCCAGCCCTGAGGCGTCAGCCTGCGCCGGCCCATTCGCGGGTAAACCCGCTCCTACAGGGGCCTGCGGGCTTGCCCGCGAGGCCCTGTAATGTCAGGAAAACTCGAACGGCGCCAACTGAAAGCCCTGTTCGTc                                                                                 <  1:611263/150‑1 (MQ=255)
                                                                    aaGGGCGCGGTAGTCAGCCAGCCCTGAGGCGTCAGCCTGCGCCGGCCCATTCGCGGGTAAACCCGCTCCTACAGGGGCCTGCGGGCTTGCCCGCGAGGCCCTGTAATGTCAGGaa                                                                                                                    >  2:2083007/1‑115 (MQ=255)
                                                                      gggggCGGTAGTCAGCCAGCCCTGAGGCGTCAGCCTGCGCCGGCCCATTCGCGGGTAAACCCGCTCCTACAGGGGCCTGCGGGCTTGCCCGCGAGGCCCTGTAATGTCAGGAAAACTCGAACGgcgc                                                                                                      <  1:185820/127‑1 (MQ=255)
                                                                      gggggCGGTAGTCAGCCAGCCCTGAGGCGTCAGCCTGCGCCGGCCCATTCGCGGGTAAACCCGCTCCTACAGGGGCCTGCGGGCTTGCCCGCGAGGCCCTGTAATGTCAGGAAAACTCGAACGgcgc                                                                                                      >  2:185820/1‑127 (MQ=255)
                                                                      gggggCGGTAGTCAGCCAGCCCTGAGGCGTCAGCCTGCGCCGGCCCATTCGCGGGTAAACCCGCTCCTACAGGGGCCTGCGGGCTTGCCCGCGAGGCCCTGTAATGTCAGGAAAACTCGAACGGCGCCAACTGAAAGCCCTGTTCGTCaa                                                                               >  2:2056070/1‑150 (MQ=255)
                                                                             gTAGTCAGCCAGCCCTGAGGCGTCAGCCTGCGCCGGCCCATTCGCGGGTAAACCCGCTCCTACAGGGGCCTGCGGGCTTGCCCGCGAGGCCCTGTAATGTCAGGAAAACTCGAACGGCGCCAACTGAAAGCCCTGTTCGTCAACCTGCaa                                                                        <  2:794214/150‑1 (MQ=255)
                                                                              tAGTCAGCCAGCCCTGAGGCGTCAGCCTGCGCCGGCCCATTCGCGGGTAAACCCGCTCCTACAGGGGCCTGCGGGCTTGCCCGCGAGGCCCTGTAATGTCAGGAAAACTCGAACGGCGCCAACTGAAAGCCCTGTTCGTCAACCTGCAAg                                                                       <  1:2434608/150‑1 (MQ=255)
                                                                                gTCAGCCAGCCCTGAGGCGTCAGGCTGCGCCGGCCCATTCGCGGGTAAACCCGCTCCTACAGGGGCCTGCGGGTTTGCCCGCGAGGCCCTGTAATGTCAGGAAAACTCGAACGGCGCCAACTGATAGCCCTGTTCGTCAACCTGCTAGGc                                                                     >  2:1192398/1‑150 (MQ=255)
                                                                                 tCAGCCAGCCCTGAGGCGTCAGCCTGCGCCGGCCCATTCGCGGGTAAACCCGCTCCTACAGGGGCCTGCGGGCTTGCCCGCGAGGCCCTGTAATGTCAGGAAAACTCGAACGGCGCCAACTGAAAGCCCTGTTCGTCAACCTGCAAGGcc                                                                    >  1:1502104/1‑150 (MQ=255)
                                                                                 tCAGCCAGCCCTGAGGCGTCAGCCTGCGCCGGCCCATTCGCGGGTAAACCCGCTCCTACAGGGGCCTGCGGGCTTGCCCGCGAGGCCCTGTAATGTCAGGAAAACTCGAACGGCGCCAACTGAAAGCCCTGTTCGTCAACCTGCAAGGcc                                                                    >  1:230884/1‑150 (MQ=255)
                                                                                  cagccagcCCTGAGGCGTCAGCCTGCGCCGGCCCATTCGCGGGTAAACCCGCTCCTACAGGGGCCTGCGGGCTTGCCCGCGAGGCCCTGTAATGTCAGGAAAACTCGAACGGCGCCAACTGAAAGCCCTGTTCGTCAACCTGCAAGGccc                                                                   <  1:1325539/150‑1 (MQ=255)
                                                                                    gccagccCTGAGGCGTCAGCCTGCGCCGGCCCATTCGCGGGTAAACCCGCGCCTACAGGGGCCTGCGGGCTTGCCCGCGAGGCCCTGTAATGTCAGGAAAACTCGAACGGCGCCAACTGAAAGCCCTGTTCGTCAACCTGCAAGGCCCaa                                                                 >  2:2233209/1‑150 (MQ=255)
                                                                                     ccagccCTGAGGCGTCAGCCTGCGCCGGCCCATTCGCGGGTAAACCCGCTCCTACAGGGGCCTGCGGGCTTGCCCGCGAGGCCCTGTAATGTCAGGAAAACTCGAACGGCGTCAACTGAAAGCCCTGTTCGTCAACCTGCAAGGCCCAAc                                                                >  2:1033259/1‑150 (MQ=255)
                                                                                     ccagccCTGAGGCGTCAGCCTGCGCCGGCCCATTCGCGGGTAAACCCGCTCCTACAGGGGCCTGCGGGCTTGCCCGCGAGGCCCTGTAATGTCAGGAAAACTCGAACGGCGCCAACTGAAAGc                                                                                           <  2:1903713/123‑1 (MQ=255)
                                                                                     ccagccCTGAGGCGTCAGCCTGCGCCGGCCCATTCGCGGGTAAACCCGCTCCTACAGGGGCCTGCGGGCTTGCCCGCGAGGCCCTGTAATGTCAGGAAAACTCGAACGGCGCCAACTGAAAGc                                                                                           >  1:1903713/1‑123 (MQ=255)
                                                                                     ccagccCTGAGGCGTCAGCCTGCGCCGGCCCATTCGCGGGTAAACCCGCTCCTACAGGGGCCTGCGGGCTTGCCCGCGAGGCCCTGTAATGTCAGGAAAACTCGAACGGCGCCAACTGAAAGCCCTGTTCGTCAACCTGCAAGGCCCAAc                                                                >  2:1700304/1‑150 (MQ=255)
                                                                                         cccTGAGGCGTCAGCCTGCGCCGGCCCATTCGCGGGTAAACCCGCTCCTACAGGGGCCTGCGGGCTTGCCGGCGAGGCCCTGTAATGTCAGGAAAACTCGAACGGGGGCAACTGAAAGCCCTGTTCGTCAACCTGCAAGGCCCAACcgcg                                                            >  1:535398/1‑150 (MQ=255)
                                                                                         cccTGAGGCGTCAGCCTGCGCCGGCCCATTCGCGGGTAAACCCGCTCCTACAGGGGCCTGCGGGCTTGCCCGCGAGGCCCTGTAATGTCAGGAAAACTCGAACGGCGTCAACTGAAAGCCCTGTTCGTCAACCTGCAAGGCCCAACcgcg                                                            <  1:1033259/150‑1 (MQ=255)
                                                                                          ccTGAGGCGTCAGCCTGCGCCGGCCCATTCGCGGGTAAACCCGCTCCTACAGGGGCCTGCGGGCTTGCCCTCGAGGCCCTGTAATGTCAGGAAAACTCGAACGGCGCCAACTGAAAGCCCTGTTCGTCAACCTGCAAGGCCCAACCgcgg                                                           >  1:1353855/1‑150 (MQ=255)
                                                                                          ccTGAGGCGTCAGCCTGCGCCGGCCCATTCGCGGGTAAACCCGCTCCTACAGGGGCCTGCGGGCTTGCCCGCGAGGCCCTGTAATGTCAGGAAAACTCGAACGGCGCCAACTGAAAGCCCTGTTCGTCACCCTGCAAGGCCCAACCgcgg                                                           <  2:950292/150‑1 (MQ=255)
                                                                                           cTGAGGCGTCAGCCTGCGCCGGCCCATTCGCGGGTAAACCCGCTCCTACAGGGGCCTGCGGGCTTGCCCGCGAGGCCCTGTAATGTCAGGAAAACTCGAACGGCGCCAACTGAAAGCCCTGTTCGTCAACCTGCAAGGCCCAACCgcggc                                                          >  1:598621/1‑150 (MQ=255)
                                                                                              aGGCGTCAGCCTGCGCCGGCCCATTCGCGGGTAAACCCGCTCCTACAGGGGCCTGCGGGCTTGCCCGCGAGGCCCTGTAATGTCAGGAAAACTCGAACGGCGCCAACTGAAAGCCCTGTTCGTCAACCTGCAAGGCCCAACCGCGGCGGt                                                       <  1:537420/150‑1 (MQ=255)
                                                                                               ggCGTCAGCCTGCGCCGGCCCATTCGCGGGTAAACCCGCTCCTACAGGGGCCTGCGGGCTTGCCCGCGAGGCCCTGTAATGTCAGGAAAACTCGAACGGCGCCAACTGAAAGCCCTGTTCGTCAACCTGCAAGGCCCAACCGCGCCGGTc                                                      >  2:1738532/1‑150 (MQ=255)
                                                                                                gCGTCAGCCTGCGCCGGCCCATTCGCGGGTAAACCCGCTCCTACAGGGGCCTGCGGGCTTGCCCGCGAGGCCCTGTAATGTCAGGAAAACTCGAACGGCGCCAACTGAAAGCCCTGTTCGTCAACCTGCAAGGCCCAACCGCGGCGGTcc                                                     <  1:831374/150‑1 (MQ=255)
                                                                                                  gTCAGCCTGCGCCGGCCCATTCGCGGGTAAACCCGCTCCTACAGGGGCCTGCGGGCTTGCCCGCGAGGCCCTGTAATGTCAGGAAAACTCGAACGGCGCCAACTGAAAGCCCTGTTCGTCAACCTGCAAGGCCCAACCGCGGCGGTNCCa                                                   <  2:1502104/150‑1 (MQ=255)
                                                                                                  gTCAGCCTGCGCCGGCCCATTCGCGGGTAAACCCGCTCCTACAGGGGCCTGCGGGCTTGCCCGCGAGGCCCTGTAATGTCAGGAAAACTCGAACGGCGCCAACTGAAAGCCCTGTTCGTCAACCTGCAAGGCCCAACCGCGGCGGTCCCa                                                   <  2:230884/150‑1 (MQ=255)
                                                                                                   tCAGCCTGCGCCGGCCCATTCGCGGGTAAACCCGCTCCTACAGGGGCCTGCGGGCTTGCCCGCGAGGCCCTGTAATGTCAGGAAAACTCGAACGGCGCCAACTGAAAGCCCTGTTCGTCAACCTGCAAGGCCCAACCGCGGCGGTCCCAg                                                  <  1:1700304/150‑1 (MQ=255)
                                                                                                   tCAGCCTGCGCCGGCCCATTCGCGGGTAAACCCGCTCCTACAGGGGCCTGCGGGCTTGCCCGCGAGGCCCTGTAATGTCAGGAAAACTCGAACGGCGCCAACTGAAAGCCCTGTTCGTCAACCTGCAAGGCCCAACCGCGGCGGCCCCac                                                  >  2:946029/1‑149 (MQ=255)
                                                                                                    cAGCCTGCGCCGGCCCATTCGCGGGTAAACCCGCTCCTACAGGGGCCTGCGGGCTTGCCCGCGAGGCCCTGTAATGTCAGGAAAACTCGAACGGCGCCAACTGAAAGCCCTGt                                                                                      <  2:1392404/113‑1 (MQ=255)
                                                                                                    cAGCCTGCGCCGGCCCATTCGCGGGTAAACCCGCTCCTACAGGGGCCTGCGGGCTTGCCCGCGAGGCCCTGTAATGTCAGGAAAACTCGAACGGCGCCAACTGAAAGCCCTGt                                                                                      >  1:1392404/1‑113 (MQ=255)
                                                                                                    cAGCCTGCGCCGGCCCATTCGCGGGTAAACCCGCTCCTACAGGGGCCTGCGGGCTTGCCCGCGAGGCCCTGTAATGTCAGGAAAACTCGAACGGCGCCAACTGAAAGCCCTGTTCGTCAACCTGCAAGGCCCAACCGCGGCGGTCCCAGt                                                 <  2:2315705/150‑1 (MQ=255)
                                                                                                        cTGCGCCGGCCCATTCGCGGGTAAACCCGCTCCTACAGGGGCCTGCGGGCTTGCCCGCGAGGCCCTGTAATGTCAGGAAAACTCGAACGGCGCCAACTGAAAGCCCTGTTCGTCAACCTGCAAGGCCCAACCgcgg                                                           >  1:1905729/1‑136 (MQ=255)
                                                                                                        cTGCGCCGGCCCATTCGCGGGTAAACCCGCTCCTACAGGGGCCTGCGGGCTTGCCCGCGAGGCCCTGTAATGTCAGGAAAACTCGAACGGCGCCAACTGAAAGCCCTGTTCGTCAACCTGCAAGGCCCAACCgcgg                                                           <  2:1905729/136‑1 (MQ=255)
                                                                                                        cTGCGCCGGCCCATTCGCGGGTAAACCCGCTCCTACAGGGGCCTGCGGGCTTGCCCGCGAGGCCCTGTAATGTCAGGAAAACTCGAACGGCGCCAACTGAAAGCCCTGTTCGTCAACCTGCAAGGCCCAACCGCGGCGGTCCCAGTCGcc                                             >  1:1825853/1‑150 (MQ=255)
                                                                                                                 cccATTCGCGGGTAAACCCGCTCCTACAGGGGCCTGCGGGCTTGCCCGCGAGGCCCTGTAATGTCAGGAAAACTCGAACTGCGCCAACTGAAAGCCCTGTTCGTCAACCTGCAAGGCCCAACCGCGGCGGTCCCAGTCGCCCAGCACGAt                                    >  1:278352/1‑150 (MQ=255)
                                                                                                                          gggTAAACCCGCTCCTACAGGGGCCTGCGGGCTTGCCCGCGAGGCCCTGTAATGTCAGGAAAACTCGAACGGCGCCAACTGAAAGCCCTGTTCGTCAACCTGCAAGGCCCAACCGCGGCGGTCCCAGTCGCCCAGCACGATACGCCGTGc                           <  1:358334/150‑1 (MQ=255)
                                                                                                                                aCCCGCTCCTACAGGGGCCTGCGGGCTTGCCCGCGAGGCCCTGTAATGTCAGGAAAACTCGAACGGCGCCAACTGAAAGCCCTGTTCGTCAACCTGCAAGGCCCAACCGCGGCGGTCCCAGTCGCCCAGCACGATACGCCGTGCCGGCTg                     <  2:1825853/150‑1 (MQ=255)
                                                                                                                                 cccGCTCCTACAGGGGCCTGCGGGCTTGCCCGCGAGGCCCTGTAATGTCAGGAAAACTCGAACGGCGCCAACTGAAAGCCCTGTTCGTCAACCTGCAAg                                                                       <  2:1640417/99‑1 (MQ=255)
                                                                                                                                 cccGCTCCTACAGGGGCCTGCGGGCTTGCCCGCGAGGCCCTGTAATGTCAGGAAAACTCGAACGGCGCCAACTGAAAGCCCTGTTCGTCAACCTGCAAg                                                                       >  1:1640417/1‑99 (MQ=255)
                                                                                                                                 cccGCTCCTACAGGGGCCTGCGGGCTTGCCCGCGAGGCCCTGTAATGTCAGGAAAACTCGAACGGCGCCAACTGAAAGCCCTGTTCGTCAACCTGCAAGGCCCAACCGCGGCGGTCCCAGTCGCCCAGCACGATACGCCGTGCCGGCTgg                    <  2:1150982/150‑1 (MQ=255)
                                                                                                                                  ccGCTCCTACAGGGGCCTGCGGGCTTGCCCGCGAGGCCCTGTAATGTCAGGAAAACTCGAACGGCGCCAACTGAAAGCCCTGTTCGTCAACCTGCAAGGCCCAACCGCGGCGGTCCCAGTCGCCCAGCACGATACGCCGTGCCGGCTGGc                   >  1:743079/1‑150 (MQ=255)
                                                                                                                                  ccGCTCCTACAGGGGCCTGCGGGCTTGCCCGCGAGGCCCTGTAATGTCAGGAAAACTCGAACGGCGCCAACTGAAAGCCCTGTTCGTCAACCTGCAAGGCCCAACCGCGGCGGTCCCAGTCGCCCAGCACGATACGCCGTGCCGGCTGGc                   >  1:621468/1‑150 (MQ=255)
                                                                                                                                   cGCTCCTACAGGGGCCTGCGGGCTTGCCCGCGAGGCCCTGTAATGTCAGGAAAACTCGAACGGCGCCAACTGAAAGCCCTGTTCGTCAACCTGCAAGGCCCAACCGCGGCGGTCCCAGTCGCCCAGCACGATACGCCGTGCCGGCTGGcc                  >  1:2155428/1‑150 (MQ=255)
                                                                                                                                    gCTCCTACAGGGGCCTGCGGGCTTGCCCGCGAGGCCCTGTAATGTCAGGAAAACTCGAACGGCGCCAACTGAAAGCCCTGTTCGTCAACCTGCAAGGCCCAACCGCGGCGGTCCCAGTCGCCCCGCACGATACGCCGTGCCGGCTGGccc                 >  2:1012868/1‑149 (MQ=255)
                                                                                                                                      tCCTACAGGGGCCTGCGGGCTTGCCCGCGAGGCCCTGTAATGTCAGGAAAACTCGAACGGCGCCAACTGAAAGCCCTGTTCGTCAACCTGCAAGGCCCAACCGCGGCGGTCCCAGTCGCCCAGCACGATACGCCGTGCCGGCTGGCCATc               <  2:1639737/150‑1 (MQ=255)
                                                                                                                                        cTACAGGGGCCTGCGGGCTTGCCCGCGAGGCCCTGTAATGTCAGGAAAACTCGAACGGCGCCAACTGAAAGCCCTGTTCGTCAACCTGCAAGGCCCAACCGCGGCGGTCCCAGTCGCCCAGCACGATACGCCGTGCCGGCTGGCCATCGa             >  2:1938347/1‑150 (MQ=255)
                                                                                                                                        attacagggCCTCGCGGGCAAGCCCGCAGGCCCTGTAATGTCAGGAAAACTCGAACGGCGCCAACTGAAAGCCCTGTTCGTCAACCTGCAAGGCCCAACCGCGGCGGTCCCAGTCGCCCAGCACGATACGCCGTGCCGGCTGGCCATCGa             >  1:1009939/7‑150 (MQ=255)
                                                                                                                                         tACAGGGGCCTGCGGGCTTGCCCGCGAGGCCCTGTAATGTCAGGAAAACTCGAACGGCGCCAACTGAAAGCCCTGTTCGTCAACCTGCAAGGCCCAACCGCGGCGGTCCCa                                                   >  2:401046/1‑111 (MQ=255)
                                                                                                                                         tACAGGGGCCTGCGGGCTTGCCCGCGAGGCCCTGTAATGTCAGGAAAACTCGAACGGCGCCAACTGAAAGCCCTGTTCGTCAACCTGCAAGGCCCAACCGCGGCGGTCCCa                                                   <  1:401046/111‑1 (MQ=255)
                                                                                                                                         tACAGGGGCCTGCGGGCTTGCCCGCGAGGCCCTGTAATGTCAGGAAAACTCGAACGGCGCCAACTGAAAGCCCTGTTCGTCAACCTGCAAGGCCCAACCGCGGCGGTCCCAGTCGCCCAGCACGATACGCCGTGCCGGCTGGCCATCGAc            >  1:1148933/1‑150 (MQ=255)
                                                                                                                                                  cTGCGGGCTTGCCCGCGAGGCCCTGTAATGTCAGGAAAACTCGAACGGCGCCAACTGAAAGCCCTGTTCGTCAACCTGCAAGGCCCAACCGCGGCGGTCCCAGTCGCCCAGCACGATACGCCGTGCCGGCTGGCCATCGACCATCAGCtt   <  2:2155428/150‑1 (MQ=255)
                                                                                                                                                   tGCGGGCTTGCCCGCGAGGCCCTGTAATGTCAGGAAAACTCGAACGGCGCCAACTGAGAGCCCTGTTCGTCAACCTGCAAGGCCCAACCGCGGCGGTCCCAGTCGCCCAGCACGATACGCCGTGCCGGCTGGCCATCGACCATCAGCTtg  <  1:2176703/150‑1 (MQ=255)
                                                                                                                                                   tGCGGGCTTGCCCGCGAGGCCCTGTAATGTCAGGAAAACTCGAACGGCGCCAACTGAAAGCCCTGTTCGTCAACCTGCAAGGCCCAACCGCGGCGGTCCCAGTCGCCCAGCACGATACGCCGTGcc                          >  1:1433020/1‑126 (MQ=255)
                                                                                                                                                   tGCGGGCTTGCCCGCGAGGCCCTGTAATGTCAGGAAAACTCGAACGGCGCCAACTGAAAGCCCTGTTCGTCAACCTGCAAGGCCCAACCGCGGCGGTCCCAGTCGCCCAGCACGATACGCCGTGcc                          <  2:1433020/126‑1 (MQ=255)
                                                                                                                                                   tGCGGGCTTGCCCGCGAGGCCCTGTAATGTCAGGAAAACTCGAACGGCGCCAACTGAAAGCCCTGTTCGTCAACCTGCAAGGCCCAACCGCGGCGGTCCCAGTCGCCCAGCACGATACGCCGTGCCGGCTGGCCATCGACCATc        <  1:230136/144‑1 (MQ=255)
                                                                                                                                                   tGCGGGCTTGCCCGCGAGGCCCTGTAATGTCAGGAAAACTCGAACGGCGCCAACTGAAAGCCCTGTTCGTCAACCTGCAAGGCCCAACCGCGGCGGTCCCAGTCGCCCAGCACGATACGCCGTGCCGGCTGGCCATCGACCATc        >  2:230136/1‑144 (MQ=255)
                                                                                                                                                   tGCGGGCTTGCCCGCGAGGCCCTGTAATGTCAGGAAAACTCGAACGGCGCCAACTGAAAGCCCTGTTCGTCAACCTGCAAGGCCCAACCGCGGCGGTCCCAGTCGCCCAGCACGATACGCCGTGCCGGCTGGCCATCGACCATCAGCTtg  <  1:2176690/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TTCCGTTTACCGAGGCGCAGATCAAGGCAGATCCTGAGTACCGGGAAGTGGTGATCAGTGAGCGGGACAAGGGGGCGGTAGTCAGCCAGCCCTGAGGCGTCAGCCTGCGCCGGCCCATTCGCGGGTAAACCCGCTCCTACAGGGGCCTGCGGGCTTGCCCGCGAGGCCCTGTAATGTCAGGAAAACTCGAACGGCGCCAACTGAAAGCCCTGTTCGTCAACCTGCAAGGCCCAACCGCGGCGGTCCCAGTCGCCCAGCACGATACGCCGTGCCGGCTGGCCATCGACCATCAGCTTG  >  NC_002947/3300051‑3300347

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