Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,791,129 T→C 5.9% I197V (ATC→GTC) ‡ hemE ← uroporphyrinogen decarboxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,791,1290TC 5.9% 198.2 / 2.3 68I197V (ATC→GTC) ‡hemEuroporphyrinogen decarboxylase
Reads supporting (aligned to +/- strand):  ref base T (30/34);  new base C (2/2);  total (32/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 = 2.89e-01

ACCGGCACCTTGCGCCCCTCGTGCTCGCGGATCAGGCCGCTGACGATCTTGCGCATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCATGGCCTTGGTCTTGCGGAAGTCCTTCGACGAGCCGCCTTCGACCATGTAGGTAGCCAGGGTCCACG  >  NC_002947/5790984‑5791275
                                                                                                                                                 |                                                                                                                                                   
aCCGGCACCTTGCGCCCCTCGTGCTCGCGGATCAGGCCGCTGACGATCTTGCGCATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCt                                                                                                                                                   <  2:117353/148‑1 (MQ=255)
aCCGGCACCTTGCGCCCCTCGTGCTCGCGGATCAGGCCGCTGACGATCTTGCGCATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCt                                                                                                                                                   <  4:158252/148‑1 (MQ=255)
aCCGGCACCTTGCGCCCCTCGTGCTCGCGGATCAGGCCGCTGACGATCTTGCGCATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTg                                                                                                                                                  >  6:311236/1‑149 (MQ=255)
 ccGGCACCTTGCGCCCCTCGTGCTCGCGGTTCAGGCCGCTGACGATCTTGCGCATGTAGGCCAGGGAGCACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTgg                                                                                                                                                 <  8:134646/149‑1 (MQ=255)
  cGGCACCTTGCGCCCCTCGTGCTCGCGGATCAGGCCGCTGACGATCTTGCGCATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGc                                                                                                                                                >  6:233726/1‑149 (MQ=255)
  cGGCACCTTGCGCCCCTCGTGCTCGCGGATCAGGCCGCTGACGATCTTGCGCATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGc                                                                                                                                                >  4:58667/1‑149 (MQ=255)
       ccTTGCGCCCCTCGTGCTCGCGGATCAGGCCGCTGACGATCTTGCGCATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGATCCTGGCCggcg                                                                                                                                           <  4:35799/149‑4 (MQ=255)
        catGCGCCCCTCGTGCTCGCGGATCAGGCCGCTGACGATCTTGCGCATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTaa                                                                                                                                          <  8:290624/147‑1 (MQ=255)
        cTTGCGCCCCTCGTGCTCGCGGATCAGGCCGCTGACGATCTTGCGCATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGTGCGCCGGCCAGGATCTGGCCGtt                                                                                                                                            >  4:11595/1‑147 (MQ=255)
        cTTGCGCCCCTCGTGCTCGCGGATCAGGCCGCTGACGATCTTGCGCATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGtt                                                                                                                                            >  3:215744/1‑147 (MQ=255)
               cccTCGTGCTCGCGGATCAGGCCGCTGACGATCTTGCGCATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCt                                                                                                                                   <  5:326884/149‑1 (MQ=255)
                ccTCGTGCTCGCGGATCAGGCCGCTGACGATCTTGCGCATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTg                                                                                                                                 <  2:81870/149‑1 (MQ=255)
                       cTCGCGGATCAGGCCGCTGACGATCTTGCGCATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGAcc                                                                                                                           <  5:311236/148‑1 (MQ=255)
                            ggATCAGGCCGCTGACGATCTTGCGCATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTgg                                                                                                                     <  6:279652/149‑1 (MQ=255)
                                  ggCCGCTGACGATCTTGCGCATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGc                                                                                                               <  1:19102/149‑1 (MQ=255)
                                     cGCTGACGATCTTGCGCATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCGAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTg                                                                                                            <  7:211644/149‑1 (MQ=255)
                                      gCTGACGATCTTGCGCATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGt                                                                                                           >  7:48063/1‑149 (MQ=255)
                                       cTGACGATCTTGCGCATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGTCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTc                                                                                                          >  3:72645/1‑149 (MQ=255)
                                       cTGACGATCTTGCGCATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTc                                                                                                          <  3:11595/149‑1 (MQ=255)
                                          aCGATCTTGCGCATGTAGGCCAGGGAGAACTCCTGATAGGTCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGa                                                                                                                         <  6:192308/131‑1 (MQ=255)
                                          aCGATCTTGCGCATGTAGGCCAGGGAGAACTCCTGATAGGTCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGa                                                                                                                         >  5:192308/1‑131 (MQ=255)
                                            gATCTTGCGCATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCcag                                                                                                       >  7:221481/1‑147 (MQ=255)
                                              tCTTGCGCATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCcagcag                                                                                                    >  8:14858/1‑148 (MQ=255)
                                                   cgcATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTg                                                                                              >  5:8709/1‑149 (MQ=255)
                                                    gcATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGc                                                                                             <  8:150851/149‑1 (MQ=255)
                                                     cATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCa                                                                                            <  2:310094/149‑1 (MQ=255)
                                                            gCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCa                                                                                                                 >  4:203952/1‑121 (MQ=255)
                                                            gCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCa                                                                                                                 <  3:203952/121‑1 (MQ=255)
                                                                 ggAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGt                                                                                >  4:107211/1‑149 (MQ=255)
                                                                 ggAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGt                                                                                >  2:285575/1‑149 (MQ=255)
                                                                       cTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCg                                                                           >  1:19780/1‑148 (MQ=255)
                                                                        tCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCg                                                                           >  2:314126/1‑147 (MQ=255)
                                                                        tCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCg                                                                           >  8:259577/1‑147 (MQ=255)
                                                                            gATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGc                                                                      >  3:111342/1‑148 (MQ=255)
                                                                                 gccgccgccGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCATGGc                                                                >  4:283864/1‑149 (MQ=255)
                                                                                 gccgccgccGACAGGTTGCCGCCCAAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCAGGTGACCGACTGGGCCAGCTTGTGCAGCAGCAGGTGGCGGGCTTGCGGGTTGTCGTAGGCCATGGc                                                                >  7:49490/1‑149 (MQ=255)
                                                                                  ccgccgccgACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACTGCCTGGGCGCCGGCCAGGATCTGGCAGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTAGTAGGCCATGGcc                                                               >  2:313313/1‑149 (MQ=255)
                                                                                  ccgccgccgACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCATg                                                                  >  8:236862/1‑146 (MQ=255)
                                                                                  ccgccgccgACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCATg                                                                  <  7:236862/146‑1 (MQ=255)
                                                                                  ccgccgccgACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCATGGcc                                                               >  6:170543/1‑149 (MQ=255)
                                                                                        ccgACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCATGGCCTTGGTc                                                         <  4:215744/149‑1 (MQ=255)
                                                                                           aCAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGTCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCATGGCCTTGGTCTTg                                                      <  4:72645/149‑1 (MQ=255)
                                                                                                   ccgccCCAGTGGGCGAAGTTCTGTTCGGCCTGGGCGGCGGCCAGGTTCTGGCCGTTATGGTTTCTGGTGACCGACTGGGCCAGCTTGTCCAGGTGCATGTGCAGGGCTTGCGGGTTGTCTTAGGTTATGGCCTTTGTCTTGCGGATGTc                                              <  8:49490/149‑1 (MQ=255)
                                                                                                      ccccAGGTGTCGAAGATCTGTACGGCTTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCATGGCCTTGGTCTTGCGGAAGTCCtt                                           >  8:194875/1‑149 (MQ=255)
                                                                                                       cccAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCATGGCCTTGGTCTTGCGGAAGTCCTTc                                          <  1:314126/149‑1 (MQ=255)
                                                                                                       cccAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCATGGCCTTGGTCTTGCGGAAGTCCTTc                                          <  7:259577/149‑1 (MQ=255)
                                                                                                       cccAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCATGGCCTTGGTCTTGCGGAAGTCCTTc                                          <  7:14858/149‑1 (MQ=255)
                                                                                                            gtgtCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCATGGCCTTGGTCTTGCGGAAGTCCTTcgac                                       >  7:97465/1‑147 (MQ=255)
                                                                                                                 gAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCATGGCCTTGGTCTTGCGGAAGTCCTTCGACGAgccgc                                <  1:193847/149‑1 (MQ=255)
                                                                                                                     aTCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCATGGCCTTGGTCTTGCGGAAGTCCTTCGACGAGCCGCCTTc                            >  6:280751/1‑149 (MQ=255)
                                                                                                                          tACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCATGGCCTTGGTCTTGCGGAAGTCCTTCGACGAGCCGCCTTCGAcc                        <  1:285575/148‑1 (MQ=255)
                                                                                                                          tACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCATGGCCTTGGTCTTGCGGAAGTCCTTCGACGAGCCGCCTTCGACCa                       >  8:23182/1‑149 (MQ=255)
                                                                                                                          tACGGCCTGGGCGCCGGCCAGGATCAGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCATGGCCTTGGTCTTGCGGAAGTCCTTCGACGAGCCGCCTTCg                           <  4:145140/145‑1 (MQ=255)
                                                                                                                           aCGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCATGGCCTTGGTCTTGCGGAAGTCCTTCGACGAGCCGCCTTc                            >  1:229742/1‑144 (MQ=255)
                                                                                                                           aCGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCATGGCCTTGGTCTTGCGGAAGTCCTTCGACGAGCCGCCTTc                            <  2:229742/143‑1 (MQ=255)
                                                                                                                            cGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCATGGCCTTGGTCTTGCGGAAGTCCTTCGACGAGCCGCCTTCCACCATg                     >  4:322496/1‑149 (MQ=255)
                                                                                                                              gCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCATGGCCTTGGTCTTGCGGAAGTCCTTCGACGAGCCGCCTTCGACCATgta                   <  5:280751/149‑1 (MQ=255)
                                                                                                                               ccTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCATGGCCTTGGTCTTGCGGAAGTCCTTCGACGAGCCGCCTTCGACCATgtag                  <  2:19780/149‑1 (MQ=255)
                                                                                                                                cTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCATGGCCTTGGTCTTGCGGAAGTCCTTCGACGAGCCGCCTTCGACCATgtagg                 <  5:233726/149‑1 (MQ=255)
                                                                                                                                cTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCATGGCCTTGGTCTTGCGGAAGTCCTTCGACGAGCCGCCTTCGACCATgtagg                 <  3:283864/149‑1 (MQ=255)
                                                                                                                                cTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCATAGGCCATGGCCTTGGTCTTGCGGAAGTCCTTCGACGAGCCGCCTTCGACCATgtagg                 <  8:48063/149‑1 (MQ=255)
                                                                                                                                   cttaacGGCCAGATCCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGTAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCATGGCCTTGGTCTTGCGGAAGTCCTTCGACGAGCCGCCTTCGACCATgtaggtag              >  6:75895/6‑149 (MQ=255)
                                                                                                                                    ttaacGGCCAGATCCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCATGGCCTTGGTCTTGCGGAAGTCCTTCGACGAGCCGCCTTCGACCATgta                   <  7:260524/139‑1 (MQ=255)
                                                                                                                                    ttaacGGCCAGATCCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCATGGCCTTGGTCTTGCGGAAGTCCTTCGACGAGCCGCCTTCGACCATgta                   >  8:260524/5‑143 (MQ=255)
                                                                                                                                    gcgcCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCATGGCCTTGGTCTTGCGGAAGTCCTTCGACGAGCCGCCTTCGACCATGTAGGTAGc             <  6:8709/149‑1 (MQ=255)
                                                                                                                                        cGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCATGGCCTTGGTCTTGCGGAAGTCCTTCGACGAGCCGCCTTCGACCATGTAGGTAGCCAgg         <  7:194875/149‑1 (MQ=255)
                                                                                                                                         ggCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCATGGCCTTGGTCTTGCGGAAGTCCTTCGACGAGCCGCCTTCGACCATGTAGGTAGCCAggg        <  1:313313/149‑1 (MQ=255)
                                                                                                                                               gATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCATGGCCTTGGTCTTGCGGAAGTCCTTCGACGAGCCGCCTTCGACCATGTAGGTAGCCAGGGTCCACg  >  5:164738/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
ACCGGCACCTTGCGCCCCTCGTGCTCGCGGATCAGGCCGCTGACGATCTTGCGCATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCATGGCCTTGGTCTTGCGGAAGTCCTTCGACGAGCCGCCTTCGACCATGTAGGTAGCCAGGGTCCACG  >  NC_002947/5790984‑5791275

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: