Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,038,225 G→C 8.6% G390G (GGC→GGG PP_0897 ← class 1 fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,038,2250GC 8.6% 206.9 / 7.4 68G390G (GGC→GGGPP_0897class 1 fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base G (30/32);  new base C (4/2);  total (34/34)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAACGGCCATCAGGTACACGGCCTTGTTGTCCTTGATCGCCTCGATGGCGGTCGGGCCACGCTCGGACTTGCCAATCATGCCCAGCAGGCCAGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTCTACCGGCAACTCTTCGCCACGGTTGAGCATTTCAACCATGCGCTTGTGCGCGGCATCACGGCCGGTCAGCATC  >  NC_002947/1038080‑1038372
                                                                                                                                                 |                                                                                                                                                     
cAACGGCCATCAGGTACACGGCCTTGTTGTCCTTGATCGCCTCGATGGCGGTCGGGCCACGCTCGGACTTGCCAATCATGCCCAGCAGGCCAGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCg                                                                                                                                                    >  4:10891/1‑149 (MQ=255)
 aaCGGCCATCAGGTACACGGCCTTGTTGTCCTTGATCGCCTCGATGGCGGTCGGGCCACGCTCGGACTTGCCAATCATGCCCAGCAGGCCAGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCgg                                                                                                                                                   >  3:147126/1‑149 (MQ=255)
         tCAGGTACACGGCCTTGTTGTCCTTGATCGCCTCGATGGCGGTCGGGCCACGCTCGGACTTGCCAATCATGCCCAGCAGGCCAGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTAGGCCCAGCCGGcccg                                                                                                                                           >  7:345298/1‑148 (MQ=255)
           aGGTACACGGCCTTGTTGTCCTTGATCGCCTCGATGGCGGTCGGGCCACGCTCGGACTTGCCAATCATGCCCAGCAGGCCAGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCC‑‑TAc                                                                                                                                       <  3:10891/149‑1 (MQ=255)
             gTACACGGCCTTGTTGTCCTTGATCGCCTCGATGGCGGTCGGGCCACGCTCGGACTTGCCAATCATGCCCAGCAGGCCAGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTAGGCCCAGCCGGCCC‑‑GACCa                                                                                                                                     <  8:345298/149‑1 (MQ=255)
                cacGGCCTTGTTGTCCTTGATCGCCTCGATGGCGGTCGGGCCACGCTCGGACTTGCCAATCATGCCCAGCAGGCCAGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCC‑‑TACCACtt                                                                                                                                  >  1:120666/1‑149 (MQ=255)
                cacGGCCTTGTTGTCCTTGATCGCCTCGATGGCGGTCGGGCCACGCTCGGACTTGCCAATCATGCCCAGCAGGCCAGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATAAGGGTAGCAGTTGTCGGGCCGGCTGGGCC‑‑TAACACt                                                                                                                                   >  6:91758/1‑148 (MQ=255)
                            gtCCTTGATCGCCTCGATGGCGGTCGGGCCACGCTCGGACTTGCCAATCATGCCCAGCAGGCCAGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCC‑‑TACCACTTCGTcacccacc                                                                                                                       <  3:185704/148‑1 (MQ=255)
                            gtCCTTGATCGCCTCGATGGCGGTCGGGCCACGCTCGGACTTGCCAATCATGCCCAGCAGGCCAGTCTGCTCAAGGATATGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCC‑‑TACCACTTCGTCACCCACCg                                                                                                                      >  2:283383/1‑149 (MQ=255)
                                      gCCTCGATGGCGGTCGGGCCACGCTCGGACTTGCCAATCATGCCCAGCAGGCCAGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCgg                                                                                                            >  6:205932/1‑149 (MQ=255)
                                      gCCTCGATGGCGGTCGGGCCACGCTCGGACTTGCCAATCATGCCCAGCAGGCCAGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCgg                                                                                                            >  2:368604/1‑149 (MQ=255)
                                       ccTCGATGGCGGTCGGGCCACGCTCGGACTTGCCAATCATGCCCAGCAGGCCAGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCggg                                                                                                           >  1:90764/1‑149 (MQ=255)
                                            aTGGCGGTCGGGCCACGCTCGGACTTGCCAATCATGCCCAGCAGGCCAGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCGGGCCGAc                                                                                                      <  7:300336/149‑1 (MQ=255)
                                            aTGGCGGTCGGGCCACGCTCGGACTTGCCAATCATGCCCAGCAGGCCAGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCGGGCCGAc                                                                                                      <  5:361445/149‑1 (MQ=255)
                                             tGGCGGTCGGGCCACTCTCGGACTTGCCAATCATGCCCAGCAGGCCAGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGTCTGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCGGGCCGAc                                                                                                      <  8:331768/148‑1 (MQ=255)
                                             tGGCGGTCGGGCCACGCTCGGACTTGCCAATCATGCCCAGCAGGCCAGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCGGGCCGACg                                                                                                     <  1:368604/149‑1 (MQ=255)
                                                cGGTCGGGCCACGCTCGGACTTGCCAATCATGCCCAGCAGGCCAGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCGGGCCGACgtag                                                                                                  <  4:173509/149‑1 (MQ=255)
                                                 ggTCGGGCCACGCTCGGACTTGCCAATCTTGCCCAGCAGGCCAGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCGGGCCGACgtag                                                                                                  >  1:207705/1‑148 (MQ=255)
                                                         cACGCTCGGACTTGCCAATCATGCCCAGCAGGCCAGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAg                                                                                         <  2:131443/149‑1 (MQ=255)
                                                         cACGCTCGGACTTGCCAATCATGCCCAGCAGGCCAGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAg                                                                                         >  4:99757/1‑149 (MQ=255)
                                                            gCTCGGACTTGCCAATCATGCCCAGCAGGCCAGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCgg                                                                                      >  1:148387/1‑149 (MQ=255)
                                                                 gACTTGCCAATCATGCCCAGCAGGCCAGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCttt                                                                                 >  8:151888/1‑149 (MQ=255)
                                                                 gACTTGCCAATCATGCCCAGCAGGCCAGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCttt                                                                                 >  6:273731/1‑149 (MQ=255)
                                                                    ttGCCAATCATGCCCAGCAGGCCAGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAg                                                                              >  6:133015/1‑149 (MQ=255)
                                                                       ccAATCATGCCCAGCAGGCCAGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTc                                                                           >  4:314657/1‑149 (MQ=255)
                                                                       ccAATCATGCCCAGCAGGCCAGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTc                                                                           >  2:314421/1‑149 (MQ=255)
                                                                         aaTCATGCCCAGCAGGCCAGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGtagtag                                                                                               <  1:29821/127‑1 (MQ=255)
                                                                         aaTCATGCCCAGCAGGCCAGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGtagtag                                                                                               >  2:29821/1‑127 (MQ=255)
                                                                          aTCATGCCCAGCAGGCCAGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTCTAc                                                                        <  1:283383/149‑1 (MQ=255)
                                                                             aTGCCCAGCAGGCCAGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTCTACCg                                                                      <  3:314657/148‑1 (MQ=255)
                                                                                 ccAGCAGGCCAGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTCTACCGGCAAc                                                                 >  8:231774/1‑149 (MQ=255)
                                                                                  cagcagGCCAGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTCTACCGGCAAc                                                                 >  5:8970/1‑148 (MQ=255)
                                                                                   agcagGCCAGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTCTACCGGCAActc                                                               <  7:151888/149‑1 (MQ=255)
                                                                                   agcagGCCAGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTCTACCGGCAActc                                                               <  5:273731/149‑1 (MQ=255)
                                                                                      agGCCAGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTCTACCGGCAACTCTTc                                                            <  2:120666/149‑1 (MQ=255)
                                                                                      agGCCAGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTCTACCGGCAACTCTTc                                                            <  3:99757/149‑1 (MQ=255)
                                                                                       ggCCAGTCTGCTCAAGGAGCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTCTACCGGCAACTCTTCg                                                           <  6:251291/149‑1 (MQ=255)
                                                                                        gCCAGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGc                                                                                        >  4:183064/1‑119 (MQ=255)
                                                                                        gCCAGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGc                                                                                        <  3:183064/119‑1 (MQ=255)
                                                                                             tCTGCTCAAGGATCTGGCGAGTGAACTTGTCGATACGGGTAGCAGTTTTCGGGCCGGCTGGGCC‑‑TACCACTTCGTCACCCACCTGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTCTACCGGCAACTCTTCGCCACgg                                                     <  2:90764/149‑1 (MQ=255)
                                                                                             tCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTCTACCGGCAACTCTTCGCCACgg                                                     <  3:78594/149‑1 (MQ=255)
                                                                                                 cTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTCTACCGGCAACTCTTc                                                            >  7:3711/1‑138 (MQ=255)
                                                                                                 cTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTCTACCGGCAACTCTTc                                                            <  8:3711/138‑1 (MQ=255)
                                                                                                 cTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTCTACCGGCAACTCTTCGCCACGGtt                                                   <  5:91758/147‑1 (MQ=255)
                                                                                                 cTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTCTACCGGCAACTCTTCGCCACGGTTGa                                                 <  4:356072/149‑1 (MQ=255)
                                                                                                 cTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTCTACCGGCAACTCTTCGCCACGGTTGa                                                 <  3:126511/149‑1 (MQ=255)
                                                                                                  tCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCCGATACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTCTACCGGCAACTCTTCGCCACGGtt                                                   >  1:290109/1‑148 (MQ=255)
                                                                                                  tCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTCTACCGGCAACTCTTCGCCACGGTTGAg                                                >  5:338952/1‑149 (MQ=255)
                                                                                                  tCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTCTACCGGCAACTCTTCGCCACGGTTGAg                                                >  1:20502/1‑149 (MQ=255)
                                                                                                  tCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTCTACCGGCAACTCTTCGCCACGGTTGAg                                                >  4:170094/1‑149 (MQ=255)
                                                                                                          cTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTCTACCGGCAACTCTTCGCCACGGTTGAGCATTTCa                                         <  1:314421/148‑1 (MQ=255)
                                                                                                             gCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCgg                                                                                      >  5:216364/1‑100 (MQ=255)
                                                                                                             gCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCgg                                                                                      <  6:216364/100‑1 (MQ=255)
                                                                                                             gCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTc                                                                                <  6:282949/106‑1 (MQ=255)
                                                                                                             gCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTCTACCGGCAACTCTTCGCCACGGTTGAGCATTTCaa                                        <  5:205932/146‑1 (MQ=255)
                                                                                                             gCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTCTACCGGCAACTCTTCGCCACGGTTGAGCATTTCAACCa                                     >  6:354366/1‑149 (MQ=255)
                                                                                                               gAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCttt                                                                                 <  2:114493/103‑1 (MQ=255)
                                                                                                               gAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCAGGCCGACGTAGTAGATGAAGCGGCCttt                                                                                 >  1:114493/1‑103 (MQ=255)
                                                                                                                  tGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTCTACCGGCAACTCTTCGCCACGGTTGAGCATTTCAACCATgcgc                                >  8:57842/1‑149 (MQ=255)
                                                                                                                     aCTTGTCCATACGGGTAGCAGTGGGCGGGCCGGCGGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTCTACCGGCAACTCTTCGCCACGGTTGAGCATTTCAACCATGCGCTtg                             <  5:302890/149‑1 (MQ=255)
                                                                                                                                cGGGTAGCAGTGGTCGGGCCGGCTGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTCTACCGGCAACTCTTCGCCACGGTTGAGCATTTCAACCATGCGCTTGTGCGCGGCATc                  <  6:338952/149‑1 (MQ=255)
                                                                                                                                    acgaaGTGGTAGGCCCAGCCGGCCC‑‑GACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTCTACCGGCAACTCTTCGCCACGGTTGAGCATTTCAACCATGCGCTTGTgcgc                        >  6:99245/5‑139 (MQ=255)
                                                                                                                                    acgaaGTGGTAGGCCCAGCCGGCCC‑‑GACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTCTACCGGCAACTCTTCGCCACGGTTGAGCATTTCAACCATGCGCTTGTgcgc                        <  5:99245/135‑1 (MQ=255)
                                                                                                                                    acgaaGTGGTAGGCCCAGCCGGCCC‑‑GACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTCTACCGGCAACTCTTCGCCACGGTTGAGCATTTCAACAATGCGATTGTGCTCGGCGTCACgg              >  6:256416/5‑149 (MQ=255)
                                                                                                                                    acgaaGTGGTAGGCCCAGCCGGCCC‑‑GACACCTTAGTCACCCACAGGGTCGACCGGGCCGACGTAGTAGATGTCGCGGGCTTTCAGGTCTACCGGGCAATCTTCGCAACGGTTGAGAATTTCACAAATACGAATGTGCGAGTCATCACg               >  6:256418/5‑148 (MQ=255)
                                                                                                                                       cAGTGGTCGGGCCGGCTGGGCC‑‑TACCCATTCGTAACACACCGTGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTCTACCGGCAACTCTTCGCCACGGTTGAGCATTTCAACCATGCGCTTGTGCGCGGCATCACGGCCg           >  7:37928/1‑149 (MQ=255)
                                                                                                                                       cAGTGGTCGGGCCGGCTGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTCTACCGGCAACTCTTCGCCACGGTTGAGCATTTCAACCATGCGCTTGTGCGCGGCATCACGGCCg           <  2:20502/149‑1 (MQ=255)
                                                                                                                                          tGGTCGGGCCGGCTGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTCTACCGGCAACTCTTCGCCACGGTTGAGCATTTCAACCATGCGCTTGTGCGCGGCATCACGGCCGGTc        <  4:147126/149‑1 (MQ=255)
                                                                                                                                                ggCCGGCTGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTCTACCGGCAACTCTTCGCCACGGTTGAGCATTTCAACCATGCGCTTGTGCGCGGCATCACGGCCGGTCAGCATc  <  3:170094/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
CAACGGCCATCAGGTACACGGCCTTGTTGTCCTTGATCGCCTCGATGGCGGTCGGGCCACGCTCGGACTTGCCAATCATGCCCAGCAGGCCAGTCTGCTCAAGGATCTGGCGAGTGAACTTGTCCATACGGGTAGCAGTGGTCGGGCCGGCTGGGCC‑‑TACCACTTCGTCACCCACCGGGTCGACCGGGCCGACGTAGTAGATGAAGCGGCCTTTCAGGTCTACCGGCAACTCTTCGCCACGGTTGAGCATTTCAACCATGCGCTTGTGCGCGGCATCACGGCCGGTCAGCATC  >  NC_002947/1038080‑1038372

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: