Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,558,505 A→T 7.2% S250T (TCC→ACC)  purA ← adenylosuccinate synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,558,5050AT 7.2% 210.3 / 5.6 69S250T (TCC→ACC) purAadenylosuccinate synthetase
Reads supporting (aligned to +/- strand):  ref base A (31/33);  new base T (2/3);  total (33/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

CGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAGTGGTGTTGGAGCTGGTGACGTACGGGTAGGTGCCGTGGTCGATGTCCAGCAGCGAGCCCTGGGCGCCTTCGAACATGATGTCCTTGCCGGCGCGACGCAGGTTGTGCAGCTCGGCGGTGACGTCGAGC  >  NC_002947/5558358‑5558647
                                                                                                                                                   |                                                                                                                                              
cGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGa                                                                                                                                                >  3:73043/1‑148 (MQ=255)
  aaCTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGaaa                                                                                                                                              >  4:232794/1‑148 (MQ=255)
    accttGGCCAAGCGTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATac                                                                                                                                           >  6:148145/5‑149 (MQ=255)
          gCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGTATTTCCACCgg                                                                                                                                     <  5:14297/149‑1 (MQ=255)
                 ttGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAgtggt                                                                                                                              >  5:295516/1‑149 (MQ=255)
                      cAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGTATTTCCACCggtt                                                                                                                                   >  4:172960/1‑137 (MQ=255)
                      cAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGTATTTCCACCggtt                                                                                                                                   <  3:172960/139‑3 (MQ=255)
                      cAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCAATTGGGGTAGGCCTTGGTGATACCGAGCTTTTTGTGGAGTGGCAACGGCCCGCCCCCGCAACCGGTGGAAATACAACCGGAAGTGTTGTTgg                                                                                                                         >  6:27184/1‑149 (MQ=255)
                               aCCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAGTGGTGTTGGAGCTGGTGa                                                                                                                >  1:7618/1‑149 (MQ=255)
                                ccGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAGTGGTGTTGGAGCTGGTGAc                                                                                                               <  6:191236/149‑1 (MQ=255)
                                   ggctcGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAGTGGTGTTGGAGCTGGTGACGTa                                                                                                            <  6:40777/147‑1 (MQ=255)
                                   gTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAGTGGTGTTGGAGCTGGTGACGTa                                                                                                            >  5:197953/1‑149 (MQ=255)
                                    tctcGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGTATTTCCACCggttcc                                                                                                                                 >  8:220811/1‑123 (MQ=255)
                                    tctcGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGTATTTCCACCggttcc                                                                                                                                 <  7:220811/127‑5 (MQ=255)
                                      tcgtcgAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAGTGGTGTTGGAGCTGGTGACGTAc                                                                                                           >  4:67501/1‑147 (MQ=255)
                                      tcgtcgAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAGTGGTGTTGGAGCTGGTGACGTAc                                                                                                           <  3:67501/147‑1 (MQ=255)
                                        gtcgAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAGTGGTGTTGGAGCTGGTGACGTACggg                                                                                                        <  8:21138/148‑1 (MQ=255)
                                        gtcgAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAGTGGTGTTGGAGCTGGTGACGTACgg                                                                                                         <  6:295516/147‑1 (MQ=255)
                                           gAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAGTGGTGTTGGAGCTGGTGACGTACGGGTAgg                                                                                                    <  5:148145/149‑1 (MQ=255)
                                            aaCAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGATGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAGTGGTGTTGGAGCTGGTGACGTAc                                                                                                           <  7:117473/141‑1 (MQ=255)
                                            aaCAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAGTGGTGTTGGAGCTGGTGACGTACGGGTAGGt                                                                                                   >  3:293400/1‑149 (MQ=255)
                                            aaCAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACAACCGGCAGTGGTGTTGGAGCTGGTGACGTAc                                                                                                           >  8:117473/1‑141 (MQ=255)
                                            aaCAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGAGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAGTGGTGTTGGAGCTGGTGACGTACGGGTAGGt                                                                                                   >  1:216195/1‑149 (MQ=255)
                                                 ttCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAGTGGTGTTGGAGCTGGTGACGTACGGGTAGGTGCCGt                                                                                              >  8:23032/1‑149 (MQ=255)
                                                    ggtcggGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAGTGGTGTTGGAGCTGGTGACGTACGGGTAGGTGCCGTGGt                                                                                           <  4:58764/149‑1 (MQ=255)
                                                           aaCGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAGTGGTGTTGGAGCTGGTGACGTACTGGTa                                                                                                      <  7:166778/131‑1 (MQ=255)
                                                           aaCGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAGTGGTGTTGGAGCTGGTGACGTACGGGTa                                                                                                      >  8:166778/1‑131 (MQ=255)
                                                                   ggAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAGTGGTGTTGGAGCTGGTGACGTACGGGTAGGTGCCGTGGTCGATGTCCagcagc                                                                             >  7:149863/1‑148 (MQ=255)
                                                                       cccACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAGTGGTGTTGGAGCTGGTg                                                                                                                 <  1:208245/108‑1 (MQ=255)
                                                                       cccACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAGTGGTGTTGGAGCTGGTg                                                                                                                 >  2:208245/1‑108 (MQ=255)
                                                                        ccACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAGTGGTGTTGGAGCTGGTGACGTACGGGTAGGTGCCGTGGTCGATGTCCAGCAGCGAGccc                                                                       >  2:191776/1‑149 (MQ=255)
                                                                          aCGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAGTGGTGTTGGAGCTGGTGACGTACGGGTAGGTGCCGTGGTCGa                                                                                        >  4:176749/1‑130 (MQ=255)
                                                                          aCGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAGTGGTGTTGGAGCTGGTGACGTACGGGTAGGTGCCGTGGTCGa                                                                                        <  3:176749/130‑1 (MQ=255)
                                                                          aCGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAGTGGTGTTGGAGCTGGTGACGTACGGGTAGGTGCCGTGGTCGATGTCCAGCAGCGAGCCCTg                                                                     >  5:77947/1‑149 (MQ=255)
                                                                                   gtgtAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAGTGGTGTTGGAGCTGGTGACGTACGGGTAGGTGCCGTGGTCGATGTCCAGCAGCGAGCCCTGGGCGCCtt                                                             <  7:23032/148‑1 (MQ=255)
                                                                                   gtgtAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAGTGGTGTTGGAGCTGGTGACGTACGGGTAGGTGCCGTGGTCGATGTCCAGCAGCGAGCCCTGGGCGCCTTc                                                            >  6:272246/1‑149 (MQ=255)
                                                                                    tgtAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAGTGGTGTTGGAGCTGGTGACGTACGGGTAGGTGCCGTGGTCGATGTCCAGCAGCGAGCCCTGGGCGCCTTc                                                            >  2:306828/1‑148 (MQ=255)
                                                                                    tgtAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCCGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAGTGGTGTTGGAGCTGGTGACGTACGGGTAGGTGCCGTGGTCGATGTCCAGCAGCGAGCCCTGGGCGCCTTCg                                                           <  5:272246/149‑1 (MQ=255)
                                                                                      tAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAGTGGTGTTGGAGCTGGTGACGTACGGGTAGGTGCCGTGGTCGATGTCcag                                                                                 >  2:37680/1‑125 (MQ=255)
                                                                                      tAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAGTGGTGTTGGAGCTGGTGACGTACGGGTAGGTGCCGTGGTCGATGTCcag                                                                                 <  1:37680/125‑1 (MQ=255)
                                                                                        ggCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAGTGGTGTTGGAGCTGGTGACGTACGGGTAGGTGCCGTGGTCGATGTCCAGCAGCGAGCCCTGGGCGCCTTCGAAc                                                        <  7:103169/148‑1 (MQ=255)
                                                                                          ccTTGGTGATACCGAGGATGTAGGCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAGTGGTGTTGGAGCTGGTGACGTACGGGTAGGTGCCGTGGTCGATGTCCAGCAGCGAGCCCTGGGCGCCTTCGAACatg                                                     <  6:50978/149‑1 (MQ=255)
                                                                                           cTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAGTGGTGTTGGAGCTGGTGACGTACGGGTAGGTGCCGTGGTCGATGTCCAGCAGCGAGCCCTGGGCGCCTTCGAACatga                                                    >  7:158605/1‑149 (MQ=255)
                                                                                               gTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAGTGGTGTTGGAGCTGGTGACGTACGGGTAGGTGCCGTGGTCGATGTCCAGCAGCGAGCCCTGGGCGCCTTCGAACATGATGt                                                 <  7:11684/148‑1 (MQ=255)
                                                                                                  aTACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAGTGGTGTTGGAGCTGGTGACGTACGGGTAGGTGCCGTGGTCGATGTCCAGCAGCGAGCCCTGGGCGCCTTCGAACATGATGTCCtt                                             <  4:293400/149‑1 (MQ=255)
                                                                                                  aTACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAGTGGTGTTGGAGCTGGTGACGTACGGGTAGGTGCCGTGGTCGATGTCCAGCAGCGAGCCCTGGGCGCCTTCGAACATGATGTCCtt                                             <  3:232794/149‑1 (MQ=255)
                                                                                                  aTACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAGTGGTGTTGGAGCTGGTGACGTACGGGTAGGTGCCGTGGTCGATGTCCAGCAGCGAGCCCTGGGCGCCTTCGAACATGATGTCCtt                                             <  2:216195/149‑1 (MQ=255)
                                                                                                      cGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAGTGGTGTTGGAGCTGGTGACGTACGGGTAGGTGCCGTGGTCGATGTCCAGCAGCGAGCCCTGGGCGCCTTCGAACATGATGTCCTTGCCg                                         >  6:102596/1‑149 (MQ=255)
                                                                                                       gAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAGTGGTGTTGGAGCTGGTGACGTACGGGTAGGTGCCGTGGTCGATGTCCAGCAGCGAGCCCTGGGCGCCTTCGaa                                                         <  4:170186/132‑1 (MQ=255)
                                                                                                       gAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAGTGGTGTTGGAGCTGGTGACGTACGGGTAGGTGCCGTGGTCGATGTCCAGCAGCGAGCCCTGGGCGCCTTCGaa                                                         >  3:170186/1‑132 (MQ=255)
                                                                                                          gATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAGTGGTGTTGGAGCTGGTGACGTACGGGTAGGTGCCGTGGTCGATGTCCAGCAGCGAGCCCTGGGCGCCTTCGAACATGATGTCCTTGCCGgcgc                                     <  6:77947/149‑1 (MQ=255)
                                                                                                            tGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAGTGGTGTTGGAGCTGGTGACGTACGGGTAGGTGCCGTGGTCGATGTCCAGCAGCGAGCCCTGGGCGCCTTCGAACATGATGTCCTTGCCGGCGCGa                                   <  1:191776/149‑1 (MQ=255)
                                                                                                             gTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAGTGGTGTTGGAGCTGGTGACGTACGGGTAGGTGCCGTGGTCGATGTCCAGCAGCGAGCCCTGGGCGCCTTCGAACATGATGTCCTTGCCGGCGCGAc                                  >  6:187425/1‑149 (MQ=255)
                                                                                                               aGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAGTGGTGTTGGAGCTGGTGACGTACGGGTAGGTGCCGTGGTCGATGTCCAGCAGCGAGCCCTGGGCGCCTTCGAACATGATGTCCTTGc                                           <  1:130573/138‑1 (MQ=255)
                                                                                                               aGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAGTGGTGTTGGAGCTGGTGACGTACGGGTAGGTGCCGTGGTCGATGTCCAGCAGCGAGCCCTGGGCGCCTTCGAACATGATGTCCTTGc                                           >  2:130573/1‑138 (MQ=255)
                                                                                                                   gAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAGTGGTGTTGGAGCTGGTGACGTACGGGTAGGTGCCGTGGTCGATGTCCAGCAGCGAGCCCTGGGCGCCTTCGAACATGATGTCCtt                                             <  6:321899/132‑1 (MQ=255)
                                                                                                                   gAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAGTGGTGTTGGAGCTGGTGACGTACGGGTAGGTGCCGTGGTCGATGTCCAGCAGCGAGCCCTGGGCGCCTTCGAACATGATGTCCtt                                             >  5:321899/1‑132 (MQ=255)
                                                                                                                   gAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAGTGGTGTTGGAGCTGGTGACGTACGGGTAGGTGCCGTGGTCGATGTCCAGCAGCGAGCCCTGGGCGCCTTCGAACATGATGTCCTTGCCGGCGCGACGCAgg                             <  8:158605/148‑1 (MQ=255)
                                                                                                                      gTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAGTGGTGTTGGAGCTGGTGACGTACGGGTAGGTGCCGTGGTCGATGTCCAGCAGCGAGCCCTGGGCGCCTTCGAACATGATGTCCTTGCCGGCGCGACGCAGGTtg                          <  2:7618/148‑1 (MQ=255)
                                                                                                                        aCATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAGTGGTGTTGGAGCTGGTGACGTACGGGTAGGTGCCGTGGTCGATGTCCAGCAGCGAGCCCTGGGCGCCTTCGAAc                                                        >  4:183226/1‑116 (MQ=255)
                                                                                                                        aCATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAGTGGTGTTGGAGCTGGTGACGTACGGGTAGGTGCCGTGGTCGATGTCCAGCAGCGAGCCCTGGGCGCCTTCGAAc                                                        <  3:183226/116‑1 (MQ=255)
                                                                                                                         cATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAGTGGTGTTGGAGCTGGTGACGTACGGGTAGGTGCCGTGGTCGATGTCCAGCAGCGAGCCCTGGGCGCCTTCGAACATGATGTCCTTGCCGGCGCGACGCAGGTTGTGCa                      <  6:197953/149‑1 (MQ=255)
                                                                                                                                 cGACGCCGGAACCGGTGGAAATACCACCGGCAGTGGTGTTGGAGCTGGTGACGTACGGGTAGGTGCCGTGGTCGATGTCCAGCAGCGAGCCCTGGGCGCCTTCGAACATGATGTCCTTGCCGGCGCGACGCCGGTTGTGCAGCTcggcg              >  8:295738/1‑149 (MQ=255)
                                                                                                                                   aCGCCGGAACCGGTGGAAATACCACCGGCAGTGGTGTTGGAGCTGGTGACGTACGGGTAGGTGCCGTGGTCGATGTCCAGCAGCGAGCCCTGGGCGCCTTCGAACATGATGTCCTTGCCGGCGCGACGCAGGTTGTGCAGCTCGGCGGt            <  8:149863/149‑1 (MQ=255)
                                                                                                                                      ccGGAACCGGTGGAAATACCACCGGCAGTGGTGTTGGAGCTGGTGACGTACGGGTAGGTGCCGTGGTCGATGTCCAGCAGCGAGCCCTGGGCGCCTTCGAACATGATGTCCTTGCCGGCGCGACGCAGGTTGTGCAGCTCGGCGGTGAc         >  6:23028/1‑149 (MQ=255)
                                                                                                                                       cGGAACCGGTGGAAATACCACCGGCAGTGGTGTTGGAGCTGGTGACGTACGGGTAGGTGCCGTGGTCGATGTCCAGCAGCGAGCCCTGGGCGCCTTCGAACATGATGTCCTTGCCGGCGCGACGCAGGTTGTGCAGCTCGGCGGTGACg        >  1:303472/1‑149 (MQ=255)
                                                                                                                                        ggAACCGGTGGAAATACCACCGGCAGTGGTGTTGGAGCTGGTGACGTACGGGTAGGTGCCGTGGTCGATGTCCAGCAGCGAGCCCTGGGCGCCTTCGAACATGATGTCCTTGCCGGCGCGa                                   >  8:130907/1‑121 (MQ=255)
                                                                                                                                        ggAACCGGTGGAAATACCACCGGCAGTGGTGTTGGAGCTGGTGACGTACGGGTAGGTGCCGTGGTCGATGTCCAGCAGCGAGCCCTGGGCGCCTTCGAACATGATGTCCTTGCCGGCGCGa                                   <  7:130907/121‑1 (MQ=255)
                                                                                                                                         gAACCGGTGGAAATACCACCGGCAGTGGTGTTGGAGCTGGTGACGTACGGGTAGGTGCCGTGGTCGATGTCAAGCAGCGAGCCCTGGGCGCATTCGAACATGATGTGCTTGCCGGCGCGACGCAGGTTGTGCAGCTCGGCGGTGACGTc      >  5:323559/1‑149 (MQ=255)
                                                                                                                                           aCCGGTGGAAATACCACCGGCAGTGGTGTTGGAGCTGGTGACGTACGGGTAGGTGCCGTGGTCGATGTCCAGCAGCGAGCCCTGGGCGCCTTCGAACATGATGTCCTTGCCGGCGCGACGCAGGTTGTGCAGCTCGGCGGTGACGTCGa    <  1:306828/149‑1 (MQ=255)
                                                                                                                                             cGGTGGAAATACCACCGGCAGTGGTGTTGGAGCTGGTGACGTACGGGTAGGTGCCGTGGTCGATGTCCAGCAGCGAGCCCTGGGCGCCTTCGAACATGATGTCCTTGCCGGCGCGACGCAGGTTGTGCAGCTCGGCGGTGACGTCGAGc  <  5:187425/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                              
CGAACTCGTGGCCACGCTTGGCCAAGGTAGCACCGGTCTCGTCGAACAGTTCGGTCGGGAACGGACCGGAACCCACGCGAGTGGTGTAGGCCTTGGTGATACCGAGGATGTAGTCGAGGTACATCGGGCCGACGCCGGAACCGGTGGAAATACCACCGGCAGTGGTGTTGGAGCTGGTGACGTACGGGTAGGTGCCGTGGTCGATGTCCAGCAGCGAGCCCTGGGCGCCTTCGAACATGATGTCCTTGCCGGCGCGACGCAGGTTGTGCAGCTCGGCGGTGACGTCGAGC  >  NC_002947/5558358‑5558647

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: