Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 239,942 G→T 7.2% E79D (GAG→GAT)  PP_0188 → uroporphyrin‑III C‑methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947239,9420GT 7.2% 217.4 / 5.6 69E79D (GAG→GAT) PP_0188uroporphyrin‑III C‑methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (31/33);  new base T (3/2);  total (34/35)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.95e-01

CCACCCCACCTGCCAAGCGCTCCGGCAGTGGCCTGGCAGCGCTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCA‑GCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGCCCAGTTGCAGGGCGACCAGCAGCGTCTCAGC  >  NC_002947/239797‑240083
                                                                                                                                                 |                                                                                                                                              
ccaccccaccTGCCAAGCGCTCCGGCAGTGGCCTGGCAGCGCTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAAGTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGgcg                                                                                                                                             >  5:310405/1‑149 (MQ=255)
 caccccaccTGCCAAGCGCTCCGGCAGTGGCCTGGCAGCGCTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGgcgc                                                                                                                                            >  6:291852/1‑149 (MQ=255)
 caccccaccTGCCAAGCGCTCCGGCAGTGGCCTGGCAGCGCTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGgcgc                                                                                                                                            >  6:270548/1‑149 (MQ=255)
 caccccaccTGCCAAGCGCTCCGGCAGTGGCCTGGCAGCGCTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGgcgc                                                                                                                                            >  1:300465/1‑149 (MQ=255)
 caccccaccTGCCAAGCGCTCCGACAGTGGCCTGGCAGCGCTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGCGAGTTCAACCAGGGCCAGCACATCGAGgcgc                                                                                                                                            >  6:256354/1‑149 (MQ=255)
  accccaccTGCCAAGCGCTCCGGCAGTGGCCTGGCAGCGCTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCt                                                                                                                                           >  2:268115/1‑149 (MQ=255)
    cccaccTGCCAAGCGCTCCGGCAGTGGCCTGGCAGCGCTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGa                                                                                                                                         <  2:300465/149‑1 (MQ=255)
            ccAAGCGCTCCGGCAGTGGCCTGGCAGCGCTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGa                                                                                                                                         >  6:178584/1‑141 (MQ=255)
            ccAAGCGCTCCGGCAGTGGCCTGGCAGCGCTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGa                                                                                                                                         <  5:178584/141‑1 (MQ=255)
             cAAGCGCTCCGGCAGTGGCCTGGCAGCGCTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCA‑gcgcg                                                                                                                               <  5:275962/149‑1 (MQ=255)
             cAAGCGCTCCGGCAGTGGCCTGGCAGCGCTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCA‑gcgcg                                                                                                                               >  6:312571/1‑149 (MQ=255)
              aaGCGCGCCAGCAGTGGCCTGGCAGCGCTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGGGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCA‑gcgcgc                                                                                                                              <  5:23738/149‑1 (MQ=255)
               aGCGCTCCGGCAGTGGCCTGGCAGCGCTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCAA‑GCGCGcc                                                                                                                             >  2:245317/1‑149 (MQ=255)
                  gcTCCGGCAGTGGCCTGGCAGCGCTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCA‑GCGCGCCGAg                                                                                                                          >  5:37251/1‑149 (MQ=255)
                      cGGCAGTGGCCTGGCAGCGCTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCA‑GCGCGCCGAGgcg                                                                                                                       <  2:225877/148‑1 (MQ=255)
                         cAGTGGCCTGGCAGCGCTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCA‑GCGCCTCGGCGCGCTgg                                                                                                                   <  3:308483/149‑2 (MQ=255)
                          aGTGGCCTGGCAGCGCTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCA‑GCGCGCCGTGGCGCTgca                                                                                                                  <  1:150427/149‑1 (MQ=255)
                              gCCTGGCAGCGCTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCA‑GCGCGCCGAGGCGCTgca                                                                                                                  <  8:298041/145‑1 (MQ=255)
                              gCCTGGCAGCGCTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCA‑GCGCGCCGAGGCGCTgca                                                                                                                  >  7:298041/1‑145 (MQ=255)
                              gCCTGGCAGCGCTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCA‑GCGCGCCGAGGCGCTGCAGCaa                                                                                                              <  3:130037/149‑1 (MQ=255)
                                   gCAGCGCTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGgcgc                                                                                                                                            <  8:327880/115‑1 (MQ=255)
                                   gCAGCGCTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGgcgc                                                                                                                                            >  7:327880/1‑115 (MQ=255)
                                   gCAGCGCTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCA‑GCGCGCCGAGGCGCTGCAGCAACGt                                                                                                           <  7:10521/147‑1 (MQ=255)
                                    cAGCGCTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCA‑GCGCGCCGAGGCGCTGCAGCAACGTGaa                                                                                                        <  5:270548/149‑1 (MQ=255)
                                    cAGCGCTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCA‑GCGCGCCGAGGCGCTGCAGCAACGTGaa                                                                                                        <  6:37251/149‑1 (MQ=255)
                                      gcgcTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGCGGCGCTGACCCA‑GCGCGCCGAGGCGCTGCAGCAACGTGAAca                                                                                                      <  1:250575/149‑1 (MQ=255)
                                        gcTGGCCTTGCTGCTTGTTGCGGCCGGGGTGGCGGTTGTTGGCTGGGGTGTCTGGCATGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCA‑GCGCGCCGAGGCTCTGCATCAACGTGAAcagc                                                                                                    <  8:246534/149‑1 (MQ=255)
                                        gcTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCA‑GCGCGCCGAGGCGCTGCAGCAACGTGAAcagc                                                                                                    <  2:42635/149‑1 (MQ=255)
                                            gCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCA‑GCGCGCCGAGGCGCTGCAGCAACGTGaa                                                                                                        >  3:103472/1‑141 (MQ=255)
                                            gCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCA‑GCGCGCCGAGGCGCTGCAGCAACGTGaa                                                                                                        <  4:103472/141‑1 (MQ=255)
                                            gCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCA‑GCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGa                                                                                                 <  1:245317/148‑1 (MQ=255)
                                               ttGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCA‑GCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAg                                                                                             >  1:317559/1‑149 (MQ=255)
                                                            cGGCCGGGGTTGCGGTAGGTGGCTGGGGTGTCTTGGAGGTGCGTCAACTTCAGGGCAGTGTGTTCAACCAGGGCCAGCACATCGAGGCGGTGACCCATCCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGc                                                                                 <  1:268115/149‑1 (MQ=255)
                                                               ccGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCA‑GCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGt                                                                             >  2:106663/1‑149 (MQ=255)
                                                                         ggTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCA‑GCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAgccg                                                                   >  2:204946/1‑149 (MQ=255)
                                                                         ggTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCA‑GCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAgccg                                                                   >  2:260671/1‑149 (MQ=255)
                                                                             ggTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGCGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCA‑GCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCa                                                                <  5:256354/148‑1 (MQ=255)
                                                                                      ggTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCA‑GCGCGCCGAGGCGCTGCAGCAACGTGAAc                                                                                                       >  4:223683/1‑100 (MQ=255)
                                                                                      ggTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCA‑GCGCGCCGAGGCGCTGCAGCAACGTGAAc                                                                                                       <  3:223683/100‑1 (MQ=255)
                                                                                        tgtCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCA‑GCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTagaag                                                    <  2:188914/149‑1 (MQ=255)
                                                                                          tCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCA‑GCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGAc                                                  >  7:25932/1‑149 (MQ=255)
                                                                                              gCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCA‑GCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGAccgc                                               <  5:312571/148‑1 (MQ=255)
                                                                                                aGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCA‑GCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCgg                                            >  1:176170/1‑149 (MQ=255)
                                                                                                     cGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCA‑GCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGttt                                       <  5:291852/149‑1 (MQ=255)
                                                                                                     cGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCA‑GCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGttt                                       >  1:227291/1‑149 (MQ=255)
                                                                                                              cAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCA‑GCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGcc                                                                        <  5:274942/107‑1 (MQ=255)
                                                                                                              cAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCA‑GCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGcc                                                                        >  6:274942/1‑107 (MQ=255)
                                                                                                                    cgacaGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCA‑GCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGCCCAGTTGCa                        <  3:281398/145‑1 (MQ=255)
                                                                                                                      tGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCA‑GCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGCCCAGTTGCAgg                      <  5:288992/149‑1 (MQ=255)
                                                                                                                         gTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCA‑GCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGCCCAGTTGCAGGGc                    >  7:110901/1‑148 (MQ=255)
                                                                                                                          ttCAACCAGGGCCAGCACATCGAGGCGCTGACCCA‑GCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGACAGCGACCTAGAAGACCGCCGGCGTTTGGTAGCCCAGTTGCAGGGCGa                  >  7:312588/1‑149 (MQ=255)
                                                                                                                            cAACCAGGGCCAGCACATCGAGGCGCTGACCCA‑GCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGCCCAGTTGCAGGGCGAc                 <  2:317559/148‑1 (MQ=255)
                                                                                                                            cAACCAGGGCCAGCACATCGAGGCGCTGACCCA‑GCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGCCCAGTTGCAGGGCGAc                 <  1:260671/148‑1 (MQ=255)
                                                                                                                              aCCAGGGCCAGCACATCGAGGCGCTGACCCA‑GCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGCCCAGTTGCAGGGCGACcag              >  8:303941/1‑149 (MQ=255)
                                                                                                                               ccAGGGCCAGCACATCGAGGCGCTGACCCA‑GCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCATTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGCCCAGTTGc                         <  7:221269/137‑1 (MQ=255)
                                                                                                                               ccAGGGCCAGCACATCGAGGCGCTGACCCA‑GCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGCCCAGTTGc                         >  8:221269/1‑137 (MQ=255)
                                                                                                                                 tgggTCAGCGCCTCGATGTGCTGGCCCA‑GCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGCCCAGTTGCAGGGCGACcagcag           >  8:116259/2‑149 (MQ=255)
                                                                                                                                   ggtcAGCGCCTCGATGTGCTGGCCCA‑GCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGc                                                          <  1:3058/97‑1 (MQ=255)
                                                                                                                                   ggtcAGCGCCTCGATGTGCTGGCCCA‑GCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGc                                                          >  2:3058/4‑100 (MQ=255)
                                                                                                                                    gtcAGCGCCTCGATGTGCTGACCCA‑GCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGCCCAGTTGCAGGGCGACCAGCAGCGt        >  7:60654/3‑149 (MQ=255)
                                                                                                                                        gCACATCGAGGCGCTGACCCA‑GCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTgg                                     <  2:4419/116‑1 (MQ=255)
                                                                                                                                        gCACATCGAGGCGCTGACCCA‑GCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTgg                                     >  1:4419/1‑116 (MQ=255)
                                                                                                                                          gcctCGATGTGCTGGCCCA‑GCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGCCCAGTTGCAGGGCGACCAGCAGCGTCTcagc  <  7:116259/146‑1 (MQ=255)
                                                                                                                                           caTCGAGGCGCTGACCCA‑GCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGCCCAGTTGCAGGGCGACCAGCAGCGTCTcagc  >  3:327343/1‑148 (MQ=255)
                                                                                                                                            aTCGAGGCGCTGACCCA‑GCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGCCCAGTTGCAgg                      >  4:328617/1‑127 (MQ=255)
                                                                                                                                            aTCGAGGCGCTGACCCA‑GCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGCCCAGTTGCAgg                      <  3:328617/127‑1 (MQ=255)
                                                                                                                                            aTCGAGGCGCTGACCCA‑GCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGCCCAGTTGCAGGGCGACcagcag           <  3:224596/138‑1 (MQ=255)
                                                                                                                                            aTCGAGGCGCTGACCCA‑GCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGCCCAGTTGCAGGGCGACcagcag           >  4:224596/1‑138 (MQ=255)
                                                                                                                                            aTCGAGGCGCTGACCCA‑GCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGCCCAGTTGCAGGGCGACCAGCAGCGTCTcagc  >  2:103482/1‑147 (MQ=255)
                                                                                                                                                 |                                                                                                                                              
CCACCCCACCTGCCAAGCGCTCCGGCAGTGGCCTGGCAGCGCTGGCCTTGCTGCTTGGTGCGGCCGGGGTGGCGGTAGGTGGCTGGGGTGTCTGGCAGGTGCGTCAACTGCAGGGCAGTGAGTTCAACCAGGGCCAGCACATCGAGGCGCTGACCCA‑GCGCGCCGAGGCGCTGCAGCAACGTGAACAGCAGATCAGTGCGCAACTGGCCAGTCTGCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGCCCAGTTGCAGGGCGACCAGCAGCGTCTCAGC  >  NC_002947/239797‑240083

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