Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 524,391 G→T 6.2% P377P (CCG→CCT)  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,3910GT 6.2% 276.7 / 6.2 80P377P (CCG→CCT) rhsDRhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor
Reads supporting (aligned to +/- strand):  ref base G (29/46);  new base T (2/3);  total (31/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.17e-01

GTATATGACCGCAGCAATACGCAGGTGCGCGCTTTCACGTATGACGCGCAGCACCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCGGGGCCGGGCTGAAACGGGTGGTGAAAAAAGAACTGGCGGACGGCAGCGTC  >  NC_000913/524245‑524532
                                                                                                                                                  |                                                                                                                                             
gTATATGACCGCAGCAATACGCAGGTGCGCGCTTTCACGTATGACGCGCAGCACCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGc                                                                                                                                             <  1:85867/149‑1 (MQ=255)
 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                                                                                                               <  3:491312/149‑1 (MQ=255)
                              gcTTTCACGTATGACGCGCAGCACCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCa                                                                                                               >  4:299453/1‑149 (MQ=255)
                              gcTTTCACGTATGACGCGCAGCACCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCa                                                                                                               >  2:248527/1‑149 (MQ=255)
                                    aCGTATGACGCGCAGCACCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCg                                                                                                         >  1:398459/1‑149 (MQ=255)
                                    aCGTATGACGCGCAGCACCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCg                                                                                                         <  2:390066/149‑1 (MQ=255)
                                      gTATGACGCGCAGCACCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCa                                                                                                       >  1:168733/1‑149 (MQ=255)
                                         tGACGCGCAGCACCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCAtca                                                                                                    >  2:275067/1‑149 (MQ=255)
                                         tGACGCGCAGCACCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCAtca                                                                                                    <  4:379031/149‑1 (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                                                                          <  4:331618/149‑1 (MQ=255)
                                                                   tggCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGt                                                                          <  1:482461/149‑1 (MQ=255)
                                                                   tggCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGt                                                                          <  3:528415/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                                                                                                  <  2:219368/117‑1 (MQ=255)
                                                                           cGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCAtcacc                                                                                                  >  1:219368/1‑117 (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                                                   <  2:139894/149‑1 (MQ=255)
                                                                                          ccGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCgg                                                   <  3:357090/149‑1 (MQ=255)
                                                                                          ccGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCgg                                                   <  4:149809/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:229594/149‑1 (MQ=255)
                                                                                                            tacGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCGGGGCCGGGCTGAAACGggt                                 <  1:479654/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                                                                       <  2:646981/78‑1 (MQ=255)
                                                                                                                                           ttaaCCCTGCCGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAg                                                                       >  1:646981/3‑80 (MQ=25)
                                                                                                                                           ttaaCCCTGCCGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCGGGGCCggg                                            >  2:233035/3‑107 (MQ=255)
                                                                                                                                           ttaaCCCTGCCGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCGGGGCCggg                                            <  1:233035/105‑1 (MQ=255)
                                                                                                                                           tGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCACATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCGGGGCCGGGCTGAAACGGGTGGTGAAAAAAGAACTGGCGGACGGCAGCGTc  <  4:529734/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                             
GTATATGACCGCAGCAATACGCAGGTGCGCGCTTTCACGTATGACGCGCAGCACCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCGGGGCCGGGCTGAAACGGGTGGTGAAAAAAGAACTGGCGGACGGCAGCGTC  >  NC_000913/524245‑524532

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