Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,087,189 C→G 6.4% G45A (GGC→GCC)  astB ← N‑succinylarginine dihydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,087,1890CG 6.4% 250.8 / 3.9 76G45A (GGC→GCC) astBN‑succinylarginine dihydrolase
Reads supporting (aligned to +/- strand):  ref base C (39/32);  new base G (2/3);  total (41/35)
Fisher's exact test for biased strand distribution p-value = 6.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.66e-01

CTCGCGGGCGGCGCGCTGGATAACTTCGGCATCGCTGCCACTGAACCCCAGGCGGCGCAGCGCCGCCACGTCCGGGCGCTCTTGTGGCGCCAGAACACCCTGCTTGAAGCCCATGTCGGCCAGCGCTTTCATCTTCGCCAGGCCCTGGCGCGCGGCCTCACGCGGGTTGGACGCCTGCTGGCTATTGCTCTGCGAAGCCACGTTACCGTAAGACAGGCCGCCATAGTTGTGGGTAGGCCCCACCAGGCCATCAAAATTCACTTCATAGGATTTCATCGGCTAGGCTCCGT  >  NC_002947/5087047‑5087336
                                                                                                                                                 |                                                                                                                                                   
cTCGCGGGCGGCGCGCTGGATAACTTCGGCATCGCTGCCACTGAACCCCAGGCGGCGCAGCGCCGCCACGTCCGGGCGCTCTTGTGGCGCCAGAACACCCTGCTTGAAGCCCATGTCGGCCAGCGCTTTCATCTTCGCCAGGCCCTGgc                                                                                                                                               <  2:282760/149‑1 (MQ=255)
   gcgggcggCGCGCTGGATAACTTCGGCATCGCTGCCACTGAACCCCAGGCGGCGCAGCGCCGCCACGTCCGGGCGCTCTTGTGGCGCCAGAACACCCTGCTTGAAGCCCATGTCGGCCAGCGCTTCATCTTCGCCAGGGCCTGgcgaa                                                                                                                                            <  4:236146/148‑3 (MQ=255)
    atcgcggcgCGCTGGATAACTTCGGCATCGCTGCCACTGACCCCCCGGCGGCGCAGCGCCGCCACGTCCGGGCGCTCTTGTGGCGCCAGAACACCCTGCTTGAAGCCCATGTCGGCCAGCGCTTTCATCTTCGCCAGGGCCTGgcgaag                                                                                                                                           <  1:125412/146‑4 (MQ=255)
      ggcggcGCGGTGGATAACTTCGGCATCGCGCCCACTCACCCCCCGGCGGCGCAGCGCCGCCACGTCCGGGCGCTCTTGTGGCGCCAGAACACCCTGCTTGAAGCCCATGTCGGCCAGCGCTTTCATCTTCGCCAGGCCCTGGCGCGCgg                                                                                                                                         <  1:165707/149‑1 (MQ=255)
         ggcgCGCTGGATAACTTCGGCATCGCTGCCACTGAACCCCAGGCGGCGCAGCGCCGCCACGTCCGGGCGCTCTTGTGGCGCCAGAACACCCTGCTTGAAGCCCATGTCGGCCAGCGCTTTCATCTTCGCCAGGGCCTGgcgaag                                                                                                                                           >  2:125412/1‑141 (MQ=255)
                  gATAACTTCGGCATCGCTGCCACTGAACCCCAGGCGGCGCAGCGCCGCCACGTCCGGGCGCTCTTGTGGCGCCAGAACACCCTGCTTGAAGCCCATGTCGGCCAGCGCTTTCATCTTCGCCAGGCCCTGGCGCGCGGCCTCACGCGGGt                                                                                                                             >  7:124272/1‑149 (MQ=255)
                   aTAACTTCGGCATCGCTGCCACTGAACCCCAGGCGGCGCAGCGCCGCCACGTCCGGGCGCTCTTGTGGCGCCAGAACACCCTGCTTGAAGCCCATGTCGGCCAGCGCTTTCATCTTCGCCAGGCCCTGGCGCGCGGCCTCACGCGGGtt                                                                                                                            <  3:182646/149‑1 (MQ=255)
                   aTAACTTCGGCATCGCTGCCACTGAACCCCAGGCGGCGCAGCGCCGCCACGTCCGGGCGCTCTTGTGGCGCCAGAACACCCTGCTTGAAGCCCATGTCGGCCAGCGCTTTCATCTTCGCCAGGCCCTGGCGCGCGGCCTCACGCGGGtt                                                                                                                            <  7:60362/149‑1 (MQ=255)
                   aTAACTTCGGCATCGCTGCCACTGAACCCCAGGCGGCGCAGCGCCGCCACGTCCGGGCGCTCTTGTGGCGCCAGAACACCCTGCTTGAAGCCCATGTCGGCCAGCGCTTTCATCTTCGACCGGCCCTGGCGCGCGGCCTCACGCGGGtt                                                                                                                            >  8:259961/1‑149 (MQ=255)
                    tAACTTCGGCATCGCTGCCACTGAACCCCAGGCGGCGCAGCGCCGCCACGTCCGGGCGCTCTTGTGGCGCCAGAACACCCTGCTTGAAGCCCATGTCGGCCAGCGCTTTCATCTTCGCCAGGCCCTGGCGCGCGGCCTCACGCGGGTTg                                                                                                                           <  3:179162/149‑1 (MQ=255)
                        ttCGGCATCGCTGCCACTGAACCCCAGGCGGCGCAGCGCCGCCACGTCCGGGCGCTCTTGTGGCGCCAGAACACCCTGCTTGAAGCCCATGTCGGCCAGCGCTTTCATCTTCGCCAGGGCCTGgcgaagat                                                                                                                                         >  1:30679/1‑126 (MQ=255)
                        ttCGGCATCGCTGCCACTGAACCCCAGGCGGCGCAGCGCCGCCACGTCCGGGCGCTCTTGTGGCGCCAGAACACCCTGCTTGAAGCCCATGTCGGCCAGCGCTTTCATCTTCGCCAGGGCCTGgcgaagat                                                                                                                                         <  2:30679/131‑6 (MQ=255)
                               tggCTGCCACTGAACCCCCGGCGGCGCAGCGCCGCCACGTCCGGGCGCTCTTGTGGCGCCAGAACACCCTGCTTGAAGCCCATGTCGGCCAGCGCTTTCATCTTCGCCAGGCCCTGGCGCGCGGCCTCACGCGGGTTGGACGCctgctg                                                                                                                <  3:123066/147‑1 (MQ=255)
                                  cTGCCACTGAACCCCAGGCGGCGCAGCGCCGCCACGTCCGGGCGCTCTTGTGGCGCCAGAACACCCTGCTTGAAGCCCATGTCGGCCAGCGCTTTCATCTTCGCCAGGCCCTGGCGCGCGGCCTCACGCGGGTTGGACGCCTGCTGGc                                                                                                              >  4:315372/1‑148 (MQ=255)
                                   tGCCACTGAACCCCAGGCGGCGCAGCGCCGCCACGTCCGGGCGCTCTTGTGGCGCCAGAACACCCTGCTTGAAGCCCATGTCGGCCAGCGCTTTCATCTTCGCCAGGCCCTGGCGCGCGGCCTCACGCGGGTTGGACGCCTGCTGGCTa                                                                                                            >  7:120635/1‑149 (MQ=255)
                                      cACTGAACCCCAGGCGGCGCAGCGCCGCCACGTCCGGGCGCTCTTGTGGCGCCAGAACACCCTGCTTGAAGCCCATGTCGGCCAGCGCTTTCATCTTCGCCAGGCCCTGGCGCGCGGCCTCACGCGGGTTGGACGCCTGCTGGCTATTg                                                                                                         <  7:205351/149‑1 (MQ=255)
                                       aCTGAACCCCAGGCGGCGCAGCGCCGCCACGTCCGGGCGCTCTTGTGGCGCCAGAACACCCTGCTTGAAGCCCATGTCGGCCAGCGCTTTCATCTTCGCCAGGCCCTGGCGCGCGGCCTCACGCGGGTTGGACGCCTGCTGGCTATTGc                                                                                                        >  7:265282/1‑149 (MQ=255)
                                        cTGAACCCCAGGCGGCGCAGCGCCGCCACGTCCGGGCGCTCTTGTGGCGCCAGAACACCCTGCTTGAAGCCCATGTCGGCCAGCGCTTTCATCTTCGCCAGGCCCTGGCGCGCGGCCTCACGCGGGTTGGACGCctgct                                                                                                                 >  3:21127/1‑139 (MQ=255)
                                           aaCCCCAGGCGGCGCAGCGCCGCCACGTCCGGGCGCTCTTGTGGCGCCAGAACACCCTGCTTGAAGCCCATGTCGGCCAGCGCTTTCATCTTCGCCAGGCCCt                                                                                                                                                  <  4:311969/103‑1 (MQ=255)
                                           aaCCCCAGGCGGCGCAGCGCCGCCACGTCCGGGCGCTCTTGTGGCGCCAGAACACCCTGCTTGAAGCCCATGTCGGCCAGCGCTTTCATCTTCGCCAGGCCCt                                                                                                                                                  >  3:311969/1‑103 (MQ=255)
                                           aaCCCCAGGCGGCGCAGCGCCGCCACGTCCGGGCGCTCATGTGGCGCCAGAACACCCTGCTTGAAGCCCATGTCGGCCAGCGCTTTCATCTTCGCCAGGCCCTGGCGCGCGGCCTCACGCGGGTTGGACGCCTGCTGGCTATTGctc                                                                                                      <  2:219851/147‑1 (MQ=255)
                                              cccAGGCGGCGCAGCGCCGCCACGTCCGGGCGCTCTTGTGGCGCCAGAACACCCTGCTTGAAGCCCATGTCGGCCAGCGCTTTCATCTTCGCCAGGCCCTGGCGCGCGGCCTCACGCGGGTTGGACGCCTGCTGGCTATTGCTCTGCg                                                                                                  >  7:275214/1‑148 (MQ=255)
                                               ccAGGCGGCGCAGCGCCGCCACGTCCGGGCGCTCTTGTGGCGCCAGAACACCCTGCTTGAAGCCCATGTCGGCCAGCGCTTTCATCTTCGCCAGGCCCTGGCGCGCGGCCTCACGCGGGTTGGACGCCTGCTGGCTATTGCTCTGCGaa                                                                                                <  8:275214/149‑1 (MQ=255)
                                                 aGGCGGCGCAGCGCCGCCACGTCCGGGCGCTCTTGTGGCGCCAGAACACCCTGCTTGAAGCCCATGTCGGCCAGCGCTTTCATCTTCGCCAGGCCCTGGCGCGCGGCCTCACGCGGGTTGGACGCCTGCTGGCTATTGCTCTGCGAAGc                                                                                              >  5:300524/1‑149 (MQ=255)
                                                   gcggcgCAGCGCCGCCACGTCCGGGCGCTCTTGTGGCGCCAGAACACCCTGCTTGAAGCCCATGTCGGCCAGCGCTTTCATCTTAGCCAGGCACTGGCGCGCGGCCTCACGCGGGTTGGACGCCTGCTGGCTATTGCTCTGCGAAGCCa                                                                                            >  2:135551/1‑149 (MQ=255)
                                                         cAGCGCCGCCACGTCCGGGCGCTCTTGTGGCGCCAGAACACCCTGCTTGAAGCCCATGTCGGCCAGCGCTTTCATCTCTCCCCCAGGCCCTGGCGCGCGGCCTCACGCGGGTTGGACGCCTGCTGGCTATTGCTCTGCGAAGCCACGtt                                                                                        <  1:261574/149‑1 (MQ=255)
                                                                  cACGTCCGGGCGCTCTTGTGGCGCCAGAACACCCTGCTTGAAGCCCATGTCGGCCAGCGCTTTCATCTTCGCCAGGCCCTGGCGCGCGGCCTCACGCGGGTTGGACGCCTGCTGGCTATTGCTCTGCGAAGCCACGTTACCGTAAGACa                                                                             <  2:6637/149‑1 (MQ=255)
                                                                       ccGGGCGCTCTTGTGGCGCCAGAACACCCTGCTTGAAGCCCATGTCGGCCAGCGCTTTCATCTTCGCCAGGCCCTGGCGCGCGGCCTCACGCGGGTTGGACGCCTGCTGGCTATTGCTCTGCGAAGCCACGTTACCg                                                                                    <  7:226146/137‑1 (MQ=255)
                                                                       ccGGGCGCTCTTGTGGCGCCAGAACACCCTGCTTGAAGCCCATGTCGGCCAGCGCTTTCATCTTCGCCAGGCCCTGGCGCGCGGCCTCACGCGGGTTGGACGCCTGCTGGCTATTGCTCTGCGAAGCCACGTTACCg                                                                                    >  8:226146/1‑137 (MQ=255)
                                                                               tctTGTGGCGCCAGAACACCCTGCTTGAAGCCCATGTCGGCCAGCGCTTTCATCTTCGCCAGGCCCTGGCGCGCGGCCTCACGCGGGTTGGACGCCTGCTGGCTATTGCTCTGCGAAGCCACGTTACCGTAAGACAGGCCGCCATAGtt                                                                >  2:228792/1‑149 (MQ=255)
                                                                               tctTGTGGCGCCAGAACACCCTGCTTGAAGCCCATGTCGGCCAGCGCTTTCATCTTCGCCAGGCCCTGGCGCGCGGCCTCACGCGGGTTGGACGCCTGCTGGCTATTGCTCTGCGAAGCCACGTTACCGTAAGACAGGCCGCAATAGtt                                                                >  3:244508/1‑149 (MQ=255)
                                                                                 ttGTGGCGCCAGAACACCCTGCTTGAAGCCCATGTCGGCCAGCGCTTTCATCTTCGCCAGGCCCTGGCGCGCGGCCTCACGCGGGTTGGACGCCTGCTGGCTATTGCTCTGCGAAGCCACGTTACCGTAAGACAGGCCGCCATAGTtgt                                                              >  5:101109/1‑149 (MQ=255)
                                                                                   gtgGCGCCAGAACACCCTGCTTGAAGCCCATGTCGGCCAGCGCTTTCATCTTCGCCAGGCCCTGGCGCGCGGCCTCACGCGGGTTGGACGCCTGCTGGCTATTGCTCTGCGAAGCCACGTTACCGTAAGACAGGCCGCCATAGTtgt                                                              <  8:120635/147‑1 (MQ=255)
                                                                                   gtgGCGCCAGAACACCCTGCTTGAAGCCCATGTCGGCCAGCGCTTTCATCTTCGCCAGGCCCTGGCGCGCGGCCTCACGCGGGTTGGACGCCTGCTGGCTATTGCTCTGCGAAGCCACGTTACCGTAAGACAGGCCGCCATAGTtgt                                                              <  2:73865/147‑1 (MQ=255)
                                                                                      gcgcCAGAACACCCTGCTTGAAGCCCATGTCGGCCAGCGCTTTCATCTTCGCCAGGCCCTGGCGCGCGGCCTCACGCGGGTTGGACGCCTGCTGGCTATTGCTCTGCGACGCCACGTTACCGTAAGACAGGCCGCCATAGTTGTggg                                                           >  7:135126/1‑147 (MQ=255)
                                                                                        gcCAGAACACCCTGCTTGAAGCCCATGTCGGCCAGCGCTTTCATCTTCGCCAGGCCCTGGCGCGCGGCCTCACGCGGGTTGGACGCCTGCTGGCTATTGCTCTGCGAAGCCACGTTACCGTAAGACAGGCCGCCATAGTTGTGGGTAgg                                                       >  1:191435/1‑149 (MQ=255)
                                                                                                  ccTGCTTGAAGCCCATGTCGGCCAGCGGCTTTCATCTTCGCCAGGCCCTGGCGCGCGGCCTCACGCGGGTTGGACGCCTGCTGGCTATTGCTCTGCGAAGCCACGTTACCGTAAGACAGGCCGCCATAGTTGTGGGTAGGCCCCACCAg                                              >  1:14481/1‑149 (MQ=255)
                                                                                                  ccTGCTTGAAGCCCATGTCGGCCAGCGCTTTCATCTTCGCCAGGCCCTGGCGCGCGGCCTCACGCGGGTTGGACGCCTGCTGGCTATTGCTCTGCGAAGCCACGt                                                                                         <  6:104599/105‑1 (MQ=255)
                                                                                                  ccTGCTTGAAGCCCATGTCGGCCAGCGCTTTCATCTTCGCCAGGCCCTGGCGCGCGGCCTCACGCGGGTTGGACGCCTGCTGGCTATTGCTCTGCGAAGCCACGt                                                                                         >  5:104599/1‑105 (MQ=255)
                                                                                                  ccTGCTTGAAGCCCATGTCGGCCAGCGCTTTCATCTTCGCCAGGCCCTGGCGCGCGGCCTCACGCGGGTTGGACGCCTGCTGGCTATTGATCTGCGAAGCCACCTTACCGTAAGACAGGCCGCCATAGTTGTGGGTATGCCCAACCAgg                                             >  4:237023/1‑149 (MQ=255)
                                                                                                    tGCTTGAAGCCCATGTCGGCCAGCGCTTTCATCTTCGCCAGGCCCTGGCGCGCGGCCTCACGCGGGTTGGACGCCTGCTGGCTATTGCTCTGCGAAGCCACGTTACCGTAAGACAGGCCGCCATAGTTGTGGGTAGGCCCCACCAGGc                                            >  3:91146/1‑148 (MQ=255)
                                                                                                         gAAGCCCATGTCGGCCAGCGCTTTCATCTTCGCCAGGCCCTGGCGCGCGGCCTCACGCGGGTTGGACGCCTGCTGGCTATTGCTCTGCGAAGCCACGTTACCGTAAGa                                                                               <  5:287492/108‑1 (MQ=255)
                                                                                                         gAAGCCCATGTCGGCCAGCGCTTTCATCTTCGCCAGGCCCTGGCGCGCGGCCTCACGCGGGTTGGACGCCTGCTGGCTATTGCTCTGCGAAGCCACGTTACCGTAAGa                                                                               >  6:287492/1‑108 (MQ=255)
                                                                                                         gAAGCCCATGTCGGCCAGCGCTTTCATCTTCGCCAGGCCCTGGCGCGCGGCCTCACGCGGGTTGGACGCCTGCTGGCCATTGCTCTGCGAAGCCACGTTACCGTAAGACAGGCCGCCATAGTTGTGGGTAGGCCCCACCAGGCCATCa                                       >  3:179973/1‑148 (MQ=255)
                                                                                                          aaGCCCATGTCGGCCAGCGCTTTCATCTTCGCCAGGCCCTGGCGCGCGGCCTCACGCGGGTTGGACGCCTGCTGGCTATTGCTCTGCGAAGCCACGTTACCGTAAGACAGGCCGCCATAGTTGTGGGTAGGCCCCACCAGGCCATCaaa                                     <  3:237023/149‑1 (MQ=255)
                                                                                                           aGCCCATGTCGGCCAGCGCTTTCATCTTCGCCAGGCCCTGGCGCGCGGCCTCACGCGGGTTGGACGCCTGCTGGCTATTGCTCTGCGAAGCCACGTTACCGTAAGACAGGCCGCCATAGTTGTGGGTAGGCCCCACCAGGCCATCaaaa                                    <  8:265282/149‑1 (MQ=255)
                                                                                                           aGCCCATGTCGGCCAGCGCTTTCATCTTCGCCAGGCCCTGGCGCGCGGCCTCACGCGGGTTGGACGCCTGCTGGCTATTGCTCTGCGAAGCCACGTTACCGTAAGACAGGCCGCCATAGTTGTGGGTAGGCCCCACCAGGCCATCaaaa                                    >  5:192248/1‑149 (MQ=255)
                                                                                                              ccATGTCGGCCAGCGCTTTCATCTTCGCCAGGCCCTGGCGCGCGGCCTCACGCGGGTTGGACGCCTGCTGGCTATTGCTCTGCGAAGCCACGTTACCGTAAGACa                                                                             >  6:19096/1‑105 (MQ=255)
                                                                                                              ccATGTCGGCCAGCGCTTTCATCTTCGCCAGGCCCTGGCGCGCGGCCTCACGCGGGTTGGACGCCTGCTGGCTATTGCTCTGCGAAGCCACGTTACCGTAAGACa                                                                             <  5:19096/105‑1 (MQ=255)
                                                                                                              ccATGTCGGCCAGCGCTTTCATCTTCGCCAGGCCCTGGCGCGCGGCCTCACGCGGGTTGGACGCCTGCTGGCTATTGCTCTGCGAAGCCACGTTACCGTAAGACAGGCCGCCATAGTTGTGGGTAGGCCCCACCAGGCCATCAAAAtt                                  >  5:195146/1‑148 (MQ=255)
                                                                                                                aTGTCGGCCAGCGCTTTCATCTTCGCCAGGCCCTGGCGCGCGGCCTCACGCGGGTTGGACGCCTGCTGGCTATTGCTCTGCGAAGCCACGTTACCGTAAGACAGGCCGCCATATTTGTGGGTAGGc                                                      >  8:285322/1‑126 (MQ=255)
                                                                                                                aTGTCGGCCAGCGCTTTCATCTTCGCCAGGCCCTGGCGCGCGGCCTCACGCGGGTTGGACGCCTGCTGGCTATTGCTCTGCGAAGCCACGTTACCGTAAGACAGGCCGCCATAGTTGTGGGTAGGc                                                      <  7:285322/126‑1 (MQ=255)
                                                                                                                 tGTCGGCCAGCGCTTTCATCTTCGCCAGGCCCTGGCGCGCGGCCTCACGCGGGTTGGACGCCTGCTGGCTATTGCTCTGCGAAGCCACGTTACCGTAAGACAGGCCGCCATAGTTGTGGGTAGGCCCCACCAGGCCATCAAAATTCACt                              >  4:101350/1‑149 (MQ=255)
                                                                                                                  gTCGGCCAGCGCTTTCATCTTCGCCAGGCCCTGGCGCGCGGCCTCACGCGGGTTGGACGCCTGCTGGCTATTGCTCTGCGAAGCCACGTTACCGTAAGACAGGCCGCCATAGTTGTGGGTAGGCCCCACCAGGCCATCAAAATTCACtt                             >  6:275529/1‑149 (MQ=255)
                                                                                                                  gTCGGCCAGCGCTTTCATCTTCGCCAGGCCCTGGCGCGCGGCCTCACGCGGGTTGGACGCCTGCTGGCTATTGCTCTGCGAAGCCACGTTACCGTAAGACAGGCCGCCATAGTTGTGGGTAGGCCCCACCAGGCCATCAAAATTCACtt                             >  3:297906/1‑149 (MQ=255)
                                                                                                                      gCCAGCGCTTTCATCTTCGCCAGGCCCTGGCGCGCGGCCTCACGCGGGTTGGACGCCTGCTGGCTATTGCTCTGCGAAGCCACGTTACCGTAAGACAGGCCGCCATAGTTGTGGGTAGGCCCCACCAGGCCATCAAAATTCACTTCa                           >  5:296102/1‑147 (MQ=255)
                                                                                                                      gCCAGCGCTTTCATCTTCGCCAGGCCCTGGCGCGCGGCCTCACGCGGGTTGGACGCCTGCTGGCTATTGCTCTGCGAAGCCACGTTACCGTAAGACAGGCCGCCATAGTTGTGGGTAGGCCCCACCAGGCCATCAAAATTCACTTCa                           >  2:65723/1‑147 (MQ=255)
                                                                                                                      gCCAGCGCTTTCATCTTCGCCAGGCCCTGGCGCGCGGCCTCACGCGGGTTGGACGCCTGCTGGCTATTGCTCTGCGAAGCCACGTTACCGTAAGACAGGCCGCCATAGTTGTGGGTAGGCCCCACCAGGCCATCAAAATTCACTTCATa                         <  8:124272/149‑1 (MQ=255)
                                                                                                                      gCCAGCGCTTTCATCTTCGCCAGGCCCTGGCGCGCGGCCTCACGCGGGTTGGACGCCTGCTGGCTATTGCTCTGCGAAGCCACGTTACCGTAAGACAGGCCGCCATAGTTGTGGGTAGGCCCCACCAGGCCATCAAAATTCACTTCATa                         <  3:101350/149‑1 (MQ=255)
                                                                                                                       ccAGCGCTTTCATCTTCGCCAGGCCCTGGCGCGCGGCCTCACGCGGGTTGGACGCCTGCTGGCTATTGCTCTGCGAAGCCACGTTACCGTAAGACAGGCCGCCATAGTTGTGGGTAGGCCCCACCAGGCCATCAAAATTCACTTCATAg                        <  1:135551/149‑1 (MQ=255)
                                                                                                                         aGCGCTTTCATCTTCGCCAGGCCCTGGCGCGCGGCCTCACGCGGGTTGGACGCCTGCTGGCTATTGCTCTGCGAAGCCACGTTACCGTAAGACAGGCCGCCATAGTTGTGGGTAGGCCCCACCAGGCCATCAAAATTCACTTCATAGGa                      <  6:300524/149‑1 (MQ=255)
                                                                                                                                 tCATCTTCGCCAGGCCCTGGCGCGCGGCCTCACGCGGGTTGGACGCCTGCTGGCTATTGCTCTGCGAAGCCACGTTACCGTAAGACAGGCCGCCATAGTTGTGGGTAGGCCCCACCAGGCCATCAAAATTCACTTCATAGGATTTCATc               >  4:168208/1‑149 (MQ=255)
                                                                                                                                 tCATCTTCGCCAGGCCCTGGCGCGCGGCCTCACGCGGGTTGGACGCCTGCTGGCTATTGCTCTGCGAAGCCACGTTACCGTAAGACAGGCCGCCATAGTTGTGGGTAGGCCCCACCAGGCCATCAAAATTCACTTCATAGGATTTCATc               >  4:135747/1‑149 (MQ=255)
                                                                                                                                   aTCTTCGCCAGGCCCTGGCGCGCGGCCTCACGCGGGTTGGACGCCTGCTGGCTATTGCTCTGCGAAGCCACGTTACCGTAAGACAGGCCGCCATAGTTGTGGGTAGGCCCCACCAGGCCATCAAAATTCACTTCATAGGATTTCa                 <  7:58468/145‑1 (MQ=255)
                                                                                                                                   aTCTTCGCCAGGCCCTGGCGCGCGGCCTCACGCGGGTTGGACGCCTGCTGGCTATTGCTCTGCGAAGCCACGTTACCGTAAGACAGGCCGCCATAGTTGTGGGTAGGCCCCACCAGGCCATCAAAATTCACTTCATAGGATTTCa                 >  8:58468/1‑145 (MQ=255)
                                                                                                                                   aTCTTCGCCAGGCCCTGGCGCGCGGCCTCACGCGGGTTGGACGCCTGCTGGCTATTGCTCTGCGAAGCCACGTTACCGTAAGACAGGCCGCCATAGTTGTGGGTAGGCCCCACCAGGCCATCAAAATTCACTTCATAGGATTTCATCgg             >  5:22709/1‑149 (MQ=255)
                                                                                                                                    tCTTCGCCAGGCCCTGGCGCGCGGCCTCACGCGGGTTGGACGCCTGCTGGCTATTGCTCTGCGAAGCCACGTTACCGTAAGACAGGCCGCCATAGTTGTGGGTAGGCCCCACCAGGCCATCAAAATTCACTTCATAGGATTTCATCGGc            >  3:219758/1‑149 (MQ=255)
                                                                                                                                          cGCCAGGCCCTGGCGCGCGGCCTCACGCGGGTTGGACGCCTGCTGGCTATTGCTCTGCGAAGCCACGTTACCGTAAGACAGGCCGCCATAGTTGTGGGTAGGCCCCACCAGGCCATCAAAATTCACTTCATAGGATTTCATCGGCTa          <  1:126012/147‑1 (MQ=255)
                                                                                                                                            ccAGGCCCTGGCGCGCGGCCTCACGCGGGTTGGACGCCTGCTGGCTATTGCTCTGCGAAGCCACGTTACCGTAAGACAGGCCGCCATAGTTGTGGGTAGGCCCCACCAGGCCATCAAAATTCACTTCATAGGa                      <  6:122017/133‑1 (MQ=255)
                                                                                                                                            ccAGGCCCTGGCGCGCGGCCTCACGCGGGTTGGACGCCTGCTGGCTATTGCTCTGCGAAGCCACGTTACCGTAAGACAGGCCGCCATAGTTGTGGGTAGGCCCCACCAGGCCATCAAAATTCACTTCATAGGa                      >  5:122017/1‑133 (MQ=255)
                                                                                                                                            ccAGGCCCTGGCGCGCGGCCTCACGCGGGTTGGACGCCTGCTGGCTATTGCTCTGCGAAGCCACGTTACCGTAAGACAGGCCGCCATAGTTGTGGGTAGGCCCCACCAGGCCATCAAAATTCACTTCATAGGATTTCATCGGCTAGGCt      >  3:139583/1‑149 (MQ=255)
                                                                                                                                            ccAGGCCCTGGCGCGCGGCCTCACGCGGGTTGGACGCCTGCTGGCCATTGCTCTGCGAAGCCACGTTACCGTAAGACAGGCCGCCATAGTTGTGGGTAGGCCCCACCAGGCCATCAAAATTCACTTCATAGGATTTCATCGGCTAGGCt      >  4:200333/1‑149 (MQ=255)
                                                                                                                                              aGGCCCTGGCGCGCGGCCTCACGCGGGTTGGACGCCTGCTGGCTATTGCTCTGCGAAGCCACGTTACCGTAAGACAGGCCGCCATAGTTGTGGGTAGGCCCCACCAGGCCATCAAAATTCACTTCATAGGATTTCATCGGCTAGGCTcc    <  6:296102/149‑1 (MQ=255)
                                                                                                                                              aGGCCCTGGCGCGCGGCCTCACGCGGGTTGGACGCCTGCTGGCTATTGCTCTGCGAAGCCACGTTACCGTAAGACAGGCCGCCATAGTTGTGGGTAGGCCCCACCAGGCCATCAAAATTCACTTCATAGGATTTCATCGGCTAGGCTcc    <  6:22709/149‑1 (MQ=255)
                                                                                                                                                gCCCTGGCGCGCGGCCTCACGCGGGTTGGACGCCTGCTGGCTATTGCTCTGCGAAGCCACGTTACCGTAAGACAGGCCGCCATAGTTGTGGGTAGGCCCCACCCGGCCATCAAAATTCACTTCATAGGATTTCATCGGCTAGGCTCCGt  <  1:228792/149‑1 (MQ=255)
                                                                                                                                                gCCCTGGCGCGCGGCCTCACGCGGGGTGGACGCCTGCTGGCTATTGCTCTGCGAAGCCACGTTACCGTAAGACAGGCCGCCATAGTTGTGGGTAGGCCCCACCAGGCCATCAAAATTCACTTCATAGGATTTCATCGGCTAGGCTCCGt  <  4:248077/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
CTCGCGGGCGGCGCGCTGGATAACTTCGGCATCGCTGCCACTGAACCCCAGGCGGCGCAGCGCCGCCACGTCCGGGCGCTCTTGTGGCGCCAGAACACCCTGCTTGAAGCCCATGTCGGCCAGCGCTTTCATCTTCGCCAGGCCCTGGCGCGCGGCCTCACGCGGGTTGGACGCCTGCTGGCTATTGCTCTGCGAAGCCACGTTACCGTAAGACAGGCCGCCATAGTTGTGGGTAGGCCCCACCAGGCCATCAAAATTCACTTCATAGGATTTCATCGGCTAGGCTCCGT  >  NC_002947/5087047‑5087336

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: