Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,037,957 C→G 5.9% D480H (GAC→CAC)  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,9570CG 5.9% 219.2 / 4.0 68D480H (GAC→CAC) PP_0897class 1 fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base C (31/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 = 8.86e-01

CGCTGTACGTGTGGGAGCGGGCTTGCCCGCGAAGCAGCCGCCGCAGTTACATCAGGTTACTTCACTTCGACGGCCAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAGGCCAGTACCTTCGACTTGCGGATGGCCTGGGCGACCAGGTAGGCAGCGCCGCCAACGGCCATCAGGTACACGGCC  >  NC_002947/1037815‑1038102
                                                                                                                                              |                                                                                                                                                 
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:202204/1‑149 (MQ=255)
                          ccGCGAAGCAGCCGCCGCAGTTACATCAGGTTACTTCACTTCGACGGCCAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACtt                                                                                                                   >  7:196897/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                                                                                       >  5:294903/1‑149 (MQ=255)
                                                      ggTTACTTCACTTCGACGGCCAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAg                                                                                       >  7:104271/1‑149 (MQ=255)
                                                       ggtACTTCACTTCGACGGCCAGGCTTTCGGCAATCTTGCACTGCCACAGGGCAGGGCCAGTGATGTGAACCGATACACCGTTGCTGTACACGGCCACGGTCACCGGCATTTCCTTGACTTCCAACTCGTAGAACGCTTCCATGCCCAg                                                                                       >  3:234050/3‑148 (MQ=255)
                                                            ttcacttcGACGGCCAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCAt                                                                                             >  3:266779/1‑137 (MQ=255)
                                                            ttcacttcGACGGCCAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCAt                                                                                             <  4:266779/137‑1 (MQ=255)
                                                                  tcGACGGCCAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAGGc                                                                           >  2:138084/1‑149 (MQ=255)
                                                                   cGACGGCCAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAGGc                                                                           <  1:138084/148‑1 (MQ=255)
                                                                      cGGCCAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAGGCCAGt                                                                       <  1:138440/149‑1 (MQ=255)
                                                                         ccAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGccc                                                                                         >  5:218037/1‑128 (MQ=255)
                                                                         ccAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGccc                                                                                         <  6:218037/128‑1 (MQ=255)
                                                                         ccAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAGGCCAGTAcc                                                                    >  2:181638/1‑149 (MQ=255)
                                                                          cAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTCGCCGTGGACAGCAACGGt                                                                                                                                       >  8:162569/1‑81 (MQ=37)
                                                                          cAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTCGCCGTGGACAGCAACGGt                                                                                                                                       <  7:162569/81‑1 (MQ=37)
                                                                            ggCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAGGCCAGTACCTTc                                                                 <  3:61383/149‑1 (MQ=255)
                                                                            ggCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAGGCCAGTACCTTc                                                                 <  2:50578/149‑1 (MQ=255)
                                                                              cggtCGGCAATCTTGCTCTGCCCCAGGGCAGGGCCAGTGATGTGTCCCGATTCACCCTTGCTGTCCACCGCGACGGTCACCGGCATCTCCTTGCCATCGAACTCCTAGATCGCTTCCATGCCCAGCTCGGCGAAGGCCAGTACCTGCGa                                                               <  8:104271/146‑1 (MQ=255)
                                                                                   ggCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAGGCCAGTACCTTCGACTTGc                                                          >  1:286961/1‑149 (MQ=255)
                                                                                   ggCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAGGCCAGTACCTTCGACTTGc                                                          <  8:242654/149‑1 (MQ=255)
                                                                                   ggCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAGGCCAGTACCTTCGACTTGc                                                          >  7:242654/1‑149 (MQ=255)
                                                                                       aTCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAGGCCAGTACCTTCGACTTGCGGAt                                                      <  6:294903/149‑1 (MQ=255)
                                                                                       aTCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAGGCCAGTACCTTCGACTTGCGGAt                                                      <  2:286961/149‑1 (MQ=255)
                                                                                        tCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAGGc                                                                           <  2:318420/127‑1 (MQ=255)
                                                                                        tCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAGGc                                                                           >  1:318420/1‑127 (MQ=255)
                                                                                         cTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAGGCCAGTACCTTCGACt                                                             >  8:303736/1‑140 (MQ=255)
                                                                                         cTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAGGCCAGTACCTTCGACt                                                             <  7:303736/140‑1 (MQ=255)
                                                                                            gCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAGGCCAGTACCTTCGACTTGCGGATGGc                                                   <  3:297778/147‑1 (MQ=255)
                                                                                                tGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAGGCCAGTACCTTCGACTTGCGGATGGCCTGGGc                                             >  4:137528/1‑149 (MQ=255)
                                                                                                  ccACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAGGCCAGTACCTTCGACTTGCGGATGGCCTGGGCGa                                           >  3:112442/1‑149 (MQ=255)
                                                                                                  ccACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAGGCCAGTACCTTCGACTTGCGGATGGCCTGGGCGa                                           <  5:293905/149‑1 (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       >  7:38374/1‑142 (MQ=255)
                                                                                                                                             tCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAGGCCAGTACCTTCGACTTGCGGATGGCCTGGGCGACCAGGTAGGCAGCGCCGCCAACGGCCATCAGGTaca       <  8:38374/142‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                 
CGCTGTACGTGTGGGAGCGGGCTTGCCCGCGAAGCAGCCGCCGCAGTTACATCAGGTTACTTCACTTCGACGGCCAGGCTTTCGGCAATCTTGCTCTGCCACAGGGCAGGGCCAGTGATGTGTACCGATTCACCGTTGCTGTCCACGGCGACGGTCACCGGCATGTCCTTGACTTCGAACTCGTAGATCGCTTCCATGCCCAGCTCGGCGAAGGCCAGTACCTTCGACTTGCGGATGGCCTGGGCGACCAGGTAGGCAGCGCCGCCAACGGCCATCAGGTACACGGCC  >  NC_002947/1037815‑1038102

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