Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,983,573 C→T 7.0% L184L (CTG→CTA) ‡ kduI ← hexuronate isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,983,5730CT 7.0% 270.5 / 8.5 86L184L (CTG→CTA) ‡kduIhexuronate isomerase
Reads supporting (aligned to +/- strand):  ref base C (37/43);  new base T (3/3);  total (40/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.40e-01

GCATCACAATATGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTGGTGAGGTTATCGCCTAACGTGACTGGAGATACTTCGTCCGGTGTGACTTTTTTGGTGGGATACGTC  >  NC_000913/2983431‑2983720
                                                                                                                                              |                                                                                                                                                   
gCATCACAATATGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTcc                                                                                                                                               <  3:856750/149‑1 (MQ=255)
  aTCACAATATGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTcc                                                                                                                                               <  3:368520/147‑1 (MQ=255)
  aTCACAATATGACGCGTTTCTTGCGGCTGCCCCATAATGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTcc                                                                                                                                               >  4:368520/1‑147 (MQ=255)
   tCACAATATGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTc                                                                                                                                            >  3:17363/1‑149 (MQ=255)
     acaATATGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCTGGCTCCgggta                                                                                                                                          <  4:493049/149‑5 (MQ=255)
     acaATATGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAg                                                                                                                                          >  3:323944/1‑149 (MQ=255)
       aaTATGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCTGGCTCCgggt                                                                                                                                           >  2:757446/1‑143 (MQ=255)
       aaTATGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCTGGCTCCgggt                                                                                                                                           <  1:757446/146‑4 (MQ=255)
       aaTATGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGcc                                                                                                                                        >  2:323704/1‑149 (MQ=255)
         tatGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGcccc                                                                                                                                      <  2:635513/149‑1 (MQ=255)
           tGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTTGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCAt                                                                                                                                    <  4:405652/149‑1 (MQ=255)
                 gTTTCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCaa                                                                                                                              <  2:387976/149‑1 (MQ=255)
                 gTTTCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCaa                                                                                                                              >  2:38175/1‑149 (MQ=255)
                 gTTTCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCaa                                                                                                                              >  1:561220/1‑149 (MQ=255)
                  tttCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAAt                                                                                                                             <  4:323944/149‑1 (MQ=255)
                      ttGCGGCTGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGTGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCCGCCCCATACCCAATTGGc                                                                                                                         <  4:457649/149‑1 (MQ=255)
                           gCTGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGtt                                                                                                                    >  3:25160/1‑149 (MQ=255)
                            cTGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGttt                                                                                                                   >  1:53264/1‑149 (MQ=255)
                             tGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTc                                                                                                                  >  3:860563/1‑149 (MQ=255)
                             tGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTc                                                                                                                  <  1:367478/149‑1 (MQ=255)
                               ccccATCATGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCa                                                                                                                <  4:193473/149‑1 (MQ=255)
                               ccccATCATGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCa                                                                                                                <  3:326048/149‑1 (MQ=255)
                                cccATCATGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAg                                                                                                               <  3:542989/149‑1 (MQ=255)
                                cccATCATGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAg                                                                                                               >  4:289100/1‑149 (MQ=255)
                                 ccATCATGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGt                                                                                                              >  3:91786/1‑149 (MQ=255)
                                  catcatGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTa                                                                                                             >  2:17688/1‑149 (MQ=255)
                                     catGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACAt                                                                                                          <  2:53264/149‑1 (MQ=255)
                                       tgtgGAAATCGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATcc                                                                                                        >  3:815014/1‑149 (MQ=255)
                                         tgGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCgg                                                                                                      >  2:609842/1‑149 (MQ=255)
                                            aaaaCGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGAc                                                                                                   >  3:545551/1‑149 (MQ=255)
                                                  cAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAata                                                                                             <  4:545551/149‑1 (MQ=255)
                                                   aGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAatat                                                                                            >  4:691461/1‑149 (MQ=255)
                                                   aGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAatat                                                                                            >  3:56935/1‑149 (MQ=255)
                                                       gTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACaaaa                                                                                               <  4:66180/142‑1 (MQ=255)
                                                       gTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACaaaa                                                                                               >  3:66180/1‑142 (MQ=255)
                                                       gTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGt                                                                                        >  1:272062/1‑149 (MQ=255)
                                                            catcCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATc                                                                                   <  4:636613/149‑1 (MQ=255)
                                                             atcCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCg                                                                                  >  2:109543/1‑149 (MQ=255)
                                                               ccATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTg                                                                                <  3:330665/149‑1 (MQ=255)
                                                                cATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACATGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCg                                                                              >  4:789653/1‑149 (MQ=255)
                                                                  tatTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGa                                                                             >  4:319405/1‑149 (MQ=255)
                                                                    ttGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACg                                                                           >  3:239919/1‑149 (MQ=255)
                                                                       aaatagaaataAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGtt                                                                        <  2:173037/149‑1 (MQ=255)
                                                                         atagaaataAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTAc                                                                      >  1:754586/1‑149 (MQ=255)
                                                                         atagaaataAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTAc                                                                      <  2:147934/149‑1 (MQ=255)
                                                                           agaaataAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCGGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTg                                                                    >  4:559955/1‑149 (MQ=255)
                                                                           agaaataAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTg                                                                    <  2:802219/149‑1 (MQ=255)
                                                                            gaaataAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTgg                                                                   <  4:815014/149‑1 (MQ=255)
                                                                            gaaataAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGATTACTgg                                                                   <  3:601858/149‑1 (MQ=255)
                                                                             gaataAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTACGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCATGCGACGGTTACTGGt                                                                  >  4:413380/2‑149 (MQ=255)
                                                                               ataAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTGGTGa                                                                <  1:561954/149‑1 (MQ=255)
                                                                                taAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTGGTGAg                                                               <  3:266366/149‑1 (MQ=255)
                                                                                taAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTGGTGAg                                                               <  4:17363/149‑1 (MQ=255)
                                                                                taAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTGGTGAg                                                               <  4:71401/149‑1 (MQ=255)
                                                                                     ttCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTGGTGAGGTTAt                                                          >  4:349234/1‑149 (MQ=255)
                                                                                     ttCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTGGTGAGGTTAt                                                          >  2:896909/1‑149 (MQ=255)
                                                                                          tCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTGGTGAGGTTATCGCCt                                                     >  3:515949/1‑149 (MQ=255)
                                                                                               cgcTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTGGTGAGGTTATCGCCTAACGt                                                <  3:109461/149‑1 (MQ=255)
                                                                                                gcTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTGGTGAGGTTATCGCCTAACgag                                               <  1:209036/149‑3 (MQ=255)
                                                                                                gcTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTGGTGAGGTTATCGCCTAACGTg                                               >  4:240869/1‑149 (MQ=255)
                                                                                                 cTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTGGTGAGGTTATCGCCTAACGTGa                                              <  2:66893/149‑1 (MQ=255)
                                                                                                  tCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTGGTGAGGTTATCGCCTAACGTGAc                                             >  4:152661/1‑149 (MQ=255)
                                                                                                  tCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTGGTGAGGTTATCGCCTAACGTGAc                                             <  2:104492/149‑1 (MQ=255)
                                                                                                    gTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTGGTGAGGTTATCGCCTAACGTGACTg                                           >  1:504987/1‑149 (MQ=255)
                                                                                                    gTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTGGTGAGGTTATCGCCTAACGTGACTg                                           <  1:17688/149‑1 (MQ=255)
                                                                                                       ggTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTGGTGAGGTTATCGCCTAACGTGACTGgag                                        >  4:156961/1‑149 (MQ=255)
                                                                                                         tgtgACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTGGTGAGGTTATCGCCTAACGTGACTGGAGAt                                      <  1:171981/149‑1 (MQ=255)
                                                                                                              cacGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTGGTGAGGTTATCGCCTAACGTGACTGGAGATACTTc                                 <  3:240869/149‑1 (MQ=255)
                                                                                                              acatgCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTGGTGAGGTTATCGCCTAACGTGACTGGAGATACTTc                                 <  3:789653/145‑1 (MQ=255)
                                                                                                                  gCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAAGTGGCAGGTTTCCAGTACATCCGGGGCAAAATATTTGTTAATCGTGCTACGGTTACTGGTGAGGTTATCGCCTAACGTGACTGGAGATACTTCGTcc                             >  2:697717/1‑149 (MQ=255)
                                                                                                                     tGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTGGTGAGGTTATCGCCTAACGTGACTGGAGATACTTCGTCCGgt                          <  4:91786/149‑1 (MQ=255)
                                                                                                                       gtgtTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTGGTGAGGTTATCGCCTAACGTGACTGGAGATACTTCGTCCGgtgt                        <  1:109543/149‑1 (MQ=255)
                                                                                                                             cacaAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTgg                                                                   >  2:10867/1‑100 (MQ=255)
                                                                                                                             cacaAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTgg                                                                   <  1:10867/100‑1 (MQ=255)
                                                                                                                             cacaAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTGGTg                                                                 <  4:133834/102‑1 (MQ=255)
                                                                                                                             cacaAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTGGTg                                                                 >  3:133834/1‑102 (MQ=255)
                                                                                                                             cacaAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTGGTGAGGTTATCGCCTAACGTGACTGGAGATACTTCGTCCGGTGTGACttt                  >  3:484794/1‑149 (MQ=255)
                                                                                                                                     tgacGGAGCTGGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTGGTGAGGTTATCGCCTAACGTGACTGGAGATACTTCGTCCGGTGTGACTTTTTtggtgg          >  3:700145/4‑149 (MQ=255)
                                                                                                                                     tgacGGAGCTGGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTGGTGAGGTTATCGCCTAACGTGACTGGAGATACTTCGTCCGGTGTGACTGTTTtggtgg          >  1:15710/4‑149 (MQ=255)
                                                                                                                                     tgacGGAGCTGGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTAGTTAATCGTGCGACGGTTACTGGTGAGGTTATCGCCTAACGTGACTGGAGATACTTCGTCCGGTGTGACTTTTTtggtgg          <  4:700145/146‑1 (MQ=255)
                                                                                                                                     aCCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTGGTGAGGTTATCGCCTAACGTGACTGGAGATACTTCGTCCGGTGTGACTTTTTtggtgg          <  3:319405/149‑1 (MQ=255)
                                                                                                                                     aCCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTGGTGAGGTTATCGCCTAACGTGACTGGAGATACTTCGTCCGGTGTGACTTTTTtggtgg          <  1:323704/149‑1 (MQ=255)
                                                                                                                                      cccGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTGGTGAGGTTATCGCCTAACGTGACTGGAGATACTTCGTCCGGTGTGACTTTTTTGGTggg         <  4:432745/149‑1 (MQ=255)
                                                                                                                                      cccGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTGGTGAGGTTATCGCCTAACGTGACTGGAGATACTTCGTCCGGTGTGACTTTTTTGGTggg         <  3:321762/149‑1 (MQ=255)
                                                                                                                                      cccGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTGGTGAGGTTATCGCCTAACGTGACTGGAGATACTTCGTCCGGTGTGACTTTTTTGGTggg         <  1:327322/149‑1 (MQ=255)
                                                                                                                                          gAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTGGTGAGGTTATCGCCTAACGTGACTGGAGATACTTCGTCCGGTGTGACTTTTTTGGTGGGATAc     <  4:141510/149‑1 (MQ=255)
                                                                                                                                             ccAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTGGTGAGGTTATCGCCTAACGTGACTGGAGATACTTCGTCCGGTGTGACTTTTTTGGTGGGATAcgtc  <  1:362133/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                   
GCATCACAATATGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTGTTCCACAAGTTACCCGGAGCCAGCTCCGTCAGCCCCATACTCAATTGGCAGGTTTCCAGTACATCCGGGACAAAATATTTGTTAATCGTGCGACGGTTACTGGTGAGGTTATCGCCTAACGTGACTGGAGATACTTCGTCCGGTGTGACTTTTTTGGTGGGATACGTC  >  NC_000913/2983431‑2983720

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: