Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,993,162 C→G 11.0% V357V (GTC→GTG bkdAA → branched‑chain alpha‑keto acid dehydrogenase complex subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,993,1620CG11.0% 173.4 / 10.8 62V357V (GTC→GTGbkdAAbranched‑chain alpha‑keto acid dehydrogenase complex subunit alpha
Reads supporting (aligned to +/- strand):  ref base C (26/29);  new base G (4/3);  total (30/32)
Fisher's exact test for biased strand distribution p-value = 7.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

TGCCGGCCCGCACTCGACGTCGGACGACCCCTCCAAGTACCGCCCTGCCGACGACTGGAGCCACTTCCCGCTGGGTGACCCGATCGCCCGCCTGAAGCAGCACCTGATCAAGATCGGCCACTGGTCCGAGGAAGAACACCAGGCCGTCACGGCCGAGCTCGAAGCTGCGGTGATTGCCGCACAGAAAGAAGCCGAGCAGTACGGCACCCTGGCCAACGGGCACATCCCGAGCGCCGCCTCGATGTTCGAGGATGTGTACAAGGAAATGCCCGACCACCTGCGCCGTCAACGCCAG  >  NC_002947/4993015‑4993309
                                                                                                                                                   |                                                                                                                                                   
tGCCGGCCCGCACTCGACGTCGGACGACCCCTCCAAGTACCGCCCTGCCGACGACTGGAGCCACTTCCCGCTGGGTGACCCGATCGCCCGCCTGAAGCAGCACCTGATCAAGATCGGCCACTGGTCCGAGGAAGAACACCAGGCCGTCAc                                                                                                                                                   <  2:607446/150‑1 (MQ=255)
   cGGCCCGCACTCGACGTCGGACGACCCCTCCAAGTACCGCCCTGCCGACGACTGGAGCCACTTCCCGCTGGGTGACCCGATCGCCCGCCTGAAGCAGCACCTGATCAAGATCGGCCACTGGTCCGAGGAAGAACACCAGGCCGTCACGGc                                                                                                                                                >  1:574692/1‑150 (MQ=255)
             tCGACGTCGGACGACCCCTCCAAGTACCGCCCTGCCGACGACTGGAGCCACTTCCCGCTGGGTGACCCGATCGCCCGCCTGAAGCAGCACCTGATCAAGATCGGCCACTGGTCCGAGGAAGAACACCAGGCCGTGACGGcctg                                                                                                                                             <  2:520482/143‑3 (MQ=255)
             tCGACGTCGGACGACCCCTCCAAGTACCGCCCTGCCGACGACTGGAGCCACTTCCCGCTGGGTGACCCGATCGCCCGCCTGAAGCAGCACCTGATCAAGATCGGCCACTGGTCCGAGGAAGAACACCAGGCCGTGACGGcctg                                                                                                                                             >  1:520482/1‑141 (MQ=255)
               gACGTCGGACGACCCCTCCAAGTACCGCCCTGCCGACGACTGGAGCCACTTCCCGCTGGGTGACCCGATCGCCCGCCTGAAGCAGCACCTGATCAAGATCGGCCACTGGTCCGAGGAAGAACACCAGGCCGTCACGGCCGAGCTCGAAGc                                                                                                                                    >  1:1382967/1‑150 (MQ=255)
                  gTCGGACGACCCCTCCAAGTACCGCCCTGCCGACGACTGGAGCCACTTCCCGTTGGGTGACCCGATCGCCCGCCTGAAGCAGCACCTGATCAAGATCGGCCACTGGTCCGAGGAAGAACACCAGGCCGTCACGGCCGAGCTCGAAGCTGc                                                                                                                                 >  1:1415731/1‑150 (MQ=255)
                      gacgacCCCTCCAAGTACCGCCCTGCCGACGACTGGAGCCACTTCCCGCTGGGTGACCCGATCGCCCGCCTGAAGCAGCACCTGATCAAGATCGGCCACTGGTCCGAGGAAGAACACCAGGCCGTCACGGCCGAGCTCGAAGCTGCGGTg                                                                                                                             >  1:128141/1‑150 (MQ=255)
                       acgacCCCTCCAAGTACCGCCCTGCCGACGACTGGAGCCACTTCCCGTTGGGTGACCCGATCGCCCGCCTGAAGCAGCACCTGATCAAGATCGGCCACTGGTCCGAGGAAGAACACCAGGCCGTCACGGCCGAGCTCGAAGCTGCGGTGa                                                                                                                            <  2:1415731/150‑1 (MQ=255)
                              cTCCAAGTACCGCCCTGCCGACGACTGGAGCCACTTCCCGCTGGGTGACCCGATCGCCCGCCTGAAGCAGCACCTGATCAAGATCGGCCACTGGTCAGAGGAAGAACACCAGGCCGTCACGGCCGAGCTCGAAGCTGCGGTGATTGCCGc                                                                                                                     >  1:137607/1‑150 (MQ=255)
                               tCCAAGTACCGCCCTGCCGACGACTGGAGCCACTTCCCGCTGGGTGACCCGATCGCCCGCCTGAAGCAGCACCTGATCAAGATCGGCCACTGGTCCGAGGAAGAACACCAGGCCGTCACGGCCGAGCTCGAAGCTGCGGTGATTGCCg                                                                                                                      <  2:435139/148‑1 (MQ=255)
                               tCCAAGTACCGCCCTGCCGACGACTGGAGCCACTTCCCGCTGGGTGACCCGATCGCCCGCCTGAAGCAGCACCTGATCAAGATCGGCCACTGGTCCGAGGAAGAACACCAGACCGTCACGGCCGAGCTCGAAGCTGCGGTGATTGCCg                                                                                                                      >  1:435139/1‑148 (MQ=255)
                                 cAAGTACCGCCCTGCCGACGACTGGAGCCACTTCCCGCTGGGTGACCCGATCGCCCGCCTGAAGCAGCACCTGATCAAGATCGGCCACTGGTCCGAGGAAGAACACCAGGCCGTCACGGCCGAGCTCGAAGCTGCGGTGATTGCCGcaca                                                                                                                  <  1:615387/150‑1 (MQ=255)
                                     tACCGCCCTGCCGACGACTGGAGCCACTTCCCGCTGGGTGACCCGATCGCCCGCCTGAAGCAGCACCTGATCAAGATCGGCCACTGGTCCGAGGAAGAACACCAGGCCGTGACGGcctggtg                                                                                                                                          <  1:558634/122‑6 (MQ=255)
                                     tACCGCCCTGCCGACGACTGGAGCCACTTCCCGCTGGGTGACCCGATCGCCCGCCTGAAGCAGCACCTGATCAAGATCGGCCACTGGTCCGAGGAAGAACACCAGGCCGTGACGGcctggtg                                                                                                                                          >  2:558634/1‑117 (MQ=255)
                                       ccGCCCTGCCGACGACTGGAGCCACTTCCCGCTGGGTGACCCGATCGCCCGCCTGAAGCAGCACCTGATCAAGATCGGCCACTGGTCCGAGGAAGAACACCAGGCCGTCACGGCCGAGCTCGAAGCTGCGGTGATTGCCGCACagaaaga                                                                                                            <  2:1382967/150‑1 (MQ=255)
                                           ccTGCCGACGACTGGAGCCACTTCCCGCTGGGTGACCCGATCGCCCGCCTGAAGCAGCACCTGATCAAGATCGGCCACTGGTCCGAGGAAGAACACCAGGCCGTCACGGCCGAGCTCGAAGCTGCGGTGATTGCCGCACAGAAAGAAGcc                                                                                                        <  2:433921/150‑1 (MQ=255)
                                             tGCCGACGACTGGAGCCACTTCCCGCTGGGTGACCCGATCGCCCGCCTGAAGCAGCACCTGATCAAGATCGGCCACTGGTCCGAGGAAGAACACCAGGCCGTCACGGCCGAGCTCGAAGCTGCGGTGATTGCCGCACAGAAAGAAGCCGa                                                                                                      <  2:27901/150‑1 (MQ=255)
                                              gCCGACGACTGGAGCCACTTCCCGCTGGGTGACCCGATCGCCCGCCTGAAGCAGCACCTGATCAAGATCGGCCACTGGTCCGAGGAAGAACACCAGGCCGTCACGGCCGAGCTCGAAGCTGCGGTGATTGCCGCACAGAAAGAAGCCGAg                                                                                                     <  2:41005/150‑1 (MQ=255)
                                                         gAGCCACTTCCCGCTGGGTGACCCGATCGCCCGCCTGAAGCAGCACCTGATCAAGATCGGCCACTGGTCCGAGGAAGAACACCAGGCCGTCACGGCCGAg                                                                                                                                            <  1:1284276/100‑1 (MQ=255)
                                                         gAGCCACTTCCCGCTGGGTGACCCGATCGCCCGCCTGAAGCAGCACCTGATCAAGATCGGCCACTGGTCCGAGGAAGAACACCAGGCCGTCACGGCCGAg                                                                                                                                            >  2:1284276/1‑100 (MQ=255)
                                                         gAGCCACTTCCCGCTGGGTGACCCGATCGCCCGCCTGAAGCAGCACCTGATCAAGATCGGCCACTGGTCCGAGGAAGAACACCAGGCCGTCACGGCCGAGATCGAAGCAAAAAAAATTGCCGCATAGAAAGAAACCGAGCAGTACTGCAc                                                                                          >  2:1374388/1‑150 (MQ=255)
                                                              aCTTCCCGCTGGGTGACCCGATCGCCCGCCTGAAGCAGCACCTGATCAAGATCGGCCACTGGTCCGAGGAAGAACACCAGGCCGTCACGGCCGAGCTCGAAGCTGCGGTGATTGCCGCACAGACAGAAGCCGAGCAGTACGGCACCCTgg                                                                                     >  1:50182/1‑150 (MQ=255)
                                                                  cccGCTGGGTGACCCGATCGCCCGCCTGAAGCAGCACCTGATCAAGATCGGCCACTGGTCCGAGGAAGAACACCAGGCCGTCACGGCCGAGCTCGAAGCTGCGGTGATTGCCGCACAGAAAGAAGCCGAGCAGTACGGCACCCTGGCCaa                                                                                 >  2:986365/1‑150 (MQ=255)
                                                                    ggcgggggGACCCGATCGGCCGCCTGAAGGAGCACGGGATGAAGAGCGGCCGCGGGTCGGAGGAAGAGCGCCAGGGCGTCACGGCCGAGCTGGAAGGTGCGGTGAGTGCCGCACAGAAAGAGGCCGAGCAGTACGGCACCCTGGCCAACg                                                                               <  2:460723/146‑1 (MQ=255)
                                                                    cGCTGGGTGACCCGATCGCCCGCCTGAAGCAGCACCTGATCAAGATCGGCCACTGGTCCGAGGAAGAACACCAGGCCGTCACGGCCGAGCTCGAAGCTGCGGTGATTGCCGCACAGAAAGAAGCCGAGCAGTACGGCAc                                                                                          >  1:845692/1‑139 (MQ=255)
                                                                    cGCTGGGTGACCCGATCGCCCGCCTGAAGCAGCACCTGATCAAGATCGGCCACTGGTCCGAGGAAGAACACCAGGCCGTCACGGCCGAGCTCGAAGCTGCGGTGATTGCCGCACAGAAAGAAGCCGAGCAGTACGGCAc                                                                                          <  2:845692/139‑1 (MQ=255)
                                                                    cGCTGGGTGACCCGATCGCCCGCCTGAAGCAGCACCTGATCAAGATCGGCCACTGGTCCGAGGAAGAACACCAGGCCGTCACGGCCGAGCTCGAAGCTGCGGTGATTGCCGCACAGAAAGAAGCCGAGCAGTACGGCACCCTGGCCAACg                                                                               >  2:249225/1‑150 (MQ=255)
                                                                     gctggctgACCCGATCGGCCGCCTAAAGCAGCACCTGATCAAGATCGGCCACTGGTCCGAGGAAGAACACCAGGCCGTCACGGCCGAGCTCGAAGCTGCGGTGATTGCCGCACAGAAAGAAGCCGAGCAGTACGGCACCCTGGCCAACgg                                                                              >  2:313896/1‑150 (MQ=255)
                                                                     cctcggTGACCCGATCCCCCGCCTGAAGCAGCACCTGATCAAGATCGGCCACTGGTCCGAGGAAGAACACCAGGCCGTCACGGCCGAGCTCGAAGCTGCGGTGATTGCCGCACAGAAAGAAGCCGAGCAGTACGGCACCCTGGCCAACgg                                                                              <  2:679156/146‑1 (MQ=255)
                                                                      cTGGGTGACCCGATCGCCCTCCTGAAGCAGCACCTGATCAAGATCGGCCACTGGTCCGAGGAAGAACACCAGGCCGTCACGGCCGAGCTCGAAGCTGCGGTGATTGCCGCACAGAAAGAAGCCGAGCAGTACGGCACCCTGGCCAACggg                                                                             <  2:933556/150‑1 (MQ=255)
                                                                            gACCCGATCGCCCGCCTGAAGCAGCACCTGATCAAGATCGGCCACTGGTCCGAGGAAGAACACCAGGCCGTCACGGACGAGCTCGAAGCTGCGGTGATTGCCGCACAGAAAGAAGCCGAGCAGTACGGCACCCTGGCCAACGGGCACATc                                                                       <  2:381854/150‑1 (MQ=255)
                                                                              cccGATCGCCCGCCTGAAGCAGCACCTGATCAAGATCGGCCACTGGTCCGAGGAAGAACACCAGGCCGTCACGGCCGAGCTCGAAGCTGCGGTGATTGCCGCACAGAAAGAAGCCGAGCAGTACGGCACCCTGGCCAACGGGCACATccc                                                                     >  2:451049/1‑150 (MQ=255)
                                                                               ccGATCGCCCGCCTGAAGCAGCACCTGATCAAGATCGGCCACTGGTCCGAGGAAGAACACCAGGCCGTCACGGCCGAGCTCGAAGCTGCGGTGATTGCCGCACAGAAAGAAGCCGAGCAGTACGGCACCCTGGCCAACGGGCACATCCCg                                                                    >  2:1352490/1‑150 (MQ=255)
                                                                                             gAAGCAGCACCTGATCAAGATCGGCCACTGGTCCGAGGAAGAACACCAGGCCGTCACGGCCGAGCTCGAAGCTGCGGTGATTGCCGCACAGAAAGAAGCCGAGCAGTACGGCACCCTGGCCAACggg                                                                             >  1:1011602/1‑127 (MQ=255)
                                                                                             gAAGCAGCACCTGATCAAGATCGGCCACTGGTCCGAGGAAGAACACCAGGCCGTCACGGCCGAGCTCGAAGCTGCGGTGATTGCCGCACAGAAAGAAGCCGAGCAGTACGGCACCCTGGCCAACggg                                                                             <  2:1011602/127‑1 (MQ=255)
                                                                                             gAAGCAGCACCTGATCAAGATCGGCCACTGGTCCGAGGAAGAACACCAGGCCGTCACGGCCGAGCTCGAAGCTGCGGTGATTGCCGCACAGAAAGAAGCCGAGCAGTACGGCACCCTGGCCAACGGGCACATCCCGAGCGCCGCCTCGAt                                                      >  1:179583/1‑150 (MQ=255)
                                                                                                 cagcaCCTGATCAAGATCGGCCACTGGTCCGAGGAAGAACACCAGGCCGTCACGGCCGAGCTCGAAGCTGCGGTGATTGCCGCACAGAAAGAAGCCGAGCAGTACGGCACCCTGGCCAACGGGCACATCCCGAGCGCCGCCTCGATGTTc                                                  <  2:179583/150‑1 (MQ=255)
                                                                                                  agcaCCTGATCAAGATCGGCCACTGGTCCGAGGAAGAACACCAGGCCGTCACGGCCGAGCTCGAAGCTGCGGTGATTGCCGCACAGAAAGAAGCCGAGCAGTACGGCACCCTGGCCAACgg                                                                              <  1:464372/121‑1 (MQ=255)
                                                                                                  agcaCCTGATCAAGATCGGCCACTGGTCCGAGGAAGAACACCAGGCCGTCACGGCCGAGCTCGAAGCTGCGGTGATTGCCGCACAGAAAGAAGCCGAGCAGTACGGCACCCTGGCCAACgg                                                                              >  2:464372/1‑121 (MQ=255)
                                                                                                  agcaCCTGATCAAGATCGGCCACTGGTCCGAGGAAGAACACCAGGCCGTCACGGCCGAGCTCGAAGCTGCGGTGATTGCCGCACAGAAAGAAGCCGAGCAGTACGGCACCCTGGCCAACGGGCACATCCCGAGCGCCGCCTCGATGTTCg                                                 >  2:195242/1‑150 (MQ=255)
                                                                                                   ccaCCTGATCAAGATCGGCCCCTGCTCCGAGGAAGAACCCCAGGCCGTCACGGCCGAGCTCGAAGCTGCGGTGATTGCCGCACAGAAAGAAGCCGAGCAGTACGGCACCCTGGCCAACGGGCACATCCCGAGCGCCGCCTCGATGTTCGa                                                <  1:654393/149‑1 (MQ=255)
                                                                                                    caCCTGATCAAGATCGGCCACTGGTCCGAGGAAGAACACCAGGCCGTCACGGCCGAGCTCGAAGCTGCGGTGATTGCCGCACAGAAAGAAGCCGAGCAGTACGGCACCCTGGCCAACGGGCACATCCCGAGCGCCGCCTCGATGTTCGAg                                               <  1:195242/150‑1 (MQ=255)
                                                                                                       cTGATCAAGATCGGCCACTGGTCCGAGGAAGAACACCAGGCCGTCACGGCCGAGCTCGAAGCTGCGGTGATTGCCGCACAGAAAGAAGCCGAGCAGTACGGCACCCTGGCCAACGGGCACATCCCGAGCGCCGCCTCGATGTTCGAGGAt                                            >  1:427561/1‑150 (MQ=255)
                                                                                                               gATCGGCCACTGGTCCGAGGAAGAACACCAGGCCGTCACGGCCGAGCTCGAAGCTGCGGTGATTGCCGCACAGAAAGAAGCCGAGCAGTACGGCACCCTGGCCAACGGGCACATCCCGAGCGCCGCCTCGATGTTc                                                  <  1:522243/136‑1 (MQ=255)
                                                                                                               gATCGGCCACTGGTCCGAGGAAGAACACCAGGCCGTCACGGCCGAGCTCGAAGCTGCGGTGATTGCCGCACAGAAAGAAGCCGAGCAGTACGGCACCCTGGCCAACGGGCACATCCCGAGCGCCGCCTCGATGTTc                                                  >  2:522243/1‑136 (MQ=255)
                                                                                                               gATCGGCCACTGGTCCGAGGAAGAACACCAGGCCGTCACGGCCGAGCTCGAAGCTGCGGTGATTGCCGCACAGAAAGAAGCCGAGCAGTACGGCACCCTGGCCAACGGGCACATCCCGAGCGCCGCCTCGATGTTCGAGGATGTGTACaa                                    <  2:574692/150‑1 (MQ=255)
                                                                                                                aTCGGCCACTGGTCCGAGGAAGAACACCAGGCCGTCACGGCCGAGCTCGAAGCTGCGGTGATTGCCGCACAGAAAGAAGCCGAGCAGTACGGCACCCTGGCCAACGGGCACATCCCGAGCGCCGCCTCGATGTTCGAGGATGTGTACAAg                                   <  1:1352490/150‑1 (MQ=255)
                                                                                                                aTCGGCCACTGGTCCGAGGAAGAACACCAGGCCGTCACGGCCGAGCTCGAAGCTGCGGTGATTGCCGCACAGAAAGAAGCCGAGCAGTACGGCACCCTGGCCAACGGGCACATCCCGAGCGCCGCCTCGATGTTCGAGGATGTGTACAAg                                   <  1:986365/150‑1 (MQ=255)
                                                                                                                   ggCCACTGGTCCGAGGAAGAACACCAGGCCGTCACGGCCGAGCTCGAAGCTGCGGTGATTGCCGCACAGAAAGAAGCCGAGCAGTACGGCACCCTGGCCAACGGGCACATCCCGAGCGCCGCCTCGATGTTCGAGGATGTGTACAAGGaa                                <  2:881772/150‑1 (MQ=255)
                                                                                                                     ccACTGGTCCGAGAAAGAACACCAGGCCGTCACGGCCGAGCTCGAAGCTGCGGTGATTGCCGCACAGAAAGAAGCCGAGCAGTACGGCACCCTGGCCAACGGGCACATCCCGAGCGCCGCCTCGATGTTCGAGGATGTGTACAAGGAAAt                              >  2:411324/1‑150 (MQ=255)
                                                                                                                         tGGTCCGAGGAAGAACACCAGGCCGTCACGGCCGAGCTCGAAGCTGCGGTGATTGCCGCACAGAAAGAAGCCGAGCAGTACGGCACCCTGGCCAACGGGCACATCCCGAGCGCCGCCTCGATGTTAGAGGATGTGTACAAGGAAATGccc                          <  2:128141/150‑1 (MQ=255)
                                                                                                                           gTCCGAGGAAGAACACCAGGCCGTCACGGCCGAGCTCGAAGCTGCGGTGATTGCCGCACAGAAAGAAGCCGAGCAGTACGGCACCCTGGCCAACGGGCACATCCCGAGCGCCGCCTCGATGTTCGAGGATGTGTACAAGGAAATGCCCGa                        <  1:68613/150‑1 (MQ=255)
                                                                                                                                     gaaCACCAGGCCGTCACGGCCGAGCTCGAAGCTGCGGTGATTGCCGCACAGAAAGAAGCCGAGCAGTACGGCACCCTGGCCAACGGGCACATCCCGAGCGCCGCCTCGATGTTCGAGGATGTGTACAAGGAAATGccc                          <  2:594234/138‑1 (MQ=255)
                                                                                                                                     gaaCACCAGGCCGTCACGGCCGAGCTCGAAGCTGCGGTGATTGCCGCACAGAAAGAAGCCGAGCAGTACGGCACCCTGGCCAACGGGCACATCCCGAGCGCCGCCTCGATGTTCGAGGATGTGTACAAGGAAATGccc                          >  1:594234/1‑138 (MQ=255)
                                                                                                                                       acacCAGGCCGTCACGGCCGAGCTCGAAGCTGCGGTGATTGCCGCACAGAAAGAAGCCGAGCAGTACGGCACCCTGGCCAACGGGCACATCCCGAGCGCCGCCTCGATGTTCGAGGATGTGTACAAGGAAATGCCCGACCACCTGCGCCg            <  2:137607/150‑1 (MQ=255)
                                                                                                                                         acCAGGCCGTCACGGCCGAGCTCGAAGCTGCGGTGATTGCCGCACAGCAAGAAGCCGAGCAGTACGGCACCCTGGCCAACGGGCACATCCCGAGCGCAGCCTCGATGTTCGAGGCTGTGTACACGGAAATGCCCGACCACCTGCGCCGTc          >  2:249475/1‑150 (MQ=255)
                                                                                                                                          ctcggCCGTGACGGCCGAGCTCGAAGCTGCGGTGATTGCCGCACAGAAAGAAGCCGAGCAGTACGGCACCCTGGCCAACGGGCACATCCCGAGCGCCGCCTCGATGTTCGa                                                <  2:944863/108‑1 (MQ=255)
                                                                                                                                          ctcggCCGTGACGGCCGAGCTCGAAGCTGCGGTGATTGCCGCACAGAAAGAAGCCGAGCAGTACGGCACCCTGGCCAACGGGCACATCCCGAGCGCCGCCTCGATGTTCGa                                                >  1:944863/4‑111 (MQ=255)
                                                                                                                                           tcggCCGTGACGGCCGAGCTCGAAGCTGCGGTGATTGCCGCACAGAAAGAAGCCGAGCAGTACGGCACCCTGGCCAACGGGCACATCCCGAGCGCCGCCTCGATGTTCGAGGATGTGTACAAGGAAATGCCCGACCACCTGCGCCGTCaa        >  1:62122/3‑150 (MQ=255)
                                                                                                                                           cAGGCCGTCACGGCCGAGCTCGAAGCTGCGGTGATTGCCGCACAGAAAGAAGCCGAGCAGTACGGCACCCTGGCCAACGGGCACATCCCGAGCGCCGCCTCGATGTTCGAGGATGTGTACAAGGAAATGCCCGACCACCTGCGCCGTCaa        >  2:559096/1‑150 (MQ=255)
                                                                                                                                              gCCGTCACGGCCGAGCTCGAAGCTGCGGTGATTGCCGCACAGAAAGAAGCCGAGCAGTACGGCACCCTGGCCAACGGGCACATCCCGAGCGCCGCCTCGATGTTCGAGGATGTGTACAAGGAAATGTCCGACCACCTGCGCCGTCAACGc     >  1:259986/1‑150 (MQ=255)
                                                                                                                                                 gTCACGGCCGAGCTCGAAGCTGCGGTGACTGCCGAACAGAAAGAAGCCGAGCAGTACGGCACCCTGGCCAACGGGCACATCCCGAGCTCAGCCTCGATGTTCGAGGATGTGTACAAGGAAATGCCCGACCACCTGCGCCGTCAACGCCAg  <  2:50182/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TGCCGGCCCGCACTCGACGTCGGACGACCCCTCCAAGTACCGCCCTGCCGACGACTGGAGCCACTTCCCGCTGGGTGACCCGATCGCCCGCCTGAAGCAGCACCTGATCAAGATCGGCCACTGGTCCGAGGAAGAACACCAGGCCGTCACGGCCGAGCTCGAAGCTGCGGTGATTGCCGCACAGAAAGAAGCCGAGCAGTACGGCACCCTGGCCAACGGGCACATCCCGAGCGCCGCCTCGATGTTCGAGGATGTGTACAAGGAAATGCCCGACCACCTGCGCCGTCAACGCCAG  >  NC_002947/4993015‑4993309

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: