Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 36,758 T→A 15.6% D356V (GAT→GTT)  caiC ← putative crotonobetaine/carnitine‑CoA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091336,7580TA15.6% 2245.1 / inf 789D356V (GAT→GTT) caiCputative crotonobetaine/carnitine‑CoA ligase
Reads supporting (aligned to +/- strand):  ref base T (350/314);  new base A (61/62);  total (412/377)
Fisher's exact test for biased strand distribution p-value = 5.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.98e-01

Only 100 of 813 total aligned reads displayed.

GGTATCGCCGGTATGCAGCCAGCCATCGGCTTCCAGCACTTTCGCAGTGGCTTGTGGGTTGAGAAAGTACTCTTTGAAGATGGTTTTCCCAGGTATGCCTTTAATGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAGGACGATCGCCGATAATGCCCACAATGGTTTCCGTCATCCCATAAGACGTCAGCAAGCGAACGCCGAAGCGTTCACAAAA  >  NC_000913/36610‑36906
                                                                                                                                                    |                                                                                                                                                    
ggtatcGCCGGTATGCAGCCAGCCATCGGCTTCCAGCACTTTCGCAGTGGCTTGTGGGTTGAGAAAGTACTCTTTGAAGATGGTTTTCCCAGGTATGCCTTTAATGCAGATTTCACCGATCTCACCAGCCGGGAGTGGGCGATTGTGAt                                                                                                                                                      <  1:2919472/149‑1 (MQ=255)
ggtatcGCCGGTATGCAGCCAGCCATCGGCTTCCAGCACTTTCGCAGTGGCTTGTGGGTTGAGAAAGTACTCTTTGAAGATGGTTTTCCCAGGTATGCCTTTAATGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGTGATTGTGAt                                                                                                                                                      <  2:2753944/149‑1 (MQ=255)
          gTATGCAGCCAGCCATCGGCTTCCAGCACTTTCGCAGTGGCTTGTGGGTTGAGAAAGTACTCTTTGAAGATGGTTTTCCCAGGTATGCCTTTAATGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGAt                                                                                                                                            <  3:1399472/149‑1 (MQ=255)
             tGCAGCCAGCCATCGGCTTCCAGCACTTTCGCAGTGGCTTGTGGGTTGAGAAAGTACTCTTTGAAGATGGTTTTCCCAGGTATGCCTTTAATGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGAtctc                                                                                                                                         <  3:271566/149‑1 (MQ=255)
                agccagccATCGGCTTCCAGCACTTTCGCAGTGGCTTGTGGGTTGAGAAAGTACTCTTTGAAGATGGTTTTCCCAGGTATGCCTTTAATGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGc                                                                                                                                      >  2:5085029/1‑149 (MQ=255)
                     gccaTCGGCTTCCAGCACTTTCGCAGTGGCTTGTGGGTTGAGAAAGTACTCTTTGAAGATGGTTTTCCCAGGTATGCCTTTAATGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGt                                                                                                                                 >  4:2547333/1‑149 (MQ=255)
                      ccaTCGGCTTCCAGCACTTTCGCAGTGGCTTGTGGGTTGAGAAAGTACTCTTTGAAGATGGTTTTCCCAGGTATGCCTTTAATGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTa                                                                                                                                <  4:2001372/149‑1 (MQ=255)
                      ccaTCGGCTTCCAGCACTTTCGCAGTGGCTTGTGGGTTGAGAAAGTACTCTTTGAAGATGGTTTTCCCAGGTATGCCTTTAATGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGCGATTGTGAACGTCGCGGATCTCCGCTTCGTa                                                                                                                                <  3:3883497/149‑1 (MQ=255)
                       caTCGGCTTCCAGCACTTTCGCAGTGGCTTGTGGGTTGAGAAAGTACTCTTTGAAGATGGTTTTCCCAGGTATGCCTTTAATGCAGATTTCACCGATCTCACCAGCCGGGAGTGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAg                                                                                                                               >  2:3264856/1‑149 (MQ=255)
                       caTCGGCTTCCAGCACTTTCGCAGTGGCTTGTGGGTTGAGAAAGTACTCTTTGAAGATGGTTTTCCCAGGTATGCCTTTAATGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAg                                                                                                                               >  1:3533561/1‑149 (MQ=255)
                        aTCGGCTTCCAGCACTTTCGCAGTGGCTTGTGGGTTGAGAAAGTACTCTTTGAAGATGGTTTTCCCAGGTATGCCTTTAATGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGc                                                                                                                              >  4:1216808/1‑149 (MQ=255)
                                 cAGCACTTTCGCAGTGGCTTGTGGGTTGAGAAAGTACTCTTTGAAGATGGTTTTCCCAGGTATGCCTTTAATGCAGATTTCACCGATCTCACCAGCCGGGAGTGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAAc                                                                                                                         >  1:2167185/1‑145 (MQ=255)
                                 cAGCACTTTCGCAGTGGCTTGTGGGTTGAGAAAGTACTCTTTGAAGATGGTTTTCCCAGGTATGCCTTTAATGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCa                                                                                                                     >  4:3036198/1‑149 (MQ=255)
                                 cAGCACTTTCGCAGTGGCTTGTGGGTTGAGAAAGTACTCTTTGAAGATGGTTTTCCCAGGTATGCCTTTAATGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCa                                                                                                                     >  4:2102519/1‑149 (MQ=255)
                                 cAGCACTTTCGCAGTGGCTTGTGGGTTGAGAAAGTACTCTTTGAAGATGGTTTTCCCAGGTATGCCTTTAATGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGCGATTGTGAACGTCGCGGATCTCCGCTTCGTAGCAAAACCCCa                                                                                                                     <  2:3671677/149‑1 (MQ=255)
                                   gCACTTTCGCAGTGGCTTGTGGGTTGAGAAAGTACTCTTTGAAGATGGTCTCCCCAGGTATGCCTTTAATGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCAcc                                                                                                                   >  1:1706664/1‑149 (MQ=255)
                                     aCTTTCGCAGTGGCTTGTGGGTTGAGAAAGTACTCTTTGAAGATGGTTTTCCCAGGTATGCCTTTAATGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCAcc                                                                                                                   <  3:4676175/147‑1 (MQ=255)
                                     aCTTTCGCAGTGGCTTGTGGGTTGAGAAAGTACTCTTTGAAGATGGTTTTCCCAGGTATGCCTTTAATGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCg                                                                                                                 >  2:741712/1‑149 (MQ=255)
                                       tttCGCAGTGGCTTGTGGGTTGAGAAAGTACTCTTTGAAGATGGTTTTCCCAGGTATGCCTTTAATGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGAc                                                                                                               <  2:4266560/149‑1 (MQ=255)
                                       tttCGCAGTGGCTTGTGGGTTGAGAAAGTACTCTTTGAAGATGGTTTTCCCAGGTATGCCTTTAATGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGAc                                                                                                               >  1:5970980/1‑149 (MQ=255)
                                           gCAGTGGCTTGTGGGTTGAGAAAGTACTCTTTGAAGATGGTTTTCCCAGGTATGCCTTTAATGCAGATTTCACCGATCTCACCAGCCGGGAGTGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAAt                                                                                                           <  4:603246/149‑1 (MQ=255)
                                              gTGGCTTGTGGGTTGAGAAAGTACTCTTTGAAGATGGTTTTCCCAGGTATGCCTTTAATGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGCTCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGa                                                                                                        >  4:4662204/1‑149 (MQ=255)
                                                ggCTTGTGGGTTGAGAAAGTACTCTTTGAAGATGGTTTTCCCAGGTATGCCTTTAATGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACg                                                                                                      <  3:2389568/149‑1 (MQ=255)
                                                 gCTTGTGGGTTGAGAAAGTACTCTTTGAAGATGGTTTTCCCAGGTATGCCTTTAATGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACgg                                                                                                     <  4:804225/149‑1 (MQ=255)
                                                   ttGTGGGTTGAGAAAGTACTCTTTGAAGATGGTTTTCCCAGGTATGCCTTTAATGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGcc                                                                                                   >  4:5697228/1‑149 (MQ=255)
                                                         gTTGAGAAAGTACTCTTTGAAGATGGTTTTCCCAGGTATGCCTTTAATGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCTCCCGCCCCACCGACGGCCAGcgtc                                                                                             <  4:6461538/149‑1 (MQ=255)
                                                         gTTGAGAAAGTACTCTTTGAAGATGGTTTTCCCAGGTATGCCTTTAATGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGcgtc                                                                                             >  2:3004914/1‑149 (MQ=255)
                                                          ttGAGAAAGTACTCTTTGAAGATGGTTTTCCCAGGTATGCCTTTAATGCAGATTTCACCGATCTCACCAGCCGGGAGTGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGcgtcg                                                                                            >  1:4653512/1‑149 (MQ=255)
                                                          ttGAGAAAGTACTCTTTGAAGATGGTTTTCCCAGGTATGCCTTTAATGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGCGATTGTGAACGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACgg                                                                                                     <  3:3246018/140‑1 (MQ=255)
                                                           tGAGAAAGTACTCTTTGAAGATGGTTTTCCCAGGTATGCCTTTAATGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGcgtcgt                                                                                           <  1:6321365/149‑1 (MQ=255)
                                                            gagaAAGTACTCTTTGAAGATGGTTTTCCCAGGTATGCCTTTAATGCAGATTTCACTGATCTCACCAGCCGGGAGTGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGtt                                                                                          >  2:215292/1‑149 (MQ=255)
                                                              gaAAGTACTCTTTGAAGATGGTTTTCCCAGGTATGCCTTTAATGCAGATTTCACCGATCTCACCAGCCGGGAGTGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTa                                                                                        >  3:3454396/1‑149 (MQ=255)
                                                                aaGTACTCTTTGAAGATGGTTTTCCCAGGTATGCCTTTAATGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGCGATTGTGAACGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATc                                                                                      <  4:1776035/149‑1 (MQ=255)
                                                                    aCTCTTTGAAGATGGTTTTCCCAGGTATGCCTTTAATGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCa                                                                                  <  3:5244800/149‑1 (MQ=255)
                                                                     ctctTTGAAGATGGTTTTCCCAGGTATGCCTTTAATGCAGATTTCACCGATCTCACCAGCCGGGAGTGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAg                                                                                 <  3:298446/149‑1 (MQ=255)
                                                                     ctctTTGAAGATGGTTTTCCCAGGTATGCCTTTAATGCAGATTTCACCGATCTCACCAGCCGGGAGTGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAg                                                                                 <  3:4970866/149‑1 (MQ=255)
                                                                        tttGAAGATGGTTTTCCCAGGTATGCCTTTAATGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAGGAc                                                                              >  3:2194928/1‑149 (MQ=255)
                                                                          tGAAGATGGTTTTCCCAGGTATGCCTTTAATGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAGGACGa                                                                            <  3:6881709/149‑1 (MQ=255)
                                                                           gAAGATGGTTTTCCCAGGTATGCCTTTAATGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAGGACGAt                                                                           <  1:5085029/149‑1 (MQ=255)
                                                                             aGATGGTTTTCCCAGGTATGCCTTTAATGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGCGATTGTGAACGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAGGACGATCg                                                                         >  1:3218412/1‑149 (MQ=255)
                                                                              gATGGTTTTCCCAGGTATGCCTTTAATGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAGGACGATCGc                                                                        <  3:4384078/149‑1 (MQ=255)
                                                                              gATGGTTTTCCCAGGTATGCCTTTAATGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGCGATTGTGAACGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACgg                                                                                                     <  1:4246749/120‑1 (MQ=255)
                                                                               aTGGTTTTCCCAGGTATGCCTTTAATGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAGGACGATCGcc                                                                       >  1:3806887/1‑149 (MQ=255)
                                                                                 ggTTTTCCCAGGTATGCCTTTAATGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAGGACGATCGCCGa                                                                     <  3:916362/149‑1 (MQ=255)
                                                                                       cccAGGTATGCCTTTAATGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGCGATTGTGAACGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAGGACGATCGCCGATAATGc                                                               <  1:1812424/149‑1 (MQ=255)
                                                                                        ccAGGTATGCCTTTAATGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGtt                                                                                          >  1:6569779/1‑121 (MQ=255)
                                                                                        ccAGGTATGCCTTTAATGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAGGACGATCGCCGATAATGcc                                                              <  2:6014994/149‑1 (MQ=255)
                                                                                          aGGTATGCCTTTAATGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAGGACGATCGCCGATAATGCCca                                                            >  1:2562504/1‑149 (MQ=255)
                                                                                           ggTATGCCTTTAATGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGCGATTGTGAACGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAGGACGATCGCCGATAATGCCcac                                                           <  1:7098405/149‑1 (MQ=255)
                                                                                             tATGCCTTTAATGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGCGATTGTGAACGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAGGACGATCGCCGATAATGCCCACaa                                                         <  3:6080399/149‑1 (MQ=255)
                                                                                              aTGCCTTTAATGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAGGACGATCGCCGATAATGCCCACAAt                                                        <  2:1267239/149‑1 (MQ=255)
                                                                                               tGCCTTTAATGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAGGACGATCGCCGATAATGCCCACAATg                                                       <  3:1792629/149‑1 (MQ=255)
                                                                                                gCCTTTAATGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAGGACGATCGCCGATAATGCCCACAATg                                                       <  4:2294402/148‑1 (MQ=255)
                                                                                                 ccTTTAATGCAGATTTCACCGATCTCACCAGCCGGGAGTGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAGGACGATCGCCGATAATGCCCACAATGGt                                                     >  2:1626646/1‑149 (MQ=255)
                                                                                                 ccTTTAATGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGCGATTGTGAACGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAGGACGATCGCCGATAATGCCCACAATGGt                                                     >  3:5533773/1‑149 (MQ=255)
                                                                                                  cTTTAATGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAGGACGATCGCCGATAATg                                                                >  3:2003795/1‑137 (MQ=255)
                                                                                                   tttAATGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAGGACGATCGCCGATAAt                                                                 >  1:1587965/1‑135 (MQ=255)
                                                                                                   tttAATGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAGGACGATCGCCGATAATGCCCACAATGGttt                                                   >  2:5584618/1‑149 (MQ=255)
                                                                                                        tGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGCGATTGTGATCATCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAGGACGATCGCCGATAATGCCCACAATgg                                                      >  4:6970952/1‑141 (MQ=255)
                                                                                                         gCAGATTTCACCGATCTCACCAGCCGGGAGTGGGCGATTGTGATCGTCGCGGTTCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAGGACGATCGCCGATAATGCCCACAATGGTTTCCGTCa                                             >  2:6640767/1‑149 (MQ=255)
                                                                                                           aGATTTCACCGATCTCACCAGCCGGGAGTGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAGGACGATCGCCGATAATGCCCACAATGGTTTCCGTCATc                                           <  2:3497135/149‑1 (MQ=255)
                                                                                                           aGATTTCACCGATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAGGACGATCGCCGATAATGCCCACAATGGTTTCCGTCATc                                           <  2:4485135/149‑1 (MQ=255)
                                                                                                           aGATTTCACCGATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAGGACGATCGCCGATAATGCCCACAATGGTTTCCGTCATc                                           <  4:4280566/149‑1 (MQ=255)
                                                                                                            gATTTCACCGATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAGGACGATCGCCGATAATGCCCACAATGGTTTCCGTCATcc                                          >  3:4142893/1‑149 (MQ=255)
                                                                                                             aTTTCACCGATCTCACCAGCCGGGAGTGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAGGACGATCGCCGATAATGCCCACAATGGTTTCCGTCATccc                                         <  1:812102/149‑1 (MQ=255)
                                                                                                             aTTTCACCGATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAGGACGATCGCCGATAATGCCCACAATGGTTTCCGTCATccc                                         <  2:6372797/149‑1 (MQ=255)
                                                                                                              tttCACCGATCTCACCAGCCGGGAGTGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAGGACGATCGCCGATAATGCCCACAATGGTTTCCGTCATCCCa                                        >  4:5045839/1‑149 (MQ=255)
                                                                                                              tttCACCGATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAGGACGATCGCCGATAATGCCCACAATGGTTTCCGTCATCCCa                                        >  3:3709726/1‑149 (MQ=255)
                                                                                                                  aCCGATCTCACCAGCCGGGAGTGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAGGACGATCGCCGATAATGCCca                                                            >  4:3828340/1‑125 (MQ=255)
                                                                                                                  aCCGATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAGGACGATCGCCGATAATGCCCACAATGGTTTCCGTCATCCCATAAg                                    <  2:5285804/149‑1 (MQ=255)
                                                                                                                  aCCGATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAGGACGATCGCCGATAATGCCCACAATGGTTTCCGTCATCCCATAAg                                    >  3:6556324/1‑149 (MQ=255)
                                                                                                                   ccGATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAGGACGATCGCCGATAATGCCCACAATGGTTTCCGTCATCCCATAAGa                                   <  1:4029014/149‑1 (MQ=255)
                                                                                                                   ccGATCTCACCAGCCGGGAGCGGGCGATTGTGAACGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAGGACGATCGCCGATAATGCCCACAATGGTTTCCGTCATCCCATAAGa                                   <  1:1412920/149‑1 (MQ=255)
                                                                                                                      aTCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAGGACGATCGCCGATAATGCCCACAATGGTTTCCGTCATCCCATAAGACGt                                >  2:4666526/1‑149 (MQ=255)
                                                                                                                       tctcACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAGGACGATCGCCGATAATGCCCACAATGGTTTCCGTCATCCCATAAGACGTc                               <  1:2371811/149‑1 (MQ=255)
                                                                                                                       tctcACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAGGACGATCGCCGATAATGCCCACAATGGTTTCCGTCATCCCATAAGACGTc                               <  2:4077114/149‑1 (MQ=255)
                                                                                                                       tctcACCAGCCGGGAGCGGGCGATTGTGAACGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAGGACGATCGCCGATAATGCCCACAATGGTTTCCGTCATCCCATAAGACGTc                               <  4:321958/149‑1 (MQ=255)
                                                                                                                        ctcACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAGGACGATCGCCGATAATGCCCACAATGGTTTCCGTCATCCCATAAGACGTCa                              <  2:6987418/149‑1 (MQ=255)
                                                                                                                         tcACCAGCCGGGAGTGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAGGACGATCGCCGATAATGCCCACAATGGTTTCCGTCATCCCATAAGACGTCAg                             >  4:3471639/1‑149 (MQ=255)
                                                                                                                         tcACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAGGACGATCGCCGATAATGCCCACAATGGTTTCCGTCATCCCATAAGACGTCAg                             >  1:6701454/1‑149 (MQ=255)
                                                                                                                         tcACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAGGACGATCGCCGATAATGCCCACAATGGTTTCCGTCATCCCATAAGACGTCAg                             >  4:1504537/1‑149 (MQ=255)
                                                                                                                         tcACCAGCCGGGAGCGGGCGATTGTGAACGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAGGACGATCGCCGATAATGCCCACAATGGTTTCCGTCATCCCATAAGACGTCAg                             >  3:6354062/1‑149 (MQ=255)
                                                                                                                          cACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAGGACGATTGCCGATAATGCCCACAATGGTTTCCGTCATCCCATAAGACGTCAGc                            >  2:1877675/1‑149 (MQ=255)
                                                                                                                            ccAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAGGACGATCGCCGATAATGCCCACAATGGTTTCCGTCATCCCATAAGACGTCAGCaa                          >  1:5563762/1‑149 (MQ=255)
                                                                                                                            ccAGCCGGGAGCGGGCGATTGTGAACGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAGGACGATCGCCGATAATGCCCACaa                                                         <  2:932226/118‑1 (MQ=255)
                                                                                                                              aGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAGGACGATCGCCGATAATGCCCACAATGGTTTCCGTCATCCCATAAGACGTCAGCAAGc                        >  1:5350295/1‑149 (MQ=255)
                                                                                                                                    gAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAGGACGATCGCCGATAATGCCCACAATGGTTTCCGTCATCCCATAAGACGTCAGCAAGCGAACGc                  <  2:3858442/149‑1 (MQ=255)
                                                                                                                                    gAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAGGACGATCGCCGATAATGCCCACAATGGTTTCCGTCATCCCATAAGACGTCAGCAAGCGAACGc                  <  1:894464/149‑1 (MQ=255)
                                                                                                                                        gggCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAGGACGATCGCCGATAATGCCCACAATGGTTTCCGTCATCCCATAAGACGTCAGCAAGCGAACGCCGaa              <  4:6556324/149‑1 (MQ=255)
                                                                                                                                         ggCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAGGACGATCGCCGATAATGCCCACAATGGTTTcc                                                 >  4:4041102/1‑113 (MQ=255)
                                                                                                                                             aTTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAGGACGATCGCCGATAATGCCCACAATGGTTTCCGTCATCCCATAAGACGTCAGCAAGCGAACGCCGAAGCGtt         <  3:7063960/149‑1 (MQ=255)
                                                                                                                                              ttGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAGGACGATCGCCGATAATGCCCACAATGGTTTCCGTCATCCCATAAGACGTCAGCAAGCGAACGCCGAAGCGTTc        >  2:5807690/1‑149 (MQ=255)
                                                                                                                                               tgtgATCGTCGCGGATCTCCGCTTCGTAGCAAAGCCCCTCCCGGCCGATCGACGGCCAGCTTCGTTTATCGCCAGGACGATCGCCGATACTGCCCACAAGGGTTTCCGTCACCCCATGAGGCGTCAGCAAGCGAACGCCGAAGCGTTca       <  4:6685828/149‑1 (MQ=255)
                                                                                                                                               tgtgATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAGGACGATCGCCGATAATGCCCACAATGGTTTCCGTCATCCCATAAGACGTCAGCAAGCGAACGCCGAAGCGTTca       >  4:2691710/1‑149 (MQ=255)
                                                                                                                                                  gATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAGGACGATCGCCGATAATGCCCACAATGGTTTCCGTCATCCCATAAGACGTCAGCAAGCGAACGCCGAAGCGTTca       >  1:5790886/1‑146 (MQ=255)
                                                                                                                                                  gATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAGGACGATCGCCGATAATGCCCACAATGGTTTCCGTCATCCCATAAGACGTCAGCAAGCGAACGCCGAAGCGTTCACaa    >  1:3443687/1‑149 (MQ=255)
                                                                                                                                                  gATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAGGACGATCGCCGATAATGCCCACAATGGTTTCCGTCATCCCATAAGACGTCAGCAAGCGAACGCCGAAGCGTTCACaa    >  2:4940653/1‑149 (MQ=255)
                                                                                                                                                    tcgtcgCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAGGACGATCGCCGATAATGCCCACAATGGTTTCCGTCATCCCATAAGACGTCAGCAAGCGAACGCCGAAGCGTTCACaaaa  <  4:5382054/149‑1 (MQ=255)
                                                                                                                                                    tcgtcgCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAGGACGATCGCCGATAATGCCCACAATGGTTTCCGTCATCCCATAAGACGTCAGCAAGCGAACGCCGAAGCGTTCACaaaa  <  2:2426902/149‑1 (MQ=255)
                                                                                                                                                    tcgtcgCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGCCGTTTATCGCCAGGACGATCGCCGATAATGCCCACAATGGTTTCCGTCATCCCATAAGACGTCATCAAGCGAACGCCGAAGCGTTCACaaaa  >  1:1024663/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GGTATCGCCGGTATGCAGCCAGCCATCGGCTTCCAGCACTTTCGCAGTGGCTTGTGGGTTGAGAAAGTACTCTTTGAAGATGGTTTTCCCAGGTATGCCTTTAATGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAGGACGATCGCCGATAATGCCCACAATGGTTTCCGTCATCCCATAAGACGTCAGCAAGCGAACGCCGAAGCGTTCACAAAA  >  NC_000913/36610‑36906

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