Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,822,673 C→T 7.3% intergenic (‑34/+35) recX ← / ← recA regulatory protein for RecA/DNA recombination and repair protein; ssDNA‑dependent ATPase; synaptase; ssDNA and dsDNA binding protein; ATP‑dependent homologous DNA strand exchanger; recombinase A; LexA autocleavage cofactor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,822,6730CT 7.3% 338.0 / 10.8 109intergenic (‑34/+35)recX/recAregulatory protein for RecA/DNA recombination and repair protein; ssDNA‑dependent ATPase; synaptase; ssDNA and dsDNA binding protein; ATP‑dependent homologous DNA strand exchanger; recombinase A; LexA autocleavage cofactor
Reads supporting (aligned to +/- strand):  ref base C (53/48);  new base T (5/3);  total (58/51)
Fisher's exact test for biased strand distribution p-value = 7.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

Only 100 of 109 total aligned reads displayed.

GCGAGTTTACGTCGCAGTTCTTGCTCACTGTGATCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAACTCACGTACTTTCTTCTCGATCTCTTTC  >  NC_000913/2822530‑2822819
                                                                                                                                               |                                                                                                                                                  
gcgAGTTTACGTCGCAGTTCTTGCTCACTGTGATCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAg                                                                                                                                               >  1:945346/1‑149 (MQ=255)
 cgAGTTTACGTCGCAGTTCTTGCTCACTGTGATCGCGCACCGCCAGAATGCGTACCGCACGATCCAATAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATTCTTACAACTTAAAAAAGCAAAAGGGCCGCAGa                                                                                                                                              <  3:93894/149‑1 (MQ=255)
      ttACGTCGCAGTTCTTGCTCACTGTGATCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGa                                                                                                                                         <  2:683332/149‑1 (MQ=255)
         cGTCGCAGTTCTTGCTCACTGTGATCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCGTATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGTCGCATCTGCGGccc                                                                                                                                      <  4:858957/149‑1 (MQ=255)
             gCAGTTCTTGCTCACTGTGATCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTtgt                                                                                                                                  >  3:1008904/1‑149 (MQ=255)
               aGTTCTTGCTCACTGTGATCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTgtgt                                                                                                                                >  3:715250/1‑149 (MQ=255)
                 ttCTTGCTCACTGTGATCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGTCGCATCTGCGGCCCTtttgctt                                                                                                                              <  1:113249/149‑4 (MQ=255)
                    ttGCTCACTTTGATCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTATAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTATATCaa                                                                                                                           >  2:660155/1‑149 (MQ=255)
                    ttGCTCACTGTGATCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCaa                                                                                                                           >  4:718209/1‑149 (MQ=255)
                        tCACTGTGATCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGTCGCATCTGCGGCCCttt                                                                                                                                   >  4:891842/1‑136 (MQ=255)
                        tCACTGTGATCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGTCGCATCTGCGGCCCttt                                                                                                                                   <  3:891842/137‑2 (MQ=255)
                        tCACTGTGATCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCaaa                                                                                                                          <  2:298540/146‑1 (MQ=255)
                        tCACTGTGATCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCaaa                                                                                                                          >  1:298540/1‑146 (MQ=255)
                        tCACTGTGATCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACaa                                                                                                                       <  1:995997/149‑1 (MQ=255)
                            tgtgATCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACg                                                                                                                   >  3:197654/1‑149 (MQ=255)
                            tgtgATCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACg                                                                                                                   >  4:422690/1‑149 (MQ=255)
                             gtgATCGCGCACCGCCAGAATGCGTACCGCTCGAGCCAGCAGGCGAGCATATGCCGAGCGACGGGAGGTTGATCCTGACTGGGCATATCCTTATAACTACAAAAAGCCGAAGGGCCGCAGATGCGAACCTTTTGTACCAAACAGGActa                                                                                                                  <  4:667743/149‑3 (MQ=255)
                             gtgATCGCGCACCGCCAGAATGCGTACCGCACGATCGAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCGTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGa                                                                                                                  >  1:551132/1‑149 (MQ=255)
                              tgATCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGAt                                                                                                                 >  4:263820/1‑149 (MQ=255)
                               gATCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGAtt                                                                                                                >  4:757803/1‑149 (MQ=255)
                                   gcgcACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTaaaa                                                                                                            <  4:526557/149‑1 (MQ=255)
                                    cgcACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTaaaaa                                                                                                           <  1:531572/149‑1 (MQ=255)
                                    cgcACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTaaaaa                                                                                                           <  4:699380/149‑1 (MQ=255)
                                     gcACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAAt                                                                                                          <  3:380520/149‑1 (MQ=255)
                                     gcACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAAt                                                                                                          >  4:698920/1‑149 (MQ=255)
                                     gcACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAAt                                                                                                          <  1:130032/149‑1 (MQ=255)
                                       aCCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCt                                                                                                        >  3:255425/1‑149 (MQ=255)
                                          gCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCg                                                                                                     >  2:534233/1‑149 (MQ=255)
                                           ccAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGt                                                                                                    <  3:585784/149‑1 (MQ=255)
                                           ccAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGt                                                                                                    <  2:19508/149‑1 (MQ=255)
                                               aaTGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGt                                                                                                <  4:508366/149‑1 (MQ=255)
                                                    gTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTa                                                                                                  >  1:305772/1‑142 (MQ=255)
                                                    gTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTa                                                                                                  <  2:305772/142‑1 (MQ=255)
                                                         gCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGTGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACg                                                                                      >  3:239169/1‑149 (MQ=255)
                                                           aCGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGcc                                                                                    >  1:392847/1‑149 (MQ=255)
                                                               tCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCg                                                                                >  1:291992/1‑149 (MQ=255)
                                                               tCCAACAGGCGAGCATATGCCGGGCGACGGGATCTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCg                                                                                <  4:715250/149‑1 (MQ=255)
                                                                 cAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCt                                                                              >  2:612613/1‑149 (MQ=255)
                                                                  aaCAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTa                                                                             >  4:614085/1‑149 (MQ=255)
                                                                   aCAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTat                                                                            <  3:468584/149‑1 (MQ=255)
                                                                    cAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTatc                                                                           <  3:165877/149‑1 (MQ=255)
                                                                       gCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTatcatc                                                                        <  3:757803/149‑1 (MQ=255)
                                                                         gAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTa                                                                      <  3:422690/149‑1 (MQ=255)
                                                                         gAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTa                                                                      >  1:888800/1‑149 (MQ=255)
                                                                          aGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTAc                                                                     <  3:614085/149‑1 (MQ=255)
                                                                          aGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGAAAAAGGGCCGCAGATGCGACCCTTGTGTATCACACAAGACGATTCAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTAc                                                                     <  4:843872/149‑1 (MQ=255)
                                                                             atatGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACaga                                                                  >  3:929242/1‑149 (MQ=255)
                                                                                 gCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGaaa                                                              >  2:153646/1‑149 (MQ=255)
                                                                                 gCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGaaa                                                              >  1:543582/1‑149 (MQ=255)
                                                                                 gCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGaaa                                                              >  2:197904/1‑149 (MQ=255)
                                                                                          cGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGtt                                                     <  3:1078966/149‑1 (MQ=255)
                                                                                            ggATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGa                                                   >  4:539059/1‑149 (MQ=255)
                                                                                            ggATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGa                                                   <  1:411933/149‑1 (MQ=255)
                                                                                             gATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTGTCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAg                                                  >  2:627303/1‑149 (MQ=255)
                                                                                              aTGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGt                                                 >  4:517395/1‑149 (MQ=255)
                                                                                                gTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTc                                               >  3:496190/1‑149 (MQ=255)
                                                                                                 ttGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCg                                              <  2:469373/149‑1 (MQ=255)
                                                                                                  tGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCgg                                             <  2:284373/149‑1 (MQ=255)
                                                                                                   gATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTatcatc                                                                        <  2:446030/121‑1 (MQ=255)
                                                                                                   gATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTatcatc                                                                        >  1:446030/1‑121 (MQ=255)
                                                                                                    aTTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGtt                                                     >  1:9849/1‑139 (MQ=255)
                                                                                                    aTTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGtt                                                     <  2:9849/139‑1 (MQ=255)
                                                                                                     ttCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACaga                                                                  <  4:351392/125‑1 (MQ=255)
                                                                                                     ttCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACaga                                                                  >  3:351392/1‑125 (MQ=255)
                                                                                                     ttCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGACATCCGGCGTTGAGTTCGGGtt                                          >  4:781394/1‑149 (MQ=255)
                                                                                                      tCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTg                                         <  2:1080457/149‑1 (MQ=255)
                                                                                                      tCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTg                                         >  2:506313/1‑149 (MQ=255)
                                                                                                        tGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCt                                       <  2:102158/149‑1 (MQ=255)
                                                                                                         gTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTc                                      <  3:1052229/149‑1 (MQ=255)
                                                                                                            aTGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTcagc                                   <  3:643337/149‑1 (MQ=255)
                                                                                                             tGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTcagca                                  <  4:211158/149‑1 (MQ=255)
                                                                                                              ggCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCgg                                                         <  3:1063387/125‑1 (MQ=255)
                                                                                                              ggCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCgg                                                         >  4:1063387/1‑125 (MQ=255)
                                                                                                                 atatCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCaa                              >  1:879264/1‑149 (MQ=255)
                                                                                                                 atatCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCaa                              >  3:1013900/1‑149 (MQ=255)
                                                                                                                       ttACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAACTCACg                        >  1:262148/1‑149 (MQ=255)
                                                                                                                       ttACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAACTCACg                        >  1:707148/1‑149 (MQ=255)
                                                                                                                       ttACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAACTCACg                        >  1:833761/1‑149 (MQ=255)
                                                                                                                          cAACTAAAAAAAGCAAAAGGGCCTCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTTCTCAGCAGCAACTCACGTAc                     <  2:551132/149‑1 (MQ=255)
                                                                                                                             cTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTa                                                                                                  <  1:804912/69‑1 (MQ=255)
                                                                                                                             cTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTa                                                                                                  >  2:804912/1‑69 (MQ=255)
                                                                                                                             cTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAACTCACGTACttt                  >  4:1024394/1‑149 (MQ=255)
                                                                                                                              ttAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAACTCACGTACTttc                 >  4:350753/1‑149 (MQ=255)
                                                                                                                              ttAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAACTCACGTACTttc                 >  3:22518/1‑149 (MQ=255)
                                                                                                                                 aaaaaGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAACTCACGTACTttcttc              >  2:383152/1‑149 (MQ=255)
                                                                                                                                 aaaaaGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAACTCACGTACTttcttc              <  1:387031/149‑1 (MQ=255)
                                                                                                                                  aaaaGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAACTCACGTACTTtcttct             >  2:765890/1‑149 (MQ=255)
                                                                                                                                   aaaGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAACTCACGTACTTTCTtctc            >  4:113951/1‑149 (MQ=255)
                                                                                                                                   aaaGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAACTCACGTACTTTCTtctc            <  1:197904/149‑1 (MQ=255)
                                                                                                                                     aGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAACTCACGTACTTTCTTCTCGa          <  1:153646/149‑1 (MQ=255)
                                                                                                                                      gCAAAAGTGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTACCATATACAGAGAAATCCGGCGTTGAGTTCGGGTCGCTCAGCAGCAACTCACGTACTTTCTTCTCGAt         >  4:114991/1‑149 (MQ=255)
                                                                                                                                         caaGGGTCGCATCTGCGGCCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAACTCACGTACTTTCTTCTCGAtctc      >  2:488276/2‑149 (MQ=255)
                                                                                                                                          aaGGGTCGCATCTGCGGCCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCaa                              >  4:726646/1‑124 (MQ=255)
                                                                                                                                           aGGGTCGCATCTGCGGCCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAACTCACGTACTTTCTTCTCGATCTCtt    >  4:337834/1‑149 (MQ=255)
                                                                                                                                            gggTCGCATCTGCGGCCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAACTCACGTACTTTCTTCTCGATCTCttt   >  1:1104120/1‑149 (MQ=255)
                                                                                                                                             ggCCGCAGATGCGACCCTTGTGTGTCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAACTCACGTACTTTCTTCTCGATCTCTTTc  <  1:627303/149‑1 (MQ=255)
                                                                                                                                             ggCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCtt                                                                                                       >  4:181845/1‑48 (MQ=255)
                                                                                                                                             ggCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCtt                                                                                                       <  3:181845/48‑1 (MQ=255)
                                                                                                                                             ggCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAACTCACGTACTTTCTTCTCGATCTCTTTc  <  4:496190/149‑1 (MQ=255)
                                                                                                                                             ggCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCGTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAACTCACGTACTTTCTTCTCGATCTCTTTc  <  1:323691/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                  
GCGAGTTTACGTCGCAGTTCTTGCTCACTGTGATCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAACTCACGTACTTTCTTCTCGATCTCTTTC  >  NC_000913/2822530‑2822819

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