Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 508,778:1 +G 6.9% coding (233/1890 nt) estP ← esterase EstP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947508,7781.G 6.9% 282.0 / 15.9 72coding (233/1890 nt)estPesterase EstP
Reads supporting (aligned to +/- strand):  ref base . (37/30);  new base G (2/3);  total (39/33)
Fisher's exact test for biased strand distribution p-value = 6.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.82e-01

ATTGGGGTTTCCTGGGGGGATGACGGTCTCGGACGTTGTCGTGATCGAGTCCAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGCTGTCGCCAAACACGATCAGGGTT  >  NC_002947/508630‑508924
                                                                                                                                                     |                                                                                                                                                  
aTTGGGGTTTCCTGGGGGGATGACGGTCTCGGACGTTGTCGTGATCGAGTCCAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCTACGGGTTAAGGTCGCCTGGCGCGACGc                                                                                                                                                     <  7:146979/149‑1 (MQ=255)
      gTTTCCTGGGGGGATGACGGTCTCGGACGTGGTCGTGATCGAGTCCAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATGGGGCTGGATACCTACGGATGTCGACGGGTTACGGTCGCCTGGCGCGCGGC‑CCAggc                                                                                                                                              >  4:198927/1‑149 (MQ=255)
      gTTTCCTGGAGGGATGACGGTCTCGGACGTTGTCGTGATCGAGTCCAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAggc                                                                                                                                              <  5:336286/149‑1 (MQ=255)
       tttCCTGGGGGGATGACGGTCTCGGACGTTGTCGTGATCGAGTCCAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGG‑GCGTCGCGCCAggcg                                                                                                                                             <  6:347880/149‑1 (MQ=255)
         tCCTGGGGGGATGACGGTCTCGGACGTTGTCGTGATCGAGTCCAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGCCGCCTGG‑GCGTCGCGCCAggcga                                                                                                                                            <  4:278517/148‑2 (MQ=255)
           cTGGGGGGATGACGGTCTCGGACGTTGTCGTGATCGAGTCCAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCa                                                                                                                                         >  7:201030/1‑149 (MQ=255)
                gggATGACGGTCTCGGACGTTGTCGTGATCGAGTCCAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAg                                                                                                                                    >  5:377340/1‑149 (MQ=255)
                 ggaggaCGGTCTCGGACGTTGTCGTGATCGACTCCAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAgg                                                                                                                                   <  3:198927/149‑1 (MQ=255)
                 ggATGACGGTCTCGGACGTTGTCGTGATCGAGTCCAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAgg                                                                                                                                   <  7:182309/149‑1 (MQ=255)
                 ggATGACGGTCTCGGACGTTGTCGTGATCGAGTCCAGGATCTGCGGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACGTCCGGATGTCGCCGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAgg                                                                                                                                   <  5:308253/149‑1 (MQ=255)
                              ggACGTTGTCGTGATCGAGTCCAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTAGTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTACGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGaca                                                                                                                      >  8:323891/1‑149 (MQ=255)
                              ggACGTTGTCGTGATCGAGTACAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGaca                                                                                                                      >  2:269902/1‑149 (MQ=255)
                                   ttGTCGTGCTCGAGTCCAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGGCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCggg                                                                                                                 <  2:343546/149‑1 (MQ=255)
                                      tCGTGATCGAGTCCAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGc                                                                                                               >  7:263786/1‑148 (MQ=255)
                                        gTGATCGAGTCCAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCg                                                                                                              >  3:113582/1‑147 (MQ=255)
                                           aTCGAGTCCAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGtt                                                                                                         >  2:108876/1‑149 (MQ=255)
                                                    aGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGc                                                                                                >  7:26751/1‑149 (MQ=255)
                                                    aGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGc                                                                                                >  5:184874/1‑149 (MQ=255)
                                                    aGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGG‑GCGTCGCGCCAggcgac                                                                                                                                           >  3:123507/1‑104 (MQ=255)
                                                    aGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGG‑GCGTCGCGCCAggcgac                                                                                                                                           <  4:123507/106‑3 (MQ=255)
                                                          tgctgGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACg                                                                                          >  7:316527/1‑149 (MQ=255)
                                                            ctgGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCa                                                                                            <  6:377340/145‑1 (MQ=255)
                                                            ctgGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACg                                                                                          <  1:365858/147‑1 (MQ=255)
                                                            ctgGGGGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTTGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGt                                                                                        <  2:331167/149‑1 (MQ=255)
                                                                gtggtgTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGt                                                                                    <  7:340376/149‑1 (MQ=255)
                                                                 tggtgTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTg                                                                                   <  8:316527/149‑1 (MQ=255)
                                                                     gtACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAAc                                                                               >  6:77767/1‑149 (MQ=255)
                                                                      tACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAAc                                                                               <  1:269902/148‑1 (MQ=255)
                                                                       aCCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAAcgc                                                                             >  2:211433/1‑149 (MQ=255)
                                                                         ccgccgACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAAcgcgc                                                                           >  5:249296/1‑149 (MQ=255)
                                                                         ccgccgACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAAcgcgc                                                                           >  5:144497/1‑149 (MQ=255)
                                                                         ccgccgACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAAcgcgc                                                                           >  5:250456/1‑149 (MQ=255)
                                                                         ccgccgACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAAcgcgc                                                                           >  2:122826/1‑149 (MQ=255)
                                                                         ccgccgACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAAcgcgc                                                                           <  8:37611/149‑1 (MQ=255)
                                                                          cgccgACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAAcgcgcg                                                                          <  8:201030/149‑1 (MQ=255)
                                                                           gccgACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACgcgcgc                                                                         >  4:88998/1‑149 (MQ=255)
                                                                                cTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTggg                                                                    <  8:12/149‑1 (MQ=255)
                                                                                 tGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTgggg                                                                   >  6:306954/1‑149 (MQ=255)
                                                                                    ccAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTAc                                                                >  6:128084/1‑149 (MQ=255)
                                                                                       gTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGc                                                             <  4:333515/149‑1 (MQ=255)
                                                                                           ttaCCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCcaa                                                          <  1:108876/148‑1 (MQ=255)
                                                                                           ttaCCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCcaa                                                          <  3:88998/148‑1 (MQ=255)
                                                                                             aCCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCcaaca                                                        <  6:144497/148‑1 (MQ=255)
                                                                                              ccATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGt                                                                                                    >  4:336975/1‑103 (MQ=255)
                                                                                              ccATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGt                                                                                                    <  3:336975/103‑1 (MQ=255)
                                                                                              ccATAGGGCTGGATACATACGGATGTCGACGGGTTAAGGTCGACTGTCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCcaacaa                                                       <  6:272033/148‑1 (MQ=255)
                                                                                               cATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCcaac                                                         >  5:234902/1‑145 (MQ=255)
                                                                                               cATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGa                                                     <  6:234902/149‑1 (MQ=255)
                                                                                                 tCGGGCTGTATACCTACGGATGTCGACGGGTTAAGTTCGCCTGGAGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGTTACCGCCAACAAGATc                                                   <  8:263786/149‑1 (MQ=255)
                                                                                                   gggCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCgg                                                 <  5:77767/149‑1 (MQ=255)
                                                                                                    ggCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCgg                                                 >  3:250170/1‑148 (MQ=255)
                                                                                                       tGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGaaa                                             >  1:135665/1‑149 (MQ=255)
                                                                                                       tGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGaaa                                             >  5:357331/1‑149 (MQ=255)
                                                                                                         gATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGaaa                                             <  4:353618/147‑1 (MQ=255)
                                                                                                         gATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACt                                           >  4:204832/1‑149 (MQ=255)
                                                                                                              cTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATc                                                   >  6:207014/1‑136 (MQ=255)
                                                                                                              cTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATc                                                   <  5:207014/136‑1 (MQ=255)
                                                                                                              cTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCg                                      >  2:39922/1‑149 (MQ=255)
                                                                                                              cTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCg                                      >  3:115196/1‑149 (MQ=255)
                                                                                                              cTACGGATGTCGACGGGATAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCg                                      >  5:27305/1‑149 (MQ=255)
                                                                                                               tACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATAGGGAAACTGCCCgg                                     >  2:267485/1‑149 (MQ=255)
                                                                                                                         gACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGa                            <  6:184874/148‑1 (MQ=255)
                                                                                                                         gACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGa                            <  5:306954/148‑1 (MQ=255)
                                                                                                                         gACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGa                            <  6:250456/148‑1 (MQ=255)
                                                                                                                              gTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGCTGt                      <  1:211433/149‑1 (MQ=255)
                                                                                                                               ttAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGCTGTc                     >  1:151622/1‑149 (MQ=255)
                                                                                                                                  aGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCCTCGGCGAGGCTGTCGc                   >  4:309625/1‑148 (MQ=255)
                                                                                                                                   ggccGCCTGG‑GCGTCGCGCCAGGCGACCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGCTGTCGcc                  >  2:323080/4‑148 (MQ=255)
                                                                                                                                    gTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGCTGTCGCCa                 >  6:75963/1‑148 (MQ=255)
                                                                                                                                         cTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGCTGTCGCCAAacac            >  7:260512/1‑148 (MQ=255)
                                                                                                                                         cTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGCTGTCGCCAAacac            >  1:255703/1‑148 (MQ=255)
                                                                                                                                                gACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGCTGTCGCCAAACACGATCAggg    <  8:260512/149‑1 (MQ=255)
                                                                                                                                                   gC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGCTGTCGCCAAACACGATCAGGGtt  >  3:380738/1‑148 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
ATTGGGGTTTCCTGGGGGGATGACGGTCTCGGACGTTGTCGTGATCGAGTCCAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGCTGTCGCCAAACACGATCAGGGTT  >  NC_002947/508630‑508924

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