Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,810,873 A→G 11.8% H35R (CAC→CGC) ‡ mprA → transcriptional repressor of microcin B17 synthesis and multidrug efflux

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,810,8730AG11.8% 215.9 / 18.0 76H35R (CAC→CGC) ‡mprAtranscriptional repressor of microcin B17 synthesis and multidrug efflux
Reads supporting (aligned to +/- strand):  ref base A (29/38);  new base G (4/5);  total (33/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.20e-01

CTGTCGTTACTATATCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTTCTGAATTAAGTTGTGCTCTTGGATCATCC  >  NC_000913/2810727‑2811021
                                                                                                                                                  |                                                                                                                                                    
cTGTCGTTACTATATCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGcattcaca                                                                                                                                                    >  1:574981/1‑149 (MQ=255)
cTGTCGTTACTATATCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGcatgcaca                                                                                                                                                    >  4:517357/1‑149 (MQ=255)
      ttACTATATCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGTGCATGCaaa                                                                                                                                              <  3:203286/149‑1 (MQ=255)
      ttACTATATCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCaaa                                                                                                                                              <  3:172654/149‑1 (MQ=255)
      ttACTATATCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCaaa                                                                                                                                              <  1:350182/149‑1 (MQ=255)
      ttACTATATCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCaaa                                                                                                                                              <  3:667468/149‑1 (MQ=255)
          tataTCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGTGCATGCaaa                                                                                                                                              <  2:271277/145‑1 (MQ=255)
          tataTCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGTGCATGCaaa                                                                                                                                              >  1:271277/1‑145 (MQ=255)
          tataTCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAaagcaa                                                                                                                                          >  3:246556/1‑149 (MQ=255)
           atatCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAaagcaag                                                                                                                                         >  3:691462/1‑149 (MQ=255)
               cGGCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGTGCATGCaaa                                                                                                                                              <  4:681270/140‑1 (MQ=255)
               cGGCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGTGCATGCaaa                                                                                                                                              >  3:681270/1‑140 (MQ=255)
                  cTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTgg                                                                                                                                  <  4:246556/149‑1 (MQ=255)
                    gAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGgag                                                                                                                                >  3:664677/1‑149 (MQ=255)
                       aTTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGTGCATGCAaaga                                                                                                                                            <  3:18490/134‑2 (MQ=255)
                       aTTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGTGCATGCAaaga                                                                                                                                            >  4:18490/1‑133 (MQ=255)
                       aTTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAAc                                                                                                                             >  3:115405/1‑149 (MQ=255)
                        ttAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAAcc                                                                                                                            <  3:679959/149‑1 (MQ=255)
                         tAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCg                                                                                                                           <  2:95944/149‑1 (MQ=255)
                                 tCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATaaaa                                                                                                                   <  3:10112/149‑1 (MQ=255)
                                 tCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATaaaa                                                                                                                   >  2:670551/1‑149 (MQ=255)
                                   aTACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATg                                                                                                                 >  3:548219/1‑149 (MQ=255)
                                        cAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGaa                                                                                                            <  1:670551/149‑1 (MQ=255)
                                        cAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGaa                                                                                                            <  4:495859/149‑1 (MQ=255)
                                               aTAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAg                                                                                                     <  4:154849/149‑1 (MQ=255)
                                               aTAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAg                                                                                                     >  2:663985/1‑149 (MQ=255)
                                                   ttCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGa                                                                                                 <  1:128380/149‑1 (MQ=255)
                                                   ttCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGa                                                                                                 >  2:628605/1‑149 (MQ=255)
                                                   ttCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGa                                                                                                 >  1:165266/1‑149 (MQ=255)
                                                   ttCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGa                                                                                                 <  2:341767/149‑1 (MQ=255)
                                                    tCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGAt                                                                                                >  4:592408/1‑149 (MQ=255)
                                                        ttACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAAc                                                                                            <  3:454530/149‑1 (MQ=255)
                                                             cccATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGAc                                                                                       <  4:664677/149‑1 (MQ=255)
                                                                aTTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACgtt                                                                                    >  4:82083/1‑149 (MQ=255)
                                                                      cAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTAt                                                                              <  4:115405/149‑1 (MQ=255)
                                                                       aaaTGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATg                                                                             >  1:313048/1‑149 (MQ=255)
                                                                       aaaTGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATg                                                                             <  3:517357/149‑1 (MQ=255)
                                                                             tAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACgtt                                                                                    >  4:438783/1‑136 (MQ=255)
                                                                             tAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACgtt                                                                                    <  3:438783/136‑1 (MQ=255)
                                                                             tAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTAt                                                                              <  1:18869/142‑1 (MQ=255)
                                                                             tAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTAt                                                                              >  2:18869/1‑142 (MQ=255)
                                                                             tAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTg                                                                       >  3:533575/1‑149 (MQ=255)
                                                                              aaaaTTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGa                                                                      <  3:343087/149‑1 (MQ=255)
                                                                               aaaTTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATTCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGAt                                                                     <  2:574981/149‑1 (MQ=255)
                                                                                        gcCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGa                                                            <  2:313048/149‑1 (MQ=255)
                                                                                        gcCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGa                                                            <  2:165266/149‑1 (MQ=255)
                                                                                         ccAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAg                                                           <  3:82083/149‑1 (MQ=255)
                                                                                              cgccACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCa                                                      <  4:691462/149‑1 (MQ=255)
                                                                                                 cACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGa                                                   >  3:626517/1‑149 (MQ=255)
                                                                                                  aCGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGaa                                                  >  1:113216/1‑149 (MQ=255)
                                                                                                    gAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGaaaa                                                >  3:187043/1‑149 (MQ=255)
                                                                                                    gAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGaaaa                                                <  3:592408/149‑1 (MQ=255)
                                                                                                     aaGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAAc                                               <  3:204636/149‑1 (MQ=255)
                                                                                                              cTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTAtt                                      >  1:615326/1‑149 (MQ=255)
                                                                                                                 aTCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAGAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGa                                                            >  4:478429/1‑124 (MQ=255)
                                                                                                                 aTCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGa                                                            <  3:478429/124‑1 (MQ=255)
                                                                                                                    aGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCt                                >  1:250410/1‑149 (MQ=255)
                                                                                                                      gagaTCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAAc                                                                                            <  4:479107/87‑1 (MQ=255)
                                                                                                                      gagaTCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAAc                                                                                            >  3:479107/1‑87 (MQ=255)
                                                                                                                       agaTCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTTCt                             <  4:33179/149‑1 (MQ=255)
                                                                                                                       agaTCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTTCt                             <  4:533575/149‑1 (MQ=255)
                                                                                                                         aTCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTTCTGa                           >  2:658538/1‑149 (MQ=255)
                                                                                                                            cTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTTCTGa                           <  3:121616/146‑1 (MQ=255)
                                                                                                                            cTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTTCTGa                           >  4:121616/1‑146 (MQ=255)
                                                                                                                             ttCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTTCTGAATTa                       >  1:401314/1‑149 (MQ=255)
                                                                                                                              tCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTTCTGAATTaa                      <  1:663985/149‑1 (MQ=255)
                                                                                                                              tCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTTCTGAATTaa                      <  2:401314/149‑1 (MQ=255)
                                                                                                                               cTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTTCTGAATTAAg                     <  1:628605/149‑1 (MQ=255)
                                                                                                                                tGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTTCTGAATTAAGt                    <  4:626517/149‑1 (MQ=255)
                                                                                                                                  aCTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAAGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTATCAGGAAAACCACAGTATTCAGCCTTCTGAATTAGGTtg                  >  4:582684/1‑149 (MQ=255)
                                                                                                                                    tCGTCTTTGTATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTTCTGAATTAAGTtgtg                >  2:189959/1‑149 (MQ=255)
                                                                                                                                    tCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTTCTGAATTAAGTtgtg                >  4:179822/1‑149 (MQ=255)
                                                                                                                                    tCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTTCTGAATTAAGTtgtg                <  2:113216/149‑1 (MQ=255)
                                                                                                                                    tCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTTCTGAATTAAGTtgtg                <  4:548219/149‑1 (MQ=255)
                                                                                                                                    tCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTAGATTACGc                                                                >  1:167300/1‑101 (MQ=255)
                                                                                                                                    tCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTAGATTACGc                                                                <  2:167300/101‑1 (MQ=255)
                                                                                                                                         tttGCATGTGCATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACg                                                                 >  2:85059/1‑95 (MQ=255)
                                                                                                                                         tttGCATGTGCATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACg                                                                 <  1:85059/95‑1 (MQ=255)
                                                                                                                                            gcatgcaCATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGaaaa                                                >  3:578939/1‑109 (MQ=255)
                                                                                                                                            gcatgcaCATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGaaaa                                                <  4:578939/109‑1 (MQ=255)
                                                                                                                                             catgcacatgcaAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTTCTGAATTAAGTTGTGCTCTTGGat       <  2:615326/149‑1 (MQ=255)
                                                                                                                                              atgcacatgcaAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCa                                    >  4:129881/1‑119 (MQ=255)
                                                                                                                                              atgcacatgcaAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCa                                    <  3:129881/119‑1 (MQ=255)
                                                                                                                                              atgcacatgcaAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTTCTGAATTAAGTTGTGCTCTTGGatc      >  2:541336/1‑149 (MQ=255)
                                                                                                                                              atgcacatgcaAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTTCTGAATTAAGTTGTGCTCTTGGatc      >  3:170727/1‑149 (MQ=255)
                                                                                                                                                  acatgcaAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTTCTGAATTAAGTTGTGCTCTTGGATCATcc  >  1:373668/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
CTGTCGTTACTATATCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGAAAACCACAGTATTCAGCCTTCTGAATTAAGTTGTGCTCTTGGATCATCC  >  NC_000913/2810727‑2811021

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