Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,037,958 C→G 5.9% V479V (GTG→GTC)  PP_0897 ← class 1 fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,037,9580CG 5.9% 205.3 / 3.9 68V479V (GTG→GTC) PP_0897class 1 fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base C (32/32);  new base G (2/2);  total (34/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.71e-01

CGCTGTACGTGTGGGAGCGGGCTTGCCCGCGAAGCAGCCGCCGCAGTTACATCAGGTTACTTCACTTCGACGGCCAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAGGCCAGTACCTTCGACTTGCGGATGGCCTGGGCGACCAGGTAGGCAGCGCCGCCAACGGCCATCAGGTACACGGCCTTGT  >  NC_002947/1037815‑1038106
                                                                                                                                               |                                                                                                                                                    
cgcTGTACGTGTGGGAGCGGGCTTGCCCGCGAAGCAGCCGCCGCAGTTACATCAGGTTACTTCACTTCGACGGCCAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACgg                                                                                                                                                  <  7:100822/148‑1 (MQ=255)
            gggAGCGGGCTTGCCCGCGAAGCAGCCGCCGCAGTTACATCAGGTTACTTCACTTCGACGGCCAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTAGCGATTCACCGTCGCCGTGGACAGCAACGGt                                                                                                                                           <  1:30054/143‑1 (MQ=255)
            gggAGCGGGCTTGCCCGCGAAGCAGCCGCCGCAGTTACATCAGGTTACTTCACTTCGACGGCCAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTCGCCGTGGACAGCAACGGt                                                                                                                                           >  2:30054/1‑143 (MQ=255)
            cttcgCGGGCAAGCCCGCGAAGCAGCCGCCGCAGTTACATCAGGTTACTTCACTTCGACGGCCAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCgg                                                                                                                                     >  1:50578/5‑149 (MQ=255)
             ggAGCGGGCTTGCCCGCGAAGCAGCCGCCGCAGTTACATCAGGTTACTTCACTTCGACGGCCAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGc                                                                                                                                                 >  8:211454/1‑136 (MQ=255)
             ggAGCGGGCTTGCCCGCGAAGCAGCCGCCGCAGTTACATCAGGTTACTTCACTTCGACGGCCAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGc                                                                                                                                                 <  7:211454/136‑1 (MQ=255)
             ggAGCGGGCTTGCCCGCGAAGCAGCCGCCGCAGTTACATCAGGTTACTTCACTTCGACGGCCAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCgg                                                                                                                                     >  6:331545/1‑148 (MQ=255)
               aGCGGGCTTGCCCGCGAAGCAGCCGCCGCAGTTACATCAGGTTACTTCACTTCGACGGCCAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCAt                                                                                                                                  <  6:117603/149‑1 (MQ=255)
                gCGGGCTTGCCCGCGAAGCAGCCGCCGCAGTTACATCAGGTTACTTCACTTCGACGGCCAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACg                                                                                                                                             <  3:34942/137‑1 (MQ=255)
                gCGGGCTTGCCCGCGAAGCAGCCGCCGCAGTTACATCAGGTTACTTCACTTCGACGGCCAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACg                                                                                                                                             >  4:34942/1‑137 (MQ=255)
                gCGGGCTTGCCCGCGAAGCAGCCGCCGCAGTTACATCAGGTTACTTCACTTCGACGGCCAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATg                                                                                                                                 >  4:293734/1‑149 (MQ=255)
                   ggCTTGCCCGCGAAGCAGCCGCCGCAGTTACATCAGGTTACTTCACTTCGACGGCCAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTcc                                                                                                                              >  4:61383/1‑149 (MQ=255)
                      ttGCCCGCGAAGCAGCCGCCGCAGTTACATCAGGTTACTTCACTTCGACGGCCAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACgg                                                                                                                                                  <  4:252224/126‑1 (MQ=255)
                      ttGCCCGCGAAGCAGCCGCCGCAGTTACATCAGGTTACTTCACTTCGACGGCCAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACgg                                                                                                                                                  >  3:252224/1‑126 (MQ=255)
                        gCCCGCGAAGCAGCCGCCGCAGTTACATCAGGTTACTTCACTTCGACGGCCAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCTGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGAc                                                                                                                         >  7:197753/1‑149 (MQ=255)
                          ccGCGAAGCAGCCGCCGCAGTTACATCAGGTTACTTCACTTCGACGGCCAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACtt                                                                                                                       >  7:196897/1‑149 (MQ=255)
                          ccGCGAAGCAGCCGCCGCAGTTACATCAGGTTACTTCACTTCGACGGCCAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACtt                                                                                                                       >  7:202204/1‑149 (MQ=255)
                           cgcgAAGCAGCCGCCGCAGTTACATCAGGTTACTTCACTTCGACGGCCAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTc                                                                                                                      <  3:148120/149‑1 (MQ=255)
                             cgAAGCAGCCGCCGCAGTTACATCAGGTTACTTCACTTCGACGGCCAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGa                                                                                                                    <  3:293734/149‑1 (MQ=255)
                                 gcagcCGCCGCAGTTACATCAGGTTACTTCACTTCGACGGCCAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCTGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTc                                                                                                                <  8:197753/149‑1 (MQ=255)
                                 gcagcCGCCGCAGTTACATCAGGTTACTTCACTTCGACGGCCAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTGACCGTTGCTGTCCACGGCGTCGGTCACCGGCATGTCCTTGACTTCGAACTc                                                                                                                <  5:331545/149‑1 (MQ=255)
                                             ggtccATCAGGTTACTGCACTTCGACGGCCAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATGCACCGGTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTc                                                                                                    <  8:202204/145‑1 (MQ=255)
                                             gTTACATCAGGTTACTTCACTTCGACGGCCAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTc                                                                                                    <  8:196897/149‑1 (MQ=255)
                                                     aGGTTACTTCACTTCGACGGCCAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGcc                                                                                              >  4:297778/1‑147 (MQ=255)
                                                     aGGTTACTTCACTTCGACGGCCAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGcc                                                                                              >  2:86971/1‑147 (MQ=255)
                                                      ggTTACTTCACTTCGACGGCCAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAg                                                                                           >  7:104271/1‑149 (MQ=255)
                                                      ggTTACTTCACTTCGACGGCCAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAg                                                                                           >  5:294903/1‑149 (MQ=255)
                                                       ggtACTTCACTTCGACGGCCAGGCTTTCGGCAATCTTGCACTGCCACAGGGCAGGGCCAGTGATGTGAACCGATACACCGTTGCTGTACACGGCCACGGTCACCGGCATTTCCTTGACTTCCAACTCGTAGAACGCTTCCATGCCCAg                                                                                           >  3:234050/3‑148 (MQ=255)
                                                            ttcacttcGACGGCCAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCAt                                                                                                 <  4:266779/137‑1 (MQ=255)
                                                            ttcacttcGACGGCCAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCAt                                                                                                 >  3:266779/1‑137 (MQ=255)
                                                                  tcGACGGCCAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAGGc                                                                               >  2:138084/1‑149 (MQ=255)
                                                                   cGACGGCCAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAGGc                                                                               <  1:138084/148‑1 (MQ=255)
                                                                      cGGCCAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAGGCCAGt                                                                           <  1:138440/149‑1 (MQ=255)
                                                                         ccAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGccc                                                                                             <  6:218037/128‑1 (MQ=255)
                                                                         ccAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGccc                                                                                             >  5:218037/1‑128 (MQ=255)
                                                                         ccAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAGGCCAGTAcc                                                                        >  2:181638/1‑149 (MQ=255)
                                                                          cAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTCGCCGTGGACAGCAACGGt                                                                                                                                           <  7:162569/81‑1 (MQ=37)
                                                                          cAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTCGCCGTGGACAGCAACGGt                                                                                                                                           >  8:162569/1‑81 (MQ=37)
                                                                            ggCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAGGCCAGTACCTTc                                                                     <  3:61383/149‑1 (MQ=255)
                                                                            ggCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAGGCCAGTACCTTc                                                                     <  2:50578/149‑1 (MQ=255)
                                                                              cggtCGGCAATCTTGCTCTGCCCCAGGGCAGGGCCAGTGATGTGTCCCGATTCACCCTTGCTGTCCACCGCGACGGTCACCGGCATCTCCTTGCCATCGAACTCCTAGATCGCTTCCATGCCCAGCTCGGCGAAGGCCAGTACCTGCGa                                                                   <  8:104271/146‑1 (MQ=255)
                                                                                   ggCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAGGCCAGTACCTTCGACTTGc                                                              >  7:242654/1‑149 (MQ=255)
                                                                                   ggCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAGGCCAGTACCTTCGACTTGc                                                              >  1:286961/1‑149 (MQ=255)
                                                                                   ggCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAGGCCAGTACCTTCGACTTGc                                                              <  8:242654/149‑1 (MQ=255)
                                                                                       aTCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAGGCCAGTACCTTCGACTTGCGGAt                                                          <  2:286961/149‑1 (MQ=255)
                                                                                       aTCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAGGCCAGTACCTTCGACTTGCGGAt                                                          <  6:294903/149‑1 (MQ=255)
                                                                                        tCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAGGc                                                                               >  1:318420/1‑127 (MQ=255)
                                                                                        tCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAGGc                                                                               <  2:318420/127‑1 (MQ=255)
                                                                                         cTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAGGCCAGTACCTTCGACt                                                                 <  7:303736/140‑1 (MQ=255)
                                                                                         cTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAGGCCAGTACCTTCGACt                                                                 >  8:303736/1‑140 (MQ=255)
                                                                                            gCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAGGCCAGTACCTTCGACTTGCGGATGGc                                                       <  3:297778/147‑1 (MQ=255)
                                                                                                tGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAGGCCAGTACCTTCGACTTGCGGATGGCCTGGGc                                                 >  4:137528/1‑149 (MQ=255)
                                                                                                  ccACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAGGCCAGTACCTTCGACTTGCGGATGGCCTGGGCGa                                               <  5:293905/149‑1 (MQ=255)
                                                                                                  ccACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAGGCCAGTACCTTCGACTTGCGGATGGCCTGGGCGa                                               >  3:112442/1‑149 (MQ=255)
                                                                                                   cacaGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAGGCCAGTACCTTCGACTTGCGGATGGCCTGGGCGAc                                              <  3:137528/149‑1 (MQ=255)
                                                                                                      agggcagggcCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAGGCCAGTACCTTCGACTTGCGGATGGCCTGGGCGACCAg                                           >  3:317627/1‑149 (MQ=255)
                                                                                                          cagggcCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAGGCCAGTACCTTCGACTTGCGGATGGCCTGGGCGACCAGGTAg                                       <  5:295103/149‑1 (MQ=255)
                                                                                                                cAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAGGCCAGTACCTTCGa                                                                   >  6:239762/1‑115 (MQ=255)
                                                                                                                cAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAGGCCAGTACCTTCGa                                                                   <  5:239762/115‑1 (MQ=255)
                                                                                                                         gtACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGATCTCGTAGATCGCTTCCATGACCAGCTCGGCGAAGCCCAGTACCTTCGACTTGCGGATGGCCTGGGCGACCAGGTAGGCAGCGCCGACAACg                        >  5:42090/1‑149 (MQ=255)
                                                                                                                          tACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAGGCCAGTACCTTCGACTTGCGGATGGCCTGGGCGACCAGGTAGGCAGCGCCGCCAACgg                       >  2:202661/1‑149 (MQ=255)
                                                                                                                          tACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAGGCCAGTACCTTCGACTTGCGGATGGCCTGGGCGACCAGGTAGGCAGCGCCGCCAACg                        <  3:158879/148‑1 (MQ=255)
                                                                                                                                ttCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAGGCCAGTACCTTCGACTTGCGGATGGCCTGGGCGACCAGGTAGGCAGCGCCGCCAACGGCCAt                   >  6:313875/1‑147 (MQ=255)
                                                                                                                                   aCCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCTGCGAAGGCCAGTACCTTCGACTTGCGGATGGCCTGGGCGACCAg                                           <  2:230277/120‑1 (MQ=255)
                                                                                                                                   aCCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAGGCCAGTACCTTCGACTTGCGGATGGCCTGGGCGACCAg                                           >  1:230277/1‑120 (MQ=255)
                                                                                                                                             tCCACGGCGACGGTCACCGGCATGTCCTTGACTTCTAACTCGTACATCGCTTCCATGCCCAGCTCGGCGAAGGCCAGTACCTTCGACTTGCGTATGGCCTGGGCGACCAGGTAGGCAGCGCCGCCAACGGCCATCAGGTACACGGcc      <  6:42090/147‑1 (MQ=255)
                                                                                                                                             tCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAGGCCAGTACCTTCGACTTGCGGATGGCCTGGGCGACCAGGTAGGCAGCGCCGCCAACGGCCATCAGGTaca           <  8:38374/142‑1 (MQ=255)
                                                                                                                                             tCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAGGCCAGTACCTTCGACTTGCGGATGGCCTGGGCGACCAGGTAGGCAGCGCCGCCAACGGCCATCAGGTaca           >  7:38374/1‑142 (MQ=255)
                                                                                                                                               cACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAGGCCAGTACCTTCGACTTGCGGATGGCCTGGGCGACCAGGTAGGCAGCGCCGCCAACGGCCATCAGGTACACGGCCttgt  >  4:338347/1‑149 (MQ=255)
                                                                                                                                               cACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAGGCCAGTACCTTCGACTTGCGGATGGCCTGGGCGACCAGGTAGGCAGCGCCGCCAACGGCCATCAGGTACACGGCCttgt  >  4:137844/1‑149 (MQ=255)
                                                                                                                                               cACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAGGCCAGTACCTTCGACTTGCGGATGGCCTGGGCGACCAGGTAGGCAGCGCCGCCAACGGCCATCAGGTACACGGCCttgt  >  1:81403/1‑149 (MQ=255)
                                                                                                                                               cACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAGGCCAGTACCTTCGACTTGCGGATGGCCTGGGCGACCAGGTAGGCAGCGCCGCCAACGGCCATCAGGTACACGGCCttgt  >  1:233641/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
CGCTGTACGTGTGGGAGCGGGCTTGCCCGCGAAGCAGCCGCCGCAGTTACATCAGGTTACTTCACTTCGACGGCCAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAGGCCAGTACCTTCGACTTGCGGATGGCCTGGGCGACCAGGTAGGCAGCGCCGCCAACGGCCATCAGGTACACGGCCTTGT  >  NC_002947/1037815‑1038106

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