Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,656,158 G→C 6.7% A11A (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,1580GC 6.7% 235.8 / 6.2 74A11A (GCC→GCG) cysHphosphoadenosine phosphosulfate reductase
Reads supporting (aligned to +/- strand):  ref base G (34/35);  new base C (3/2);  total (37/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.60e-01

CCTGCTTGTTCAGCTTCCAGGCCATGTCGACCAGCACCACATCCTCGGCGCCGCTGAAGGAGATCCACAGGTCATCACCGAAGTGCTCGAAGGCAAGCTTGAGGATGTCCTGCGGGGACTTGTTGGCA‑TAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTGCGCTCGTAGTAGGGGGATGGTAACAAAAAAATCGTGGTCTAGATGGGCTCACAAGCCCTGCAGGGTCTGCATCAG  >  NC_002947/2656022‑2656297
                                                                                                                                         |                                                                                                                                           
ccTGCTTGTTCAGCTTCCAGGCCATGTCGACCAGCACCACATCCTCGGCGCCGCTGAAGGAGATCCACAGGTCATCACCGAAGTGCTCGAAGGCAAGCTTGAGGATGTCCTGCGGGGACTTGTTGGCA‑TAGGTCGCGGCCAGggcggc                                                                                                                                  <  5:111192/148‑1 (MQ=255)
     ttGTTCAGCTTCCAGGCCATGTCGACCAGCACCACATCCTCGGCGCCGCTGAAGGAGATCCACAGGTCATCACCGAAGTGCTCGAAGGCAAGCTTGAGGATGTCCTGCGGGGACTTGTTGGCA‑TAGGTCGCGGCCAGGGCGGCGACGTc                                                                                                                            >  4:116661/1‑149 (MQ=255)
              ttCCAGGCCATGTCGACCAGCACCACATCCTCGGCGCCGCTGAAGGAGATCCACAGGTCATCACCGAAGTGCTCGAAGGCAAGCTTGAGGATGTCCTGCGGGGACTTGTTGGCA‑TAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTg                                                                                                                   >  5:28019/1‑149 (MQ=255)
                ccAGGCCATGTCGACCAGCACCACATCCTCGGCGCCGCTGAAGGAGATCCACAGGTCATCACCGAAGTGCTCGAAGGCAAGCTTGAGGATGTCCTGCGGGGACTTGTTGGCA‑TAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTgg                                                                                                                  <  3:116661/148‑1 (MQ=255)
                  aGGCCATGTCGACCAGCACCACATCCTCGGCGCCGCTGAAGGAGATCCACAGGTCATCACCGAAGTGCTCGAAGGCAAGCTTGAGGATGTCCTGCGGGGACTTGTTGGCA‑TAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGtt                                                                                                                    <  5:290936/144‑1 (MQ=255)
                  aGGCCATGTCGACCAGCACCACATCCTCGGCGCCGCTGAAGGAGATCCACAGGTCATCACCGAAGTGCTCGAAGGCAAGCTTGAGGATGTCCTGCGGGGACTTGTTGGCA‑TAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCt                                                                                                                >  5:128582/1‑148 (MQ=255)
                  aGGCCATGTCGACCAGCACCACATCCTCGGCGCCGCTGAAGGAGATCCACAGGTCATCACCGAAGTGCTCGAAGGCAAGCTTGAGGATGTCCTGCGGGGACTTGTTGGCA‑TAGGTCGCCGCCCTGGCCGCGAc                                                                                                                               <  6:289900/133‑1 (MQ=255)
                  aGGCCATGTCGACCAGCACCACATCCTCGGCGCCGCTGAAGGAGATCCACAGGTCATCACCGAAGTGCTCGAAGGCAAGCTTGAGGATGTCCTGCGGGGACTTGTTGGCA‑TAGGTCGCCGCCCTGGCCGCGAc                                                                                                                               >  5:289900/1‑133 (MQ=255)
                  aGGCCATGTCGACCAGCACCACATCCTCGGCGCCGCTGAAGGAGATCCACAGGTCATCACCGAAGTGCTCGAAGGCAAGCTAGAGGATGTCCTGCGGGGACTTGTTGGCA‑TAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGtt                                                                                                                    >  6:290936/1‑144 (MQ=255)
                      cATGTCGACCAGCACCACATCCTCGGCGCCGCTGAAGGAGATCCACAGGTCATCACCGAAGTGCTCGAAGGCAAGCTTGAGGATGTCCTGCGGGGACTTGTTGGCACTAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCtcatc                                                                                                            <  1:68549/149‑1 (MQ=255)
                            gACCAGCACCACATCCTCGGCGCCGCTGAAGGAGATCCACAGGTCATCACCGAAGTGCTCGAAGGCAAGCTTGAGGATGTCCTGCGGGGACTTGTTGGCA‑TAGGTCGCGGCCAGggcggc                                                                                                                                  <  4:320695/120‑1 (MQ=255)
                                    ccacATCCTCGGCGCCGCTGAAGGAGATCCACAGGTCATCACCGAAGTGCTCGAAGGCAAGCTTGAGGATGTCCTGCGGGGACTTGTTGGCA‑TAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCt                                                                                             >  2:119080/1‑149 (MQ=255)
                                     cacaTCCTCGGCGCCGCTGAAGGAGATCCACAGGTCATCACCGAAGTGCTCGAAGGCAAGCTTGAGGATGTCCTGCGGGGACTTGTTGGCA‑TAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCtcatc                                                                                                            <  3:131643/133‑1 (MQ=255)
                                     cacaTCCTCGGCGCCGCTGAAGGAGATCCACAGGTCATCACCGAAGTGCTCGAAGGCAAGCTTGAGGATGTCCTGCGGGGACTTGTTGGCA‑TAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCtcatc                                                                                                            >  4:131643/1‑133 (MQ=255)
                                     cacaTCCTCGGCGCCGCTGAAGGAGATCCACAGGTCATCACCGAAGTGCTCGAAGGCAAGCTTGAGGATGTCCTGCGGGGACTTGTTGGCA‑TAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCTCATCAGGCTCTTTCCATCtt                                                                                            >  6:6070/1‑149 (MQ=255)
                                     cacaTCCTCGGCGCCGCTGAAGGAGATCCACAGGTCATCACCGAAGTGCTCGAAGGCAAGCTTGAGGATGTCCTGCGGGGACTTGTTGGCA‑TAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCtt                                                                                            >  3:265903/1‑149 (MQ=255)
                                            tCGGCGCCGCTGAAGGAGATCCACAGGTCATCACCGAAGTGCTCGAAGGCAAGCTTGAGGATGTCCTGCGGGGACTTGTTGGCA‑TAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCtgtg                                                                                     >  8:241975/1‑149 (MQ=255)
                                                     cTGAAGGAGATCCACAGGTCATCACCGAAGTGCTCGAAGGCAAGCTTGAGGATGTCCTGCGGGGACTTGTTGGCA‑TAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTgc                                                                            >  5:316834/1‑149 (MQ=255)
                                                       gAAGGAGATCCACAGGTCATCACCGAAGTGCTCGAAGGCAAGCTTGAGGATGTCCTGCGGGGACTTGTTGGCA‑TAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTgcgc                                                                          >  3:11295/1‑149 (MQ=255)
                                                         aGGAGATCCACAGGTCATCACCGAAGTGCTCGAAGGCAAGCTTGAGGATGTCCTGCGGGGACTTGTTGGCA‑TAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTGCGCTc                                                                        <  6:128582/149‑1 (MQ=255)
                                                            agaTCCACAGGTCATCACCGAAGTGCTCGAAGGCAAGCTTGAGGATGTCCTGCGGGGACTTGTTGGCA‑TAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTGCGCTCgta                                                                     <  4:265903/149‑1 (MQ=255)
                                                                 cacaGGTCATCACCGAAGTGCTCGAAGGCAAGCTTGAGGATGTCCTGCGGGGACTTGTTGGCA‑TAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTGCGCTCGTAGTAgg                                                                <  7:241975/149‑1 (MQ=255)
                                                                 cacaGGTCATCACCGAAGTGCTCGAAGGCAAGCTTGAGGATGTCCTGCGGGGACTTGTTGGCA‑TAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTGCGCTCGTAGTAgg                                                                >  7:195477/1‑149 (MQ=255)
                                                                  acaGGTCATCACCGAAGTGCTCGAAGGCAAGCTTGAGGATGTCCTGCGGGGACTTGTTGGCA‑TAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTGCGCTCGTAGTAggg                                                               >  1:303195/1‑149 (MQ=255)
                                                                      gTCATCACCGAAGTGCTCGAAGGCAAGCTTGAGGATGTCCTGCGGGGACTTGTTGGCA‑TAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTGCGCTCGTAGTAGGGGGAt                                                           >  5:255810/1‑149 (MQ=255)
                                                                         atcaCCGAAGTGCTCGAAGGCAAGCTTGAGGATGTCCTGCGGGGACTTGTTGGCA‑TAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTGCGCTCGTAGTAggg                                                               <  8:276723/142‑1 (MQ=255)
                                                                         atcaCCGAAGTGCTCGAAGGCAAGCTTGAGGATGTCCTGCGGGGACTTGTTGGCA‑TAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTGCGCTCGTAGTAggg                                                               >  7:276723/1‑142 (MQ=255)
                                                                             ccGAAGTGCTCGAAGGCAAGCTTGAGGATGTCCTGCGGGGACTTGTTGGCA‑TAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTGCGCTCGTAGTAgg                                                                >  3:37333/1‑137 (MQ=255)
                                                                             ccGAAGTGCTCGAAGGCAAGCTTGAGGATGTCCTGCGGGGACTTGTTGGCA‑TAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTGCGCTCGTAGTAgg                                                                <  4:37333/137‑1 (MQ=255)
                                                                                     cTCGAAGGCAAGCTTGAGGATGTCCTGCGGGGACTTGTTGGCA‑TAGGTCGCGGCCAGGGTGGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTGCGCTCGTAGTAGGGGGATGGTAACAAAAAAATc                                            >  5:306997/1‑149 (MQ=255)
                                                                                     cTCGAAGGCAAGCTTGAGGATGTCCTGCGGGGACTTGTTGGCA‑TAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTGCGCTCGTAGTAGGGGGATGGTAACAAAAAAAt                                             <  1:119080/148‑1 (MQ=255)
                                                                                     cTCGAAGGCAAGCTTGAGGATGTCCTGCGGGGACTTGTTGGCA‑TAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTGCGCTCGTAGTAGGGGGATGGTAACAAAAAAATc                                            >  4:223729/1‑149 (MQ=255)
                                                                                      tCGAAGGCAAGCTTGAGGATGTCCTGCGGGGACTTGTTGGCA‑TAGGTCGCGGCCAgg                                                                                                                                       >  5:288638/1‑57 (MQ=255)
                                                                                      tCGAAGGCAAGCTTGAGGATGTCCTGCGGGGACTTGTTGGCA‑TAGGTCGCGGCCAgg                                                                                                                                       <  6:288638/57‑1 (MQ=255)
                                                                                      tCGAAGGCAAGCTTGAGGATGTCCTGCGGGGACTTGTTGGCA‑TAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTGCGCTCGTAGTAGGGGGATGGTAACAAAAAAATCg                                           >  6:349375/1‑149 (MQ=255)
                                                                                      tCGAAGGCAAGCTTGAGGATGTCCTGCGGGGACTTGTTGGCA‑TAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTGCGCTCGTAGTAGGGGGATGGTAACAAAAAAATCg                                           >  8:210765/1‑149 (MQ=255)
                                                                                        gAAGGCTATCTTGAGGATGTCCTGCGGGGACTTGTTGTCA‑TAGGTCGCGGCCAGGGCGGCGACGTAGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTGCGCTCGTAGTAGGGGGATGGTAACAAAAAAATCGt                                          <  5:6070/148‑1 (MQ=255)
                                                                                        gAAGGCAAGCTTGAGGATGTCCTGCGGGGACTTGTTGGCA‑TAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTGCGCTCGTAGTAGGGGGATGGTAACAAAAAAATCGTg                                         >  8:64538/1‑149 (MQ=255)
                                                                                             cAAGCTTGAGGATGTCCTGCGGGGACTTGTTGGCA‑TAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTGCGCTCGTAGTAGGGGGATGGTAACAAAAAAATCGTGGTCTa                                    <  8:195477/149‑1 (MQ=255)
                                                                                                gCTTGAGGATGTCCTGCGGGGACTTGTTGGCA‑TAGGTCGCGGCCAGGGTGGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTGCGCTCGTAGTAGGGGGATGGTAACAAAAAAATCGTGGTCTAGAt                                 <  6:306997/149‑1 (MQ=255)
                                                                                                gCTTGAGGATGTCCTGCGGGGACTTGTTGGCA‑TAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTGCGCTCGTAGTAGGGGGATGGTAACAAAAAAATCGTGGTCTAGAt                                 >  1:174215/1‑149 (MQ=255)
                                                                                                gCTTGAGGATGTCCTGCGGGGACTTGTTGGCA‑TAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTGCGCTCGTAGTAGGGGGATGGTAACAAAAAAATCGTGGTCTAGAt                                 <  5:349375/149‑1 (MQ=255)
                                                                                                  ttGAGGATGTCCTGCGGGGACTTGTTGGCA‑TAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTGCGCTCGTAGTAGGGGGGTGGTAACAAAAAAATCGTGGTCTAGATgg                               >  1:65101/1‑149 (MQ=255)
                                                                                                     aGGATGTCCTGCGGGGACTTGTTGGCA‑TAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTGCGCTCGTAGTAGGGGGATGGTAACAAAAAAATCGTGGTCTAGATGGGCt                            >  1:327126/1‑149 (MQ=255)
                                                                                                             cTGCGGGGACTTGTTGGCA‑TAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTGCGCTCGTAGTAGGGGGATGGTAACAAAAAAATCGTGGTCTAGATGGGCTCACAAGc                     <  6:316834/148‑1 (MQ=255)
                                                                                                               gCGGGGACTTGTTGGCA‑TAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTGCGCTCGTAGTAGGGGGATGGTAACaaaaaa                                               >  8:70706/1‑120 (MQ=255)
                                                                                                               gCGGGGACTTGTTGGCA‑TAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTGCGCTCGTAGTAGGGGGATGGTAACaaaaaa                                               <  7:70706/120‑1 (MQ=255)
                                                                                                               gCGGGGACTTGTTGGCA‑TAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTGCGCTCGTAGTAGGGGGATGGTAACAAAAAAATCGTGGTCTAGATGGGCTCACAAGCCCt                  >  2:41917/1‑149 (MQ=255)
                                                                                                                cGGGGACTTGTTGGCA‑TAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTGCGCTCGTAGTAGGGGGATGGTAACAAAAAAATCGTGGTCTAGATGGGCTCACAAGccc                   <  7:210765/147‑1 (MQ=255)
                                                                                                                cGGGGACTTGTTGGCA‑TAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTGCGCTCGTAGTAGGGGGATGGTAACAAAAAAATCGTGGTCTAGATGGGCTCACAAGCCCt                  <  6:255810/148‑1 (MQ=255)
                                                                                                                cGGGGACTTGTTGGCA‑TAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTGCGCTCGTAGTAGGGGGATGGTAACAAAAAAATCGTGGTCTAGATGGGCTCACAAGCCCt                  <  3:223729/148‑1 (MQ=255)
                                                                                                                  gggACTTGTTGGCA‑TAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTGCGCTCGTAGTAGGGGGGTGGTAACAAAAAAATCGTGGTCTAGATGGGCTCACAAGCCCTGCa               <  2:65101/149‑1 (MQ=255)
                                                                                                                  gggACTTGTTGGCA‑TAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTGCGCTCGTAGTAGGGGGATGGTAACAAAAAAATCGTGGTCTAGATGGGCTCACAAGCCCTg                 >  5:112310/1‑147 (MQ=255)
                                                                                                                   ggACTTGTTGGCA‑TAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTGCGCTCGTAGTAGGGGGATGGTAACAAAAAAATCGTGGTCTAGATGGGCTCACAAGCCCTGCAg              <  1:41917/149‑1 (MQ=255)
                                                                                                                      cTTTTTGGCA‑TAGTTCTCGGCGAGTGCGGCGACGGCGAAGGGTGGGCTCATCAGGCGGTTTCCATCTTGGCTGTGTCGCTCTGCGCGCGTAGTAGGGGGTTGGTAACAAAAAAATCGTGGTCTAGATGGGCTCACAAGCCCTGCAGGGt           <  6:326762/149‑1 (MQ=255)
                                                                                                                         gttgGCA‑TAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTGCGCTCGTTGTAGGGGGATGGTAACAAAAAAATCGTGGTCTAGATGGGCTCACAAGCCCTGCAGGGTCTg        >  1:136997/1‑149 (MQ=255)
                                                                                                                         gttgGCA‑TAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTGCGCTCGTAGTAGGGGGATGGTAACAAAAAAATCGTGGTCTAGATGGGCTCACAAGCCCTGCAGGGTCt         <  2:303195/148‑1 (MQ=255)
                                                                                                                         gttgGCA‑TAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTGCGCTCGTAGTAGGGGGATGGTAACAAAAAAATCGTGGTCTAGATGGGCTCACAAGCCCTGCAGGGTCTg        >  1:305331/1‑149 (MQ=255)
                                                                                                                         gttgGCA‑TAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTGCGCTCGTAGTAGGGGGATGGTAACAAAAAAATCGTGGTCTAGATGGGCTCACAAGCCCTGCAGGGTCTg        >  5:20492/1‑149 (MQ=255)
                                                                                                                         gttgGCA‑TAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTGCGCTCGTAGTAGGGGGATGGTAACAAAAAAATCGTGGTCTAGATGGGCTCACAAGCCCTGCAGGGTCTg        >  4:323800/1‑149 (MQ=255)
                                                                                                                          ttgGCA‑TAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTGCGCTCGTAGTAGGGGGATGGTAACAAAAAAATCGTGGTCTAGATGGGCTCACAAGCCCTGCAGGGTCTGc       <  2:174215/149‑1 (MQ=255)
                                                                                                                          ttgGCA‑TAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTGCGCTCGTAGTAGGGGGATGGTAACAAAAAAATCGTGGTCTAGATGGGCTCACAAGCCCTGCAGGGTCTGc       <  7:64538/149‑1 (MQ=255)
                                                                                                                             gCA‑TAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTGCGCTCGTAGTAGGGGGATGGTAACAAAAAAATCGTGGTCTAGATGGGCTCACAAGCCCTGCAGGGTCTGCATc    <  4:11295/149‑1 (MQ=255)
                                                                                                                              cA‑TAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTGCGCTCGTAGTAGGGGGATGGTAACAAAAAAATCGTGGTCTAGATGGGCTCACAAGCCCTGCAGGGTCTGCATCa   <  3:323800/149‑1 (MQ=255)
                                                                                                                              cA‑TAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTGCGCTCGTAGTAGGGGGATGGTAACAAAAAAATCGTGGTCTAGATGGGCTCACAAGCCCTGCAGGGTCTGCATCa   <  2:305331/149‑1 (MQ=255)
                                                                                                                              cA‑TAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTGCGCTCGTAGTAGGGGGATGGTAACAAAAAAATCGTGGTCTAGATGGGCTCACAAGCCCTGCAGGGTCTGCATCa   <  6:20492/149‑1 (MQ=255)
                                                                                                                               tcgacgTCGCCGCCCTGGCCGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTGCGCTCGTAGTAGGGGGATGGTAACAAAAAAATCGTGGTCTAGATGGGCTCACAAGCCCTGCAGGGTCTGCAt     >  4:326196/6‑147 (MQ=255)
                                                                                                                               a‑TAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTGCGCTCGTAGTAGGGGGATGGTAACAAAAAAATCGTGGTCTAGAt                                 <  5:344538/118‑1 (MQ=255)
                                                                                                                               a‑TAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTGCGCTCGTAGTAGGGGGATGGTAACAAAAAAATCGTGGTCTAGAt                                 >  6:344538/1‑118 (MQ=255)
                                                                                                                               a‑TAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTGCGCTCGTAGTAGGGGGATGGTAACAAAAAAATCGTGGTCTAGATgg                               <  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)
                                                                                                                                         |                                                                                                                                           
CCTGCTTGTTCAGCTTCCAGGCCATGTCGACCAGCACCACATCCTCGGCGCCGCTGAAGGAGATCCACAGGTCATCACCGAAGTGCTCGAAGGCAAGCTTGAGGATGTCCTGCGGGGACTTGTTGGCA‑TAGGTCGCGGCCAGGGCGGCGACGTCGAAGGGTTGGCTCATCAGGCGGTTTCCATCTTGGCTGTGGCGCTGTGCGCTCGTAGTAGGGGGATGGTAACAAAAAAATCGTGGTCTAGATGGGCTCACAAGCCCTGCAGGGTCTGCATCAG  >  NC_002947/2656022‑2656297

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: ‑