Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 92,186 T→C 5.0% Q65R (CAG→CGG) ‡ trpA ← tryptophan synthase subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00294792,1860TC 5.0% 301.2 / 4.0 99Q65R (CAG→CGG) ‡trpAtryptophan synthase subunit alpha
Reads supporting (aligned to +/- strand):  ref base T (51/43);  new base C (2/3);  total (53/46)
Fisher's exact test for biased strand distribution p-value = 6.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.89e-01

GTGTCCACGCCAAAGCGGTGGATAGGGTTGTAGTAGCCCATCAGCACGATCGGCGTGGTCTGGTTATCCACACGGA‑‑ATTCGCGAACCATCTGCAGGGTCTTGACCAAGGTCTGGCCAGCTTCAAGGGCGCGCAGGGTGGCCAGCTGGATGGCCACGCCATCGGCCATCGGGTCGGTGAACGGCATGCCCAGCTCGATCACGTCGGCGCCAGCAGCCGGCAGGCCCTTGAGGATCTGCAACGAGGCGTCGTAGCCCGGGTCGCCAGCTGTGACGAAGGTGACCAGTGCCGAGC  >  NC_002947/92042‑92333
                                                                                                                                                   |                                                                                                                                                    
gtgtCCACGCCAAAGCGGTGGATAGGGTTGTAGTAGCCCATCAGCACGATCGGCGTGGTCTGGTTATCCACACGGA‑‑ATTCGCGAACCATCTGCAGGGTCTTGACCAAGGTCTGGCCAGCTTCAAGGGCGCGCAGGGTGGCCAGCTGGAt                                                                                                                                                  >  7:143505/1‑149 (MQ=255)
     cACGCCAAAGCGGTGGATAGGGTTGTAGTAGCCCATCAGCACGATCGGCGTGGTCTGGTTATCCACACGGA‑‑ATTCGCGAACCATCTGCAGGGTCTTGACCAAGGTCTGGCCAGCTTCAAGGGCGCGCAGGGTGGCCATCCAGCTGGCCa                                                                                                                                             <  2:136280/149‑1 (MQ=255)
     cACGCCAAAGCGGTGGATAGGGTTGTAGTAGCCCATCAGCACGAACGGCGTGGTCTGGTTATCCACACGGA‑‑ATTCGCGAACCATCTGCAGGGTCTTGACCAAGGTCTGGCCAGCTTCAAGGGCGCGCAGGGTGGCCAGCTGGATGGCCa                                                                                                                                             >  1:89800/1‑149 (MQ=255)
       cGCCAAAGCGGTGGATAGGGTTGTAGTAGCCCATCAGCACGATCGGCGTGGTCTGGTTATCCACACGGA‑‑ATTCGCGAACCATCTGCAGGGTCTTGACCAAGGTCTGGCCAGCTTCAAGGGCGCGCAGGGTGGCCATCCAGCTGGCCAcc                                                                                                                                           <  6:137608/149‑2 (MQ=255)
        gCCAAAGCGGTGGATAGGGTTGTAGTAGACCATCAGCACGATCGGCGTGGTCTGGTTATCCACACGGA‑‑ATTCGCGAACCATCTGCAGGGTCTTGACCAAGGTCTGGCCAGCTTCAAGGGCGCGCAGGGTGGCCAGCTGGATGGCCACgc                                                                                                                                          >  6:348909/1‑149 (MQ=255)
         ccAAAGCGGTGGATAGGGTTGTAGTAGCCCATCAGCACGATCGGCGTGGTCTGGTTATCCACACGGA‑‑ATTCGCGAACCATCTGCAGGGTCTTGACCAAGGTCTGGCCAGCTTCAAGGGCGCGCAGGGTGGCCAGCTGGATGGCCACgcc                                                                                                                                         <  3:90481/149‑1 (MQ=255)
                         ggTTGTAGTAGCCCATCAGCACGATCGGCGTGGTCTGGTTATCCACACGGA‑‑ATTCGCGAACCATCTGCAGGGTCTTGACCAAGGTCTGGCCAGCTTCAAGGGCGCGCAGGGTGGCCAGCTGGATGGCCACGCCATCGGCCATCGGGTc                                                                                                                          <  2:212982/148‑1 (MQ=255)
                         ggTTGTAGTAGCCCATCAGCACGATCGGCGTGGTCTGGTTATCCACACGGA‑‑ATTCGCGAACCATCTGCAGGGTCTTGACCAAGGTCTGGCCAGCTTCAAGGGCGCGCAGGGTGGCCAGCTGGATGGCCACGCCATCGGCCATCGGGTc                                                                                                                          >  1:212982/1‑148 (MQ=255)
                         ggTTGTAGTAGCCCATCAGCACGATCGGCGTGGTCTGGTTATCCACACGGA‑‑ATTCGCGAACCATCTGCAGGGTCTTGACCAAGGTCTGGCCAGCTTCAAGGGCGCGCAGGGTGGCCAGCTGGATGGCCACGCCATCGGCCATCGGGTCg                                                                                                                         >  3:354472/1‑149 (MQ=255)
                         ggTTGTAGTAGCCCATCAGCACGATCGGCGTGGTCTGGTTATCCACACGGA‑‑ATTCGCGAACCATCTGCAGGGTCTTGACCAAGGTCTGGCCAGCTTCAAGGGCGCGCAGGGTGGCCAGCTGGATGGCCACGCCATCGGCCATCGGGTCg                                                                                                                         <  2:330815/149‑1 (MQ=255)
                          gTTGTAGTAGCCCATCAGCACGATCGGCGTGGTCTGGTTATCCACACGGA‑‑ATTAGCGAACCATCTGCAGGGTCTTGACCAAGGTCTGGCCAGCTTCAAGGGCGCGCAGGGTGGCCAGCTGGATGGCCACGCCATCGGCCATCGGGTCgg                                                                                                                       >  8:19099/1‑149 (MQ=255)
                             gtagtaGCCCATCAGCACGATCGGCGTGGTCTGGTTATCCACACGGA‑‑ATTCGCGAACCATCTGCAGGGTCTTGACCAAGGTCTGGCCAGCTTCAAGGGCGCGCAGGGTGGCCAGCTGGATGGCCACGCCATCGGCCATCGGGTCGGTg                                                                                                                     <  2:320184/148‑1 (MQ=255)
                              tagtagCCCATCAGCACGATCGGCGTGGTCTGGTTATCCACCTCTTAATATTCGCGAACCATCTGCAGGGTCTTGACCAAGGTCTGGCCAGCTTCAAGGGCGCGCAGGGTGGCCAGCTGGATGGCCACGCCATCGGCCATCGGGTCgg                                                                                                                       <  1:364536/148‑1 (MQ=255)
                                   gCCCATCAGCACGATCGGCGTGGTCTGGTTATCCACACGGA‑‑ATTCGCGAACCATCTGCAGGGTCTTGACCAAGGTCTGGCCAGCTTCAAGGGCGCGCAGGGTGGCCAGCTGGATGGCCACGCCATCGGCCATCGGGTCGGTGAACGGc                                                                                                               <  1:321591/148‑1 (MQ=255)
                                   gCCCATCAGCACGATCGGCGTGGTCTGGTTATCCACACGGA‑‑ATTCGCGAACCATCTGCAGGGTCTTGACCAAGGTCTGGCCAGCTTCAAGGGCGCGCAGGGTGGCCAGCTGGATGGCCACGCCATCGGCCATCGGGTCGGTGAACGGc                                                                                                               <  5:121547/148‑1 (MQ=255)
                                    cccATCAGCACGATCGGCGTGGTCTGGTTATCCACACGGA‑‑ATTCGCGAACCATCTGCAGGGTCTTGACCAAGGTCTGGCCAGCTTCAAGGGCGCGCAGGGTGGCCAGCTGGATGGCCACgcc                                                                                                                                         <  3:326509/122‑1 (MQ=255)
                                    cccATCAGCACGATCGGCGTGGTCTGGTTATCCACACGGA‑‑ATTCGCGAACCATCTGCAGGGTCTTGACCAAGGTCTGGCCAGCTTCAAGGGCGCGCAGGGTGGCCAGCTGGATGGCCACgcc                                                                                                                                         >  4:326509/1‑122 (MQ=255)
                                    cccATCAGCACGATCGGCGTGGTCTGGTTATCCACACGGA‑‑ATTCGCGAACCATCTGCAGGGTCTTGACCAAGGTCTGGCCAGCTTCAAGGGCGCGCAGGGTGGCCAGCTGGATGGCCACGCCATCGGGCATCGGGTCGGTGAACGGCAt                                                                                                             >  4:293802/1‑149 (MQ=255)
                                     ccATCAGCACGATCGGCGTGGTCTGGTTATCCACACGGA‑‑ATTCGCGAACCATCTGCAGGGTCTTGACCAAGGTCTGGCCAGCTTCAAGGGCGCGCAGGGTGGCCAGCTGGATGGCCACGCCATCGGCCATCGGGTCGGTGAACGGCATg                                                                                                            >  1:259222/1‑149 (MQ=255)
                                       aTCAGCACGATCGGCGTGGTCTGGTTATCCACACGGA‑‑ATTCTCGAACCATCTGCAGGGTCTTGACCAAGGTCTGGCCAGCTTCAAGGGCGCTCAGGGTGGCCAGCTGGATGGCCACGCCATCGGCCATCGGGTCGGTGAAc                                                                                                                  <  8:219863/141‑1 (MQ=255)
                                       aTCAGCACGATCGGCGTGGTCTGGTTATCCACACGGA‑‑ATTCGCGAACCATCTGCAGGGTCTTGACCAAGGTCTGGCCAGCTTCAAGGGCGCGCAGGGTGGCCAGCTGGATGGCCACGCCATCGGCCATCGGGTCGGTGAAc                                                                                                                  >  7:219863/1‑141 (MQ=255)
                                       aTCAGCACGATCGGCGTGGTCTGGTTATCCACACGGA‑‑ATTCGCGAACCATCTGCAGGGTCTTGACCAAGGTCTGGCCAGCTTCAAGGGCGCGCAGGGTGGCCAGCTGGATGGCCACGCCATCGGCCATCGGGTCGGTGAACGGCATGcc                                                                                                          <  8:412580/149‑1 (MQ=255)
                                       aTCAGCACGATCGGCGTGGTCTGGTTATCCACACGGA‑‑ATTCGCGAACCATCTGCAGGGTCTTGACCAAGGTCTGGCCAGCTTCAAGGGCGCGCAGGGTGGCCAGCTGGATGGCCACGCCATCGGCCATCGGGTCGGTGAACGGCATGcc                                                                                                          >  1:215968/1‑149 (MQ=255)
                                           gCACGATCGGCGTGGTCTGGTTATCCACACGGA‑‑ATTCGCGAACCATCTGCAGGGTCTTGACCAAGGTCTGGCCAGCTTCAAGGGCGCGCAGGGTGGCCAGCTGGATGGCCACGCCATCGGCCATCGGGTCGGTGAACGGCATGCCCAGc                                                                                                      <  8:143505/149‑1 (MQ=255)
                                            cACGATCGGCGTGGTCTGGTTATCCACACGGA‑‑ATTCGCGAACCATCTGCAGGGTCTTGACCAAGGTCTGGCCAGCTTCAAGGCCGCGCAGGGTGGCCAGCTGGATGGCCACGCCATCGGCCATCGGGACGGTGAACGGCATGCCCAGCt                                                                                                     >  4:37591/1‑149 (MQ=255)
                                              cGATCGGCGTGGTCTGGTTATCCACACGGA‑‑ATTCGCGAACCATCTGCAGGGTCTTGACCAAGGTCTGGCCAGCTTCAAGGGCGCGCAGGGTGGCCAGCTGGATGGCCACGCCATCGGCCATCGGGTCGGTGAACGGCATGCCCAGCTCg                                                                                                   <  3:282517/149‑1 (MQ=255)
                                                aTCGGCGTGGTCTGGTTATCCACACGGA‑‑ATTCGCGAACCATCTGCAGGGTCTTGACCAAGGTCTGGCCAGCTTCAAGGGCGCGCAGGGTGGCCAGCTGGATGGCCACGCCATCGGCCATCGGGTCGGTGAACGGCATGCCCAGCTCGAt                                                                                                 <  2:89800/149‑1 (MQ=255)
                                                   ggCGTGGTCTGGTTATCCACACGGA‑‑ATTCGCGAACCATCTGCAGGGTCTTGACCAAGGTCTGGCCAGCTTCAAGGGCGCGCAGGGTGGCCATCCAGCTGGCCAccct                                                                                                                                         >  8:66070/1‑104 (MQ=255)
                                                   ggCGTGGTCTGGTTATCCACACGGA‑‑ATTCGCGAACCATCTGCAGGGTCTTGACCAAGGTCTGGCCAGCTTCAAGGGCGCGCAGGGTGGCCAGCTGGATGGCCACgccatcggccatc                                                                                                                               >  2:158255/1‑117 (MQ=255)
                                                   ggCGTGGTCTGGTTATCCACACGGA‑‑ATTCGCGAACCATCTGCAGGGTCTTGACCAAGGTCTGGCCAGCTTCAAGGGCGCGCAGGGTGGCCAGCTGGATGGCCACgccatcggccatc                                                                                                                               <  1:158255/117‑1 (MQ=255)
                                                   ggCGTGGTCTGGTTATCCACACGGA‑‑ATTCGCGAACCATCTGCAGGGTCTTGACCAAGGTCTGGCCAGCTTCAAGGGCGCGCAGGCTGGCCATCCAGCTGGCCAccct                                                                                                                                         <  7:66070/107‑4 (MQ=255)
                                                    gCGTGGTCTGGTTATCCACACGGA‑‑ATTCGCGAACCATCTGCAGGGTCTTGACCAAGGTCTGGCCAGCTTCAAGGGCGCGCAGGGTGGCCAGCTGGATGGCCACGCCATCGGCCATCGGGTCGGTGAACGGCATGCCCAGCTCGATCACg                                                                                             >  8:311003/1‑149 (MQ=255)
                                                     cGTGGTCTGGTTATCCACACGGA‑‑ATTCGCGAACCATCTGCAGGGTCTTGACCAAGGTCTGGCCAGCTTCAAGGGCGCGCAGGGTGGCCAGCTGGATGGCCACGCCATCGGCCATCGGGTCGGTGAACGGc                                                                                                               <  5:421216/130‑1 (MQ=255)
                                                     cGTGGTCTGGTTATCCACACGGA‑‑ATTCGCGAACCATCTGCAGGGTCTTGACCAAGGTCTGGCCAGCTTCAAGGGCGCGCAGGGTGGCCAGCTGGATGGCCACGCCATCGGCCATCGGGTCGGTGAACGGc                                                                                                               >  6:421216/1‑130 (MQ=255)
                                                      gTGGTCTGGTTATCCACACGGA‑‑ATTCGCGAACCATCTGCAGGGTCTTGACCAAGGTCTGGCCAGCTTCAAGGGCGCGCAGGGTGGCCAGCTGGATGGCCACGCCATCGGCCATCGGGTCGGTGAACGGCATGCCCAGCTCGATCACGTc                                                                                           <  2:76547/149‑1 (MQ=255)
                                                       tGGTCTGGTTATCCACACGGA‑‑ATTCGCGAACCATCTGCAGGGTCTTGACCAAGGTCTGGCCAGCTTCAAGGGCGCGCAGGGTGGCCAGCTGGATGGCCACGCCATCGGCCATCGGGTCGGTGAACGGCATGCCCAGCTCGATCACGTCg                                                                                          >  7:49919/1‑149 (MQ=255)
                                                       tGGTCTGGTTATCCACACGGA‑‑ATTCGCGAACCATCTGCAGGGTCTTGACCAAGGTCTGGCCAGCTTCAAGGGCGCGCAGGGTGGCCAGCTGGATGGCCACGCCATCGGCCATCGGGTCGGTGAACGGCATGCCCAGCTCGATCACGTCg                                                                                          <  8:248922/149‑1 (MQ=255)
                                                        ggTCTGGTTATCCACACGGA‑‑ATTCGCGAACCATCTGCAGGGTCTTGACCAAGGTCTGGCCAGCTTCAAGGGCGCGCAGGGTGGCCAGCTGGATGGCCACGCCATCGGCCATCGGGTCGGTGAACGGCATGCCCAGCTCGATCACGTCgg                                                                                         >  3:79830/1‑149 (MQ=255)
                                                         gTCTGGTTATCCACACGGA‑‑ATTCGCGAACCATCTGCAGGGTCTTGACCAAGGTCTGGCCAGCTTCAAGGGCGCGCAGGGTGGCCAGCTGGATGGCCACGCCATCGGCCATCGGGTCGGTGAACGGCATGCCCAGCTCGATCACGTCgg                                                                                         <  4:79830/148‑1 (MQ=255)
                                                          tCTGGTTATCCACACGGA‑‑ATTCGCGAACCATCTGCAGGGTCTTGACCAAGGTCTGGCCAGCTTCAAGGGCGCGCAGGGTGGCCAGCTGGATGGCCACGCCATCGGCCATCGGGTCGGTGAACGGCATGCCCAGCTCGATCACGTCGgcg                                                                                       >  6:264683/1‑149 (MQ=255)
                                                           cTGGTTATCCACACGGA‑‑ATTCGCGAACCATCTGCAGGGTCTTGACCAAGGTCTGGCCAGCTTCAAGGGCGCGCAGGGTGGCCAGCTGGATGGCCACGCCATCGGCCATCGGGTCGGTGAACGGCATGCCCAGCTCGATCACGTCGgcgc                                                                                      >  6:286793/1‑149 (MQ=255)
                                                              gTTATCCACACGGA‑‑ATTCGCGAACCATCTGCAGGGTCTTGACCAAGGTCTGGCCAGCTTCAAGGGCGCGCAGGGTGGCCAGCTGGATGGCCACGCCATCGGCCATCGGGTCGGTGAACGGCATGCCCAGCTCGATCACGTCGgcgc                                                                                      <  4:126404/146‑1 (MQ=255)
                                                              gTTATCCACACGGA‑‑ATTCGCGAACCATCTGCAGGGTCTTGACCAAGGTCTGGCCAGCTTCAAGGGCGCGCAGGGTGGCCAGCTGGATGGCCACGCCATCGGCCATCGGGTCGGTGAACGGCATGCCCAGCTCGATCACGTCGGCGCcag                                                                                   >  7:278291/1‑149 (MQ=255)
                                                                    cacaCGGA‑‑ATTCGCGAACCATCTGCAGGGTCTTGACCAAGGTCTGGACAGCTTCAAGGGCGCGCAGGGTGGCCAGCTGGATGGCCACGCCATCGGCCATCGGGTCGGTGAACGGCATGCCCAGCTCGATCACGTCGGCGCCAGCAGCCg                                                                             >  4:150849/1‑149 (MQ=255)
                                                                     acacGGA‑‑ATTCGCGAACCATCTGCAGGGTCTTGACCAAGGTCTGGCCAGCTTCAAGGGCGCGCAGGGTGGCCAGCTGGATGGCCACGCCATCGGCCATCGGGTCGGTGAACGGCATGCCCAGCTCGATCAc                                                                                              <  3:70460/131‑1 (MQ=255)
                                                                     acacGGA‑‑ATTCGCGAACCATCTGCAGGGTCTTGACCAAGGTCTGGCCAGCTTCAAGGGCGCGCAGGGTGGCCAGCTGGATGGCCACGCCATCGGCCATCGGGTCGGTGAACGGCATGCCCAGCTCGATCAc                                                                                              >  4:70460/1‑131 (MQ=255)
                                                                           gA‑‑ATTCGCGAACCATCTGCAGGGTCTTGACCAAGGTCTGGCCAGCTTCAAGGGCGCGCAGGGTGGCCAGCTGGATGGCCACGCCATCGGCCATCGGGTCGGTGAACGGCATGCCCAGCTCGATCACGTCGGCGCCAGCAGCCGGCAGGc                                                                       <  5:348909/149‑1 (MQ=255)
                                                                                ttCGCGAACCATCTGCAGGGTCTTGACCAAGGTCTGGCCAGCTTCAAGGGCGCGCAGGGTGGCCAGCTGGATGGCCACGCCATCGGCCATCGGGTCGGTGAACGGCATGCCCAGCTCGATCACGTCGGCGCCAGCAGCCGGCAGGCCCt                                                                    >  1:402732/1‑149 (MQ=255)
                                                                                   gcgAACCATCTGCAGGGTCTTGACCAAGGTCTGGCCAGCTTCAAGGGCGCGCAGGGTGGCCAGCTGGATGGCGACGCCATCGGGCATCGGGTCGGTGAACGGCATGCCCAGCTCGATCACGTCGGCGCCAGCAGCCGGCAGGCCGttta                                                                 <  6:2075/149‑3 (MQ=255)
                                                                                   gcgAACCATCTGCAGGGTCTTGACCAAGGTCTGGCCAGCTTCAAGGGCGCGCAGGGTGGCCAGCTGGATGGCCACGCCATCGGCCATCGGGTCGGTGAACGGCATGCCCAGCTCGATCACGTCGGCGCCagcagc                                                                               >  3:173040/1‑135 (MQ=255)
                                                                                   gcgAACCATCTGCAGGGTCTTGACCAAGGTCTGGCCAGCTTCAAGGGCGCGCAGGGTGGCCAGCTGGATGGCCACGCCATCGGCCATCGGGTCGGTGAACGGCATGCCCAGCTCGATCACGTCGGCGCCagcagc                                                                               <  4:173040/135‑1 (MQ=255)
                                                                                       aCCATCTGCAGGGTCTTGACCAAGGTCTGGCCAGCTTCAAGGGCGCGCAGGGTGGCCAGCTGGATGGCCACGCCATCGGCCATCGGGTCGGTGAACGGCATGCCCAGCTCGATCACGTCGGCGCCAGCAGCCGGCAGGCCCTTGAGGAt                                                             <  7:311003/149‑1 (MQ=255)
                                                                                             tGCAGGGTCTTGACCAAGGTCTGGCCAGCTTCAAGGGCGCGCAGGGTGGCCAGCTGGATGGCCACGCCATCGGCCATCGGGTCGGTGAACGGCATGCCCAGCTCGATCACGTCGGCGCCAGCAGCCGGCAGGCCCTTGAGGCTCTGcca                                                       >  6:301054/1‑147 (MQ=255)
                                                                                               cAGGGTCTTGACCAAGGTCTGGCCAGCTTCAAGGGCGCGCAGGGTGGCCAGCTGGATGGCCACGCCATCGGCCATCGGGTCGGTGAACGGCATGCCCAGCTCGATCACGTCGGCGCCAGCAGCCGGCAGGCCCTTGAGGATCTGc                                                         <  3:331549/145‑1 (MQ=255)
                                                                                               cAGGGTCTTGACCAAGGTCTGGCCAGCTTCAAGGGCGCGCAGGGTGGCCAGCTGGATGGCCACGCCATCGGCCATCGGGTCGGTGAACGGCATGCCCAGCTCGATCACGTCGGCGCCAGCAGCCGGCAGGCCCTTGAGGATCTGCaa                                                       <  8:49919/147‑1 (MQ=255)
                                                                                                 gtgTCTTGACCAAGGTCTGGCCAGCTTCAAGGGCGCGCAGGGTGGCCAGCTGGATGGCCACGCCAACGGCCATCGGGGCGGTGCACGGCATGCCCAGCGCGAACACGACGGCGCCCGCAGCCGCCAGGCCGGGGAGGCACTGCAACgcg                                                   >  8:272514/3‑147 (MQ=255)
                                                                                                 gggTCTTGACCAAGGTCTGGCCAGCTTCAAGGGCGCGCAGGGTGGCCAGCTGGATGGCCACGCCATCGGCCATCGGGTCGGTGAACGGCATGCCCAGCTCGATCCCGTCGGCGCCAGCAGCCGGCAGGCCCTTGAGGATCTGCAACGAg                                                   >  7:102742/1‑149 (MQ=255)
                                                                                                 gggTCTTGACCAAGGTCTGGCCAGCTTCAAGGGCGCGCAGGGTGGCCAGCTGGATGGCCACGCCATCGGCCATCGGGTCGGTGAACGGCATGCCCAGCTCGATCACGTCGGCGCCAGCAGCCGGCAGGccc                                                                     >  7:188833/1‑131 (MQ=255)
                                                                                                 gggTCTTGACCAAGGTCTGGCCAGCTTCAAGGGCGCGCAGGGTGGCCAGCTGGATGGCCACGCCATCGGCCATCGGGTCGGTGAACGGCATGCCCAGCTCGATCACGTCGGCGCCAGCAGCCGGCAGGccc                                                                     <  8:188833/131‑1 (MQ=255)
                                                                                                 gggTCTTGACCAAGGTCTGGCCAGCTTCAAGGGCGCGCAGGGTGGCCAGCTGGATGGCCACGCCATCGGCCATCGGGTCGGTGAACGGCATGCCCAGCTCGATCACGTCGGCGCCAGCAGCCGGCAGGCCCTTGAGGATCTGCAACGAg                                                   <  2:259222/149‑1 (MQ=255)
                                                                                                  ggTCTTGACCAAGGTCTGGCCAGCTTCAAGGGCGCGCAGGGTGGCCAGCTGGATGGCCACGCCATCGGCCATCGGGTCGGTGAACGGCATGCCCAGCTCGATCACGTCGGCGCCAGCAGCCGGCAGGCCCTTGAGGATCTGCAACGAgg                                                  >  4:130335/1‑149 (MQ=255)
                                                                                                    tCTTGACCAAGGTCTGGCCAGCTTCAAGGGCGCGCAGGGTGGCCAGCTGGATGGCCACGCCATCGGCCATCGGGTCGGTGAACGGCATGCCCAGCTCGATCACGTCGGCGCCAGCAGCCGGCAGGCCCTTGAGGATCTGCAACGAGGcg                                                >  8:352904/1‑149 (MQ=255)
                                                                                                      ttGACCAAGGTCTGGCCAGCTTCAAGGGCGCGCAGGGTGGCCAGCTGGATGGCCACGCCATCGGCCATCGGGTCGGTGAACGGCATGCCCAGCTCGATCACGTCGGCGCCAGCAGCCGGCAGGCCCTTGAGGATCTGCAACGAGGcgtc                                              >  5:360180/1‑149 (MQ=255)
                                                                                                        gACCAAGGTCTGGCCAGCTTCAAGGGCGCGCAGGGTGGCCAGCTGGATGGCCACGCCATCGGCCATCGGGTCGGTGAACGGCATGCCCAGCTCGATCACGTCGGCGCCAGCAGCCGGCAGGCCCTTGAGGATCTGCAACGAGGcgtc                                              <  6:265856/147‑1 (MQ=255)
                                                                                                           cAAGGTCTGGCCAGCTTCAAGGGCGCGCAGGGTGGCCAGCTGGATGGCCACGCCATCGGCCATCGGGTCGGTGAACGGCATGCCCAGCTCGATCACGTCGGCGCCAGCAGCCGGCAGGCCCTTGAGGATCTGCAACGAGGCGTCGTAGc                                         >  1:100961/1‑149 (MQ=255)
                                                                                                           cAAGGTCTGGCCAGCTTCAAGGGCGCGCAGGGTGGCCAGCTGGATGGCCACGCCATCGGCCATCGGGTCGGTGAACGGCATGCCCAGCTCGATCACGTCGGCGCCAGCAGCCGGCAGGCCCTTGAGGATCTGCAACGAGGCGTCGTAGc                                         <  5:301054/149‑1 (MQ=255)
                                                                                                           cAAGGTCTGGCCAGCTTCAAGGGCGCGCAGGGTGGCCAGCTGGATGGCCACGCCATCGGCCATCGGGTCGGTGAACGGCATGCCCAGCTCGATCACGTCGGCGCCAGCAGCCGGCAGGCCCTTGAGGATCTGCAACGAGGCGTCGTAGc                                         >  7:25263/1‑149 (MQ=255)
                                                                                                               gTCTGGCCAGCTTCAAGGGCGCGCAGGGTGGCCAGCTGGATGGCCACGCCATCGGCCATCGGGTCGGTGAACGGCATGCCCAGCTCGATCACGTCGGCGCCAGCAGCCGGCAGGCCCTTGAGGATCTGCAACGAGGCGTCGTAGCCCgg                                     >  3:117497/1‑149 (MQ=255)
                                                                                                                tCTGGCCAGCTTCAAGGGCGCGCAGGGTGGCCAGCTGGATGGCCACGCCATCGGCCATCGGGTCGCGTGAACGGCATGCCCAGCTCGATCACGTCGGCGCCAGCAGCCGGCAGGCCCTTGAGGATCTGCAACGAGGCGTCGTAGCCCgg                                     >  1:131890/1‑149 (MQ=255)
                                                                                                                  tGGCCAGCTTCAAGGGCGCGCAGGGTGGCCAGCTGGATGGCCACGCCATCGGCCATCGGGTCGGTGAACGGCATGCCCAGCTCGATCACGTCGGCGCCAGCAGCCGGCAGGCCCTTGAGGATCTGCAACGAGGCGTCGTAGCCCGGGTc                                  >  1:382171/1‑149 (MQ=255)
                                                                                                                  tGGCCAGCTTCAAGGGCGCGCAGGGTGGCCAGCTGGATGGCCACGCCATCGGCCATCGGGTCGGTGAACGGCATGCCCAGCTCGATCACGTCGGCGCCAGCAGCCGGCAGGCCCTTGAGGATCTGCAACGAGGCGTCGTAGCCCGGGTc                                  >  3:247407/1‑149 (MQ=255)
                                                                                                                    gtcAGCTTCAAGGGCGCGCAGGGTGGCCAGCTGGATGGCCACGCCATCGGCCATCGGGTCGGTGAACGGCATGCCCAGCTCGATCACGTCGGCGCCAGCAGCCGGCAGGCCCTTGAGGATCTGCAACGAGGCGTCGTAGCCCGGGTc                                  <  3:150849/145‑1 (MQ=255)
                                                                                                                    gCCAGCTTCAAGGGCGCGCAGGGTGGCCAGCTGGATGGCCACGCCATCGGCCATCGGGTCGGTGAACGGCATGCCCAGCTCGATCACGTCGGCGCCAGCAGCCGGCAGGCCCTTGAGGATCTGCAACGa                                                    >  1:295171/1‑129 (MQ=255)
                                                                                                                    gCCAGCTTCAAGGGCGCGCAGGGTGGCCAGCTGGATGGCCACGCCATCGGCCATCGGGTCGGTGAACGGCATGCCCAGCTCGATCACGTCGGCGCCAGCAGCCGGCAGGCCCTTGAGGATCTGCAACGa                                                    <  2:295171/129‑1 (MQ=255)
                                                                                                                    gCCAGCTTCAAGGGCGCGCAGGGTGGCCAGCTGGATGGCCACGCCATCGGCCATCGGGTCGGTGAACGGCATGCCCAGCTCGATCACGTCGGCGCCAGCAGCCGGCAGGCCCTTGAGGATCTGCAACGAGGCGTCGTAGCCCGGGTc                                  <  6:360180/147‑1 (MQ=255)
                                                                                                                    gCCAGCTTCAAGGGCGCGCAGGGTGGCCAGCTGGATGGCCACGCCATCGGCCATCGGGTCGGTGAACGGCATGCCCAGCTCGATCACGTCGGCGCCAGCAGCCGGCAGGCCCTTGAGGATCTGCAACGAGGCGTCGTAGCCCGGGTCGc                                >  8:257929/1‑149 (MQ=255)
                                                                                                                      cAGCTTCAAGGGCGCGCAGGGTGGCCAGCTGGATGGCCACGCCATCGGCCATCGGGTCGGTGAACGGCATGCCCAGCTCGATCGCGTCGGCGCCAGCAGCCGGCAGGCCCTTGAGGATCTGCAACGAGGCGTCGTAGCCCGGGTCGCCa                              <  4:358244/149‑1 (MQ=255)
                                                                                                                         cTTCAAGGGCGCGCAGGGTGGCCAGCTGGATGGCCACGCCATCGGCCATCGGGTCGGTGAACGGCATGCCCAGCTCGATCACGTCGGCGCCagcagc                                                                               >  1:379142/1‑97 (MQ=255)
                                                                                                                         cTTCAAGGGCGCGCAGGGTGGCCAGCTGGATGGCCACGCCATCGGCCATCGGGTCGGTGAACGGCATGCCCAGCTCGATCACGTCGGCGCCagcagc                                                                               <  2:379142/97‑1 (MQ=255)
                                                                                                                          ttCAAGGGCGCGCAGGGTGGCCAGCTGGATGGCCACGCCATCGGCCATCGGGTCGGTGAACGGCATGTCCAGCTCGATCACGTCGGCGCCAGCAGCCGGCAGGCCCTTGAGGATCTGCAACGAGGCGTCGTAGCCCGGGTCGCCAGCtg                          >  8:29497/1‑149 (MQ=255)
                                                                                                                           tCAAGGGCGCGCAGGGTGGCCAGCTGGATGGCCACGCCATCGGCCATCGGGTCGGTGAACGGCATGCCCAGCTCGATCACGTCGGCGCCAGCAGCCGGCAGGCCCTTGAGGATCTGCAACGAGGCGTCGTAGCCCGGGTCGCCAGCtgt                         >  4:297225/1‑149 (MQ=255)
                                                                                                                            cAAGGGCGCGCAGGGTGGCCAGCTGGATGGCCACGCCATCGGCCATCGGTTCGGTGAACGGCATGCCCAGCTCGATCACGTCGGCGCCAGCAGCCGGCAGGCCCTTGAGGATCTGCAACGAGGcgtcgt                                            <  3:18418/129‑1 (MQ=255)
                                                                                                                            cAAGGGCGCGCAGGGTGGCCAGCTGGATGGCCACGCCATCGGCCATCGGGTCGGTGAACGGCATGCCCAGCTCGATCACGTCGGCGCCAGCAGCCGGCAGGCCCTTGAGGATCTGCAACGAGGcgtcgt                                            >  4:18418/1‑129 (MQ=255)
                                                                                                                            cAAGGGCGCGCAGGGTGGCCAGCTGGATGGCCACGCCATCGGCCATCGGGTCGGTGAACGGCATGCCCAGCTCGATCACGTCGGCGCCAGCAGCCGGCAGGCCCTTGAGGATCTGCAACGAGGCGTCGTAGCCCGGGTCGCCAGCtgt                         <  5:286793/148‑1 (MQ=255)
                                                                                                                             aaGGGCGCGCAGGGTGGCCAGCTGGATGGCCACGCCATCGGCCATCGGGTCGGTGAACGGCATGCCCAGCTCGATCACGTCGGCGCCAGCAGCCGGCAGGCCCTTGAGGATCTGCAACGAGGCGTCGTAGCCCGGGTCGCCAGCTGTGa                       >  1:106004/1‑149 (MQ=255)
                                                                                                                               gggCGCGCAGGGTGGCCAGCTGGATGGCCACGCCATCGGCCATCGGGTCGGTGAACGGCATGCCCAGCTCGATCACGTCGGCGCCAGCAGCCGGCAGGCCCTTGAGGATCTGCAACGAGGCGTCGTAGCCCGGGTCGCCAGCTGTGAc                      >  4:308487/1‑148 (MQ=255)
                                                                                                                                 gcgcgcAGGGTGGCCAGCTGGATGGCCACGCCATCGGCCATCGGGTCGGTGAACGGCATGCCCAGCTCGATCACGTCGGCGCCAGCAGCCGGCAGGCCCTTGAGGATCTGCAACGAGGCGTCGTAGCCCGGGTCGCCAGCTGTGACGaa                   <  2:100961/149‑1 (MQ=255)
                                                                                                                                  cgcgcAGGGTGGCCAGCTGGATGGCCACGCCATCGGCCATCGGGTCGGTGAACGGCATGCCCAGCTCGATCACGTCGGCGCCAGCAGCCGGCAGGCCCTTGAGGATCTGCAACGAGGCGTCGTAGCCCGGGTCGCCAGCTGTGACGAAg                  <  4:354472/149‑1 (MQ=255)
                                                                                                                                  cgcgcAGGGTGGCCAGCTGGATGGCCACGCCATCGGCCATCGGGTCGGTGAACGGCATGCCCAGCTCGATCACGTCGGCGCCAGCAGCCGGCAGGCCCTTGAGGATCTGCAACGAGGCGTCGTAGCCCGGGTCGCCAGCTGTGACGAAg                  <  2:402732/149‑1 (MQ=255)
                                                                                                                                    gatggcgTGGCCATCCAGCTGGCCACCCCATCGGCCATCGGGTCGGTGAACGGCATGCCCAGCTCGATCACGTCGGCGCCAGCAGCCGGCAGGCCCTTGAGGATCTGCAACGAGGCGTCGTAGCCCGGGTCGCCAGCTGTGACGAAGGt                >  1:101107/7‑149 (MQ=255)
                                                                                                                                       aGGGTGGCCAGCTGGATGGCCACGCCATCGGCCATCGGGTCGGTGAATGGCATGCCCAGCTCGATCACGTCGGCGCCAGCAGCCGGCAGGCCCTTGAGGATCTGCAACGAGGCGTCGTAGCCCGGGTCGCCAGCTGTGACGAAGGTGAc             >  5:276621/1‑149 (MQ=255)
                                                                                                                                         ggTGGCCAGCTGGATGGCCACGCCATCGGCCATCGGGTCGGTGAACGGCATGCCCAGCTCGATCACGTCGGCGCCAGCAGCCGGCAGGCCCTTGAGGATCTGCAACGAGGCGTCGTAGCCCGGGTCGCCAGCTGTGACGAAGGTGACCa           <  5:176375/149‑1 (MQ=255)
                                                                                                                                             gCCAGCTGGATGGCCACGCCATCGGCCATCGGGTCGGTGAATGGCATGCCCAGCTCGATCACGTCGGCGCCAGCAGCCGGCAGGCCCTTGAGGATCTGCAACGAGGCGTCGTAGCCCGGGTCGCCAGCTGTGACGAAGGTGACCAGTGc       <  6:276621/149‑1 (MQ=255)
                                                                                                                                             gCCAGCTGGATGGCCACGCCATCGGCCATCGGGTCGGTGAACGGCATGCCCAGCTCGATCACGTCGGCGCCAGCAGCCGGCAGGCCCTTGAGGATCTGCAACGAGGCGTCGTAGCCCGGGTCGCCAGCTGTGACGAAGGTGACCAGTGc       >  4:131507/1‑149 (MQ=255)
                                                                                                                                              ccAGCTGGATGGCCACGCCATCGGCCATCGGGTCGGTGAACGGCATGCCCAGCTCGATCACGTCGGCGCCAGCAGCCGGCAGGCCCTTGAGGATCTGCAACGAGGCGTCGTAGCCCGGGTCGCCAGCTGTGACGAAGGTGACCAGTGcc      >  1:128563/1‑149 (MQ=255)
                                                                                                                                              ccAGCTGGATGGCCACGCCATCGGCCATCGGGTCGGTGAACGGCATGCCCAGCTCGATCACGTCGGCGCCAGCAGCCGGCAGGCCCTTGAGGATCTGCAACGAGGAGTCGTAGCCCGGGTCGCCAGCTGTGACGAAGGTGACCAGTgat      >  7:94952/1‑147 (MQ=255)
                                                                                                                                                aGCTGGATGGCCACGCCATCGGCCATCGGGTCGGTGAACGGCATGCCCAGCTCGATCACGTCGGCGCCAGCAGCCGGCAGGCCCTTGAGGATCTGCAACGAGGCGTCGTAGCCCGGGTCGCCAGCTGTGACGAAGGTGACCAGTGCCGa    <  2:215968/149‑1 (MQ=255)
                                                                                                                                                 gCTGGATGGCCACGCCATCGGCCATCGGGTCGGTGAACGGCATGCCCAGCTCGATCACGTCGGCGCCAGCAGCCGGCAGGCCCTTGAGGATCTGCAACGAGGCGTCGTAGCCCGGGTCGCCAGCTGTGACGAAGGTGACCAGTGCCGAg   >  7:326472/1‑149 (MQ=255)
                                                                                                                                                  cTGGATGGCCACGCCATCGGCCATCGGGTCGGTGAACGGCATGCCCAGCTCGATCACGTCGGCGCCAGCAGCCGGCAGGCCCTTGAGGATCTGCAACGAGGCGTCGTAGCCCGGGTCGCCAGCTGTGACGAAGGTGACCAGTGCCGAGc  >  1:198335/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GTGTCCACGCCAAAGCGGTGGATAGGGTTGTAGTAGCCCATCAGCACGATCGGCGTGGTCTGGTTATCCACACGGA‑‑ATTCGCGAACCATCTGCAGGGTCTTGACCAAGGTCTGGCCAGCTTCAAGGGCGCGCAGGGTGGCCAGCTGGATGGCCACGCCATCGGCCATCGGGTCGGTGAACGGCATGCCCAGCTCGATCACGTCGGCGCCAGCAGCCGGCAGGCCCTTGAGGATCTGCAACGAGGCGTCGTAGCCCGGGTCGCCAGCTGTGACGAAGGTGACCAGTGCCGAGC  >  NC_002947/92042‑92333

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