Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,621,090 C→T R633R (CGC→CGT rhsB → 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_0009133,621,0900CT96.5% 723.8 / 24.3 395R633R (CGC→CGTrhsBRhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor
Reads supporting (aligned to +/- strand):  ref base C (6/8);  new base T (201/180);  total (207/188)
Fisher's exact test for biased strand distribution p-value = 5.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.06e-01

Only 100 of 399 total aligned reads displayed.

TGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACCACGACCTG  >  NC_000913/3620942‑3621234
                                                                                                                                                    |                                                                                                                                                
ggAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAg                                                                                                                                                <  3:781510/150‑1 (MQ=14)
     gACACGCAGGGCCATGAAACTCGGTATGAATACAACATCTCCGGTGACCTGACCTCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTAAGATGCGTGTGGAAAGGCCGTCCGTACCACTCAGGGAGGTCTGACGCGTAGTATgg                                                                                                                                           <  2:8903/151‑1 (MQ=11)
     gACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGCCGCGTAGTATgg                                                                                                                                           >  3:1058696/1‑151 (MQ=14)
     gACACGCAGGGCCATCAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATgg                                                                                                                                           <  8:542706/151‑1 (MQ=14)
        acGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAAt                                                                                                                                        >  8:1286650/1‑151 (MQ=14)
         cGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAAt                                                                                                                                        >  7:931568/1‑150 (MQ=14)
          gCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATAc                                                                                                                                      >  3:622309/1‑151 (MQ=14)
           cAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACg                                                                                                                                     >  4:1633950/1‑151 (MQ=14)
            aGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACg                                                                                                                                     >  1:1224882/1‑150 (MQ=14)
                ccATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGAtgct                                                                                                                                <  6:1242212/151‑1 (MQ=14)
                 cATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGAtgctg                                                                                                                               <  6:1481036/151‑1 (MQ=14)
                  aTGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCAGTATGGAATACGAtgctgc                                                                                                                              <  8:1245437/151‑1 (MQ=31)
                   tGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGcc                                                                                                                             >  6:895512/1‑151 (MQ=14)
                      aaCGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGCCGCGTAGTATGGAATACGATGCTGACgg                                                                                                                           >  4:394027/1‑150 (MQ=9)
                           ggTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAACCGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGt                                                                                                                     <  7:169669/151‑1 (MQ=14)
                             tATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTc                                                                                                                    >  7:1141108/1‑150 (MQ=14)
                              aTGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCAt                                                                                                                  <  6:432600/151‑1 (MQ=14)
                                gAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATcc                                                                                                                >  1:884700/1‑151 (MQ=14)
                                gAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATcc                                                                                                                <  2:593346/151‑1 (MQ=14)
                                gAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATcc                                                                                                                <  4:1567443/151‑1 (MQ=14)
                                gAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATcc                                                                                                                <  8:1668333/151‑1 (MQ=14)
                                 aaTACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCg                                                                                                               >  5:1294864/1‑151 (MQ=14)
                                 aaTACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCg                                                                                                               <  1:176032/151‑1 (MQ=14)
                                  aTACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGc                                                                                                              >  8:367349/1‑151 (MQ=14)
                                   tACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGcc                                                                                                             >  5:229722/1‑151 (MQ=14)
                                    acaacaTCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCt                                                                                                            >  1:1053628/1‑151 (MQ=14)
                                     caacaTCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGCCGCGTAGTATGGAATACGATGCTGCCGGACGGGTCCTCCGCCTg                                                                                                           >  4:1073641/1‑151 (MQ=9)
                                         aTCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGAcc                                                                                                        >  1:459058/1‑150 (MQ=14)
                                           cGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGt                                                                                                     <  4:476919/151‑1 (MQ=14)
                                             ccGGTGACCTGACCGCCGTCATTGCCCCTGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGa                                                                                                   <  3:892800/151‑1 (MQ=14)
                                             ccGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGa                                                                                                   <  5:1006551/151‑1 (MQ=14)
                                              cGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGa                                                                                                   >  2:290428/1‑150 (MQ=14)
                                                gTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGaaaa                                                                                                <  5:1634421/151‑1 (MQ=14)
                                                gTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGACCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGaaaa                                                                                                <  1:671621/151‑1 (MQ=14)
                                                  gacctgaccGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACg                                                                                              <  6:1224335/151‑1 (MQ=14)
                                                  gacctgaccGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACg                                                                                              <  8:1510372/151‑1 (MQ=31)
                                                     ctgaccGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGc                                                                                            >  3:239644/1‑150 (MQ=14)
                                                      tgaccGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGc                                                                                            >  6:553238/1‑149 (MQ=14)
                                                      tgaccGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGc                                                                                            >  6:1036455/1‑149 (MQ=14)
                                                       gaccGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGCCCAGTGAAAACGGCAGc                                                                                         >  8:945720/1‑151 (MQ=14)
                                                       gaccGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGc                                                                                         <  4:1182938/151‑1 (MQ=14)
                                                        accGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGCCGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGcc                                                                                        >  7:1261187/1‑151 (MQ=14)
                                                          cgccgTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCcac                                                                                      >  2:1579069/1‑151 (MQ=14)
                                                             cgTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCAcacc                                                                                   <  8:1002288/151‑1 (MQ=14)
                                                              gTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGCCGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCAcaccc                                                                                  >  4:702512/1‑150 (MQ=14)
                                                              gTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGCCGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCAcacca                                                                                  <  8:1212247/151‑1 (MQ=14)
                                                                 aTTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCt                                                                               <  8:851486/151‑1 (MQ=14)
                                                                   tGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCtt                                                                              <  2:1491329/150‑1 (MQ=14)
                                                                   tGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGGCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCtt                                                                              <  4:826326/150‑1 (MQ=14)
                                                                    gCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGTGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTcc                                                                            <  6:354537/151‑1 (MQ=14)
                                                                    gCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTcc                                                                            <  7:1075359/151‑1 (MQ=14)
                                                                     ccccGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCg                                                                           >  3:985720/1‑151 (MQ=14)
                                                                          gACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACg                                                                      <  1:1532603/151‑1 (MQ=14)
                                                                           aCGGCAGCAGAAACGTGACACAGTACGATGCGTGGGGAAAGGCGGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGa                                                                     <  1:326029/151‑1 (MQ=9)
                                                                           aCGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACg                                                                      >  8:257221/1‑150 (MQ=14)
                                                                                   agAAACGGGACACAGTACGATGCGGGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTg                                                             <  1:566735/151‑1 (MQ=14)
                                                                                    gAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTg                                                             >  5:884312/1‑150 (MQ=14)
                                                                                      aaCGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGAcc                                                          <  4:546124/151‑1 (MQ=14)
                                                                                      aaCGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGAcc                                                          <  7:602059/151‑1 (MQ=14)
                                                                                             cacaGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGAt                                                   <  5:412306/151‑1 (MQ=14)
                                                                                             cacaGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGAt                                                   >  3:315454/1‑151 (MQ=14)
                                                                                             cacaGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGAt                                                   <  3:4403/151‑1 (MQ=14)
                                                                                              acaGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATa                                                  <  2:1047363/151‑1 (MQ=14)
                                                                                              acaGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATa                                                  <  4:1304757/151‑1 (MQ=14)
                                                                                                aGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATAc                                                 >  2:1330892/1‑150 (MQ=14)
                                                                                                 gTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATAc                                                 >  7:642239/1‑149 (MQ=14)
                                                                                                  tACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAgg                                              >  4:1345229/1‑151 (MQ=14)
                                                                                                   aCGATGCGTGGGGAAACGCCGTCCGGACCACGCAGGGCGGGCTGACGCGTAGGATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAgg                                              <  2:931566/150‑1 (MQ=9)
                                                                                                    cGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGaa                                            <  6:250399/151‑1 (MQ=14)
                                                                                                    cGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGa                                             >  8:1416195/1‑150 (MQ=14)
                                                                                                     gATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGaa                                            >  2:41337/1‑150 (MQ=14)
                                                                                                       tGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAAcc                                         >  5:1196244/1‑151 (MQ=14)
                                                                                                        gCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAAcc                                         <  7:391485/150‑1 (MQ=14)
                                                                                                        gCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAAcc                                         <  3:152459/150‑1 (MQ=14)
                                                                                                        gCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCg                                        <  6:1605071/151‑1 (MQ=14)
                                                                                                          gTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGc                                      >  3:853088/1‑151 (MQ=14)
                                                                                                           nngggAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCt                                     >  2:233044/3‑151 (MQ=14)
                                                                                                           gggggAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCt                                     <  5:1536238/150‑1 (MQ=14)
                                                                                                               gAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGCAACCGGCTTTgc                                 >  1:1584391/1‑150 (MQ=14)
                                                                                                                aaaGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGAc                                <  3:1635064/151‑1 (MQ=14)
                                                                                                                aaaGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGAc                                <  5:533546/151‑1 (MQ=14)
                                                                                                                     ccgtccgtACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCg                           <  2:787014/151‑1 (MQ=14)
                                                                                                                      cgtccgtACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGc                          <  5:960816/151‑1 (MQ=14)
                                                                                                                        tccgtACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGcac                        >  2:176855/1‑151 (MQ=14)
                                                                                                                         ccgtACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGcaca                       <  3:1200455/151‑1 (MQ=14)
                                                                                                                           gtACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGcacac                      >  6:1682815/1‑150 (MQ=14)
                                                                                                                           gtACCACGCAGGGCGGGCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGcacaca                     <  2:1212544/151‑1 (MQ=31)
                                                                                                                             aCCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAg                    <  4:262058/150‑1 (MQ=14)
                                                                                                                               cACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCCTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTATGACCGCCGCACACAGcct                 >  3:1367861/1‑149 (MQ=255)
                                                                                                                                aCGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGtt                >  5:1471063/1‑151 (MQ=14)
                                                                                                                                    aGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATca            <  1:918203/151‑1 (MQ=14)
                                                                                                                                    aGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATc             >  5:690956/1‑150 (MQ=14)
                                                                                                                                     gggcgggcTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATcac           >  6:1069971/1‑151 (MQ=14)
                                                                                                                                     gggcgggcTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATcac           >  6:1137970/1‑151 (MQ=14)
                                                                                                                                     gggcgggcTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATcac           >  6:1405250/1‑151 (MQ=31)
                                                                                                                                       gcgggcTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATcacca         <  6:824964/151‑1 (MQ=14)
                                                                                                                                         gggcTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACCACg       >  4:174200/1‑151 (MQ=14)
                                                                                                                                              gACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACCACGAcc    <  7:288188/149‑1 (MQ=11)
                                                                                                                                              gACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACCACGAcc    <  7:1554068/149‑1 (MQ=11)
                                                                                                                                              gACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACCACGACCTg  >  7:759347/1‑151 (MQ=11)
                                                                                                                                                    |                                                                                                                                                
TGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACCACGACCTG  >  NC_000913/3620942‑3621234

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: