Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,822,684 A→G 6.6% intergenic (‑45/+24) 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,6840AG 6.6% 461.4 / 18.4 169intergenic (‑45/+24)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 A (72/78);  new base G (5/6);  total (84/85)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.93e-01

Only 100 of 170 total aligned reads displayed.

TTACGTCGCAGTTCTTGCTCACTGTGATCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAACTCACGTACTTTCTTCTCGATCTCTTTCGCGGTTTCCGGGT  >  NC_000913/2822536‑2822832
                                                                                                                                                    |                                                                                                                                                    
ttACGTCGCAGTTCTTGCTCACTGTGATCGCGCAACGCCAGAATGAGTACCGCACGATCCAACAGGAGAGCATATGCCGGGCGACGGGATGTTGATTCTGGCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGAc                                                                                                                                                     >  1:1228385/1‑150 (MQ=255)
 tACGTCGCAGTTCTTGCTCACTGTGATCGCGCACCGCCAGAATGCGTACCGCACGATCCATCAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAACGCAAAAGGGCCGCAGATGCGccc                                                                                                                                                    >  2:280834/1‑150 (MQ=255)
   cGTCGCAGTTCTTGCTCACTTTGAGCGCGCACCGCCAGAATGCGAACCGCACGATCCAACAGGCGAGCATATGCCGTGCGACTGGATGTTGATTCTGTCATGGAATATCCTTACAACTTAAAAAAGCAAAAGGGTCACATCTGCGGCCCt                                                                                                                                                  <  2:1871469/150‑1 (MQ=255)
     tCGCAGTTCTTGCTCACTGTGATCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTtg                                                                                                                                                >  1:178567/1‑150 (MQ=255)
       gCAGTTCTTGCTCACTGTGATCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTtgtg                                                                                                                                              <  1:1090040/150‑1 (MQ=255)
       gCAGTTCTTGCTCACTGTGATCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTtgtg                                                                                                                                              <  1:1858503/150‑1 (MQ=255)
       gCAGTTCTTGCTCACTGTGATCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTtgtg                                                                                                                                              >  1:1506420/1‑150 (MQ=255)
         aGTTCTTGCTCTCTGTGCTCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACATCTTAAAAAAGCAAAAGGGCCGCATATGCGACCCTTGTGTa                                                                                                                                            <  1:70660/150‑1 (MQ=255)
         aGTTCTTGCTCACTGTGATCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGTCGCATCTGCGGCCCTtttgct                                                                                                                                            <  2:1960851/150‑3 (MQ=255)
         aGTTCTTGCTCACTGTGATCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTa                                                                                                                                            <  2:669765/150‑1 (MQ=255)
         aGTTCTTGCTCACTGTGATCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTa                                                                                                                                            <  2:859111/150‑1 (MQ=255)
          gTTCATGCTCACTGTGATCGCGCACCTCCAGAATGCGTACCGCACGATCCAACAGGCTAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACCACTTAAAAAAGCAAAAGGGTCGCATCTGCGGCCCTtttgctt                                                                                                                                           <  1:1203004/150‑4 (MQ=255)
              ttGCTCACTGTGATCGCGCACCGCCAGAATGCCTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGTCGCATCTGCGGCCCTtttgctttttt                                                                                                                                       <  2:1092184/150‑8 (MQ=255)
               tGCTCACTGTGATCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGATGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAAc                                                                                                                                      >  2:1173188/1‑150 (MQ=255)
                 cTCACTGTGATCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACaa                                                                                                                                    <  1:401121/150‑1 (MQ=255)
                     cTGTGATCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACg                                                                                                                                >  1:1716155/1‑150 (MQ=255)
                       gtgATCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGAt                                                                                                                              >  1:907571/1‑150 (MQ=255)
                       gtgATCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGAt                                                                                                                              <  2:46270/150‑1 (MQ=255)
                          aTCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCAAGTGCCGGGCGACGGGAGGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTaa                                                                                                                           >  2:2013559/1‑150 (MQ=255)
                            cgcgCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGCCGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTaaaa                                                                                                                         >  1:107437/1‑150 (MQ=255)
                            cgcgCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGCCCCTTGTGTATCAAACAAGACGATTaaaa                                                                                                                         >  2:827214/1‑150 (MQ=255)
                                cACCGCCAGAATGCGTACCGCACGATCCCACAGGCGAGCATATGCCGGGCGACGGGATGTTGAGTCTGTCATGGCATATCCTGACAACTTAACAAAGCAAAAGGGCCGCAGCTGCGCCCCTTGTGTATCAAACAAGACGAGTAAAAATcg                                                                                                                     >  1:285191/1‑149 (MQ=255)
                                cACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCCGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCt                                                                                                                     >  2:1225374/1‑150 (MQ=255)
                                  ccGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTc                                                                                                                   >  1:1140712/1‑150 (MQ=255)
                                  ccGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTc                                                                                                                   <  2:1343742/150‑1 (MQ=255)
                                   cGCCAGAATGCGTACCGCCCGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTAAAAAAACAAAAAGGGCCGACGGGGCGACCCTGGGGTTACAAACATGACGATTAAAAATCTTCg                                                                                                                  >  2:1730036/1‑150 (MQ=255)
                                      cAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTa                                                                                                               >  2:664320/1‑150 (MQ=255)
                                      cAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTa                                                                                                               <  1:1106685/150‑1 (MQ=255)
                                         aaTGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGtt                                                                                                            <  1:590170/150‑1 (MQ=255)
                                             cGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGAAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTg                                                                                                        >  2:69829/1‑150 (MQ=255)
                                              gTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCtac                                                                                                       >  2:1670674/1‑148 (MQ=255)
                                              gTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGc                                                                                                       <  2:1109877/150‑1 (MQ=255)
                                              gTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGc                                                                                                       >  1:1145743/1‑150 (MQ=255)
                                               tACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCt                                                                                                      >  2:633441/1‑150 (MQ=255)
                                                       gATCCAACAGCCGAGCATATGCCGGGCGACGGGATGTTTATTCTGTCATGGCATATGCTTACAAATTAAAAAAGCAAAATAGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAAAATTCGGTAGTTTATGCTACGCCtac                                                                                              >  2:1438683/1‑148 (MQ=255)
                                                         tCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGc                                                                                            <  1:194897/150‑1 (MQ=255)
                                                         tCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGc                                                                                            <  1:722247/150‑1 (MQ=255)
                                                         tCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGc                                                                                            <  1:1399035/150‑1 (MQ=255)
                                                          acAACAGGCTAGCATATGCCGGGCGACGGGATGTTGATTGTGCCATGTCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCt                                                                                           <  1:1289204/149‑1 (MQ=255)
                                                           cAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTTAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTa                                                                                          >  1:1020734/1‑150 (MQ=255)
                                                           cAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTa                                                                                          <  2:1563393/150‑1 (MQ=255)
                                                           cAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTa                                                                                          >  1:1464661/1‑150 (MQ=255)
                                                           cAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTa                                                                                          <  1:1583681/150‑1 (MQ=255)
                                                           cAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTATATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCtt                                                                                          >  1:1884915/1‑149 (MQ=255)
                                                             aaacgCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCACAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTatc                                                                                        <  2:527017/146‑1 (MQ=255)
                                                             aCAGGCGAGCATATGCCGGGCGACGGGATGTTTATTCTGTCATGGCATATCCTTACAACTAAAAAAAGCAAAAGGGCCGCATATGCGACCCTTGTGTATCAAACAAGCCGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTatc                                                                                        <  2:1061679/150‑1 (MQ=255)
                                                             aCAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTatc                                                                                        <  2:884093/150‑1 (MQ=255)
                                                                ggCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTatcata                                                                                     >  1:1208165/1‑149 (MQ=255)
                                                                 gctaGCATATGCCGGGCGACGGGATGTTGTTTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAGGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGTCGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCt                                                                                    <  1:1902635/147‑1 (MQ=255)
                                                                   gatcATATGCCGGTCGACGGGATGTTGATTCTGTCATGGCTTATCCTTACAACTTAAATAAGCAAAACGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTAc                                                                                  <  1:327062/147‑1 (MQ=255)
                                                                   gAGCATATGCCTGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAATAAAGCATAAGGGCCGCTGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTAc                                                                                  <  1:1163349/150‑1 (MQ=255)
                                                                    aGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACa                                                                                 >  2:1082101/1‑150 (MQ=255)
                                                                     gCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTatcatc                                                                                     <  1:1430491/145‑1 (MQ=255)
                                                                       atatGACGGGCGACGGGAGGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACagag                                                                              <  1:1009665/150‑1 (MQ=255)
                                                                           gCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGCTGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAAt                                                                          >  2:1810403/1‑150 (MQ=255)
                                                                             cGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATcc                                                                        >  1:667029/1‑150 (MQ=255)
                                                                                    cGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTg                                                                 <  2:1473268/150‑1 (MQ=255)
                                                                                      ggATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGCCCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAg                                                               >  1:1089902/1‑150 (MQ=255)
                                                                                            tGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCTGCGTTGAGTTCggg                                                         >  2:1394298/1‑150 (MQ=255)
                                                                                            tGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCggg                                                         >  2:1552021/1‑150 (MQ=255)
                                                                                            tGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCggg                                                         <  1:1764690/150‑1 (MQ=255)
                                                                                            tGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCggg                                                         <  2:178567/150‑1 (MQ=255)
                                                                                             gATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGt                                                        >  2:1451823/1‑150 (MQ=255)
                                                                                             gATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGt                                                        <  1:1417498/150‑1 (MQ=255)
                                                                                             gATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGCTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGt                                                        <  1:198397/150‑1 (MQ=255)
                                                                                              aTTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGtt                                                       <  2:1228385/150‑1 (MQ=255)
                                                                                              aTTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGtt                                                       >  2:1968905/1‑150 (MQ=255)
                                                                                                  tGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTc                                                   <  1:1387277/150‑1 (MQ=255)
                                                                                                   gTCATGGCATATCCATACAACTTAAAAAAGCAAAAGGGCCGCAGCTGCGACCCTTGTGTATCAAACAAGACGAATAAAAATCTACGTTAGTTACTGCTACGCCTTCGCTTTCATCTACAGAGAAATCCGGCGTTGAGATCGGGTTGCTca                                                  >  1:469539/1‑150 (MQ=255)
                                                                                                    tCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTcag                                                 >  2:1548605/1‑150 (MQ=255)
                                                                                                     cATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTcagc                                                <  2:521824/150‑1 (MQ=255)
                                                                                                     cATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTcagc                                                >  2:1416728/1‑150 (MQ=255)
                                                                                                     cATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTcagc                                                <  1:195011/150‑1 (MQ=255)
                                                                                                     cATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTcagc                                                >  2:619946/1‑150 (MQ=255)
                                                                                                       tGGCATATCCTTACAACGTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCGATCATCTACAGAGAAAGCCGGCGGTGAGTTCGGGTTGCTcagcag                                              >  1:339862/1‑150 (MQ=255)
                                                                                                         gCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACCGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAgcagca                                            <  2:920813/150‑1 (MQ=255)
                                                                                                              tCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAACTCa                                       <  2:390347/150‑1 (MQ=255)
                                                                                                               ccTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAACTCAc                                      >  2:373/1‑150 (MQ=255)
                                                                                                               ccTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAACTCAc                                      >  1:1444213/1‑150 (MQ=255)
                                                                                                                cTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAACTCACg                                     >  2:943822/1‑150 (MQ=255)
                                                                                                                cTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAACTCACg                                     >  1:278169/1‑150 (MQ=255)
                                                                                                                 ttACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAACTCACGt                                    <  1:581247/150‑1 (MQ=255)
                                                                                                                       cTTAAAAAAGCAAAAGGGCCGCAGATGCGCCCCTTGTGTATCAAACAAGACGATTAAAACTCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAACTCACGTACttta                              >  2:1312236/1‑149 (MQ=255)
                                                                                                                       cTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAACTCACGTACttta                              >  1:938199/1‑149 (MQ=255)
                                                                                                                       cTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAACTCACGTACttta                              >  1:1325704/1‑149 (MQ=255)
                                                                                                                       cTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAACTCACGTACTttc                              <  2:885619/150‑1 (MQ=255)
                                                                                                                             aaaGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAACTCACGTACTTTCTTCTCg                        <  2:813722/150‑1 (MQ=255)
                                                                                                                               tacacaaGGGTCGCATCTGCGGCCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAACTCCCGTACTTTCTTCTCGAt                      >  2:858519/6‑150 (MQ=255)
                                                                                                                                acacaaGGGTCGCATCTGCGGCCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCGACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAACTCACGTACTTTCTTCTCGAtg                     >  1:1754438/5‑149 (MQ=255)
                                                                                                                                       ggCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAACTCACGTACTTTCTTCTCGATCTCTTTcg              <  1:1648342/150‑1 (MQ=255)
                                                                                                                                        gCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAACTCACGTACTTTCTTCTCGATCTCTTTcgc             >  2:1252393/1‑150 (MQ=255)
                                                                                                                                        gCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAACTCACGTACTTTCTTCTCGATCTCTTTcgc             >  1:550274/1‑150 (MQ=255)
                                                                                                                                        gCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAACTCACGTACTTTCTTCTCGATCTCTTTcgc             >  2:777071/1‑150 (MQ=255)
                                                                                                                                        cccccAGAGGCGACCCTTGTTTATAAAACAAGACTATAAAAAATCTTCTTTAGTTTCTGCACCGCCTCCGCTAACATCACCCCACAAACCCGGCGTTTCGTGCGGGTTGCCCAGCACAAACTCACGTACTTTCTTCTCGATCTCTTgcgc             <  2:1716155/146‑1 (MQ=255)
                                                                                                                                         ccGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAACTCACGTACTTTCTTCTCGATCTCTTTcgcg            <  2:1667754/150‑1 (MQ=255)
                                                                                                                                          cGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAACTCACGTACTTTCTTCTCGATCTCTTTCGCgg           >  1:1707811/1‑150 (MQ=255)
                                                                                                                                          cGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAACTCACGTACTTTCTTCTCGATCTCTTTCGCgg           <  2:893991/150‑1 (MQ=255)
                                                                                                                                          cGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAACTCACGTACTTTCTTCTCGATCTCTTTCGCgg           <  1:1173188/150‑1 (MQ=255)
                                                                                                                                               aTGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAACTCACGTACTTTCTTCTCGATCTCTTTCGCGGTTTcg      >  1:1281324/1‑149 (MQ=255)
                                                                                                                                                   gACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAACTCACGTACTTTCTTCTCGATCTCTTTCGCGGTTTCCGGGt  >  1:880636/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TTACGTCGCAGTTCTTGCTCACTGTGATCGCGCACCGCCAGAATGCGTACCGCACGATCCAACAGGCGAGCATATGCCGGGCGACGGGATGTTGATTCTGTCATGGCATATCCTTACAACTTAAAAAAGCAAAAGGGCCGCAGATGCGACCCTTGTGTATCAAACAAGACGATTAAAAATCTTCGTTAGTTTCTGCTACGCCTTCGCTATCATCTACAGAGAAATCCGGCGTTGAGTTCGGGTTGCTCAGCAGCAACTCACGTACTTTCTTCTCGATCTCTTTCGCGGTTTCCGGGT  >  NC_000913/2822536‑2822832

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