Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,318,095 C→G 100% L594F (TTG→TTC rcsC ← hybrid sensory kinase in two‑component regulatory system with RcsB and YojN

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,318,0950CG98.7% 267.5 / ‑4.7 75L594F (TTG→TTCrcsChybrid sensory kinase in two‑component regulatory system with RcsB and YojN
Reads supporting (aligned to +/- strand):  ref base C (0/1);  major base G (37/36);  minor base A (0/1);  total (37/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CTGGAAGAAGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACTGCTTGCGTACCACCAGCGGTAAATAGTTGGCGGTGATGTGGTTCATCACTTCAC  >  NC_000913/2317951‑2318240
                                                                                                                                                |                                                                                                                                                 
cTGGAAGAAGGGATCAAACAAGCGCACCACTTCTTTCGGCGGTATCCCCACGCCGGTATCGCGAACACCGATAGAGGGGTTATTGCCATCCGCGCCAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACa                                                                                                                                               <  4:439608/149‑1 (MQ=255)
 tGGAAGAAGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAg                                                                                                                                              >  3:390864/1‑149 (MQ=255)
    aagaagGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGtt                                                                                                                                           <  3:116936/149‑1 (MQ=255)
    aagaagGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGtt                                                                                                                                           <  1:359118/149‑1 (MQ=255)
     agaagGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAGCTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTg                                                                                                                                          <  1:356133/149‑1 (MQ=255)
     agaagGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTg                                                                                                                                          >  1:575190/1‑149 (MQ=255)
     agaagGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTg                                                                                                                                          <  3:799491/149‑1 (MQ=255)
       aagGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGga                                                                                                                                        <  1:10164/149‑1 (MQ=255)
       aagGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGga                                                                                                                                        <  1:813740/149‑1 (MQ=255)
        agGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGgag                                                                                                                                       >  4:61748/1‑149 (MQ=255)
        agGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGgag                                                                                                                                       <  1:710502/149‑1 (MQ=255)
         gggATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGTCGTTACTGAACAGGTTGgaga                                                                                                                                      >  3:506972/1‑149 (MQ=255)
           gATCGAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATg                                                                                                                                    >  2:797289/1‑149 (MQ=255)
           gATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATg                                                                                                                                    >  3:839851/1‑149 (MQ=255)
           gATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATg                                                                                                                                    >  1:132042/1‑149 (MQ=255)
           gATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATg                                                                                                                                    >  1:619809/1‑149 (MQ=255)
           gATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATg                                                                                                                                    >  2:462831/1‑149 (MQ=255)
             tCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGAc                                                                                                                                  <  3:3971/149‑1 (MQ=255)
             tCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCAACAGGTTGGAGATGAc                                                                                                                                  <  4:91114/149‑1 (MQ=255)
                aaCAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACctg                                                                                                                               >  4:780259/1‑149 (MQ=255)
                 aCAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACctgc                                                                                                                              >  4:223732/1‑149 (MQ=255)
                       gcACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTaaa                                                                                                                        <  4:805273/149‑1 (MQ=255)
                       gcACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTaaa                                                                                                                        >  2:109665/1‑149 (MQ=255)
                        caccacTTCTGTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAAc                                                                                                                       <  2:63846/149‑1 (MQ=255)
                           cacTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATTGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTaaa                                                                                                                        >  3:493989/1‑145 (MQ=255)
                           cacTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTaaa                                                                                                                        <  4:493989/145‑1 (MQ=255)
                              ttcttTCGCCGGTATCCCCGCGCCGGTATCGCGAACACGGATAGAGAGATAAGCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTAAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCg                                                                                                                 >  4:733549/1‑149 (MQ=255)
                              ttcttTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCg                                                                                                                 >  1:207343/1‑149 (MQ=255)
                                 tttCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGt                                                                                                              >  2:499785/1‑149 (MQ=255)
                                 tttCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGt                                                                                                              >  4:268387/1‑149 (MQ=255)
                                  ttCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTc                                                                                                             <  3:223732/149‑1 (MQ=255)
                                  ttCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTc                                                                                                             <  1:669062/149‑1 (MQ=255)
                                    cGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGc                                                                                                           <  4:376296/149‑1 (MQ=255)
                                              tccACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGAGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGc                                                                                                 >  2:744505/2‑149 (MQ=255)
                                              cccACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGc                                                                                                 <  3:443244/149‑1 (MQ=255)
                                              cccACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTAAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGc                                                                                                 <  4:74773/149‑1 (MQ=255)
                                               ccACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGcc                                                                                                >  4:655152/1‑149 (MQ=255)
                                                    ccGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTgg                                                                                           <  3:780259/149‑1 (MQ=255)
                                                       gTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGcac                                                                                        >  1:100231/1‑149 (MQ=255)
                                                          tCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATc                                                                                     <  1:599477/149‑1 (MQ=255)
                                                          tCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATc                                                                                     <  3:562799/149‑1 (MQ=255)
                                                           cgcgAACACGGATAGACAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATcc                                                                                    <  3:333219/149‑1 (MQ=255)
                                                            gcgAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCg                                                                                   <  1:273236/149‑1 (MQ=255)
                                                                  acGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCaa                                                                             <  1:499785/149‑1 (MQ=255)
                                                                         gagagaTAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCa                                                                      <  1:240858/149‑1 (MQ=255)
                                                                          agagaTAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAg                                                                     >  3:430147/1‑149 (MQ=255)
                                                                           gagaTAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGt                                                                    <  2:100231/149‑1 (MQ=255)
                                                                             gaTAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTAc                                                                  >  2:50288/1‑149 (MQ=255)
                                                                               tAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACaa                                                                <  1:315562/149‑1 (MQ=255)
                                                                               tAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACaa                                                                <  1:441467/149‑1 (MQ=255)
                                                                                 aTCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGc                                                              >  3:838217/1‑149 (MQ=255)
                                                                                      cATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAAc                                                         >  3:178219/1‑149 (MQ=255)
                                                                                        tCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACTg                                                       >  3:102412/1‑149 (MQ=255)
                                                                                            cgcgAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACTGCTTg                                                   <  3:773098/149‑1 (MQ=255)
                                                                                                   aTGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACTGCTTGCGTacca                                            <  3:217286/149‑1 (MQ=255)
                                                                                                    tGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACTGCTTGCGTaccac                                           >  1:159035/1‑149 (MQ=255)
                                                                                                     gCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACTGCTTGCGTaccacc                                          >  2:96435/1‑149 (MQ=255)
                                                                                                            tataCAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCaa                                                                             >  1:687074/1‑107 (MQ=255)
                                                                                                            tataCAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCaa                                                                             <  2:687074/107‑1 (MQ=255)
                                                                                                            tataCAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACTGCTTGCGTACCACCAGCGGTa                                   <  1:797289/149‑1 (MQ=255)
                                                                                                            tataCAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACTGCTTGCGTACCACCAGCGGTa                                   >  4:749634/1‑149 (MQ=255)
                                                                                                            tataCAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACTGCTTGCGTACCACCAGCGGTa                                   >  2:782698/1‑149 (MQ=255)
                                                                                                              taCAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACTGCTTGCGTACCACCAGTGGTaaa                                 <  3:800491/149‑1 (MQ=255)
                                                                                                                 aGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACTGCTTGCGTACCACCAGCGGTAAATAg                              >  4:181877/1‑149 (MQ=255)
                                                                                                                    cGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACCGCTTGCGTACCACCAGCGGTAAATAGTTg                           <  2:575190/149‑1 (MQ=255)
                                                                                                                     ggTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACTGCTTGCGTACCACCAGCGGTAAATAGTTgg                          <  1:70981/149‑1 (MQ=255)
                                                                                                                      gTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGCAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACTGCTTGCGTACCACCAGCGGTAAATAGTTGGc                         <  2:531166/149‑1 (MQ=255)
                                                                                                                            gTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACTGCTTGCGTACCACCAGCGGTAAATAGTTGGCGGTGAt                   >  3:174118/1‑149 (MQ=255)
                                                                                                                                 ttttATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACTGCTTGCGTACCACCAGCGGTAAATAGTTGGCGGTGATGTGGt              >  1:569/1‑149 (MQ=255)
                                                                                                                                 ttttATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACTGCTTGCGTACCACCAGCGGTAAATAGTTGGCGGTGATGTGGt              >  4:45265/1‑149 (MQ=255)
                                                                                                                                 ttttATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACTGCTTGCGTACCACCAGCGGTAAATAGTTGGCGGTGATGTGGt              >  1:854497/1‑149 (MQ=255)
                                                                                                                                   ttATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACTGCTTGCGTACCACCAGCGGTAAATAGTTGGCGGTGATGTGGTtc            <  2:159035/149‑1 (MQ=255)
                                                                                                                                           ttACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACTGCTTGCGTACCACCAGCGGTAAATAGTTGGCGGTGATGTGGTTCATCACTTc    >  2:244681/1‑149 (MQ=255)
                                                                                                                                            tACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACTGCTTGCGTACCACCAGCGGTAAATAGTTGGCGGTGATGTGGTTCATCACTTCa   >  3:29662/1‑149 (MQ=255)
                                                                                                                                             aCTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACTGCTTGCGTACCACCAGCGGTAAATAGTTGGCGGTGATGTGGTTCATCACTTCAc  <  4:331041/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
CTGGAAGAAGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACTGCTTGCGTACCACCAGCGGTAAATAGTTGGCGGTGATGTGGTTCATCACTTCAC  >  NC_000913/2317951‑2318240

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