Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,791,128 A→T 6.8% I197N (ATC→AAC) ‡ hemE ← uroporphyrinogen decarboxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,791,1280AT 6.8% 203.7 / 3.2 68I197N (ATC→AAC) ‡hemEuroporphyrinogen decarboxylase
Reads supporting (aligned to +/- strand):  ref base A (30/33);  new base T (2/3);  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 = 8.14e-01

GATCACCGGCACCTTGCGCCCCTCGTGCTCGCGGATCAGGCCGCTGACGATCTTGCGCATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCATGGCCTTGGTCTTGCGGAAGTCCTTCGACGAGCCGCCTTCGACCATGTAGGTAGCCAGGGTCCACG  >  NC_002947/5790980‑5791275
                                                                                                                                                    |                                                                                                                                                    
gATCACCGGCACCTTGCGCCCCTCGTGCTCGCGGATCAGGCCGCTGACGATCTTGCGCATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGa                                                                                                                                                      >  2:193847/1‑149 (MQ=255)
    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                                                                                                                 <  3:203952/121‑1 (MQ=255)
                                                                gCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCa                                                                                                                 >  4:203952/1‑121 (MQ=255)
                                                                     ggAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGt                                                                                >  2:285575/1‑149 (MQ=255)
                                                                     ggAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGt                                                                                >  4:107211/1‑149 (MQ=255)
                                                                           cTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCg                                                                           >  1:19780/1‑148 (MQ=255)
                                                                            tCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCg                                                                           >  8:259577/1‑147 (MQ=255)
                                                                            tCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCg                                                                           >  2:314126/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                                                                  <  7:236862/146‑1 (MQ=255)
                                                                                      ccgccgccgACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCATg                                                                  >  8:236862/1‑146 (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                                          <  7:14858/149‑1 (MQ=255)
                                                                                                           cccAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCATGGCCTTGGTCTTGCGGAAGTCCTTc                                          <  7:259577/149‑1 (MQ=255)
                                                                                                           cccAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCATGGCCTTGGTCTTGCGGAAGTCCTTc                                          <  1:314126/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                 <  3:283864/149‑1 (MQ=255)
                                                                                                                                    cTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCATGGCCTTGGTCTTGCGGAAGTCCTTCGACGAGCCGCCTTCGACCATgtagg                 <  5:233726/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)
                                                                                                                                                    |                                                                                                                                                    
GATCACCGGCACCTTGCGCCCCTCGTGCTCGCGGATCAGGCCGCTGACGATCTTGCGCATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCATGGCCTTGGTCTTGCGGAAGTCCTTCGACGAGCCGCCTTCGACCATGTAGGTAGCCAGGGTCCACG  >  NC_002947/5790980‑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: