Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 163,467 T→G 5.2% K123Q (AAG→CAG)  scpC ← propionyl‑CoA:succinate CoA transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947163,4670TG 5.2% 519.3 / 13.9 153K123Q (AAG→CAG) scpCpropionyl‑CoA:succinate CoA transferase
Reads supporting (aligned to +/- strand):  ref base T (69/76);  new base G (3/5);  total (72/81)
Fisher's exact test for biased strand distribution p-value = 7.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.55e-01

Only 100 of 159 total aligned reads displayed.

GAGGTTGGTGTTGTGCGAGAGGTTGATCTCCACAATCACCTGCTTGGCGAAGATCGCGAAGCTGGCCGAGTTGCCCACCGAAGTGGTTGGCACGATGTGGCCTTGCTCGGTGATGGCCACGGCTTCGATGACCGCGATGTCCGGCAGCTTCAGCTGCTGGTTGCGCAGTTGCTCGACGGTTTCCGACAGGTGCTGGTCAATGAACATCACCTGGCCGTCATTGATGGCTTTGCGCAGGGTGCTGTCGACCTGGAACGGCATGCGGCGGGCCAACACACCGGCCTCGGTCAGTTGTTTGT  >  NC_002947/163318‑163616
                                                                                                                                                     |                                                                                                                                                     
gAGGTTGGTGTTGTGCGAGAGGTTGATCTCCACAATCACCTGCTTGGCGAAGATCGCGAAGCTGGACGAGTTGCCCACCGAAGTGGTTGGCACGATGTGGCCTTGCTCGGTGATGGCCACGGCTTCGATGACCGCGATGTCCGGCAGCtt                                                                                                                                                       >  2:5019610/1‑150 (MQ=255)
  ggTTGGTGTTGTGCGAGAGGTTGATCTCCACAATCACCTGCTTGGCGAAGATCGCGAAGCTGGCCGAGTTGCCCACCGAAGTGGTTGGCACGATGTGGCCTTGCTCGGTGATGGCCACGGCTTCGATGACCGCGATGTCCGGCAGCTTCa                                                                                                                                                     <  1:1986235/150‑1 (MQ=255)
       gtgtTGTGCGAGAGGTTGATCTCCACAATCACCTGCTTGGCGAAGATCGCGAAGCTGGCCGAGTTGCCCACCGAAGTGGTTGGCACGATGTGGCCTTGCTCGGTGATGGCCACGGCTTCGATGACCGCGATGTCCGGCAGCTTCAgctgc                                                                                                                                                <  1:1950641/150‑1 (MQ=255)
        tgttgtGCGAGAGGTTGATCTCCACAATCACCTGCTTGGCGAAGATCGCGAAGCTGGCCGAGTTGCCCACCGAAGTGGTTGGCACGATGTGGCCTTGCTCGGTGATGGCCACGGCTTCGATGACCGCGATGTCCGGCAGCTTCAgctgct                                                                                                                                               >  2:4644365/1‑150 (MQ=255)
         gttgtGCGAGAGGTTGATCTCCACAATCACCTGCTTGGCGAAGATCGCGAAGCTGGCCGAGTTGCCCACCGAAGTGGTTGGCACGATGTGGCCTTGCTCGGTGATGGCCACGGCTTCGATGACCGCGATGTCCGGCAGCTTCAGctgctg                                                                                                                                              >  1:3890264/1‑150 (MQ=255)
         gttgtGCGAGAGGTTGATCTCCACAATCACCTGCTTGGCGAAGATCGCGAAGCTGGCCGAGTTGCCCACCGAAGTGGTTGGCACGATGTGGCCTTGCTCGGTGATGGCCACGGCTTCGATGACCGCGATGTCCGGCAGCTTCAGctgctg                                                                                                                                              >  1:1529324/1‑150 (MQ=255)
         gttgtGCGAGAGGTTGATCTCCACAATCACCTGCTTGGCGAAGATCGCGAAGCTGGCCGAGTTGCCCACCGAAGTGGTTGGCACGATGTGGCCTTGCTCGGTGATGGCCACGGCTTCGATGACCGCGATGTCCAGCAGCTGAAgctgccg                                                                                                                                              <  1:2478311/150‑3 (MQ=255)
         gttgtGCGAGAGGTTGATCTCCACAATCACCTGCTTGGCGAAGATCGCGAAGCTGGCCGAGTTGCCCACCGAAGTGGTTGGCACGATGTGGCCTTGCTCGGTGATGGCCACGGCTTCGATGACCGCGATGTCCAGCAGCTGAAgctgccg                                                                                                                                              <  1:3719528/150‑3 (MQ=255)
            gtgCGAGAGGTTGATCTCCACAATCACCTGCTTGGCGAAGATCGCGAAGCTGGCCGAGTTGCCCACCGAAGTGGTTGGCACGATGTGGCCTTGCTCGGTGATGGCCACGGCTTCGATGACCGCGATGTCCGGCAGCTTCAGCTGCTGGtt                                                                                                                                           <  1:4861127/150‑1 (MQ=255)
             tgCGAGAGGTTGATCTCCACAATCACCTGCTTGGCGAAGATCGCGAAGCTGGCCGAGTTGCCCACCGAAGTGGTTGGCACGATGTGGCCTTGCTCGGTGATGGCCACGGCTTCGATGACCGCGATGTCCGGCAGCTTCAGCTGCTGGTTg                                                                                                                                          <  2:4975870/150‑1 (MQ=255)
             tgCGAGAGGTTGATCTCCACAATCACCTGCTTGGCGAAGATCGCGAAGCTGGCCGAGTTGCCCACCGAAGTGGTTGGCACGATGTGGCCTTGCTCGGTGATGGCCACGGCTTCGATGACCGCGATGTCCGGCAGCTTCAGCTGCTGGTTg                                                                                                                                          <  2:4975753/150‑1 (MQ=255)
                     gTTGATCTCCACAATCACCTGCTTGGCGAAGATCGCGAAGCTGGCCGAGTTGCCCACCGAAGTGGTTGGCACGATGTGGCCTTGCTCGGTGATGGCCACGGCTTCGATGACCGCGATGTCCGGCAGCTTCAGCTGCTGGTTGCGCAGTTg                                                                                                                                  >  2:2285717/1‑150 (MQ=255)
                     gTTGATCTCCACAATCACCTGCTTGGCGAAGATCGCGAAGCTGGCCGAGTTGCCCACCGAAGTGGTTGGCACGATGTGGCCTTGCTCGGTGATGGCCACGGCTTCGATGACCGCGATGTCCGGCAGCTTCAGCTGCTGGTTGCGCAGTTg                                                                                                                                  >  1:2042545/1‑150 (MQ=255)
                        gATCTCCACAATCACCTGCTTGGCGAAGATCGCGAAGCTGGCCGAGTTGCCCACCGAAGTGGTTGGCACGATGTGGCCTTGCTCGGTGATGGCCACGGCTTCGATGACCGCGATGTCCGGCAGCTTCAGCTGCTGGTTGCGCAGTTGCTc                                                                                                                               >  2:4672733/1‑150 (MQ=255)
                        gATCTCCACAATCACCTGCTTGGCGAAGATCGCGAAGCTGGCCGAGTTGCCCACCGAAGTGGTTGGCACGATGTGGCCTTGCTCGGTGATGGCCACGGCTTCGATGACCGCGATGTCCGGCAGCTTCAGCTGCTGGTTGCGCAGTTGCTc                                                                                                                               <  2:1033194/150‑1 (MQ=255)
                          tctcCACAATCACCTGCTTGGCGAAGATCGCGAAGCTGGCCGAGTTGCCCACCGAAGTGGTTGGCACGATGTGGCCTTGCTCGGTGATGTCCACGGCTTCGATGACCGCGATGTCCGGCAGCTTCAGCTGCTGGTTGCGCAGTTGCTCGa                                                                                                                             >  1:945666/1‑150 (MQ=255)
                           ctcCACAATCACCTGCTTGGCGAAGATCGCGAAGCTGGCCGAGTTGCCCACCGAAGTGGTTGGCACGATGTGGCCTTGCTCGGTGATGGCCACGGCTTCGATGACCGCGATGTCCGGCAGCTTCAGCTGCTGGTTGCGCAGTTGCTCGAc                                                                                                                            >  2:1730413/1‑150 (MQ=255)
                            tcCACAATCACCTGCTTGGCGAAGATCGCGAAGCTGGCCGAGTTGCCCACCGAAGTGGTTGGCACGATGTGGCCTTGCTCGGTGATGGCCACGGCTTCGATGACCGCGATGTCCGGCAGCTTCAGCTGCTGGTTGCGCAGTTGCTCGAc                                                                                                                            >  1:3983433/1‑149 (MQ=255)
                            tcCACAATCACCTGCTTGGCGAAGATCGCGAAGCTGGCCGAGTTGCCCACCGAAGTGGTTGGCACGATGTGGCCTTGCTCGGTGATGGCCACGGCTTCGATGACCGCGATGTCCGGCAGCTTCAGCTGCTGGTTGCGCAGTTGCTCGAc                                                                                                                            <  2:3983433/149‑1 (MQ=255)
                              cacaATCACCTGCTTGGCGAAGATCGCGAAGCTGGCCGAGTTGCCCACCGAAGTGGTTGGCACGATGTGGCCTTGCTCGGTGATGGCCACGGCTTCGATGACCGCGATGTCCGGCAGCTTCAGCTGCTGGTTGCGCAGTTGCTCGACGGt                                                                                                                         >  2:4374964/1‑150 (MQ=255)
                               acaATCACCTGCTTGGCGAAGATCGCGAAGCTGGCCGAGTTGCCCACCGAAGTGGTTGGCACGATCTGGCCTTGCTCGGTGATGGCCACGGCTTCGATGACCGCGATGTCCGGCAGCTTCAGCTGCTGGTTGCGCAGTTGCTCGACGGtt                                                                                                                        >  2:1708500/1‑150 (MQ=255)
                                 aaTCACCTGCTTGGCGAAGATCGCGAAGCTGGCCGAGTTGCCCACCGAAGTGGTTGGCACGATGTGGCCTTGCTCGGTGATGGCCACGGCTTCGATGACCGCGATGTCCGGCAGCTTCAGCTGCTGGTTGCGCAGTTGCTCGACGGTTTc                                                                                                                      >  2:3187822/1‑150 (MQ=255)
                                 aaTCACCTGCTTGGCGAAGATCGCGAAGCTGGCCGAGTTGCCCACCGAAGTGGTTGGCACGATGTGGCCTTGCTCGGTGATGGCCACGGCTTCGATGACCGCGATGTCCGGCAGCTTCAGCTGCTGGTTGCGCAGTTGCTCGACGGTTTc                                                                                                                      >  2:3187795/1‑150 (MQ=255)
                                  aTCACCTGCTTGGCGAAGATCGCGAAGCTGGCCGAGTTGCCCACCGAAGTGGTTGGCACGATGTGGCCTTGCTCGGTGATGGCCACGGCTTCGATGACCGCGATGTCCGGCAGCTTCAGCTGCTGGTTGCGCAGTTGCTCGACGGTTTcc                                                                                                                     <  2:4245083/150‑1 (MQ=255)
                                   tCACCTGCTTGGCGAAGATCGCGAAGCTGGCCGAGTTGCCCACCGAAGTGGTTGGCACGATGTGGCCTTGCTCGGTGATGGCCACGGCTTCGATGACCGCGATGTCCGGCAGCTTCAGCTGCTGGTTGCGCAGTTGCTCGACGGTTTCCg                                                                                                                    >  1:2917842/1‑150 (MQ=255)
                                    cACCTGCTTGGCGAAGATCGCGAAGCTGGCCGAGTTGCCCACCGAAGTGGTTGGCACGATGTGGCCTTGCTCGGTGATGGCCACGGCTTCGATGACCGCGATGTCCGGCAGCTTCAGCTGCTGGTTGCGCAGTTGCTCGACGGTTTCCGa                                                                                                                   <  2:4088576/150‑1 (MQ=255)
                                      ccTGCTTGGCGAAGATCGCGAAGCTGGCCGAGTTGCCCACCGAAGTGGTTGGCACGATGTGGCCTTGCTCGGTGATGGCCACGGCTTCGATGACCGCGATGTCCGGCAGCTTCAGCTGCTGGTTGCGCAGTTGCTCGACGGTTTCCGACa                                                                                                                 <  1:56034/150‑1 (MQ=255)
                                       cTGCTTGGCGAAGATCGCGAAGCTGGCCGAGTTGCCCACCGAAGTGGTTGGCACGATGTGGCCTTGCTCGGTGATGGCCACGGCTTCGATGACCGCGATGTCCGGCAGCTTCAGCTGCTGGTTGCGCAGTTGCTCGACGGTTTCCGACAg                                                                                                                >  1:4972155/1‑150 (MQ=255)
                                       cTGCTTGGCGAAGATCGCGAAGCTGGCCGAGTTGCCCACCGAAGTGGTTGGCACGATGTGGCCTTGCTCGGTGATGGCCACGGCTTCGATGACCGCGATGTCCGGCAGCTTCAGCTGCTGGTTGCGCAGTTGCTCGACGGTTTCCGACAg                                                                                                                <  2:1220905/150‑1 (MQ=255)
                                       cTGCTTGGCGAAGATCGCGAAGCTGGCCGAGTTGCCCACCGAAGTGGTTGGCACGATGTGGCCTTGCTCGGTGATGGCCACGGCTTCGATGACCGCGATGTCCGGCAGCTTCAGCTGCTGGTTGCGCAGTTGCTCGACGGTTTCCGACAg                                                                                                                <  2:4292141/150‑1 (MQ=255)
                                           ttGGCGAAGATCGCGAAGCTGGCCGAGTTGCCCACCGAAGTGGTTGGCACGATGTGGCCTTGCTCGGTGATGGCCACGGCTTCGATGACCGCGATGTCCGGCAGCTTCAGCTGCTGGTTGCGCAGTTGCTCGACGGTTTCCGACAGGTGc                                                                                                            <  1:4895493/150‑1 (MQ=255)
                                             ggCGAAGATCGCGAAGCTGGCCGAGTTGCCCACCGAAGTGGTTGGCACGATGTGGCCTTGCTCGGTGATGGCCACGGCTTCGATGACCGCGATGTCCGGCAGCTTCAGCTGCTGGTTGCGCAGTTGCTCGACGGTTTCCGACAGGTGCTg                                                                                                          <  1:862300/150‑1 (MQ=255)
                                             ggCGAAGATCGCGAAGCTGGCCGAGTTGCCCACCGAAGTGGTTGGCACGATGTGGCCTTGCTCGGTGATGGCCACGGCTTCGATGACCGCGATGTCCGGCAGCTTCAGCTGCTGGTTGCGCAGTTGCTCGACGGTTTCCGACAGGTGCTg                                                                                                          <  2:5117038/150‑1 (MQ=255)
                                              gCGAAGATCGCGAAGCTGGCCGAGTTGCCCACCGAAGTGGTTGGCACGATGTGGCCTTGCTCGGTGATGGCCACGGCTTCGATGACCGCGATGTCCGGCAGCTTCAGCTGCTGGTTGCGCAGTTGCTCGACGGTTTCCGACAGGTGCTgg                                                                                                         >  1:4702776/1‑150 (MQ=255)
                                              gCGAAGATCGCGAAGCTGGCCGAGTTGCCCACCGAAGTGGTTGGCACGATGTGGCCTTGCTCGGTGATGGCCACGGCTTCGATGACCGCGATGTCCGGCAGCTTCAGCTGCTGGTTGCGCAGTTGCTCGACGGTTTCCGACAGGTGCTgg                                                                                                         <  1:435110/150‑1 (MQ=255)
                                                gAAGATCGCGAAGCTGGCCGAGTTGCCCACCGAAGTGGTTGGCACGATGTGGCCTTGCTCGGTGATGGCCACGGCTTCGATGACCGCGATGTCCGGCAGCTTCAGCTGCTGGTTGCGCAGTTGCTCGACGGTTTCCGACAGGTGCTGGTc                                                                                                       >  2:725555/1‑150 (MQ=255)
                                                gAAGATCGCGAAGCTGGCCGAGTTGCCCACCGAAGTGGTTGGCACGATGTGGCCTTGCTCGGTGATGGCCACGGCTTCGATGACCGCGATGTCCGGCAGCTTCAGCTGCTGGTTGCGCAGTTGCTCGACGGTTTCCGACAGGTGCTGGTc                                                                                                       >  1:1755209/1‑150 (MQ=255)
                                                 aaGATCGCGAAGCTGGCCGAGTTGCCCACCGAAGTGGTTGGCACGATGTGGCCTTGCTCGGTGATGGCCACGGCTTCGATGACCGCGATGTCCGGCAGCTTCAGCTGCTGGTTGCGCAGTTGCTCGACGGTTTCCGACAGGTGCTGGTCa                                                                                                      <  1:3427973/150‑1 (MQ=255)
                                                  aGATCGCGAAGCTGGCCGAGTTGCCCACCGAAGTGGTTGGCACGATGTGGCCTTGCTCGGTGATGGCCACGGCTTCGATGACCGCGATGTCCGGCAGCTTCAGCTGCTGGTTGCGCAGTTGCTCGACGGTTTCCGACAGGTGCTGGTCaa                                                                                                     <  1:1835486/150‑1 (MQ=255)
                                                   gATCGCGAAGCTGGCCGAGTTGCCCACCGAAGTGGTTGGCACGATGTGGCCTTGCTCGGTGATGGCCACGGCTTCGATGACCGCGATGTCCGGCAGCTTCAGCTGCTGGTTGCGCAGTTGCTCGACGGTTTCCGACAGGTGCTGGTCAAt                                                                                                    <  2:1956154/150‑1 (MQ=255)
                                                   gATCGCGAAGCTGGCCGAGTTGCCCACCGAAGTGGTTGGCACGATGTGGCCTTGCTCGGTGATGGCCACGGCTTCGATGACCGCGATGTCCGGCAGCTTCAGCTGCTGGTTGCGCAGTTGCTCGACGGTTTCCGACAGGTGCTGGTCAAt                                                                                                    <  1:2491275/150‑1 (MQ=255)
                                                    aTCGCGAAGCTGGCCGAGTTGCCCACCGAAGTGGTTGGCACGATGTGGCCTTGCTCGGTGATGGCCACGGCTTCGATGACCGCGATGTCCGGCAGCTTCAGCTGCTGGTTGCGCAGTTGCTCGACGGTTTCCGACAGGTGCTGGTCAATg                                                                                                   >  2:2448860/1‑150 (MQ=255)
                                                    aTCGCGAAGCTGGCCGAGTTGCCCACCGAAGTGGTTGGCACGATGTGGCCTTGCTCGGTGATGGCCACGGCTTCGATGACCGCGATGTCCGGCAGCTTCAGCTGCTGGTTGCGCAGTTGCTCGACGGTTTCCGACAGGTGCTGGTCAATg                                                                                                   >  1:4725964/1‑150 (MQ=255)
                                                     tCGCGAAGCTGGCCGAGTTGCCCACCGAAGTGGTTGGCACGATGTGGCCTTGCTCGGTGATGGCCACGGCTTCGATGACCGCGATGTCCGGCAGCTTCAGCTGCTGGTTGCGCAGTTGCTCGACGGTTTCCGACAGGTGCTGGTCAATGa                                                                                                  <  1:3820073/150‑1 (MQ=255)
                                                      cgcgAAGCTGGCCGAGTTGCCCACCGAAGTGGTTGGCACGATGTGGCCTTGCTCGGTGATGGCCACGGCTTCGATGACCGCGATGTCCGGCAGCTTCAGCTGCTGGTTGCGCAGTTGCTCGACGGTTTCCGACAGGTGCTGGTCAATGaa                                                                                                 <  1:2774129/150‑1 (MQ=255)
                                                            gCTGGCCGAGTTGCCCACCGAAGTGGTTGGCACGATGTGGCCTTGCTCGGTGATGGCCACGGCTTCGATGACCGCGATGTCCGGCAGCTTCAGCTGCTGGTTGCGCAGTTGCTCGACGGTTTCCGACAGGTGCTGGTCAATGAACATCAc                                                                                           >  2:3409954/1‑150 (MQ=255)
                                                                  cGAGTTGCCCACCGAAGTGGTTGTCACGATGTGGCCTTGCTCGGTGATGGCCACGGCTTCGATGACCGCGATGTCCGGCATCTTCAGCTGCTGGTTGCGCAGTTGCTCGACGGTTTCCGACAGGTGCTGGTCAATGATCATCACCTGGcc                                                                                     <  1:143580/150‑1 (MQ=255)
                                                                          ccACCGAAGTGGTTGGCACGATGTGGCCTTGCTCGGTGATGGCCACGGCTTCGATGACCGCGATGTCCGGCAGCTTCAGCTGCTGGTTGCGCAGTTGCTCGACGGTTTCCGACAGGTGCTGGTCAATGAACATCACCTGGCCGTCATTGa                                                                             <  1:4465165/150‑1 (MQ=255)
                                                                          ccACCGAAGTGGTTGGCACGATGTGGCCTTGCTCGGTGATGGCCACGGCTTCGATGACCGCGATGTCCGGCAGCTTCAGCTGCTGGTTGCGCAGTTGCTCGACGGTTTCCGACAGGTGCTGGTCAATGAACATCACCTGGCCGTCATTGa                                                                             >  1:4534214/1‑150 (MQ=255)
                                                                          ccACCGAAGTGGTTGGCACGATGTGGCCTTGCTCGGTGATGGCCACGGCTTCGATGACCGCGATGTCCGGCAGCTTCAGCTGCTGGTTGCGCAGTTGCTCGACGGTTTCCGACAGGTGCTGGTCAATGAACATCACCTGGCCGTCATTGa                                                                             <  2:1336852/150‑1 (MQ=255)
                                                                          ccACCGAAGTGGTTGGCACGATGTGGCCTTGCTCGGTGATGGCCACGGCTTCGATGACCGCGATGTCCGGCAGCTTCAGCTGCTGGTTGCGCAGTTGCTCGACGGTTTCCGACAGGTGCTGGTCAATGAACATCACCTGGCCGTCATTGa                                                                             >  1:1493212/1‑150 (MQ=255)
                                                                           cACCGAAGTGGTTGGCACGATGTGGCCTTGCTCGGTGATGTCCACGGCTTCGATGACCGCGATGTCCGGCAGCTTCAGCTGCTGGTTGCGCAGTTGCTCGACGGTTTCCGACAGGTGCTGGTCAATGAACATCACCTGGCCGTCATTGAt                                                                            <  2:945666/150‑1 (MQ=255)
                                                                           cACCGAAGTGGTTGGCACGATGTGGCCTTGCTCGGTGATGGCCACGGCTTCGATGACCGCGATGTCCGGCAGCTTCAGCTGCTGGTTGCGCAGTTGCTCGACGGTTTCCGACAGGTGCTGGTCAATGAACATCACCTGGCCGTCATTGAt                                                                            >  2:4358634/1‑150 (MQ=255)
                                                                           cACCGAAGTGGTTGGCACGATGTGGCCTTGCTCGGTGATGGCCACGGCTTCGATGACCGCGATGTCCGGCAGCTTCAGCTGCTGGTTGCGCAGTTGCTCGACGGTTTCCGACAGGTGCTGGTCAATGAACATCACCTGGCCGTCATTGAt                                                                            >  2:3158350/1‑150 (MQ=255)
                                                                           cACCGAAGTGGTTGGCACGATGTGGCCTTGCTCGGTGATGGCCACGGCTTCGATGACCGCGATGTCCGGCAGCTTCAGCTGCTGGTTGCGCAGTTGCTCGACGGTTTCCGACAGGTGCTGGTCAATGAACATCACCTGGCCGTCATTGAt                                                                            >  2:1549214/1‑150 (MQ=255)
                                                                            aCCGAAGTGGTTGGCACGATGTGGCCTTGCTCGGTGATGGCCACGGCTTCGATGACCGCGATGTCCGGCAGCTTCAGCTGCTGGTTGCGCGGTTGCTCGACGGTTTCCGACAGGTGCTGGTCAATGAACATCACCTGGCCGTCATTGATg                                                                           >  2:3680748/1‑150 (MQ=255)
                                                                            aCCGAAGTGGTTGGCACGATGTGGCCTTGCTCGGTGATGGCCACGGCTTCGATGACCGCGATGTCCGGCAGCTTCAGCTGCTGGTTGCGCAGTTGCTCGACGGTTTCCGACAGGTGCTGGTCAATGAACATCACCTGGCCGTCATTGATg                                                                           >  1:1935339/1‑150 (MQ=255)
                                                                            aCCGAAGTGGTTGGCACGATGTGGCCTTGCTCGGTGATGGCCACGGCTTCGATGACCGCGATGTCCGGCAGCTTCAGCTGCTGGTTGCGCAGTTGCTCGACGGTTTCCGACAGGTGCTGGTCAATGAACATCACCTGGCCGTCATTGATg                                                                           >  2:4072243/1‑150 (MQ=255)
                                                                                     gTTGGCACGATGTGGCCTTGCTCGGTGATGGCCACGGCTTCGATGACCGCGATGTCCGGCAGCTTCAGCTGCTGGTTGCGCAGTTGCTCGACGGTTTCCGACAGGTGCTGGTCAATGAACATCACCTGGCCGTCATTGATGGCTTTgcgc                                                                  >  1:2651612/1‑150 (MQ=255)
                                                                                      ttGGCACGATGTGGCCTTGCTCGGTGATGGCCACGGCTTCGATGACCGCGATGTCCGGCAGCTTCAGCTGCTGGTTGCGCAGTTGCTCGACGGTTTCCGACAGGTGCTGGTCAATGAACATCACCTGGCCGTCATTGATGGCTTTGCGCa                                                                 >  1:1364844/1‑150 (MQ=255)
                                                                                       tGGCACGATGTGGCCTTGCTCGGTGATGGCCACGGCTTCGATGACCGCGATGTCCGGCAGCTTCAGCTGCTGGTTGCGCAGTTGCTCGACGGTTTCCGACAGGTGCTGGTCAATGAACATCACCTGGCCGTCATTGATGGCTTTGCGCAg                                                                >  1:900893/1‑150 (MQ=255)
                                                                                       tGGCACGATGTGGCCTTGCTCGGTGATGGCCACGGCTTCGATGACCGCGATGTCCGGCAGCTTCAGCTGCTGGTTGCGCAGTTGCTCGACGGTTTCCGACAGGTGCTGGTCAATGAACATCACCTGGCCGTCATTGATGGCTTTGCGCAg                                                                <  1:4085241/150‑1 (MQ=255)
                                                                                          cACGATGTGGCCTTGCTCGGTGATGGCCACGGCTTCGATGACCGCGATGTCCGGCAGCTTCAGCTGCTGGTTGCGCAGTTGCTCGACGGTTTCCGACAGGTGCTGGTCAATGAACATCACCTGGCCGTCATTGATGGCTTTGCGCAGGGt                                                             <  1:3004850/150‑1 (MQ=255)
                                                                                           aCGATGTGGCCTTGCTCGGTGATGGCCACGGCTTCGATGACCGCGATGTCCGGCAGCTTCAGCTGCTGGTTGCGCAGTTGCTCGACGGTTTCCGACAGGTGCTGGTCAATGAACATCACCTGGCCGTCATTGATGGCTTTGCGCAGGGTg                                                            >  1:1036021/1‑150 (MQ=255)
                                                                                              aTGTGGCCTTGCTCGGTGATGGCCACGGCTTCGATGACCGCGATGTCCGGCAGCTTCAGCTGCTGGTTGCGCAGTTGCTCGACGGTTTCCGACAGGTGCTGGTCAATGAACATCACCTGGCCGTCATTGATGGCTTTGCGCAGGGTGCTg                                                         <  1:2402222/150‑1 (MQ=255)
                                                                                               tgtgGCCTTGCTCGGTGATGGCCACGGCTTCGATGACCGCGATGTCCGGCAGCTTCAGCTGCTGGTTGCGCAGTTGCTCGACGGTTTCCGACAGGTGCTGGTCAATGAACATCACCTGGCCGTCATTGATGGCTTTGCGCAGGGTGCTGt                                                        >  2:2759561/1‑150 (MQ=255)
                                                                                                gtgGCCTTGCTCGGTGATGGCCACGGCTTCGATGACCGCGATGTCCGGCAGCTTCAGCTGCTGGTTGCGCAGTTGCTCGACGGTTTCCGACAGGTGCTGGTCAATGAACATCACCTGGCCGTCATTGATGGCTTTGCGCAGGGTGCTGTc                                                       <  1:4266480/150‑1 (MQ=255)
                                                                                                gtgGCCTTGCTCGGTGATGGCCACGGCTTCGATGACCGCGATGTCCGGCAGCTTCAGCTGCTGGTTGCGCAGTTGCTCGACGGTTTCCGACAGGTGCTGGTCAATGAACATCACCTGGCCGTCATTGATGGCTTTGCGCAGGGTGCTGTc                                                       <  2:2181281/150‑1 (MQ=255)
                                                                                                     cTTGCTCGGTGATGGCCACGGCTTCGATGACCGCGATGTCCGGCAGCTTCTGCTTCTGGTTGCGCAGTTGCTCGACGGTTTCCGACAGGTGCTGGTCAATGAACATCACCTGGCCGTCATTGATGGCTTTGCGCAGGGTGCTGTCGACCt                                                  <  2:1852201/150‑1 (MQ=255)
                                                                                                     cTTGCTCGGTGATGGCCACGGCTTCGATGACCGCGATGTCCGGCAGCTTCAGCTGCTGGTTGCGCAGTTGCTCGACGGTTTCCGACAGGTGCTGGTCAATGAACATCACCTGGCCGTCATTGATGGCTTTGCGCAGGGTGCTGTCGACCt                                                  <  1:4910736/150‑1 (MQ=255)
                                                                                                         cTCGGTGATGGCCACGGCTTCGATGACCGCGATGTCCGGCAGCTTCAGCTGCTGGTTGCGCAGTTGCTCGACGGTTTCCGACAGGTGCTGGTCAATGAACATCACCTGGCCGTCATTGATGGCTTTGCGCAGGGTGCTGTCGACCTGGaa                                              >  2:1394242/1‑150 (MQ=255)
                                                                                                         cTCGGTGATGGCCACGGCTTCGATGACCGCGATGTCCGGCAGCTTCAGCTGCTGGTTGCGCAGTTGCTCGACGGTTTCCGACAGGTGCTGGTCAATGAACATCACCTGGCCGTCATTGATGGCTTTGCGCAGGGTGCTGTCGACCTGGaa                                              >  2:578588/1‑150 (MQ=255)
                                                                                                          tCGGTGATGGCCACGGCTTCGATGACCGCGATGTCCGGCAGCTTCAGCTGCTGGTTGCGCAGTTGCTCGACGGTTTCCGACAGGTGCTGGTCAATGAACATCACCTGGCCGTCATTGATGGCTTTGCGCAGGGTGCTGTCGACCTGGAAc                                             <  1:4369530/150‑1 (MQ=255)
                                                                                                             gTGATGGCCACGGCTTCGATGACCGCGATGTCCGGCAGCTTCAGCTGCTGGTTGCGCAGTTGCTCGACGGTTTCCGACAGGTGCTGGTCAATGAACATCACCTGGCCGTCATTGATGGCTTTGCGCAGGGTGCTGTCGACCTGGAACGGc                                          <  2:4387585/150‑1 (MQ=255)
                                                                                                                     cACGGCTTCGATGACCGCGATGTCCGGCAGCTTCAGCTGCTGGTTGCGCAGTTGCTCGACGGTTTCCGACAGGTGCTGGTCAATGAACATCACCTGGCCGTCATTGATGGCTTTGCGCAGGGTGCTGTCGACCTGGAACGGCATgcggcg                                  <  2:2837514/150‑1 (MQ=255)
                                                                                                                     cACGGCTTCGATGACCGCGATGTCCGGCAGCTTCAGCTGCTGGTTGCGCAGTTGCTCGACGGTTTCCGACAGGTGCTGGTCAATGAACATCACCTGGCCGTCATTGATGGCTTTGCGCAGGGTGCTGTCGACCTGGAACGGCATgcggcg                                  >  1:2513350/1‑150 (MQ=255)
                                                                                                                      aCGGCTTCGATGACCGCGATGTCCGGCAGCTTCAGCTGCTGGTTGCGCAGTTGCTCGACGGTTTCCGACAGGTGCTGGTCAATGAACATCACCTGGCCGTCATTGATGGCTTTGCGCAGAGTGCTGTCGACCTGGAACGGCATGcggcgg                                 <  2:2041150/150‑1 (MQ=255)
                                                                                                                       cGGCTTCGATGACCGCGATGTCCGGCAGCTTCAGCTGCTGGTTGCGCAGTTGCTCGACGGTTTCCGACAGGTGCTGGTCAATGAACATCACCTGGCCGTCATTGATGGCTTTGCGCAGGGTGCTGTCGACCTGGAACGGCATGCGGCggg                                <  2:3890264/150‑1 (MQ=255)
                                                                                                                       cGGCTTCGATGACCGCGATGTCCGGCAGCTTCAGCTGCTGGTTGCGCAGTTGCTCGACGGTTTCCGACAGGTGCTGGTCAATGAACATCACCTGGCCGTCATTGATGGCTTTGCGCAGGGTGCTGTCGACCTGGAACGGCATGCGGCggg                                <  1:3806626/150‑1 (MQ=255)
                                                                                                                       cGGCTTCGATGACCGCGATGTCCGGCAGCTTCAGCTGCTGGTTGCGCAGTTGCTCGACGGTTTCCGACAGGTGCTGGTCAATGAACATCACCTGGCCGTCATTGATGGCTTTGCGCAGGGTGCTGTCGACCTGGAACGGCATGCGGCggg                                <  2:843720/150‑1 (MQ=255)
                                                                                                                           ttCGATGACCGCGATGTCCGGCAGCTTCAGCTGCTGGTTGCGCAGTTGCTCGACGGTTTCCGACAGGTGCTGGTCAATGAACATCACCTGGCCGTCATTGATGGCTTTGCGCAGGGTGCTGTCGACCTGGAACGGCATGCGGCGGGCCaa                            <  1:1199109/150‑1 (MQ=255)
                                                                                                                            tCGATGACCGCGATGTCCGGCAGCTTCAGCTGCTGGTTGCGCAGTTGCTCGACGGTTTCCGACAGGTGCTGGTCAATGAACATCACCTGGCCGTCATTGATGGCTTTGCGCAGGGTGCTGTCGACCTGGAACGGCATGCGGCGGGCCAac                           <  1:2316914/150‑1 (MQ=255)
                                                                                                                                  aCCGCGATGTCCGGCAGCTTCAGCTGCTGGTTGCGCAGTTGCTCGACGGTTTCCGACAGGTGCTGGTCAATGAACATCACCTGGCCGTCATTGATGGCTTTGCGCAGGGTGCTGTCGACCTGGAACGGCATGCGGCGGGCCAACACACCg                     <  1:5064129/150‑1 (MQ=255)
                                                                                                                                    cgcgATGTCCGGCAGCTTCAGCTGCTGGTTGCGCAGTTGCTCGACGGTTTCCGACAGGTGCTGGTCAATGAACATCACCTGGCCGTCATTGATGGCTTTGCGCAGGGTGCTGTCGACCTGGAACGGCATGCGGCGGGCCAACACACCGGc                   >  1:497956/1‑150 (MQ=255)
                                                                                                                                    cgcgATGTCCGGCAGCTTCAGCTGCTGGTTGCGCAGTTGCTCGACGGTTTCCGACAGGTGCGGGTCAATGAACATCACCTGGCCGTCATTGATGGCTTTGCGCAGGGTGCTGTCGACCTGGAACGGCATGCGGCGGGCCAACACACCGGc                   <  2:4752551/150‑1 (MQ=255)
                                                                                                                                           tCCGGCAGCTTCAGCTGCTGGTTGCGCAGTTGCTCGACGGTTTCCGACAGGTGCTGGTCAATGAACATCACCTGGCCGTCATTGATGGCTTTGCGCAGGGTGCTGTCGACCTGGAACGGCATGCGGCGGGCCAACACACCGGCCTCGGTc            <  2:3703538/150‑1 (MQ=255)
                                                                                                                                           accagCAGCTGAAGCTGCCGGTTGCGCAGTTGCTCGACGGTTTCCGACAGGTGCTGGTCAATGAACATCACCTGGCCGTCATTGATGGCTTTGCGCAGGGTGCTGTCGACCTGGAACGGCATGCGGCGGGCCAACACACCGGCCTCGGTc            >  1:100285/5‑150 (MQ=255)
                                                                                                                                            ccGGCAGCTTCAGCTGCTGGTTGCGCAGTTGCTCGACGGTTTCCGACAGGTGCTGGTCAATGAACATCACCTGGCCGTCATTGATGGCTTTGCGCAGGGTGCTGTCGACCTGGAACGGCATGCGGCGGGCCAACACACCGGCCTCGGTCa           >  1:3780620/1‑150 (MQ=255)
                                                                                                                                              ggCAGCTTCAGCTGCTGGTTGCGCAGTTGCTCGACGGTTTCCGACAGGTTCTGGTCAATGAACATCACCTGGCCGTCATTGATGGCTTTGCGCAGGGTGCTGTCGACCTGGAACGGCATGCGGCGGGCCAACACACCGGCCTCGGTCAgt         <  1:5019610/150‑1 (MQ=255)
                                                                                                                                               gCAGCTTCAGCTGCTGGTTGCGCAGTTGCTCGACGGTTTCCGACAGGTGCTGGTCAATGAACATCACCTGGCCGTCATTGATGGCTTTGCGCAGGGTGCTGTCGACCTGGAACGGCATGCGGCGGGCCAACACACCGGCCTCGGTCAgtt        <  2:3647456/150‑1 (MQ=255)
                                                                                                                                               gCAGCTTCAGCTGCTGGTTGCGCAGTTGCTCGACGGTTTCCGACAGGTGCTGGTCAATGAACATCACCTGGCCGTCATTGATGGCTTTGCGCAGGGTGCTGTCGACCTGGAACGGCATGCGGCGGGCCAACACACCGGCCTCGGTCAgtt        <  1:1912265/150‑1 (MQ=255)
                                                                                                                                               gCAGCTTCAGCTGCTGGTTGCGCAGTTGCTCGACGGTTTCCGACAGGTGCTGGTCAATGAACATCACCTGGCCGTCATTGATGGCTTTGCGCAGGGTGCTGTCGACCTGGAACGGCATGCGGCGGGCCAACACACCGGCCTCGGTCAgtt        <  2:2042545/150‑1 (MQ=255)
                                                                                                                                               gCAGCTTCAGCTGCTGGTTGCGCAGTTGCTCGACGGTTTCCGACAGGTGCTGGTCAATGAACATCACCTGGCCGTCATTGATGGCTTTGCGCAGGGTGCTGTCGACCTGGAACGGCATGCGGCGGGCCAACACACCGGCCTCGGTCAgtt        <  1:4236454/150‑1 (MQ=255)
                                                                                                                                                 aGCTTCAGCTGCTGGTTGCGCAGTTGCTCGACGGTTTCCGACAGGTGCTGGTCAATGAACATCACCTGGCCGTCATTGATGGCTTTGCGCAGGGTGCTGTCGACCTGGAACGGCATGCGGCGGGCCAACACACCGGCCTCGGTCAgttgt      <  2:1803125/150‑1 (MQ=255)
                                                                                                                                                 aGCTTCAGCTGCTGGTTGCGCAGTTGCTCGACGGTTTCCGACAGGTGCTGGTCAATGAACATCACCTGGCCGTCATTGATGGCTTTGCGCAGGGTGCTGTCGACCTGGAACGGCATGCGGCGGGCCAACACACCGGCCTCGGTCAgttgt      <  2:1529324/150‑1 (MQ=255)
                                                                                                                                                  gCTTCAGCTGCTGGTTGCGCAGTTGCTCGACGGTTTCCGACAGGTGCTGGTCAATGAACATCACCTGGCCGTCATTGATGGCTTTGCGCAGGGTGCTGTCGACCTGGAACGGCATGCGGCGGGCCAACACACCGGCCTCGGTCAgttgtt     <  1:3731398/150‑1 (MQ=255)
                                                                                                                                                   cTTCAGCTGCTGGTTGCGCAGTTGCTCGACGGTTTCCGACAGGTGCTGGTCAATGAACATCACCTGGCCGTCATTGATGGCTTTGCGCAGGGTGCTGTCGACCTGGAACGGCATGCGGCGGGCCAACACACCGGCCTCGGTCAGttgttt    <  1:1708500/150‑1 (MQ=255)
                                                                                                                                                    ttCAGCTGCTGGTTGCGCAGTTGCTCGACGGTTTCCGACAGGTGCTGGTCAATGAACATCACCTGGCCGTCATTGATGGCTTTGCGCAGGGTGCTGTCGACCTGGAACGGCATGCGGCGGGCCAACACACCGGCCTCGGTCAGttgtttg   >  1:3060708/1‑150 (MQ=255)
                                                                                                                                                     tCAGCTGCTGGTTGCGCAGTTGCTCGACGGTTTCCGACAGGTGCTGGTCAATGAACATCACCTGGCCGTCATTGATGGCTTTGCGCAGGGTGCTGTCGACCTGGAACGGCATGCGGCGGGCCAACACACCGGCCTCGGTCAGttgtttgt  >  1:2030938/1‑150 (MQ=255)
                                                                                                                                                     tCAGCTGCTGGTTGCGCAGTTGCTCGACGGTTTCCGACAGGTGCTGGTCAATGAACATCACCTGGCCGTCATTGATGGCTTTGCGCAGGGTGCTGTCGACCTGGAACGGCATGCGGCGGGCCAACACACCGGCCTCGGTCAGgtgtttgt  >  2:3594199/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GAGGTTGGTGTTGTGCGAGAGGTTGATCTCCACAATCACCTGCTTGGCGAAGATCGCGAAGCTGGCCGAGTTGCCCACCGAAGTGGTTGGCACGATGTGGCCTTGCTCGGTGATGGCCACGGCTTCGATGACCGCGATGTCCGGCAGCTTCAGCTGCTGGTTGCGCAGTTGCTCGACGGTTTCCGACAGGTGCTGGTCAATGAACATCACCTGGCCGTCATTGATGGCTTTGCGCAGGGTGCTGTCGACCTGGAACGGCATGCGGCGGGCCAACACACCGGCCTCGGTCAGTTGTTTGT  >  NC_002947/163318‑163616

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