Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 459,904 A→C 5.0% intergenic (‑50/+29) pqqB ← / ← pqqA coenzyme PQQ synthesis protein B/coenzyme PQQ synthesis protein A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947459,9040AC 5.0% 314.7 / 6.4 100intergenic (‑50/+29)pqqB/pqqAcoenzyme PQQ synthesis protein B/coenzyme PQQ synthesis protein A
Reads supporting (aligned to +/- strand):  ref base A (42/53);  new base C (3/2);  total (45/55)
Fisher's exact test for biased strand distribution p-value = 6.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.38e-01

Only 100 of 101 total aligned reads displayed.

TCAGGGTGCCATCACGGTAGCCCTTGCAGTTGACGCAGTTGCAGTTCCACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAAC‑‑‑GAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAGGGTTGGGTTGTTCAGCTA  >  NC_002947/459758‑460044
                                                                                                                                                     |                                                                                                                                            
tCAGGGTGCCATCACGGTAGCCCTTGCAGTTGACGCAGTTGCAGTTCCACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAAC‑‑‑GAAAACGCCCCGGCa                                                                                                                                              <  4:403975/147‑1 (MQ=255)
   gggTGCCATCACGGTAGCCCTTGCAGTTGACGCAGTTGCAGTTCCACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAAC‑‑‑GAAAACGCCCCGGCAGGcc                                                                                                                                          >  4:426338/1‑148 (MQ=255)
    ggTGCCATCACGGTAGCCCTTGCAGTTGACGCAGTTGCAGTTCCACTGCTGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCATATACCGGATCGAATAAAC‑‑‑GAAAACGCCCCGGCAGGca                                                                                                                                          >  6:389943/1‑146 (MQ=255)
        ccATCACGGTAGCCCTTGCAGTCGACGCAGTTGCAGTTCCACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAAC‑‑‑GAAAACGCCCCGGCAGGCCGAGGCg                                                                                                                                    >  7:66358/1‑149 (MQ=255)
         cATCACGGTAGCCCTTGCAGTTGACGCAGTTGCAGTTCCACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAAC‑‑‑GAAAACGCCCCGGCAGGCCGAGGCGt                                                                                                                                   <  6:263582/149‑1 (MQ=255)
             aCGGTAGCCCTTGCAGTTGACGCAGTTGCAGTTCCACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAAC‑‑‑GAAAACGCCTCGGCCTGCCGGGGCGttttc                                                                                                                               <  1:144018/149‑5 (MQ=255)
             aCGGTAGCCCTTGCAGTTGACGCAGTTGCAGTTCCACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAAC‑‑‑GAAAACGCCCCGGCAGGCCGAGGCGTGaaa                                                                                                                               >  7:255811/1‑149 (MQ=255)
             aCGGTAGCCCTTGCAGTTGACGCAGTTGCAGTTCCACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAAC‑‑‑GAAAACGCCCCGGCAGGCCGAGGCGTGaaa                                                                                                                               <  6:443706/149‑1 (MQ=255)
              cGGTAGCCCTTGCAGTTGACGCAGTTGCAGTTCCACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAAC‑‑‑GAAAACGCCCCGGCAGGCCGAGGCGTGAAAc                                                                                                                              <  1:134376/149‑1 (MQ=255)
                gTAGCCCTTGCAGTTGACGCAGTTGCAGTTCCACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAAC‑‑‑GAAAACGCCTCGGCCTGCCGGGGCGttttcgtt                                                                                                                            <  2:140727/149‑8 (MQ=255)
                   gCCCTTGCAGTTGACGCAGTTGCAGTTCCACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAAC‑‑‑GAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTaa                                                                                                                         <  1:319488/149‑1 (MQ=255)
                   gCCCTTGCAGTTGACGCAGTTGCAGTTCCACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAAC‑‑‑GAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTaa                                                                                                                         <  4:41110/149‑1 (MQ=255)
                   gCCCTTGCAGTTGACGCAGTTGCAGTTCCACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAAC‑‑‑GAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTaa                                                                                                                         <  8:263717/149‑1 (MQ=255)
                     ccTTGCAGTTGACGCAGTTGCAGTTCCACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAAC‑‑‑GAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGc                                                                                                                       >  5:271060/1‑149 (MQ=255)
                     ccTTGCAGTTGACGCAGTTGCAGTTCCACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAAC‑‑‑GAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGc                                                                                                                       >  4:249741/1‑149 (MQ=255)
                     ccTTGCAGTTGACGCAGGTGCAGTTCCACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGCATCGAAAAAAC‑‑‑GAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGc                                                                                                                       <  5:389943/149‑1 (MQ=255)
                           aGTTGACGCAGTTGCAGTTCCACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAAC‑‑‑GAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCt                                                                                                                 <  6:114466/149‑1 (MQ=255)
                               gACGCAGTTGCAGTTCCACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAAC‑‑‑GAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAg                                                                                                              <  5:318762/148‑1 (MQ=255)
                                aCGCAGTTGCAGTTCCACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAAC‑‑‑GAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCg                                                                                                            <  7:113019/149‑1 (MQ=255)
                                aCGCAGTTGCAGTTCCACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAAC‑‑‑GAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCg                                                                                                            <  3:143267/149‑1 (MQ=255)
                                aCGCAGTTGCAGTTCCACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAAC‑‑‑GAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCg                                                                                                            >  3:339921/1‑149 (MQ=255)
                                aCGCAGTTGCAGTTCCACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAAC‑‑‑GAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCAGGCTTAgg                                                                                                             <  6:436938/148‑2 (MQ=255)
                                  gcagttgcagttCCACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACTGGGGGAGAGCCCCGGCAGCCCGAGGCGAGAAACCCTAACCCTTGCTAAGc                                                                                                             >  1:364312/1‑149 (MQ=255)
                                  gcagttgcagttCCACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAAC‑‑‑GAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGt                                                                                                          >  8:389896/1‑149 (MQ=255)
                                   cagttgcagttCCACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAAC‑‑‑GAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGtt                                                                                                         <  3:110952/149‑1 (MQ=255)
                                    agttgcagttCCACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAAC‑‑‑GAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTa                                                                                                        <  1:178738/149‑1 (MQ=255)
                                      ttgcagttCCACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAAC‑‑‑GAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGc                                                                                                      >  5:199321/1‑149 (MQ=255)
                                        gcagttCCACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAAC‑‑‑GAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCg                                                                                                     >  4:219275/1‑148 (MQ=255)
                                        gcagttCCACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAAC‑‑‑GAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCg                                                                                                     <  3:426338/148‑1 (MQ=255)
                                        gcagttCCACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAAC‑‑‑GAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCg                                                                                                     >  3:38052/1‑148 (MQ=255)
                                         cagttCCTCTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAAC‑‑‑GAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGaa                                                                                                   <  8:168335/149‑1 (MQ=255)
                                         cagttCCACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAAC‑‑‑GAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGaa                                                                                                   <  4:44740/149‑1 (MQ=255)
                                         cagttCCACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAAC‑‑‑GAAAACGCCCCGGCAGGACGACGCTTCAAACCCTAAGCCTGGCTTAGCGGTAAGCg                                                                                                     >  1:329831/1‑147 (MQ=255)
                                            ttCCACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAAC‑‑‑GAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTa                                                                                                >  7:116405/1‑149 (MQ=255)
                                                aCTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAAC‑‑‑GAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTattgc                                                                                            >  7:147403/1‑145 (MQ=255)
                                                aCTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAAC‑‑‑GAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTCCATg                                                                                            <  3:249741/149‑1 (MQ=255)
                                                aCTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAAC‑‑‑GAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATg                                                                                            >  6:172967/1‑149 (MQ=255)
                                                aCTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAAC‑‑‑GAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATg                                                                                            >  5:166319/1‑149 (MQ=255)
                                                aCTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAAC‑‑‑GAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATg                                                                                            >  7:410198/1‑149 (MQ=255)
                                                aCTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAAC‑‑‑GAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATg                                                                                            >  1:297960/1‑149 (MQ=255)
                                                 cTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAAC‑‑‑GAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATgg                                                                                           >  3:270132/1‑149 (MQ=255)
                                                    cGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAAC‑‑‑GAAAACGTCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGt                                                                                          <  3:314271/147‑1 (MQ=255)
                                                       gAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAAC‑‑‑GAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACtt                                                                                     >  3:262864/1‑149 (MQ=255)
                                                        aaCCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCGCTCTCTATGCAGAAACCGGATCGAAAAACC‑‑‑CAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTc                                                                                    <  2:215573/149‑1 (MQ=255)
                                                          cccGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAAC‑‑‑GAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGa                                                                                  <  8:410198/149‑1 (MQ=255)
                                                           ccGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAAC‑‑‑GAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCg                                                                                   >  3:347039/1‑147 (MQ=255)
                                                             gCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAAC‑‑‑GAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGc                                                                               <  1:75936/149‑1 (MQ=255)
                                                               cACCGGCGGCGGAGCCGAGAACCTGGATGTAGATGGCCACTCTCTATTCAGAAACCGGATCGAAAATAG‑‑‑TAAAACGCCCCGGCAGGGCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTTCATTGTGAGTTCGATGc                                                                               <  2:78485/147‑1 (MQ=255)
                                                                   ggcggcggAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAAC‑‑‑GAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATAc                                                                         >  1:300352/1‑149 (MQ=255)
                                                                   ggcggcggAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAAC‑‑‑GAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATAc                                                                         <  6:199321/149‑1 (MQ=255)
                                                                     cggcggAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAAC‑‑‑GAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGc                                                                       <  3:322760/149‑1 (MQ=255)
                                                                      ggcggAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAAC‑‑‑GAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACg                                                                        >  8:281530/1‑147 (MQ=255)
                                                                        cggAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAAC‑‑‑GAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAgg                                                                    >  1:206435/1‑149 (MQ=255)
                                                                            gCCGAGAACCTGGATGTACATGGCCACTCTCTATGCGGAAACCGGATCGAAAAAAC‑‑‑GAAAAAGCCCCCGCAGGCCGAGGCGTGAAACCCAAAGCCAGGCGTAGCGGTTTGCGAAGTAAATGGTGACTTCGAAGCCGATACGCAGGTCAg                                                                >  4:217186/1‑149 (MQ=255)
                                                                               gagaACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAAC‑‑‑GAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAgtg                                                              >  7:33970/1‑148 (MQ=255)
                                                                                agaACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAAC‑‑‑GAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACg                                                                        >  7:202022/1‑137 (MQ=255)
                                                                                agaACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAAC‑‑‑GAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTa                                                            <  8:66358/149‑1 (MQ=255)
                                                                                agaACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAAC‑‑‑GAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTa                                                            <  2:206435/149‑1 (MQ=255)
                                                                                agaAACTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAAC‑‑‑GAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACg                                                                        <  8:202022/137‑1 (MQ=255)
                                                                                 gaACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAAC‑‑‑GAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTAt                                                           <  7:389896/149‑1 (MQ=255)
                                                                                          tGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAAC‑‑‑GAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTg                                                  <  3:406414/149‑1 (MQ=255)
                                                                                            tACATGGCCACTCTCTATGCGGAAACCGGATCGAAAAAAC‑‑‑GAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGt                                                <  3:217186/149‑1 (MQ=255)
                                                                                            tACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAAC‑‑‑GAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTgg                                                 <  4:347039/148‑1 (MQ=255)
                                                                                             aCATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAAC‑‑‑GAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTc                                               >  7:56882/1‑149 (MQ=255)
                                                                                              cATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAAC‑‑‑GAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTc                                               <  5:149947/148‑1 (MQ=255)
                                                                                                 ggCCACGCTCAATGCATAACCCGGACCGAAAAACC‑‑‑CAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCcaca                                           <  2:216557/149‑1 (MQ=255)
                                                                                                  gcaaCTCTCTATGAATAAACCGGATCGCAAAAAC‑‑‑GAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACAt                                          <  6:233040/146‑1 (MQ=255)
                                                                                                  gCCACTCTCTATGCAGAAACCGGATCGAAAAAAC‑‑‑GAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTgg                                                                                                                   >  2:249472/1‑76 (MQ=255)
                                                                                                  gCCACTCTCTATGCAGAAACCGGATCGAAAAAAC‑‑‑GAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTgg                                                                                                                   <  1:249472/76‑1 (MQ=255)
                                                                                                  gCCACTCTCTATGCAGAAACCGGATCGAAAAAAC‑‑‑GAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCcaca                                           <  1:388359/148‑1 (MQ=255)
                                                                                                  gCCACTCTCTATGCAGAAACCGGATCGAAAAAAC‑‑‑GAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCcaca                                           <  8:116405/148‑1 (MQ=255)
                                                                                                  gCCACTCTCTATGCAGAAACCGGATCGAAAAAAC‑‑‑GAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCcaca                                           <  4:38052/148‑1 (MQ=255)
                                                                                                    cACTCTCTATGCAGAAACCGGATCGAAAAAAC‑‑‑GAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAgt                                                               >  3:295066/1‑126 (MQ=255)
                                                                                                    cACTCTCTATGCAGAAACCGGATCGAAAAAAC‑‑‑GAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAgt                                                               <  4:295066/126‑1 (MQ=255)
                                                                                                    cACTCTCTATGCAGAAACCGGATCGAAAAAAC‑‑‑GAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATgg                                        <  8:33970/149‑1 (MQ=255)
                                                                                                    cACTCTCTATGCAGAAACCGGATCGAAAAAAC‑‑‑GAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATg                                         >  5:71745/1‑148 (MQ=255)
                                                                                                     aCTCTCTATGCAGAAACCGGATCGAAAAAAC‑‑‑GAAAACGCCCCCGCACGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGAAGGTCAGTGTATTAAGTTTTCTTCCAAATgg                                        >  6:102992/1‑148 (MQ=255)
                                                                                                         tctATGCAGAAACCGGATCGAAAAAAC‑‑‑GAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAAt                                     <  8:56882/147‑1 (MQ=255)
                                                                                                           tATGCAGAAACCGGATCGAAAAAAC‑‑‑GAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACt                                 <  3:219275/149‑1 (MQ=255)
                                                                                                             tGCAGAAACCGGATCGAAAAAAC‑‑‑GAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTcc                               >  2:191799/1‑149 (MQ=255)
                                                                                                               cAGAAACCGGATCGAAAAAAC‑‑‑GAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCtt                             >  7:68472/1‑149 (MQ=255)
                                                                                                               cAGAAACCGGATCGAAAAAAC‑‑‑GAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCtt                             <  7:73043/149‑1 (MQ=255)
                                                                                                                  aaaCCGGATCGAAAAAAC‑‑‑GAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTa                                   <  2:191457/140‑1 (MQ=255)
                                                                                                                  aaaCCGGATCGAAAAAAC‑‑‑GAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTa                                   >  1:191457/1‑140 (MQ=255)
                                                                                                                     ccGGATCGAAAAAAC‑‑‑GAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAAt                       >  7:342070/1‑149 (MQ=255)
                                                                                                                         aTCGAAAAAAC‑‑‑GAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAgg                   <  8:255811/149‑1 (MQ=255)
                                                                                                                            gAAAAAAC‑‑‑GAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAgggtt                <  2:297960/149‑1 (MQ=255)
                                                                                                                            gAAAAAAC‑‑‑GAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAgggtt                <  6:71745/149‑1 (MQ=255)
                                                                                                                             aaaaaaC‑‑‑GAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAgggttg               >  7:411681/1‑149 (MQ=255)
                                                                                                                             aaaaaaC‑‑‑GAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAgggttg               >  8:307159/1‑149 (MQ=255)
                                                                                                                              aaaaaC‑‑‑GAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAgggttgg              <  7:281530/149‑1 (MQ=255)
                                                                                                                               aaaac‑‑‑gaacacgCACAGTCAGTCCGATGCGTGAAACACTAAGCCTGGATTAGCTGTTAGAGAAGTACATGTTGACTTCGAAGCCGATACGCAGTTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAgtgttggg             <  8:342070/149‑1 (MQ=255)
                                                                                                                               aaaac‑‑‑gaaaacgCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAgggttggg             <  4:339921/149‑1 (MQ=255)
                                                                                                                                 aac‑‑‑gaaaacgCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAgggttgggtt           <  5:172967/149‑1 (MQ=255)
                                                                                                                                       gaaaaagCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAGGGTTGGgttgtt        >  1:87128/1‑149 (MQ=255)
                                                                                                                                           acgCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAGGGTTGGGTTGTTCAGc    <  4:270132/149‑1 (MQ=255)
                                                                                                                                            cgCCTCGGCCTGCCGGGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAGGGTTGGGTTGTTCAGc    >  2:135125/1‑148 (MQ=255)
                                                                                                                                             gCCTCGGCCTGCCGGGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAGGGTTGGGTTGTTCAGcta  >  4:60406/1‑149 (MQ=255)
                                                                                                                                             gCCTCGGCCTGCCGGGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAGGGTTGGGTTGTTCAGcta  >  8:32756/1‑149 (MQ=255)
                                                                                                                                             gCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAGGGTTGGGTTGTTCAGcta  >  7:18295/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                            
TCAGGGTGCCATCACGGTAGCCCTTGCAGTTGACGCAGTTGCAGTTCCACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAAC‑‑‑GAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAGGGTTGGGTTGTTCAGCTA  >  NC_002947/459758‑460044

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: