Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,915,982 A→G 5.7% N309S (AAT→AGT)  barA → hybrid sensory histidine kinase, in two‑component regulatory system with UvrY

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,915,9820AG 5.7% 293.1 / 5.8 89N309S (AAT→AGT) barAhybrid sensory histidine kinase, in two‑component regulatory system with UvrY
Reads supporting (aligned to +/- strand):  ref base A (43/40);  new base G (3/2);  total (47/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.47e-01

AGGCGACGTCCGATCTGCGTGAAACGCTGGAGCAGATGGAAATTCAGAACGTTGAGTTAGATCTGGCGAAAAAGCGCGCCCAGGAAGCGGCGCGTATTAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAACACCAACGCAGCGCGATCACCTGAATACGATTGAACGTTCGGCAAATAATTTGCTGGCAATTATTAATGATGTTCTCGACTTCTCGAAACTGGAAGCAGGT  >  NC_000913/2915835‑2916127
                                                                                                                                                   |                                                                                                                                                 
aGGCGACGTCCGATCTGCGTGAAACGCTGGAGCAGATGGAAATTCAGAACGTTGAGTTAGATCTGGCGAAAAAGCGCGCCCAGGAAGCGGCGCGTATTAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAAt                                                                                                                                                  <  2:22168/149‑1 (MQ=255)
   cGACGTCCGATCTGCGTGAAACGCTGGAGCAGATGGAAATTCAGAACGTTGAGTTAGATCTGGCGAAAAAGCGCGCCCAGGAAGCGGCGCGTATTAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGgt                                                                                                                                               <  3:148274/149‑1 (MQ=255)
    gACGTCCGATCTGCGTGAAACGCTGGAGCAGATGGAAATTCAGAACGTTGAGTTAGATCTGGCGAAAAAGCGCGCCCAGGAAGCGGCGCGTATTAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGgtg                                                                                                                                              <  1:869185/149‑1 (MQ=255)
         ccGATCTGCGTGAAACGCTGGAGCAGATGGAAATTCAGAACGTTGAGTTAGATCTGGCGAAAAAGCGCGCCCAGGAAGCGGCGCGTATTAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTAtt                                                                                                                                         <  3:110259/149‑1 (MQ=255)
           gATCTGCGTGAAACGCTGGAGCAGATGGAAATTCAGAACGTTGAGTTAGATCTGGCGAAAAAGCGCGCCCAGGAAGCGGCGCGTATTAAATGCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTgg                                                                                                                                       <  3:400341/149‑1 (MQ=255)
           gATCTGCGTGAAACGCTGGAGCAGATGGAAATTCAGAACGTTGAGTTAGATCTGGCGAAAAAGCGCGCCCAGGAAGCGGCGCGTATTAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTa                                                                                                                                           >  2:568371/1‑145 (MQ=255)
           gATCTGCGTGAAACGCTGGAGCAGATGGAAATTCAGAACGTTGAGTTAGATCTGGCGAAAAAGCGCGCCCAGGAAGCGGCGCGTATTAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTgg                                                                                                                                       <  1:458126/149‑1 (MQ=255)
           gATCTGCGTGAAACGCTGGAGCAGATGGAAATTCAGAACGTTGAGTTAGATCTGGCGAAAAAGCGCGCCCAGGAAGCGGCGCGTATTAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTACGTACACCACTGAATGGTGTTa                                                                                                                                           <  1:568371/145‑1 (MQ=255)
            aTCTGCGTGAAACGCTGGAGCAGATGGAAATTCAGAACGTTGAGTTAGATCTGGCGAAAAAGCGCGCCCAGGAAGCGGCGCGTATTAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGgt                                                                                                                                               >  1:667578/1‑140 (MQ=255)
            aTCTGCGTGAAACGCTGGAGCAGATGGAAATTCAGAACGTTGAGTTAGATCTGGCGAAAAAGCGCGCCCAGGAAGCGGCGCGTATTAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGgt                                                                                                                                               <  2:667578/140‑1 (MQ=255)
            aTCTGCGTGAAACGCTGGAGCAGATGGAAATTCAGAACGTTGAGTTAGATCTGGCGAAAAAGCGCGCCCAGGAAGCGGCGCGTATTAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGc                                                                                                                                      <  4:48526/149‑1 (MQ=255)
            aTCTGCGTGAAACGCTGGAGCAGATGGAAATTCAGAACGTTGAGTTAGATCTGGCGAAAAAGCGCGCCCAGGAAGCGGCGCGTATTAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGc                                                                                                                                      >  2:886191/1‑149 (MQ=255)
              cTGCGTGAAACGCTGGAGCAGATGGAAATTCAGAACGTTGAGTTAGATCTGGCGAAAAAGCGCGCCCAGGAAGCGGCGCGTATTAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCtt                                                                                                                                    <  1:465817/149‑1 (MQ=255)
                   tGAAACGCTGGAGCAGATGGAAATTCAGAACGTTGAGTTAGATCTGGCGAAAAAGCGCGCCCAGGAAGCGGCGCGTATTAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTAccc                                                                                                                               >  1:112286/1‑149 (MQ=255)
                   tGAAACGCTGGAGCAGATGGAAATTCAGAACGTTGAGTTAGATCTGGCGAAAAAGCGCGCCCAGGAAGCGGCGCGTATTAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTAccc                                                                                                                               >  1:811765/1‑149 (MQ=255)
                     aaaCGCTGGAGCAGATGGAAATTCAGAACGTTGAGTTAGATCTGGCGAAAAAGCGCGCCCAGGAAGCGGCGCGTATTAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGc                                                                                                                                      >  1:292975/1‑140 (MQ=255)
                     aaaCGCTGGAGCAGATGGAAATTCAGAACGTTGAGTTAGATCTGGCGAAAAAGCGCGCCCAGGAAGCGGCGCGTATTAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGc                                                                                                                                      <  2:292975/140‑1 (MQ=255)
                     aaaCGCTGGAGCAGATGGAAATTCAGAACGTTGAGTTAGATCTGGCGAAAAAGCGCGCCCAGGAAGCGGCGCGTATTAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGc                                                                                                                             <  1:213007/149‑1 (MQ=255)
                     aaaCGCTGGAGCAGATGGAAATTCAGAACGCTGAGTTAGATCTGGCGAAAAAGCGCGCCCAGGAAGCGGCGCGTATTAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGc                                                                                                                             >  2:627654/1‑149 (MQ=255)
                      aaCGCTGGAGCAGATGGAAATTCAGAACGTTGAGTTAGATCTGGCGAAAAAGCGCGCCCAGGAAGCGGCGCGTATTAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGcc                                                                                                                            >  2:57421/1‑149 (MQ=255)
                      aaCGCTGGAGCAGATGGAAATTCAGAACGTTGAGTTAGATCTGGCGAAAAAGCGCGCCCAGGAAGCGGCGCGTATTAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGcc                                                                                                                            >  4:370767/1‑149 (MQ=255)
                                  gATGGAAATTCAGAACGTTGAGTTAGATCTGGCGAAAAAGCGCGCCCAGGAAGCGGCGCGTATTAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGaaaa                                                                                                                <  1:224613/149‑1 (MQ=255)
                                    tGGAAATTCAGAACGTTGAGTTAGATCTGGCGAAAAAGCGCGCCCAGGAAGCGGCGCGTATTAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACa                                                                                                              >  4:355561/1‑149 (MQ=255)
                                     ggAAATTCAGAACGTTGAGTTAGATCTGGCGAAAAAGCGCGCCCAGGAAGCGGCGCGTATTAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAg                                                                                                             >  3:8362/1‑149 (MQ=255)
                                       aaaTTCAGGACGTTGAGTTAGATCTGGGGAAAAAGCGCGCCCAGGAAGCGGCGCGTATTAAATCCGAGTTTCTGGCAAATATGTCACCCGAGCTGCGTACACCACTGAATGGTGGTATTGGCTTTACCCGCCTGACGCTGAAAACAGaa                                                                                                           <  1:685278/149‑1 (MQ=255)
                                       aaaTTCAGAACGTTGAGTTAGATCTGGCGAAAAAGCGCGCCCAGGAAGCGGCGCGTATTAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGaa                                                                                                           <  3:863556/149‑1 (MQ=255)
                                        aaTTCAGAACGTTGAGTTAGATCTGGCGAAAAAGCGCGCCCAGGAAGCGGCGCGTATTAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCCTTACCCGCCTGACGCTGAAAACAGAAt                                                                                                          <  3:724903/149‑1 (MQ=255)
                                         aTTCAGAACGTTGAGTTAGATCTGGCGAAAAAGCGCGCCCAGGAAGCGGCGCGTATTAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAAtt                                                                                                         <  2:464565/149‑1 (MQ=255)
                                          ttCAGAACGTTGAGTTAGATCTGGCGAAAAAGCGCGCCCAGGAAGCGGCGCGTATTAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTgg                                                                                                                                       <  1:932907/118‑1 (MQ=255)
                                          ttCAGAACGTTGAGTTAGATCTGGCGAAAAAGCGCGCCCAGGAAGCGGCGCGTATTAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTgg                                                                                                                                       >  2:932907/1‑118 (MQ=255)
                                          ttCAGAACGTTGAGTTAGATCTGGCGAAAAAGCGCGCCCAGGAAGCGGCGCGTATTAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTa                                                                                                        >  2:405744/1‑149 (MQ=255)
                                              gAACGTTGAGTTAGATCTGGCGAAAAAGCGCGCCCAGGAAGCGGCGCGTATTAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAacac                                                                                                    >  3:173334/1‑149 (MQ=255)
                                              gAACGTTGAGTTAGATCTGGCGAAAAAGCGCGCCCAGGAAGCGGCGCGTATTAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAacac                                                                                                    >  3:322436/1‑149 (MQ=255)
                                                aCGTTGAGTTAGATCTGGCGAAAAAGCGCGCCCAGGAAGCGGCGCGTATTAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAACACCa                                                                                                  >  1:751143/1‑149 (MQ=255)
                                                   ttGAGTTAGATCTGGCGAAAAAGCGCGCCCAGGAAGCGGCGCGTATTAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAACACCAACg                                                                                               >  3:896015/1‑149 (MQ=255)
                                                      aGTTAGATCTGGCGAAAAAGCGCGCCCAGGAAGCGGCGCGTATTAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAACACCAACGCAg                                                                                            >  1:445514/1‑149 (MQ=255)
                                                           gATCTGGCGAAAAAGCGCGCCCAGGAAGCGGCGCGTATTAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAACACCAACGCAGCGCGa                                                                                       <  2:112286/149‑1 (MQ=255)
                                                                ggCGAAAAAGCACGCCCAGGAAGCGGCGCGTATTAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAACACCAACGCAGCGCGATCAcc                                                                                  >  1:490745/1‑149 (MQ=255)
                                                                   gAAAAAGCGCGCCCAGGAAGCGGCGCGTATTAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAACACCAACGCAGCGCGATCACCTGa                                                                               >  3:540142/1‑149 (MQ=255)
                                                                   caaaaaGCGCGGCCAGGAAGCGGCGCGTATTAAATTCGAGGTTCTGGCAAATATGTCACACGATACGCGTACACCACTGATTGGTGTTATTGGCTTTACCCGCCGGACGCTGAAGACAGAATTAACACCAACGCAGCGCGATCACCTGa                                                                               >  2:84470/2‑149 (MQ=255)
                                                                    aaaaaGCGCGCCCAGGAAGCGGCGCGTATTAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAACACCAACGCAGCGCGATCACCTGaa                                                                              >  2:687323/1‑149 (MQ=255)
                                                                         gcgcgcCCAGGAAGCGGCGCGTATTAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAACACCAACGCAGCGCGATCACCTGAATACGa                                                                         <  3:249947/149‑1 (MQ=255)
                                                                           gcgcCCAGGAAGCGGCGCGTATTAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAACACCAACGCAGCGCGATCACCTGAATACGAtt                                                                       <  1:356196/149‑1 (MQ=255)
                                                                           gcgcCCAGGAAGCGGCGCGTATTAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAACACCAACGCAGCGCGATCACCTGAATACGAtt                                                                       <  1:311884/149‑1 (MQ=255)
                                                                             gcCCAGGAAGCGGCGCGTATTAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAACACCAACGCAGCGCGATCACCTGAATACGATTGa                                                                     <  1:250247/149‑1 (MQ=255)
                                                                                  ggAAGCGGCGCGTATTAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAACACCAACGCAGCGCGATCACCTGAATACGATTGAACGtt                                                                <  4:67781/149‑1 (MQ=255)
                                                                                   gAAGCGGCGCGTATTAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAACACCAACGCAGCGCGATCACCTGAATACGATTGAACGTTc                                                               <  3:233545/149‑1 (MQ=255)
                                                                                    aaGCGGCGCGTATTAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAACACCAACGCAGCGCGATCACCTGAATACGATTGAACGTTCg                                                              <  3:38477/149‑1 (MQ=255)
                                                                                         gcgcGTATTAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAACACCAACGCAGCGCGATCACCTGAATACGATTGAACGTTCGGCaaa                                                         >  3:225264/1‑149 (MQ=255)
                                                                                              tATTAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAACACCAACGCAGCGCGATCACCTGAATACGATTGAACGTTCGGCAAATAAtt                                                    <  3:355561/149‑1 (MQ=255)
                                                                                              tATTAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAACACCAACGCAGCGCGATCACCTGAATACGATTGAACGTTCGGCAAATAAtt                                                    >  4:246761/1‑149 (MQ=255)
                                                                                               aTTAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAACACCAACGCAGCGCGATCACCTGAATACGATTGAACGTTCGGCAAATAAttt                                                   >  1:615208/1‑149 (MQ=255)
                                                                                                 tAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAACACCAACGCAGCGCGATCACCTGAATGCGATTGAACGTTCGGCAAATAATTTGc                                                 >  1:265406/1‑149 (MQ=255)
                                                                                                  aaaTCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAACACCAACGCAGCGCGATCACCTGAATACGATTGAACGTTCGGCAAATAATTTGCt                                                <  4:482240/149‑1 (MQ=255)
                                                                                                     tCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAACACCAACGCAGCGCGATCACCTGAATACGATTGAACGTTCGGCAAATAATTTGCTGGc                                             >  2:20840/1‑149 (MQ=255)
                                                                                                     tCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAACACCAACGCAGCGCGATCACCTGAATACGATTGAACGTTCGGCAAATAATTTGCTGGc                                             >  4:618919/1‑149 (MQ=255)
                                                                                                     tCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAACACCAACGCAGCGCGATCACCTGAATACGATTGAACGTTCGGCAAATAATTTGCTGGc                                             >  2:589246/1‑149 (MQ=255)
                                                                                                        gAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAACACCAACGCAGCGCGATCACCTGAATACGATTGAACGTTCGGCAAATAATTTGCTGGCAat                                          >  4:361064/1‑149 (MQ=255)
                                                                                                           tttCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAACACCAACGCAGCGCGATCACCTGAATACGATTGAACGTTCGGCAAATAATTTGCTGGCAattat                                       <  1:57421/149‑1 (MQ=255)
                                                                                                              cTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAACACCAACGCAGCGCGATCACCTGAATACGATTGAACGTTCGGCAAATAATTTGCTGGCAATTATTaa                                    <  2:811765/149‑1 (MQ=255)
                                                                                                               tGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAACACCAACGCAGCGCGATCACCTGAATACGATTGAACGTTCGGCAAATAATTTGCTGGCAATTATTAat                                   <  2:751143/149‑1 (MQ=255)
                                                                                                               tGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAACACCAACGCAGCGCGATCACCTGAATACGATTGAACGTTCGGCAAATAATTTGCTGGCAATTATTAat                                   >  1:242254/1‑149 (MQ=255)
                                                                                                                 gCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAACACCAACGCAGCGCGATCACCTGAATACGATTGAACGTTCGGCAAATAATTTGCTGGCAATTATTAatga                                 >  1:373563/1‑149 (MQ=255)
                                                                                                                   aaaTATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAACACCAACGCAGCGCGATCACCTGAATACGATTGAACGTTCGGCAAATAATTTGCTGGCAATTATTAatgatg                               >  4:305116/1‑149 (MQ=255)
                                                                                                                   aaaTATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAACACCAACGCAGCGCGATCACCTGAATACGATTGAACGTTCGGCAAATAATTTGCTGGCAATTATTAatgatg                               >  3:531034/1‑149 (MQ=255)
                                                                                                                   aaaTATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAACACCAACGCAGCGCGATCACCTGAATACGATTGAACGTTCGGCAAATAATTTGCTGGCAATTATTAatgatg                               >  3:575401/1‑149 (MQ=255)
                                                                                                                    aaTATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAACACCAACGCAGCGCGATCACCTGAATACGATTGAACGTTCGGCAAATAATTTGCTGGCAATTATTAATGATGt                              >  4:691722/1‑149 (MQ=255)
                                                                                                                      tatGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAACACCAACGCAGCGCGATCACCTGAATACGATTGAACGTTCGGCAAATAATTTGCTGGCAATTATTaa                                    >  3:126651/1‑141 (MQ=255)
                                                                                                                      tatGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAACACCAACGCAGCGCGATCACCTGAATACGATTGAACGTTCGGCAAATAATTTGCTGGCAATTATTaa                                    <  4:126651/141‑1 (MQ=255)
                                                                                                                      tatGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAACACCAACGCAGCGCGATCACCTGAATACGATTGAACGTTCGGCAAATAATTTGCTGGCAATTATTAATGATGTtc                            <  3:354402/149‑1 (MQ=255)
                                                                                                                        tGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAACACCAACGCAGCGCGATCACCTGAATACGATTGAACGTTCGGCAAATAATTTGCTGGCAATTATTAATGATGTtctc                          <  4:322436/149‑1 (MQ=255)
                                                                                                                        tGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAACACCAACGCAGCGCGATCACCTGAATACGATTGAACGTTCGGCAAATAATTTGCTGGCAATTATTAATGATGTtctc                          <  2:615208/149‑1 (MQ=255)
                                                                                                                         gTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAACACCAACGCAGCGCGATCACCTGAATACGATTGAACGTTCGGCAAATAATTTGCTGGCAATTATTAATGATGTTCTCg                         >  2:713648/1‑149 (MQ=255)
                                                                                                                         gTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAACACCAACGCAGCGCGATCACCTGAATACGATTGAACGTTCGGCAAATAATTTGCTGGCAATTATTAATGATGTTCTCg                         >  1:791296/1‑149 (MQ=255)
                                                                                                                          tCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAACACCAACGCAGCGCGATCACCTGAATACGATTGAACGTTCGGCAAATAATTTGCTGGCAATTATTAATGATGTTCTCGa                        <  4:188390/149‑1 (MQ=255)
                                                                                                                                    caataacacCATTCAGTGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAACACCAACGCAGCGCGATCACCTGAATACGATTGAACGTTCGGCAAATAATTTGCTGGCAATTATTAATGATGTtctc                          >  3:669539/6‑137 (MQ=255)
                                                                                                                                    caataacacCATTCAGTGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAACACCAACGCAGCGCGATCACCTGAATACGATTGAACGTTCGGCAAATAATTTGCTGGCAATTATTAATGATGTtctc                          <  4:669539/132‑1 (MQ=255)
                                                                                                                                     aataacacCATTCAGTGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAACACCAACGCAGCGCGATCACCTGAATACGAtt                                                                       <  1:478808/87‑1 (MQ=255)
                                                                                                                                     aataacacCATTCAGTGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAACACCAACGCAGCGCGATCACCTGAATACGAtt                                                                       >  2:478808/5‑91 (MQ=255)
                                                                                                                                     aataacacCATTCAGTGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAACACCAACGCAGCGCGATCACCTGAATACGATTGAACGTTCGGCAAATAATTTGCTGGCAATTATTAATGATGTTCTCGACTTCTCGAAAc             >  2:290294/5‑149 (MQ=255)
                                                                                                                                        tACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAACACCAACGCAGCGCGATCACCTGAATACGATTGAACGTTCGGCAAATAATTTGCTGGCAATTATTAATGATGTTCTCGACTTCTCGAAACTgg          <  4:225264/149‑1 (MQ=255)
                                                                                                                                          caccacTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAACACCAACGCAGCGCGATCACCTGAATACGATTGAACGTTCGGCAAATAATTTGCTGGCAATTATTAATGATGTTCTCGACTTCTCGAAACTGGaa        >  1:263686/1‑149 (MQ=255)
                                                                                                                                           accacTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAACACCAACGCAGCGCGATCACCTGAATACGATTGAACGTTCGGCAAATAATTTGCTGGCAATTATTAATGATGTTCTCGACTTCTCGAAACTGGAAg       >  4:563996/1‑149 (MQ=255)
                                                                                                                                             cacTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAACACCAACGCAGCGCGATCACCTGAATACGATTGAACGTTCGGCAAATAATTTGCTGGCAATTATTAATGATGTTCTCGACTTCTCGAAACTGGAAGCa     >  4:541319/1‑149 (MQ=255)
                                                                                                                                             cacTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAACACCAACGCAGCGCGATCACCTGAATACGATTGAACGTTCGGCAAATAATTTGCTGGCAATTATTAATGATGTTCTCGACTTCTCGAAACTGGAAGCa     >  1:466073/1‑149 (MQ=255)
                                                                                                                                              acTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAACACCAACGCAGCGCGATCACCTGAATACGATTGAACGTTCGGCAAATAATTTGCTGGCAATTATTAATGATGTTCTCGACTTCTCGAAACTGGaa        >  4:457507/1‑145 (MQ=255)
                                                                                                                                              acTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAACACCAACGCAGCGCGATCACCTGAATACGATTGAACGTTCGGCAAATAATTTGCTGGCAATTATTAATGATGTTCTCGACTTCTCGAAACTGGaa        <  3:457507/145‑1 (MQ=255)
                                                                                                                                              acTGAATGGTGATATTGGCATCACACGCCTGACGCTGAAAAAAGAATTAACACCAACGCAGCGCCATCACCTAAATACGATTGCACGCTCGGCAAATAATTTGCTGGCAATTATTAATGATGTTCTCGACTTCTCGAAACTGGAAGCAg    <  2:788792/149‑1 (MQ=255)
                                                                                                                                                tGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAACACCAACGCAGCGCGATCACCTGAATACGATTGAACGTTCGGCAAATAATTTGCTGGCAATTATTAATGATGTTCTCGACTTCTCGAAACTGGAAGCAGGt  <  2:801386/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
AGGCGACGTCCGATCTGCGTGAAACGCTGGAGCAGATGGAAATTCAGAACGTTGAGTTAGATCTGGCGAAAAAGCGCGCCCAGGAAGCGGCGCGTATTAAATCCGAGTTTCTGGCAAATATGTCACACGAGCTGCGTACACCACTGAATGGTGTTATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAACACCAACGCAGCGCGATCACCTGAATACGATTGAACGTTCGGCAAATAATTTGCTGGCAATTATTAATGATGTTCTCGACTTCTCGAAACTGGAAGCAGGT  >  NC_000913/2915835‑2916127

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