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% 201.8 / 2.3 67I197V (ATC→GTC) ‡hemEuroporphyrinogen decarboxylase
Reads supporting (aligned to +/- strand):  ref base T (30/33);  new base C (2/2);  total (32/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.86e-01

ACCGGCACCTTGCGCCCCTCGTGCTCGCGGATCAGGCCGCTGACGATCTTGCGCATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCATGGCCTTGGTCTTGCGGAAGTCCTTCGACGAGCCGCCTTCGACCATGTAGGTAGCCAGGGTC  >  NC_002947/5790984‑5791271
                                                                                                                                                 |                                                                                                                                              
aCCGGCACCTTGCGCCCCTCGTGCTCGCGGATCAGGCCGCTGACGATCTTGCGCATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTg                                                                                                                                             <  7:26880/149‑1 (MQ=255)
  cGGCACCTTGCGCCCCTCGTGCTCGCGGATCAGGCCGCTGACGATCTTGCGCATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGc                                                                                                                                           <  4:312069/149‑1 (MQ=255)
    gCACCTTGCGCCCCTCGTGCTCGCGGATCAGGCCGCTGACGATCTTGCGCATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGATCCTGGCCg                                                                                                                                         <  8:192855/149‑1 (MQ=255)
           gcgcCCCTCGTGCTCGCGGATCAGGCCGCTGACGATCTTGCGCATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAgg                                                                                                                                   <  7:104657/148‑1 (MQ=255)
            cgcCCCTCGTGCTCGCGGATCAGGCCGCTGACGATCTTGCGCATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGt                                                                                                                                  <  7:348617/148‑1 (MQ=255)
                ccTCGTGCTCGCGGATCAGGCCGCTGACGATCTTGCGCATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTg                                                                                                                             <  2:181039/149‑1 (MQ=255)
                    gTGCTCGCGGATCAGGCCGCTGACGATCTTGCGCATGTAGTCCAGGGATAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCTCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTg                                                                                                                          <  8:39622/148‑1 (MQ=255)
                    gTGCTCGCGGATCAGGCCGCTGACGATCTTGCGCATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGt                                                                                                                           <  3:58232/147‑1 (MQ=255)
                    gTGCTCGCGGATCAGGCCGCTGACGATCTTGCGCATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCAGCCAGGATCTGGCCGTTAAGGTAGCTGGt                                                                                                                           <  5:67177/147‑1 (MQ=255)
                     tGCTCGCGGATCAGGCCGCTGACGATCTTGCGCATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGAc                                                                                                                        <  3:117323/149‑1 (MQ=255)
                        tCGCGGATCAGGCCGCTGACGATCTTGCGCATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCAGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCAGGCCGTTTAGGTAGCGGGTGACCGa                                                                                                                     >  8:276075/1‑149 (MQ=255)
                         cgcgGATCAGGCCGCTGACGATCTTGCGCATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGAc                                                                                                                    >  7:348295/1‑149 (MQ=255)
                         cgcgGATCAGGCCGCTGACGATCTTGCGCATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGATCCTGGCCgg                                                                                                                                        >  1:23661/1‑128 (MQ=255)
                         cgcgGATCAGGCCGCTGACGATCTTGCGCATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGATCCTGGCCgg                                                                                                                                        <  2:23661/129‑2 (MQ=255)
                          gcgGATCAGGCCGCTGACGATCTTGCGCATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGa                                                                                                                     <  8:167045/147‑1 (MQ=255)
                          gcgGATCAGGCCGCTGACGATCTTGCGCATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGAc                                                                                                                    <  5:146847/148‑1 (MQ=255)
                                 aGGCCGCTGACGATCTTGCGCATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCGTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAg                                                                                                            <  6:303028/149‑1 (MQ=255)
                                     cGCTGACGATCTTGCGCATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTg                                                                                                        >  8:142591/1‑149 (MQ=255)
                                     cGCTGACGATCTTGCGCATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTg                                                                                                        <  8:184306/149‑1 (MQ=255)
                                      gCTGACGATCTTGCGCATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGt                                                                                                       >  8:167795/1‑149 (MQ=255)
                                       cTGACGATCTTGCGCATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTgg                                                                                                                 <  3:129223/138‑1 (MQ=255)
                                       cTGACGATCTTGCGCATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTgg                                                                                                                 >  4:129223/1‑138 (MQ=255)
                                        tGACGATCTTGCGCATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTgg                                                                                                                            <  8:101635/126‑1 (MQ=255)
                                        tGACGATCTTGCGCATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTgg                                                                                                                            >  7:101635/1‑126 (MQ=255)
                                            gATCTTGCGCATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCcag                                                                                                   >  6:182402/1‑147 (MQ=255)
                                                  gcgcATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCcagcagcag                                                                                             <  8:324665/147‑1 (MQ=255)
                                                  gcgcATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAgg                                                                                            <  8:348295/148‑1 (MQ=255)
                                                   cgcATGGAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGTATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTg                                                                                          <  8:120442/149‑1 (MQ=255)
                                                        gTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCca                                                                                                    >  6:139962/1‑134 (MQ=255)
                                                        gTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCca                                                                                                    <  5:139962/134‑1 (MQ=255)
                                                             ccAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGc                                                                                >  6:144446/1‑149 (MQ=255)
                                                               aGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCgg                                                                              >  6:384678/1‑149 (MQ=255)
                                                                 ggAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCtt                                                                                  <  8:198693/143‑1 (MQ=255)
                                                                 ggAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCtt                                                                                  >  7:198693/1‑143 (MQ=255)
                                                                  gagaACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGtt                                                                           <  2:176282/149‑1 (MQ=255)
                                                                        tCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGAGCGACTGGGCCAGCTTGTGCTGGTGGTGGTGGTGGGGTTGGGGGTTGTGGtt                                                                     >  2:149816/1‑148 (MQ=255)
                                                                        tCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCggg                                                                             >  4:258587/1‑141 (MQ=255)
                                                                        tCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCggg                                                                             <  3:258587/141‑1 (MQ=255)
                                                                          cTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAgg                                                                   >  4:177549/1‑149 (MQ=255)
                                                                          cTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCAGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGATAAGGTAGCTGGTGACCGACTCGGCCAGCTTGTCCAGCAGCCGGTGCACGGCTTGCGGGTTGTCGTAgg                                                                   >  7:62898/1‑149 (MQ=255)
                                                                           tGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGc                                                                  >  3:123879/1‑149 (MQ=255)
                                                                                ggCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCATgg                                                             <  7:167795/149‑1 (MQ=255)
                                                                                  ccgccgccgACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCATGGcc                                                           >  5:100256/1‑149 (MQ=255)
                                                                                    gccgccgACAGGTTGCCGCTCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCa                                                                <  1:146043/142‑1 (MQ=255)
                                                                                    gccgccgACAGGTTGCCGCTCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCa                                                                >  2:146043/1‑142 (MQ=255)
                                                                                     ccgccgACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCATGGCCTTg                                                        <  7:211112/149‑1 (MQ=255)
                                                                                           aCAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCATGGCCTTGGTCTTg                                                  >  1:338496/1‑149 (MQ=255)
                                                                                            cAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCATGGCCTTGGTCTTGc                                                 <  7:132911/149‑1 (MQ=255)
                                                                                               gTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCATGGCCTTGGTCTTGCgg                                               >  5:105545/1‑148 (MQ=255)
                                                                                                ttGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGGTTGTCCAGCAGCAGGTGCAGGGCTTTCGGGTTGTCGTAGGCCATGGCCTTGGTCTTGCGGaa                                             <  2:340315/149‑1 (MQ=255)
                                                                                                     gccCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCggg                                                                             >  3:245491/1‑112 (MQ=255)
                                                                                                     gccCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCggg                                                                             <  4:245491/112‑1 (MQ=255)
                                                                                                       cccTGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTACGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGTTGCAGGGCTTGCGGGTTGTCGTAGGCCATGGCCTTGGTCTTGCGGAAGTCCTTc                                      <  3:222954/149‑1 (MQ=255)
                                                                                                       cccAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCATGGCCTTGGTCTTGCGGAAGTCCTTc                                      <  7:276075/149‑1 (MQ=255)
                                                                                                         cAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCATGGCCTTGGTCTTGCGGAAGTCCTTcga                                    <  4:123879/149‑1 (MQ=255)
                                                                                                                 gAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCCGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTAGTAGGCCATGGCCTTGGTCTTGAGGAAGTCCTTCGACGAgccgc                            >  4:49448/1‑149 (MQ=255)
                                                                                                                 gAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCATGGCCTTGGTCTTGCGGAAGTCCTTCGACGAgccgc                            >  5:100162/1‑149 (MQ=255)
                                                                                                                 gAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCATGGCCTTGGTCTTGCGGAAGTCCTTCGACGAgccgc                            >  7:13809/1‑149 (MQ=255)
                                                                                                                       cTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCATGGCCTTGGTCTTGCGGAAGTCCTTCGACGAGCCGCCTTCGa                      >  3:20795/1‑149 (MQ=255)
                                                                                                                              gCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCATGGCCTTGGTCTTGCGGAAGTCCTTCGACGAGCCGCCTTCGACCATgta               <  5:144446/149‑1 (MQ=255)
                                                                                                                               ccTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCATGGCCTTGGTCTTGCGGAAGTCCTTCGACGAGCCGCCTTCGACCATgtag              <  6:105545/149‑1 (MQ=255)
                                                                                                                               ccTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCATGGCCTTGGTCTTGCGGAAGTCCTTCGACGAGCCGCCTTCGACCATgtag              >  6:356398/1‑149 (MQ=255)
                                                                                                                                 tGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCATGGCCTTGGTCTTGCGGAAGTCCTTCGACGAGCCGCCTTCGACCATgtagg             >  2:359629/1‑148 (MQ=255)
                                                                                                                                       acGGCCAGATCCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCATGGCCTTGGTCTTGCGGAAGTCCTTCGACGAGCCGCCTTCGACCATGTAGGTAGCCAg      >  7:105568/2‑149 (MQ=255)
                                                                                                                                         ggCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCATGGCCTAGGTCTTGCGGAAGTCCTTCGACGAGCCGCCTTCGACCATGTAGGTAGCCAggg    >  2:145944/1‑149 (MQ=255)
                                                                                                                                            cAGGATCTGGCCGTTAAGGTAGCTGGTGCCCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCATGGCCTTGGTCTTGCGGAAGTCCTTCGACGAGCCGCCTTCGACCATGTAGGTAGCCAGGGTc  >  8:185038/1‑148 (MQ=255)
                                                                                                                                               gATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCATGGCCTTGGTCTTGCGGAAGTCCTTcg                                     >  6:321007/1‑110 (MQ=255)
                                                                                                                                                 |                                                                                                                                              
ACCGGCACCTTGCGCCCCTCGTGCTCGCGGATCAGGCCGCTGACGATCTTGCGCATGTAGGCCAGGGAGAACTCCTGATAGGCCGCCGCCGACAGGTTGCCGCCCCAGGTGTCGAAGATCTGTACGGCCTGGGCGCCGGCCAGGATCTGGCCGTTAAGGTAGCTGGTGACCGACTGGGCCAGCTTGTCCAGCAGCAGGTGCAGGGCTTGCGGGTTGTCGTAGGCCATGGCCTTGGTCTTGCGGAAGTCCTTCGACGAGCCGCCTTCGACCATGTAGGTAGCCAGGGTC  >  NC_002947/5790984‑5791271

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