Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,039,350 G→A 8.0% A16A (GCC→GCT)  PP_0897 ← class 1 fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,039,3500GA 8.0% 173.5 / 3.8 62A16A (GCC→GCT) PP_0897class 1 fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base G (28/29);  new base A (2/3);  total (30/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

GACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTCTTGCTTGATCACGGTCATGCAGCGCGCTCCTCTTAAAGACGGGAACATTCAAAAAGACGCTCGACATAGCCGACCAGCACGTCGAAGCGCCTTGCAAGGCGCCGGCAGGT  >  NC_002947/1039208‑1039489
                                                                                                                                              |                                                                                                                                           
gACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCgg                                                                                                                                      >  2:1135406/1‑150 (MQ=255)
   gATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGa                                                                                                                                   >  2:555863/1‑150 (MQ=255)
         cGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGctct                                                                                                                             <  1:813479/150‑1 (MQ=255)
                gcACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGAGCGTCGCCGACGCCCTgcaattc                                                                                                                      <  2:1140728/150‑7 (MQ=255)
                gcACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGAGCGTCGCCGACGCCCTgcaattc                                                                                                                      >  2:251227/1‑144 (MQ=255)
                        atctgTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTATTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTctt                                                                                                              <  1:30549/146‑1 (MQ=255)
                          gAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTcttgc                                                                                                            >  2:1203022/1‑150 (MQ=255)
                                 tCAGGATCTGGGCGATGGAATCGCGTGCAGCAGTCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTCTTGCTTGATCa                                                                                                     >  1:515448/1‑150 (MQ=255)
                                    ggATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTCTTGCTTGATCACgg                                                                                                  >  2:1166643/1‑150 (MQ=255)
                                    ggATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTCTTGCTTGATCACgg                                                                                                  <  2:606296/150‑1 (MQ=255)
                                    ggATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTCTTGCTTGATCACgg                                                                                                  >  2:349922/1‑150 (MQ=255)
                                         tGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTCTTGCTTGATCACGGTCATg                                                                                             <  1:1135406/150‑1 (MQ=255)
                                         tGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTCTTGCTTGATCACGGTCATg                                                                                             <  1:555863/150‑1 (MQ=255)
                                                 ggAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGAGCGTCGCCGAc                                                                                                                                  >  2:592280/1‑105 (MQ=255)
                                                 ggAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGAGCGTCGCCGAc                                                                                                                                  <  1:592280/105‑1 (MQ=255)
                                                   aaTCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTCTTGCTTGATCACGGTCATGCAGCGCGctc                                                                                   <  2:411960/150‑1 (MQ=255)
                                                   aaTCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTCTTGCTTGATCACGGTCATGCAGCGCGctc                                                                                   >  2:442546/1‑150 (MQ=255)
                                                     tCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGa                                                                                                                           >  2:276017/1‑108 (MQ=255)
                                                     tCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGa                                                                                                                           <  1:276017/108‑1 (MQ=255)
                                                     tCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTCTTGCTTGATCACGGTCATGCAGCGCGctcct                                                                                 >  2:9884/1‑150 (MQ=255)
                                                        cgTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGc                                                                                                                                           >  2:953026/1‑89 (MQ=255)
                                                        cgTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGc                                                                                                                                           <  1:953026/89‑1 (MQ=255)
                                                              gcagGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCGTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTCTTGCTTGATCACGGTCATGCAGCGCGCTACTCTTAAAGAc                                                                        >  2:481881/1‑150 (MQ=255)
                                                               cagGGGAGTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAGTCGACGGGGTTGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTCTTGCTTGATCACGGTCATGCAGCGCGCTCCTCTTAAAGACg                                                                       <  2:522247/150‑1 (MQ=255)
                                                                agGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTCTTGCTTGATCACGGTCATGCAGCGCGCTCCTCTTAAAGACgg                                                                      >  2:847260/1‑150 (MQ=255)
                                                                agGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTCTTGCTTGATCACGGTCATGCAGCGCGCTCCTCTTAAAGACgg                                                                      <  1:841253/150‑1 (MQ=255)
                                                                    gATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTCTTGCTTGATCACGGTCATGCAGCGCGCTCCTCTTAAAGACGGGAAc                                                                  >  1:45504/1‑150 (MQ=255)
                                                                     aTTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTCTTGCTTGATCACGGTCATGCAGCGCGCTCCTCTTAAAGACGGGAACa                                                                 >  1:723206/1‑150 (MQ=255)
                                                                       tCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTCTTGCTTGATCACGGTCATGCAGCGCGCTCCTCTTAAAGACGGGAACAtt                                                               <  2:723206/150‑1 (MQ=255)
                                                                        ccTCACGCAGATAGGCCTCGTGCATGGCCTGGCTGACACAGTCGGACTACTATTACGATATTAATTGCAGGGTGTAGGCGACGCTCTGAATCAGGTCGTCTTGCTTGATCACGGTCATGCAGCGCGCTCCTCTTAAAGACGGGAACATTc                                                              <  1:570584/150‑1 (MQ=255)
                                                                          tCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTcttgct                                                                                                           >  1:607030/1‑103 (MQ=255)
                                                                          tCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTcttgct                                                                                                           <  2:607030/103‑1 (MQ=255)
                                                                          tCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTCTTGCTTGATCACGGTCATGCAGCGCGCTCCTCTTAAAGACGGGAACATTCaa                                                            >  1:89088/1‑150 (MQ=255)
                                                                             cGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTCTTGCTTGATCACGGTCATGCAGCGCGCTCCTCTTAAAGACGGGAACATTCaaaaa                                                         <  1:1166643/150‑1 (MQ=255)
                                                                              gCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAgg                                                                                                                                             >  2:53680/1‑65 (MQ=255)
                                                                              gCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAgg                                                                                                                                             <  1:53680/65‑1 (MQ=255)
                                                                               cAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTCTTGCTTGATCACGGTCATGCAGCGCGCTCCTCTTAAAGACGGGAACATTCAAAAAGa                                                       <  2:515448/150‑1 (MQ=255)
                                                                                 gATAGGCCGCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTCTTGCTTGATCACGGTCATGCAGCGCGCTCCTCTTAAAGACGGGAACATTCAAAAAGACg                                                     <  1:596195/150‑1 (MQ=255)
                                                                                      gCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTCTTGCTTGATCACGGTCATGCAGCGCGCTCCTCtt                                                                              <  2:299437/120‑1 (MQ=255)
                                                                                      gCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTCTTGCTTGATCACGGTCATGCAGCGCGCTCCTCtt                                                                              >  1:299437/1‑120 (MQ=255)
                                                                                       ccTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTCTTGCTTGATCACGGTCATGCAGCGCGCTCCTCTTAAAGACGGGAACATTCAAAAAGACGCTCGAc                                               >  1:71820/1‑150 (MQ=255)
                                                                                           gTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTCTTGCTTGATCACGGTCATGCAGCGCGCTCCTCTTAAAGACGGGAACATTCAAAAAGACGCTCGACATAg                                           <  1:442546/150‑1 (MQ=255)
                                                                                           gTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTCTTGCTTGATCACGGTCATGCAGCGCGCTCCTCTTAAAGACGGGAACATTCAAAAAGACGCTCGACATAg                                           <  1:349922/150‑1 (MQ=255)
                                                                                                 ggCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTCTTGCTTGATCACGGTCATGCAGCGCGCTCCTCTTAAAGACGGGAACATTCAAAAAGACGCTCGACATAGCCGAcc                                     <  1:9884/150‑1 (MQ=255)
                                                                                                   ccTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTCTTGCTTGATCACGGTCATGCAGCGCGCTCCTCTTAAAGACGGGAACATTCAAAAAGACGCTCGACATAGCCGACCAg                                   >  1:658116/1‑150 (MQ=255)
                                                                                                    cTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTCTTGCTTGATCACGGTCATGCAGCGCGCTACTCTTAAAGACGGGAACATTCAAAAAGACGCTCGACATAGCCGACCAGc                                  <  1:481881/150‑1 (MQ=255)
                                                                                                    cTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGAGCGTCGCCGACGccct                                                                                                                             <  1:1182121/59‑1 (MQ=37)
                                                                                                        aTGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTCTTGCTTGATCACGGTCATGCAGCGCGCTCCTCTTAAAGACGGGAACATTCAAAAAGACGCTCGACATAGCCGACCAGCACGt                              >  1:1006644/1‑150 (MQ=255)
                                                                                                            aaTCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTCTTGCTTGATCACGGTCATGCAGCGCGCTCCTCTTAAAGACGGGAACATTCAAAAAGACGCTCGACATAGCCGACCAGCACGTCGaa                          <  1:1203022/150‑1 (MQ=255)
                                                                                                             aTCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTCTTGCTTGATCACGGTCATGCAGCGCGCTCCTCTTAAAGACGGGAACATTCAAAAAGACGCTCGACATAGTCGACTAGCGCGTCGAAg                         <  2:45504/150‑1 (MQ=255)
                                                                                                             aTCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTCTTGCTTGATCACGGTCATGCAGCGCGCTCCTCTTAAAGACGGGAACATTCAAAAAGACGCTCGACATAGCCGACCAGCACGTCGAAg                         >  2:747494/1‑150 (MQ=255)
                                                                                                             aTCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTCTTGCTTGATCACGGTCATGCAGCGCGCTCCTCTTAAAGACGGGAACATTCAAAAAGACGCTCGACATAGCCGACCAGCACGTCGAAg                         >  2:98882/1‑150 (MQ=255)
                                                                                                                gACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTCTTGCTTGATCACGGTCATGCAGCGCGCTCCTCTTAAAGACGGGAACATTCAAAAAGACg                                                     <  2:1017539/119‑1 (MQ=255)
                                                                                                                gACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTCTTGCTTGATCACGGTCATGCAGCGCGCTCCTCTTAAAGACGGGAACATTCAAAAAGACg                                                     >  1:1017539/1‑119 (MQ=255)
                                                                                                                        ggtAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTCTTGCTTGATCACGGTCATGCAGCGCGCTCCTCTTAAAGACGGGAACATTCAAAAAGACGCTCGACATAGCCGACCAGCACGTCGAAGCGCCTTGCAAg              <  2:658116/150‑1 (MQ=255)
                                                                                                                         gtagtaCGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTCTTGCTTGATCACGGTCATGCAGCGCGCTCCTCTTAAAGACGGGAACATTCAAAAAGACGCTCGACATAGCCGACCAGCACGTCGAAGCGCCTTGCAAgg             <  1:847260/150‑1 (MQ=255)
                                                                                                                             taCGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTCTTGCTTGATCACGGTCATGCAGCGCGCTCCTCTTAAAGACGGGAACATTCAAAAAGACGCTCGACATAGCCGACCAGCACGTCGAAGCGCCTTGCAAGGCGcc         >  1:57735/1‑150 (MQ=255)
                                                                                                                             taCGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTCTTGCTTGATCACGGTCATGCAGCGCGCTCCTCTTAAAGACGGGAACATTCAAAAAGACGCTCGACATAGCCGACCAGCACGTCGAAGCGCCTTGCAAGGCGcc         >  1:204976/1‑150 (MQ=255)
                                                                                                                                gagaTGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTCTTGCTTGATCACGGTCATGCAgcg                                                                                        <  2:741157/68‑1 (MQ=255)
                                                                                                                                gagaTGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTCTTGCTTGATCACGGTCATGCAgcg                                                                                        >  1:741157/1‑68 (MQ=255)
                                                                                                                                 agaTGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTCTTGCTTGATCACGGTCATGCAGCGCGCTCCTCTTAAAGACGGGAACATTCAAAAAGACGCTCGACATAGCCGACCAGCACGTCGAAgcgc                      >  2:611113/1‑133 (MQ=255)
                                                                                                                                 agaTGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTCTTGCTTGATCACGGTCATGCAGCGCGCTCCTCTTAAAGACGGGAACATTCAAAAAGACGCTCGACATAGCCGACCAGCACGTCGAAgcgc                      <  1:611113/133‑1 (MQ=255)
                                                                                                                                 agaTGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTCTTGCTTGATCACGGTCATGCAGCGCGCTCCTCTTAAAGACGGGAACATTCAAAAAGACGCTCGACATAGCCGACCAGCACGTCGAAGCGCCTTGCAAGGCGCCGGCa     <  2:204976/150‑1 (MQ=255)
                                                                                                                                    tGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTCTTGCTTGATCACGGTCATGCAGCGCGCTCCTCTTAAAGACGGGAACATTCAAAAAGACGCTCGACATAGCCGACCAGCACGTCGAAGCGCCTTGCAAGGCGCCGGCAGGt  >  2:882375/1‑150 (MQ=255)
                                                                                                                                              |                                                                                                                                           
GACGATGGCCGGTGGCGCACATGCGCGAGTTGATCAGGATCTGGGCGATGGAATCGCGTGCAGCAGGCGATTCCTCACGCAGATAGGCCTCGTGCATGGCCTGGATGAAATCGACGGGGTGGTAGTACGAGATGAATTGCAGGGCGTCGGCGACGCTCTGAATCAGGTCGTCTTGCTTGATCACGGTCATGCAGCGCGCTCCTCTTAAAGACGGGAACATTCAAAAAGACGCTCGACATAGCCGACCAGCACGTCGAAGCGCCTTGCAAGGCGCCGGCAGGT  >  NC_002947/1039208‑1039489

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑