Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 92,184 G→T 5.1% L66M (CTG→ATG)  trpA ← tryptophan synthase subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00294792,1840GT 5.1% 319.5 / 5.3 97L66M (CTG→ATG) trpAtryptophan synthase subunit alpha
Reads supporting (aligned to +/- strand):  ref base G (49/43);  new base T (2/3);  total (51/46)
Fisher's exact test for biased strand distribution p-value = 6.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.15e-01

GTGTCCACGCCAAAGCGGTGGATAGGGTTGTAGTAGCCCATCAGCACGATCGGCGTGGTCTGGTTATCCACACGGA‑‑ATTCGCGAACCATCTGCAGGGTCTTGACCAAGGTCTGGCCAGCTTCAAGGGCGCGCAGGGTGGCCAGCTGGATGGCCACGCCATCGGCCATCGGGTCGGTGAACGGCATGCCCAGCTCGATCACGTCGGCGCCAGCAGCCGGCAGGCCCTTGAGGATCTGCAACGAGGCGTCGTAGCCCGGGTCGCCAGCTGTGACGAAGGTGACCAGTGCCGAG  >  NC_002947/92042‑92332
                                                                                                                                                 |                                                                                                                                                     
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                                                                                                                        <  2:330815/149‑1 (MQ=255)
                         ggTTGTAGTAGCCCATCAGCACGATCGGCGTGGTCTGGTTATCCACACGGA‑‑ATTCGCGAACCATCTGCAGGGTCTTGACCAAGGTCTGGCCAGCTTCAAGGGCGCGCAGGGTGGCCAGCTGGATGGCCACGCCATCGGCCATCGGGTCg                                                                                                                        >  3:354472/1‑149 (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                                                                                                                              <  1:158255/117‑1 (MQ=255)
                                                   ggCGTGGTCTGGTTATCCACACGGA‑‑ATTCGCGAACCATCTGCAGGGTCTTGACCAAGGTCTGGCCAGCTTCAAGGGCGCGCAGGGTGGCCAGCTGGATGGCCACgccatcggccatc                                                                                                                              >  2:158255/1‑117 (MQ=255)
                                                   ggCGTGGTCTGGTTATCCACACGGA‑‑ATTCGCGAACCATCTGCAGGGTCTTGACCAAGGTCTGGCCAGCTTCAAGGGCGCGCAGGCTGGCCATCCAGCTGGCCAccct                                                                                                                                        <  7:66070/107‑4 (MQ=255)
                                                    gCGTGGTCTGGTTATCCACACGGA‑‑ATTCGCGAACCATCTGCAGGGTCTTGACCAAGGTCTGGCCAGCTTCAAGGGCGCGCAGGGTGGCCAGCTGGATGGCCACGCCATCGGCCATCGGGTCGGTGAACGGCATGCCCAGCTCGATCACg                                                                                            >  8:311003/1‑149 (MQ=255)
                                                     cGTGGTCTGGTTATCCACACGGA‑‑ATTCGCGAACCATCTGCAGGGTCTTGACCAAGGTCTGGCCAGCTTCAAGGGCGCGCAGGGTGGCCAGCTGGATGGCCACGCCATCGGCCATCGGGTCGGTGAACGGc                                                                                                              >  6:421216/1‑130 (MQ=255)
                                                     cGTGGTCTGGTTATCCACACGGA‑‑ATTCGCGAACCATCTGCAGGGTCTTGACCAAGGTCTGGCCAGCTTCAAGGGCGCGCAGGGTGGCCAGCTGGATGGCCACGCCATCGGCCATCGGGTCGGTGAACGGc                                                                                                              <  5:421216/130‑1 (MQ=255)
                                                      gTGGTCTGGTTATCCACACGGA‑‑ATTCGCGAACCATCTGCAGGGTCTTGACCAAGGTCTGGCCAGCTTCAAGGGCGCGCAGGGTGGCCAGCTGGATGGCCACGCCATCGGCCATCGGGTCGGTGAACGGCATGCCCAGCTCGATCACGTc                                                                                          <  2:76547/149‑1 (MQ=255)
                                                       tGGTCTGGTTATCCACACGGA‑‑ATTCGCGAACCATCTGCAGGGTCTTGACCAAGGTCTGGCCAGCTTCAAGGGCGCGCAGGGTGGCCAGCTGGATGGCCACGCCATCGGCCATCGGGTCGGTGAACGGCATGCCCAGCTCGATCACGTCg                                                                                         <  8:248922/149‑1 (MQ=255)
                                                       tGGTCTGGTTATCCACACGGA‑‑ATTCGCGAACCATCTGCAGGGTCTTGACCAAGGTCTGGCCAGCTTCAAGGGCGCGCAGGGTGGCCAGCTGGATGGCCACGCCATCGGCCATCGGGTCGGTGAACGGCATGCCCAGCTCGATCACGTCg                                                                                         >  7:49919/1‑149 (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                                                                              <  4:173040/135‑1 (MQ=255)
                                                                                   gcgAACCATCTGCAGGGTCTTGACCAAGGTCTGGCCAGCTTCAAGGGCGCGCAGGGTGGCCAGCTGGATGGCCACGCCATCGGCCATCGGGTCGGTGAACGGCATGCCCAGCTCGATCACGTCGGCGCCagcagc                                                                              >  3:173040/1‑135 (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                                        >  7:25263/1‑149 (MQ=255)
                                                                                                           cAAGGTCTGGCCAGCTTCAAGGGCGCGCAGGGTGGCCAGCTGGATGGCCACGCCATCGGCCATCGGGTCGGTGAACGGCATGCCCAGCTCGATCACGTCGGCGCCAGCAGCCGGCAGGCCCTTGAGGATCTGCAACGAGGCGTCGTAGc                                        >  1:100961/1‑149 (MQ=255)
                                                                                                           cAAGGTCTGGCCAGCTTCAAGGGCGCGCAGGGTGGCCAGCTGGATGGCCACGCCATCGGCCATCGGGTCGGTGAACGGCATGCCCAGCTCGATCACGTCGGCGCCAGCAGCCGGCAGGCCCTTGAGGATCTGCAACGAGGCGTCGTAGc                                        <  5:301054/149‑1 (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)
                                                                                                                                                 |                                                                                                                                                     
GTGTCCACGCCAAAGCGGTGGATAGGGTTGTAGTAGCCCATCAGCACGATCGGCGTGGTCTGGTTATCCACACGGA‑‑ATTCGCGAACCATCTGCAGGGTCTTGACCAAGGTCTGGCCAGCTTCAAGGGCGCGCAGGGTGGCCAGCTGGATGGCCACGCCATCGGCCATCGGGTCGGTGAACGGCATGCCCAGCTCGATCACGTCGGCGCCAGCAGCCGGCAGGCCCTTGAGGATCTGCAACGAGGCGTCGTAGCCCGGGTCGCCAGCTGTGACGAAGGTGACCAGTGCCGAG  >  NC_002947/92042‑92332

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: