Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,656,167 G→C 6.9% A8A (GCC→GCG cysH ← phosphoadenosine phosphosulfate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,656,1670GC 6.9% 219.2 / 7.0 72A8A (GCC→GCGcysHphosphoadenosine phosphosulfate reductase
Reads supporting (aligned to +/- strand):  ref base G (34/33);  new base C (3/2);  total (37/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.78e-01

CCTGCTTGTTCAGCTTCCAGGCCATGTCGACCAGCACCACATCCTCGGCGCCGCTGAAGGAGATCCACAGGTCATCACCGAAGTGCTCGAAGGCAAGCTTGAGGATGTCCTGCGGGGACTTGTTGGCATAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTGCGCTCGTAGTAGGGGGATGGTAACAAAAAAATCGTGGTCTAGATGGGCTCACAAGCCCTGCAGGGTCTGCATCAGCACCCGCACTTTGGCGAT  >  NC_002947/2656022‑2656315
                                                                                                                                                  |                                                                                                                                                    
ccTGCTTGTTCAGCTTCCAGGCCATGTCGACCAGCACCACATCCTCGGCGCCGCTGAAGGAGATCCACAGGTCATCACCGAAGTGCTCGAAGGCAAGCTTGAGGATGTCCTGCGGGGACTTGTTGGCATAGGTCGCGGCCAGggcggc                                                                                                                                                    <  5:111192/148‑1 (MQ=255)
     ttGTTCAGCTTCCAGGCCATGTCGACCAGCACCACATCCTCGGCGCCGCTGAAGGAGATCCACAGGTCATCACCGAAGTGCTCGAAGGCAAGCTTGAGGATGTCCTGCGGGGACTTGTTGGCATAGGTCGCGGCCAGGGCGGCGACGTc                                                                                                                                              >  4:116661/1‑149 (MQ=255)
              ttCCAGGCCATGTCGACCAGCACCACATCCTCGGCGCCGCTGAAGGAGATCCACAGGTCATCACCGAAGTGCTCGAAGGCAAGCTTGAGGATGTCCTGCGGGGACTTGTTGGCATAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTg                                                                                                                                     >  5:28019/1‑149 (MQ=255)
                ccAGGCCATGTCGACCAGCACCACATCCTCGGCGCCGCTGAAGGAGATCCACAGGTCATCACCGAAGTGCTCGAAGGCAAGCTTGAGGATGTCCTGCGGGGACTTGTTGGCATAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTgg                                                                                                                                    <  3:116661/148‑1 (MQ=255)
                  aGGCCATGTCGACCAGCACCACATCCTCGGCGCCGCTGAAGGAGATCCACAGGTCATCACCGAAGTGCTCGAAGGCAAGCTTGAGGATGTCCTGCGGGGACTTGTTGGCATAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGtt                                                                                                                                      <  5:290936/144‑1 (MQ=255)
                  aGGCCATGTCGACCAGCACCACATCCTCGGCGCCGCTGAAGGAGATCCACAGGTCATCACCGAAGTGCTCGAAGGCAAGCTTGAGGATGTCCTGCGGGGACTTGTTGGCATAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCt                                                                                                                                  >  5:128582/1‑148 (MQ=255)
                  aGGCCATGTCGACCAGCACCACATCCTCGGCGCCGCTGAAGGAGATCCACAGGTCATCACCGAAGTGCTCGAAGGCAAGCTTGAGGATGTCCTGCGGGGACTTGTTGGCATAGGTCGCCGCCCTGGCCGCGAc                                                                                                                                                 <  6:289900/133‑1 (MQ=255)
                  aGGCCATGTCGACCAGCACCACATCCTCGGCGCCGCTGAAGGAGATCCACAGGTCATCACCGAAGTGCTCGAAGGCAAGCTTGAGGATGTCCTGCGGGGACTTGTTGGCATAGGTCGCCGCCCTGGCCGCGAc                                                                                                                                                 >  5:289900/1‑133 (MQ=255)
                  aGGCCATGTCGACCAGCACCACATCCTCGGCGCCGCTGAAGGAGATCCACAGGTCATCACCGAAGTGCTCGAAGGCAAGCTAGAGGATGTCCTGCGGGGACTTGTTGGCATAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGtt                                                                                                                                      >  6:290936/1‑144 (MQ=255)
                      cATGTCGACCAGCACCACATCCTCGGCGCCGCTGAAGGAGATCCACAGGTCATCACCGAAGTGCTCGAAGGCAAGCTTGAGGATGTCCTGCGGGGACTTGTTGGCACTAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCtcatc                                                                                                                              <  1:68549/149‑1 (MQ=255)
                            gACCAGCACCACATCCTCGGCGCCGCTGAAGGAGATCCACAGGTCATCACCGAAGTGCTCGAAGGCAAGCTTGAGGATGTCCTGCGGGGACTTGTTGGCATAGGTCGCGGCCAGggcggc                                                                                                                                                    <  4:320695/120‑1 (MQ=255)
                                    ccacATCCTCGGCGCCGCTGAAGGAGATCCACAGGTCATCACCGAAGTGCTCGAAGGCAAGCTTGAGGATGTCCTGCGGGGACTTGTTGGCATAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCt                                                                                                               >  2:119080/1‑149 (MQ=255)
                                     cacaTCCTCGGCGCCGCTGAAGGAGATCCACAGGTCATCACCGAAGTGCTCGAAGGCAAGCTTGAGGATGTCCTGCGGGGACTTGTTGGCATAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCtcatc                                                                                                                              >  4:131643/1‑133 (MQ=255)
                                     cacaTCCTCGGCGCCGCTGAAGGAGATCCACAGGTCATCACCGAAGTGCTCGAAGGCAAGCTTGAGGATGTCCTGCGGGGACTTGTTGGCATAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCtcatc                                                                                                                              <  3:131643/133‑1 (MQ=255)
                                     cacaTCCTCGGCGCCGCTGAAGGAGATCCACAGGTCATCACCGAAGTGCTCGAAGGCAAGCTTGAGGATGTCCTGCGGGGACTTGTTGGCATAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCTCATCAGGCTCTTTCCATCtt                                                                                                              >  6:6070/1‑149 (MQ=255)
                                     cacaTCCTCGGCGCCGCTGAAGGAGATCCACAGGTCATCACCGAAGTGCTCGAAGGCAAGCTTGAGGATGTCCTGCGGGGACTTGTTGGCATAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCtt                                                                                                              >  3:265903/1‑149 (MQ=255)
                                            tCGGCGCCGCTGAAGGAGATCCACAGGTCATCACCGAAGTGCTCGAAGGCAAGCTTGAGGATGTCCTGCGGGGACTTGTTGGCATAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCtgtg                                                                                                       >  8:241975/1‑149 (MQ=255)
                                                     cTGAAGGAGATCCACAGGTCATCACCGAAGTGCTCGAAGGCAAGCTTGAGGATGTCCTGCGGGGACTTGTTGGCATAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTgc                                                                                              >  5:316834/1‑149 (MQ=255)
                                                       gAAGGAGATCCACAGGTCATCACCGAAGTGCTCGAAGGCAAGCTTGAGGATGTCCTGCGGGGACTTGTTGGCATAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTgcgc                                                                                            >  3:11295/1‑149 (MQ=255)
                                                         aGGAGATCCACAGGTCATCACCGAAGTGCTCGAAGGCAAGCTTGAGGATGTCCTGCGGGGACTTGTTGGCATAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTGCGCTc                                                                                          <  6:128582/149‑1 (MQ=255)
                                                            agaTCCACAGGTCATCACCGAAGTGCTCGAAGGCAAGCTTGAGGATGTCCTGCGGGGACTTGTTGGCATAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTGCGCTCgta                                                                                       <  4:265903/149‑1 (MQ=255)
                                                                 cacaGGTCATCACCGAAGTGCTCGAAGGCAAGCTTGAGGATGTCCTGCGGGGACTTGTTGGCATAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTGCGCTCGTAGTAgg                                                                                  <  7:241975/149‑1 (MQ=255)
                                                                 cacaGGTCATCACCGAAGTGCTCGAAGGCAAGCTTGAGGATGTCCTGCGGGGACTTGTTGGCATAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTGCGCTCGTAGTAgg                                                                                  >  7:195477/1‑149 (MQ=255)
                                                                  acaGGTCATCACCGAAGTGCTCGAAGGCAAGCTTGAGGATGTCCTGCGGGGACTTGTTGGCATAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTGCGCTCGTAGTAggg                                                                                 >  1:303195/1‑149 (MQ=255)
                                                                      gTCATCACCGAAGTGCTCGAAGGCAAGCTTGAGGATGTCCTGCGGGGACTTGTTGGCATAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTGCGCTCGTAGTAGGGGGAt                                                                             >  5:255810/1‑149 (MQ=255)
                                                                         atcaCCGAAGTGCTCGAAGGCAAGCTTGAGGATGTCCTGCGGGGACTTGTTGGCATAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTGCGCTCGTAGTAggg                                                                                 >  7:276723/1‑142 (MQ=255)
                                                                         atcaCCGAAGTGCTCGAAGGCAAGCTTGAGGATGTCCTGCGGGGACTTGTTGGCATAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTGCGCTCGTAGTAggg                                                                                 <  8:276723/142‑1 (MQ=255)
                                                                             ccGAAGTGCTCGAAGGCAAGCTTGAGGATGTCCTGCGGGGACTTGTTGGCATAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTGCGCTCGTAGTAgg                                                                                  <  4:37333/137‑1 (MQ=255)
                                                                             ccGAAGTGCTCGAAGGCAAGCTTGAGGATGTCCTGCGGGGACTTGTTGGCATAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTGCGCTCGTAGTAgg                                                                                  >  3:37333/1‑137 (MQ=255)
                                                                                     cTCGAAGGCAAGCTTGAGGATGTCCTGCGGGGACTTGTTGGCATAGGTCGCGGCCAGGGTGGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTGCGCTCGTAGTAGGGGGATGGTAACAAAAAAATc                                                              >  5:306997/1‑149 (MQ=255)
                                                                                     cTCGAAGGCAAGCTTGAGGATGTCCTGCGGGGACTTGTTGGCATAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTGCGCTCGTAGTAGGGGGATGGTAACAAAAAAAt                                                               <  1:119080/148‑1 (MQ=255)
                                                                                     cTCGAAGGCAAGCTTGAGGATGTCCTGCGGGGACTTGTTGGCATAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTGCGCTCGTAGTAGGGGGATGGTAACAAAAAAATc                                                              >  4:223729/1‑149 (MQ=255)
                                                                                      tCGAAGGCAAGCTTGAGGATGTCCTGCGGGGACTTGTTGGCATAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTGCGCTCGTAGTAGGGGGATGGTAACAAAAAAATCg                                                             >  8:210765/1‑149 (MQ=255)
                                                                                      tCGAAGGCAAGCTTGAGGATGTCCTGCGGGGACTTGTTGGCATAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTGCGCTCGTAGTAGGGGGATGGTAACAAAAAAATCg                                                             >  6:349375/1‑149 (MQ=255)
                                                                                        gAAGGCTATCTTGAGGATGTCCTGCGGGGACTTGTTGTCATAGGTCGCGGCCAGGGCGGCGACGTAGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTGCGCTCGTAGTAGGGGGATGGTAACAAAAAAATCGt                                                            <  5:6070/148‑1 (MQ=255)
                                                                                        gAAGGCAAGCTTGAGGATGTCCTGCGGGGACTTGTTGGCATAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTGCGCTCGTAGTAGGGGGATGGTAACAAAAAAATCGTg                                                           >  8:64538/1‑149 (MQ=255)
                                                                                             cAAGCTTGAGGATGTCCTGCGGGGACTTGTTGGCATAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTGCGCTCGTAGTAGGGGGATGGTAACAAAAAAATCGTGGTCTa                                                      <  8:195477/149‑1 (MQ=255)
                                                                                                gCTTGAGGATGTCCTGCGGGGACTTGTTGGCATAGGTCGCGGCCAGGGTGGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTGCGCTCGTAGTAGGGGGATGGTAACAAAAAAATCGTGGTCTAGAt                                                   <  6:306997/149‑1 (MQ=255)
                                                                                                gCTTGAGGATGTCCTGCGGGGACTTGTTGGCATAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTGCGCTCGTAGTAGGGGGATGGTAACAAAAAAATCGTGGTCTAGAt                                                   >  1:174215/1‑149 (MQ=255)
                                                                                                gCTTGAGGATGTCCTGCGGGGACTTGTTGGCATAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTGCGCTCGTAGTAGGGGGATGGTAACAAAAAAATCGTGGTCTAGAt                                                   <  5:349375/149‑1 (MQ=255)
                                                                                                  ttGAGGATGTCCTGCGGGGACTTGTTGGCATAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTGCGCTCGTAGTAGGGGGGTGGTAACAAAAAAATCGTGGTCTAGATgg                                                 >  1:65101/1‑149 (MQ=255)
                                                                                                     aGGATGTCCTGCGGGGACTTGTTGGCATAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTGCGCTCGTAGTAGGGGGATGGTAACAAAAAAATCGTGGTCTAGATGGGCt                                              >  1:327126/1‑149 (MQ=255)
                                                                                                             cTGCGGGGACTTGTTGGCATAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTGCGCTCGTAGTAGGGGGATGGTAACAAAAAAATCGTGGTCTAGATGGGCTCACAAGc                                       <  6:316834/148‑1 (MQ=255)
                                                                                                               gCGGGGACTTGTTGGCATAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTGCGCTCGTAGTAGGGGGATGGTAACaaaaaa                                                                 >  8:70706/1‑120 (MQ=255)
                                                                                                               gCGGGGACTTGTTGGCATAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTGCGCTCGTAGTAGGGGGATGGTAACaaaaaa                                                                 <  7:70706/120‑1 (MQ=255)
                                                                                                               gCGGGGACTTGTTGGCATAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTGCGCTCGTAGTAGGGGGATGGTAACAAAAAAATCGTGGTCTAGATGGGCTCACAAGCCCt                                    >  2:41917/1‑149 (MQ=255)
                                                                                                                cGGGGACTTGTTGGCATAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTGCGCTCGTAGTAGGGGGATGGTAACAAAAAAATCGTGGTCTAGATGGGCTCACAAGccc                                     <  7:210765/147‑1 (MQ=255)
                                                                                                                cGGGGACTTGTTGGCATAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTGCGCTCGTAGTAGGGGGATGGTAACAAAAAAATCGTGGTCTAGATGGGCTCACAAGCCCt                                    <  6:255810/148‑1 (MQ=255)
                                                                                                                cGGGGACTTGTTGGCATAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTGCGCTCGTAGTAGGGGGATGGTAACAAAAAAATCGTGGTCTAGATGGGCTCACAAGCCCt                                    <  3:223729/148‑1 (MQ=255)
                                                                                                                  gggACTTGTTGGCATAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTGCGCTCGTAGTAGGGGGGTGGTAACAAAAAAATCGTGGTCTAGATGGGCTCACAAGCCCTGCa                                 <  2:65101/149‑1 (MQ=255)
                                                                                                                  gggACTTGTTGGCATAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTGCGCTCGTAGTAGGGGGATGGTAACAAAAAAATCGTGGTCTAGATGGGCTCACAAGCCCTg                                   >  5:112310/1‑147 (MQ=255)
                                                                                                                   ggACTTGTTGGCATAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTGCGCTCGTAGTAGGGGGATGGTAACAAAAAAATCGTGGTCTAGATGGGCTCACAAGCCCTGCAg                                <  1:41917/149‑1 (MQ=255)
                                                                                                                      cTTTTTGGCATAGTTCTCGGCGAGTGCGGCGACGGCGAAGGGTGGGCTCATCAGGCGGTTTCCATCTTGGCTGTGTCGCTCTGCGCGCGTAGTAGGGGGTTGGTAACAAAAAAATCGTGGTCTAGATGGGCTCACAAGCCCTGCAGGGt                             <  6:326762/149‑1 (MQ=255)
                                                                                                                         gttgGCATAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTGCGCTCGTTGTAGGGGGATGGTAACAAAAAAATCGTGGTCTAGATGGGCTCACAAGCCCTGCAGGGTCTg                          >  1:136997/1‑149 (MQ=255)
                                                                                                                         gttgGCATAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTGCGCTCGTAGTAGGGGGATGGTAACAAAAAAATCGTGGTCTAGATGGGCTCACAAGCCCTGCAGGGTCt                           <  2:303195/148‑1 (MQ=255)
                                                                                                                         gttgGCATAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTGCGCTCGTAGTAGGGGGATGGTAACAAAAAAATCGTGGTCTAGATGGGCTCACAAGCCCTGCAGGGTCTg                          >  1:305331/1‑149 (MQ=255)
                                                                                                                         gttgGCATAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTGCGCTCGTAGTAGGGGGATGGTAACAAAAAAATCGTGGTCTAGATGGGCTCACAAGCCCTGCAGGGTCTg                          >  4:323800/1‑149 (MQ=255)
                                                                                                                         gttgGCATAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTGCGCTCGTAGTAGGGGGATGGTAACAAAAAAATCGTGGTCTAGATGGGCTCACAAGCCCTGCAGGGTCTg                          >  5:20492/1‑149 (MQ=255)
                                                                                                                          ttgGCATAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTGCGCTCGTAGTAGGGGGATGGTAACAAAAAAATCGTGGTCTAGATGGGCTCACAAGCCCTGCAGGGTCTGc                         <  2:174215/149‑1 (MQ=255)
                                                                                                                          ttgGCATAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTGCGCTCGTAGTAGGGGGATGGTAACAAAAAAATCGTGGTCTAGATGGGCTCACAAGCCCTGCAGGGTCTGc                         <  7:64538/149‑1 (MQ=255)
                                                                                                                             gCATAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTGCGCTCGTAGTAGGGGGATGGTAACAAAAAAATCGTGGTCTAGATGGGCTCACAAGCCCTGCAGGGTCTGCATc                      <  4:11295/149‑1 (MQ=255)
                                                                                                                              cATAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTGCGCTCGTAGTAGGGGGATGGTAACAAAAAAATCGTGGTCTAGATGGGCTCACAAGCCCTGCAGGGTCTGCATCa                     <  6:20492/149‑1 (MQ=255)
                                                                                                                              cATAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTGCGCTCGTAGTAGGGGGATGGTAACAAAAAAATCGTGGTCTAGATGGGCTCACAAGCCCTGCAGGGTCTGCATCa                     <  2:305331/149‑1 (MQ=255)
                                                                                                                              cATAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTGCGCTCGTAGTAGGGGGATGGTAACAAAAAAATCGTGGTCTAGATGGGCTCACAAGCCCTGCAGGGTCTGCATCa                     <  3:323800/149‑1 (MQ=255)
                                                                                                                               tcgacgTCGCCGCCCTGGCCGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTGCGCTCGTAGTAGGGGGATGGTAACAAAAAAATCGTGGTCTAGATGGGCTCACAAGCCCTGCAGGGTCTGCAt                       >  4:326196/6‑147 (MQ=255)
                                                                                                                               aTAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTGCGCTCGTAGTAGGGGGATGGTAACAAAAAAATCGTGGTCTAGAt                                                   >  6:344538/1‑118 (MQ=255)
                                                                                                                               aTAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTGCGCTCGTAGTAGGGGGATGGTAACAAAAAAATCGTGGTCTAGAt                                                   <  5:344538/118‑1 (MQ=255)
                                                                                                                               aTAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTGCGCTCGTAGTAGGGGGATGGTAACAAAAAAATCGTGGTCTAGATgg                                                 <  6:347255/120‑1 (MQ=255)
                                                                                                                                cgacgTCGCCGCCCTGGCCGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTGCGCTCGTAGTAGGGGGATGGTAACAAAAAAATCGTGGTCTAGATGGGCTCACAAGCCCTGCAGGGTCTGCATCa                     <  3:326196/144‑1 (MQ=255)
                                                                                                                                 gacgTCGCCGCCCTGGCCGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTGCGCTCGTAGTAGGGGGATGGTAACAAAAAAATCGTGGTCTAGATGGGCTCACAAGCCCTGCAGGGTCTGCATCAg                    >  5:299012/4‑148 (MQ=255)
                                                                                                                                  aGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTGCGCTCGTAGTAGGGGGATGGTAACAAAAAAATCGTGGTCTAGATGGGCTCACAAGCCCTGCAGGGTCTGCAt                       <  5:251298/144‑1 (MQ=255)
                                                                                                                                  aGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTGCGCTCGTAGTAGGGGGATGGTAACAAAAAAATCGTGGTCTAGATGGGCTCACAAGCCCTGCAGGGTCTGCAt                       >  6:251298/1‑144 (MQ=255)
                                                                                                                                          gCCAGGGCGGCGACGTTGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTGCGCTCGTAGTAGGGGGATGGTAACAAAAAAATCGTGGTCTAGATg                                                  <  2:271366/109‑1 (MQ=255)
                                                                                                                                          gCCAGGGCGGCGACGTTGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTGCGCTCGTAGTAGGGGGATGGTAACAAAAAAATCGTGGTCTAGATg                                                  >  1:271366/1‑109 (MQ=255)
                                                                                                                                           ccctggccgcGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTGCGCTCGTAGTAGGGGGATGGTAACAAAAAAATCGTGGTCTAGATGGGCTCACAAGCCCTGCAGGGTCTGCATCAGCACCCGCACtt         <  6:299012/145‑1 (MQ=255)
                                                                                                                                                  ggcgACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTGCGCTCGTAGTAGGGGGATGGTAACAAAAAAATCGTGGTCTAGATGGGCTCACAAGCCCTGCAGGGTCTGCATCAGCACCCGCACTTTGGCGAt  <  6:112310/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
CCTGCTTGTTCAGCTTCCAGGCCATGTCGACCAGCACCACATCCTCGGCGCCGCTGAAGGAGATCCACAGGTCATCACCGAAGTGCTCGAAGGCAAGCTTGAGGATGTCCTGCGGGGACTTGTTGGCATAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTGCGCTCGTAGTAGGGGGATGGTAACAAAAAAATCGTGGTCTAGATGGGCTCACAAGCCCTGCAGGGTCTGCATCAGCACCCGCACTTTGGCGAT  >  NC_002947/2656022‑2656315

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: