Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 524,394 A→C 6.4% A378A (GCA→GCC)  rhsD → Rhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913524,3940AC 6.4% 261.3 / 8.9 78A378A (GCA→GCC) rhsDRhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor
Reads supporting (aligned to +/- strand):  ref base A (29/44);  new base C (2/3);  total (31/47)
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.39e-01

TATATGACCGCAGCAATACGCAGGTGCGCGCTTTCACGTATGACGCGCAGCACCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCGGGGCCGGGCTGAAACGGGTGGTGAAAAAAGAACTGGCGGACGGCAGCGTCACGCGCAGC  >  NC_000913/524246‑524541
                                                                                                                                                    |                                                                                                                                                   
tataTGACCGCAGCAATACGCAGGTGCGCGCTTTCACGTATGACGCGCAGCACCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCa                                                                                                                                                     >  4:199863/1‑149 (MQ=255)
tataTGACCGCAGCAATACGCAGGTGCGCGCTTTCACGTATGACGCGCAGCACCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCa                                                                                                                                                     <  4:291184/149‑1 (MQ=255)
  tatGACCGCAGCAATACGCAGGTGCGCGCTTTCACGTATGACGCGCAGCATCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAgg                                                                                                                                                   <  4:356131/149‑1 (MQ=255)
  tatGACCGCAGCAATACGCAGGTGCGCGCTTTCACGTATGACGCGCAGCACCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAgg                                                                                                                                                   >  1:131383/1‑149 (MQ=255)
   atGACCGCAGCAATACGCAGGTGCGCGCTTTCACGTATGACGCGCAGCACCCGGGCCGGATGGTGGCGCACTGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAggg                                                                                                                                                  <  3:463747/149‑1 (MQ=255)
      aCCGCAGCAATACGCAGGTGCGCGCTTTCACGTATGACGCGCAGCACCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTa                                                                                                                                               <  3:658120/149‑1 (MQ=255)
       ccGCAGCAATACGCAGGTGCGCGCTTTCACGTATGACGCGCAGCACCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTaa                                                                                                                                              <  1:56312/149‑1 (MQ=255)
        cGCAGCAATACGCAGGTGCGCGCTTTCACGTATGACGCGCAGCACCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCTGCCGGGTTCAg                                                                                                                                             <  4:255773/149‑1 (MQ=255)
           agcaATACGCAGGTGCGCGCTTTCACGTATGACGCGCAGCACCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTa                                                                                                                                          >  4:491312/1‑149 (MQ=255)
            gcaATACGCAGGTGCGCGCTTTCACGTATGACGCGCAGCACCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTAc                                                                                                                                         <  2:394068/149‑1 (MQ=255)
                     aGGTGCGCGCTTTCACGTATGACGCGCAGCACCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCtt                                                                                                                                <  3:34256/149‑1 (MQ=255)
                       gTGCGCGCTTTCACGTATGACGCGCAGCACCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTa                                                                                                                              <  2:597865/149‑1 (MQ=255)
                        tGCGCGCTTTCACGTATGACGCGCAGCACCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTAt                                                                                                                             <  1:449174/149‑1 (MQ=255)
                          cgcgcTTTCACGTATGACGCGCAGCACCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGa                                                                                                                           <  4:673098/149‑1 (MQ=255)
                           gcgcTTTCACGTATGACGCGCAGCACCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTAt                                                                                                                             >  2:61097/1‑146 (MQ=255)
                           gcgcTTTCACGTATGACGCGCAGCACCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTAt                                                                                                                             <  1:61097/146‑1 (MQ=255)
                           gcgcTTTCACGTATGACGCGCAGCACCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAg                                                                                                                          >  3:646975/1‑149 (MQ=255)
                             gcTTTCACGTATGACGCGCAGCACCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCa                                                                                                                        >  4:299453/1‑149 (MQ=255)
                             gcTTTCACGTATGACGCGCAGCACCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCa                                                                                                                        >  2:248527/1‑149 (MQ=255)
                             gcTTTCACGTATGACGCGCAGCACCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCa                                                                                                                        <  3:491312/149‑1 (MQ=255)
                                   aCGTATGACGCGCAGCACCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCg                                                                                                                  >  1:398459/1‑149 (MQ=255)
                                   aCGTATGACGCGCAGCACCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCg                                                                                                                  <  2:390066/149‑1 (MQ=255)
                                     gTATGACGCGCAGCACCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCa                                                                                                                >  1:168733/1‑149 (MQ=255)
                                        tGACGCGCAGCACCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCAtca                                                                                                             <  4:379031/149‑1 (MQ=255)
                                        tGACGCGCAGCACCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCAtca                                                                                                             >  2:275067/1‑149 (MQ=255)
                                         gACGCGCAGCACCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAggg                                                                                                                                                  >  1:581914/1‑111 (MQ=255)
                                         gACGCGCAGCACCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAggg                                                                                                                                                  <  2:581914/111‑1 (MQ=255)
                                            gcgcAGCACCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCAtcaccgt                                                                                                         <  4:122374/149‑1 (MQ=255)
                                            gcgcAGCACCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCAtcaccgt                                                                                                         >  4:293841/1‑149 (MQ=255)
                                                    cccGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACa                                                                                                 <  3:27816/149‑1 (MQ=255)
                                                      cGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGc                                                                                               >  3:331618/1‑149 (MQ=255)
                                                       gggCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGcc                                                                                              >  4:619398/1‑149 (MQ=255)
                                                           cGGATGGTGGTGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGaa                                                                                          >  2:229594/1‑149 (MQ=255)
                                                           cGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGCCGGTACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGaa                                                                                          >  1:458159/1‑149 (MQ=255)
                                                            ggATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAAc                                                                                         >  4:594609/1‑149 (MQ=255)
                                                                ggtggCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGc                                                                                     <  2:398459/149‑1 (MQ=255)
                                                                  tggCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGt                                                                                   <  1:482461/149‑1 (MQ=255)
                                                                  tggCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGt                                                                                   <  3:528415/149‑1 (MQ=255)
                                                                  tggCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGt                                                                                   <  4:331618/149‑1 (MQ=255)
                                                                     cgcACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAg                                                                                <  1:670336/149‑1 (MQ=255)
                                                                      gcACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAgg                                                                               >  3:48474/1‑149 (MQ=255)
                                                                         ccGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGtgc                                                                            <  4:332210/149‑1 (MQ=255)
                                                                         ccGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGtgc                                                                            >  1:534196/1‑149 (MQ=255)
                                                                          cGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCAtcacc                                                                                                           >  1:219368/1‑117 (MQ=255)
                                                                          cGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCAtcacc                                                                                                           <  2:219368/117‑1 (MQ=255)
                                                                           gTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGtgctg                                                                          >  3:468721/1‑149 (MQ=255)
                                                                            ttACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGtgctgc                                                                         <  3:299453/149‑1 (MQ=255)
                                                                            ttACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGtgctgc                                                                         <  1:248527/149‑1 (MQ=255)
                                                                                   ggAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGa                                                                  >  4:633500/1‑149 (MQ=255)
                                                                                    gAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGaa                                                                 <  3:460721/149‑1 (MQ=255)
                                                                                         ccGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCgg                                                            <  4:149809/149‑1 (MQ=255)
                                                                                         ccGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCgg                                                            <  3:357090/149‑1 (MQ=255)
                                                                                         ccGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCgg                                                            <  2:139894/149‑1 (MQ=255)
                                                                                          cGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCggg                                                           >  1:129070/1‑149 (MQ=255)
                                                                                               aTGCGCTACCGCTACGACGTTACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGTTGCATACAGAAGGCGGGGCCgg                                                      <  1:425454/149‑1 (MQ=255)
                                                                                               aTGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCGGGGCCgg                                                      >  3:117021/1‑149 (MQ=255)
                                                                                                  cgctaccgctacGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGGGGGGCCGGGCt                                                   <  4:468721/149‑1 (MQ=255)
                                                                                                  cgctaccgctacGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCGGGGCCGGGCt                                                   <  3:648184/149‑1 (MQ=255)
                                                                                                    ctaccgctacGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCGGGGCCGGGCTGa                                                 >  1:45901/1‑149 (MQ=255)
                                                                                                       ccgctacGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCGGGGCCGGGCTGAAAc                                              <  4:68245/149‑1 (MQ=255)
                                                                                                           tacGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCGGGGCCGGGCTGAAACGggt                                          <  1:479654/149‑1 (MQ=255)
                                                                                                           tacGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCGGGGCCGGGCTGAAACGggt                                          <  1:229594/149‑1 (MQ=255)
                                                                                                            acgacgATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCGGGGCCGGGCTGGAACGggtg                                         <  4:571881/149‑1 (MQ=255)
                                                                                                              gacgaTACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCGGGGCCGGGCTGAAACGggtggt                                       >  4:434872/1‑149 (MQ=255)
                                                                                                              gacgaTACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCGGGGCCGGGCTGAAACGggtggt                                       <  3:160249/149‑1 (MQ=255)
                                                                                                                 gaTACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCGGGGCCGGGCTGAAACGGGTGGTGaa                                    >  1:283599/1‑149 (MQ=255)
                                                                                                                     cGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAAGACCGCATCACCGTCACCGACAGGCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCGGGGCCGGGCTGAAACGGGTGGTGaaaaaa                                <  3:653981/149‑1 (MQ=255)
                                                                                                                        ggCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCGGGGCCGGGCTGAAACGGGTGGTGAAAAAAGaa                             <  1:275067/149‑1 (MQ=255)
                                                                                                                           gggTGGTGGAGTAACAGAACCGGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTTAACCGGCGTGAGGTGCGTCATACAGAAGGCGGGGCCGTGCTGAAATGGGTGGGGAAAAAAGAACTg                          >  4:663355/1‑149 (MQ=255)
                                                                                                                            ggtggtGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCGGGGCCGGGCTGAAACGGGTGGTGAAAAAAGAACTgg                         <  4:295478/149‑1 (MQ=255)
                                                                                                                             gtggtgGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCGGGGCCGGGCTGAAACGGGTGGTGAAAAAAGAACTGGc                        <  1:280654/149‑1 (MQ=255)
                                                                                                                                 tggAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCGGGGCCGGGCTGAAACGGGTGGTGAAAAAAGAACTGGCGGAc                    >  3:236089/1‑149 (MQ=255)
                                                                                                                                     gCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCGGGGCCGGGCTGAAACGGGTGGTGAAAAAAGAACTGGCGGACGGCa                <  4:48474/149‑1 (MQ=255)
                                                                                                                                     gCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCGGGGCCGGGCTGAAACGGGTGGTGAAAAAAGAACTGGCGGACGGCa                <  4:646975/149‑1 (MQ=255)
                                                                                                                                          ttaaCCCTGCCGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAg                                                                                >  1:646981/3‑80 (MQ=25)
                                                                                                                                          ttaaCCCTGCCGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAg                                                                                <  2:646981/78‑1 (MQ=255)
                                                                                                                                          ttaaCCCTGCCGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCGGGGCCggg                                                     <  1:233035/105‑1 (MQ=255)
                                                                                                                                          ttaaCCCTGCCGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCGGGGCCggg                                                     >  2:233035/3‑107 (MQ=255)
                                                                                                                                          tGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCACATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCGGGGCCGGGCTGAAACGGGTGGTGAAAAAAGAACTGGCGGACGGCAGCGTc           <  4:529734/149‑1 (MQ=255)
                                                                                                                                                   cAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCGGGGCCGGGCTGAAACGGGTGGTGAAAAAAGAACTGGCGGACGGCAGCGTCACGCGCAGc  >  3:139030/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TATATGACCGCAGCAATACGCAGGTGCGCGCTTTCACGTATGACGCGCAGCACCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCGGGGCCGGGCTGAAACGGGTGGTGAAAAAAGAACTGGCGGACGGCAGCGTCACGCGCAGC  >  NC_000913/524246‑524541

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: ‑