Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,952,107 T→C 5.9% I238T (ATT→ACT)  yecE → UPF0759 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,952,1070TC 5.9% 229.3 / 4.0 68I238T (ATT→ACT) yecEUPF0759 family protein
Reads supporting (aligned to +/- strand):  ref base T (31/33);  new base C (2/2);  total (33/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.33e-01

ATCGCGGTTTACATCAGCGCGGCGTTAATCGGGTGATTTTAGACAGCCGCCCGGTTCATGCAGCACGTCCACACAGTGAAGCTATTCGCGACGCTCAACGAAAAAAACCTAAAGTTCCGGTACATGCTGTACTGACGGCGACAAATCCGCTGATCCGTTTTATCGGTAGTGATGATATGACGCAAAACCGGGAATTATTTCAGGTCTGGTTACAAAAATTAGCGCAGTGGCATCAGACCACTACGCCTTATCTTTTTTTACATACGCCAGACATTGCCCAGGCACCGGAACTGGTACATACCCTGTGGGAAGACTTACGTAAAACGCTTCCAGAGATCGGAGCAGTTCCGGCTATTCCACAGCAATCTTCTCTTTTCTGAATTTGCCACCTATCATAGACAGGTGCCATCGGCCATTTTAAAGGGAGTTTGTATGGTAAGCGCGCTGTACGCCGTTTTAAGTGCGTTGTTATTAATGAAGTTCTCTTTTGATGTCGTTCGCCTGCGAATGCAGTACCGCGTTGCCTATGGCGACGGCGGTTTTAGCGAACTGCAAAGCGC  >  NC_000913/1951834‑1952393
                                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                                              
aTCGCGGTTTACATCAGCGCGGCGTTAATCGGGTGATTTTAGACAGCCGCCCGGTTCATGCAGCACGTCCACACAGTGAAGC‑ATTCGCGACGCTCAACGAAAAAAACCTAAAGTTCCGGTACATGCTGTACTGACGGCGACAAATCCGCTGATCCGTTTTATCGGTAGTGATGATATGACGCAAAACCGGGAATTATTTCAGGTCTGGTTACAAAAATTAGCGCAGTGGCATCAGACCACTACGCCTTATCTTTTTTTACATACGCCAGACATTGCCCAGGCACCGGaa                                                                                                                                                                                                                                                                                <  1:677176/289‑1 (MQ=255)
 tCGCGGTTTACATCAGCGCGGCGTTAATCGGGTGATTTTAGACAGCCGCCCGGTTCATGCAGCACGTCCACACAGTGAAGCTATTCGCGACGCTCAACGAAAAAAACCTAAAGTTCCGGTACATGCTGTACTGACGGCGACAAATCCGCTGATCCGTTTTATCGGTAGTGATGATATGACGCAAAACCGGGAATTATTTCGGGTCTGGTTACAAAAATTAGCGCAGTGGCATCAGACCACTACGCCTTATCTTTTTTTACATACGCCAGACATTGCCCAGGCACCGGaa                                                                                                                                                                                                                                                                                <  1:332907/289‑1 (MQ=255)
 tCGCGGTTTACATCAGCGCGGCGTTAATCGGGTGATTTTAGACAGCCGCCCGGTTCATGCAGCACGTCCACACAGTGAAGCTATTCGCGACGCTCAACGAAAAAAACCTAAAGTTCCGGTACATGCTGTACTGACGGCGACAAATCCGCTGATCCGTTTTATCGGTAGTGATGATATGACGCAAAACCGGGAATTATTTCAGGTCTGGTTACAAAAATTAGCGCAGTGGCATCAGACCACTACGCCTTATCTTTTTTTACATACGCCAGACATTGCCCAGGCACCGGaa                                                                                                                                                                                                                                                                                <  2:461735/289‑1 (MQ=255)
 tCGCGGTTTACATCAGCGCGGCGTTAATCGGGTGATTTTAGACAGCCGCCCGGTTCATGCAGCACGTCCACACAGTGAAGCTATTCGCGACGCTCAACGAAAAAAACCTAAAGTTCCGGTACATGCTGTACTGACGGCGACAAATCCGCTGATCCGTTTTATCGGTAGTGATGATATGACGCAAAACCGGGAATTATTTCAGGTCTGGTTACAAAAATTAGCGCAGTGGCATCAGACCACTACGCCTTATCTTTTTTTACATACGCCAGACATTGCCCAGGCACCGGaa                                                                                                                                                                                                                                                                                <  1:30870/289‑1 (MQ=255)
   gcgGTTTACATCAGCGCGGCGTTAATCGGGTGATTTTAGACAGCCGCCCGGTTCATGCAGCACGTCCACACAGTGAAGCTATTCGCGACGCTCAACGAAAAAAACCTAAAGTTCCGGTACATGCTGTACTGACGGCGACAAATCCGCTGATCCGTTTTATCGGTAGTGATGATATGACGCAAAACCGGGAATTATTTCAGGTCTGGTTACAAAAATTAGCGCAGTGGCATCAGACCACTACGCCTTATCTTTTTTTACATACGCCAGACATTGCCCAGGCACCGGaa                                                                                                                                                                                                                                                                                <  2:703488/287‑1 (MQ=255)
   gcgGTTTACATCAGCGCGGCGTTAATCGGGTGATTTTAGACAGCCGCCCGGTTCATGCAGCACGTCCACACAGTGAAGCTATTCGCGACGCTCAACGAAAAAAACCTAAAGTTCCGGTACATGCTGTACTGACGGCGACAAATCCGCTGATCCGTTTTATCGGTAGTGATGATATGACGCAAAACCGGGAATTATTTCAGGTCTGGTTACAAAAATTAGCGCAGTGGCATCAGACCACTACGCCTTATCTTTTTTTACATACGCCAGACATTGCCCAGGCACCGGAACt                                                                                                                                                                                                                                                                              <  1:343771/289‑1 (MQ=255)
        ttACATCAGCGCGGCGTTAATCGGGTGATTTTAGACAGCCGCCCGGTTCATGCAGCACGTCCACACAGTGAAGCTATTCGCGACGCTCAACGAAAAAAACCTAAAGTTCCGGTACATGCTGTACTGACGGCGACAAATCCGCTGATCCGTTTTATCGGTAGTGATGATATGACGCAAAACCGGGAATTATTTCAGGTCTGGTTACAAAAATTAGCGCAGTGGCATCAGACCACTACGCCTTATCTTTTTTTACATACGCCAGACATTGCCCAGGCACCGGAACTGGTAc                                                                                                                                                                                                                                                                         <  1:31109/289‑1 (MQ=255)
                 cgcgGCGTTAATCGGGTGATTTTAGACAGCCGCCCGGTTCATGCAGCACGTCCACACAGTGAAGCTATTCGCGACGCTCAACGAAAAAAACCTAAGGTTCAGGTACATGCTGTACTGAAGGCGACAAATCCGGTGATCCGTGTTATCGGTAGTGATGATATGACGCAAACCCGGGAATTATTTCAGGTCTGGTTACAAAAATTAGTGCAGTGGCATCAGACCACTACGCATTATCTTTTTGTACATACGCCAGACATTGCATAGGCACCGGAACTGGTaagaactcggt                                                                                                                                                                                                                                                                >  2:190780/1‑279 (MQ=255)
                 cgcgGCGTTAATCGGGTGATTTTAGACAGCCGCCCGGTTCATGCAGCACGTCCACACAGTGAAGCTATTCGCGACGCTCAACGAAAAAAACCTAAAGTTCCGGTACATGCTGTACTGACGGCGACAAATCCGCTGATCCGTTTTATCGGTAGTGATGATATGACGCAAAACCGGGAATTATTTCAGGTCTGGTTACAAAAATTAGCGCAGTGGCATCAGACCACTACGCCTTATCTTTTTTTACATACGCCAGACATTGCCCAGGCACCGGAACTGGTACATACCCtgt                                                                                                                                                                                                                                                                >  1:621509/1‑289 (MQ=255)
                 cgcgGCGTTAATCGGGTGATTTTAGACAGCCGCCCGGTTCATGCAGCACGTCCACACAGTGAAGCTATTCGCGACGCTCAACGAAAAAAACCTAAAGTTCCGGTACATGCTGTACTGACGGCGACAAATCCGCTGATCCGTTTTATCGGTAGTGATGATATGACGCAAAACCGGGAATTATTTCAGGTCTGGTTACAAAAATTAGCGCAGTGGCATCAGACCACTACGCCTTATCTTTTTTTACATACGCCAGACACTGCCCAGGCACCGGAACTGGTACATACCCtgt                                                                                                                                                                                                                                                                >  1:551899/1‑289 (MQ=255)
                 cgcgGCGTTAATCGGGTGATTTTAGACAGCCGCCCGGTTCATGCAGCACGTCCACACAGTGAAGCTATTCGCGACGCTCAACGAAAAAAACCTAAAGTTCCGGTACATGCTGTACTGACGGCGACAAATCCGCTGATCCGTTTTATCGGTAGTGATGATATGACGCAAAACCGGGAATTATTTCAGGTCTGGTTACAAAAATTAGCGCAGTGGCATCAGACCACTACGCCTTATCTTTTTTTACATACGCCAGACACTGCCCAGGCACCGGAACTGGTACATACCCtgt                                                                                                                                                                                                                                                                >  1:551878/1‑289 (MQ=255)
                   cggcgTTAATCGGGTGATTTTAGACAGCCGCCCGGTTCATGCAGCACGTCCACACAGTGAAGCTATTCGCGACGCTCAACGAAAAAAACCTAAAGTTCCGGTACATGCTGTACTGACGGCGACAAATCCGCTGATCCGTTTTATCGGTAGTGATGATATGACGCAAAACCGGGAATTATTTCAGGTCTGGTTACAAAAATTAGCGCAGTGGCATCAGACCACTACGCCTTATCTTTTTTTACATACGCCAGACATTGCCCAGGCACCGGAACTGGTACATACCCTGTgg                                                                                                                                                                                                                                                              >  1:343892/1‑289 (MQ=255)
                                gTGATTTTAGACAGCCGCCCGGTTCATGCAGCACGTCCACACAGTGAAGCTATTCGCGACGCTCAACGAAAAAAACCTAAAGTTCCGGTACATGCTGTACTGACGGCGACAAATCCGCTGATCCGTTTTATCGGTAGTGATGATATGACGCAAAACCGGGAATTATTTCAGGTCTGGTTACAAAAATTAGCGCAGTGGCATCAGAACACTACGCCTTATCTTTTTTTACATACGACAGACATTGCCCAGGCACCGGAACTGGTACATCCCCTGGGGGAAGACTTACGTa                                                                                                                                                                                                                                                 >  2:633784/1‑289 (MQ=255)
                                              ccgccCGGTTCATGCAGCACGTCCACACAGTGAAGCTATTCGCGACGCTCAACGAAAAAAACCTAAAGTTCCGGTACATGCTGTACTGACGGCGACAAATCCGCTGATCCGTTTTATCGGTAGTGATGATATGACGCAAAACCGGGAATTATTTCAGGTCTGGTTACAAAAATTAGCGCAGTGGCATCAGACCACTACGCCTTATCTTTTTTTACATACGCCAGACATTGCCCAGGCACCGGAAc                                                                                                                                                                                                                                                                               >  2:398682/1‑245 (MQ=255)
                                              ccgccCGGTTCATGCAGCACGTCCACACAGTGAAGCTATTCGCGACGCTCAACGAAAAAAACCTAAAGTTCCGGTACATGCTGTACTGACGGCGACAAATCCGCTGATCCGTTTTATCGGTAGTGATGATATGACGCAAAACCGGGAATTATTTCAGGTCTGGTTACAAAAATTAGCGCAGTGGCATCAGACCACTACGCCTTATCTTTTTTTACATACGCCAGACATTGCCCAGGCACCGGAAc                                                                                                                                                                                                                                                                               <  1:398682/245‑1 (MQ=255)
                                              ccgccCGGTTCATGCAGCACGTCCACACAGTGAAGCTATTCGCGACGCTCAACGAAAAAAACCTAAAGTTCCGGTACATGCTGTACTGACGGCGACAAATCCGCTGATCCGTTTTATCGGTAGTGATGATATGACGCAAAACCGGGAATTATTTCAGGTCTGGTTACAAAAATTAGCGCAGTGGCATCAGACCACTACGCCTTATCTTTTTTTACATACGCCAGACATTGCCCAGGCACCGGAACTGGTACATACCCTGTGGGAAGACTTACGTAAAACGCTTCCagag                                                                                                                                                                                                                                   >  1:574903/1‑289 (MQ=255)
                                               cgccAGGTTCATGCAGCACGGCCTCACAGTGAAGCTATTCGCGACGCGCAACGAAAAAAACCTAAAGTTCCGGTACATGCTGTACTGACGGCGACAAATCCGCTGATCCGGTTTATCGGTAGTGATGATATGACGCAAAACCGGGAATTATTTCATGTCTGGTTACAAAAAGTAGCGCAGTGGCATCAGACCACTACGCCTTATCTTTTTTTACATACGCCAGACATTGCCCAGGCACCGGAACTGGTACATACCCTGTGGGAAGACTTACGTAAAACGCTTCCa                                                                                                                                                                                                                                      <  2:553593/285‑1 (MQ=255)
                                                gccCGGTTCATGCAGCACGTCCACACAGTGAAGCTATTCGCGACGCTCAACGAAAAAAACCTAAAGTTCCGGTACATGCTGTACTGACGGCGACAAATCCGCTGATCCGTTTTATCGGTAGTGATGATATGACGCAAAACCGGGAATTATTTCAGGTCTGGTTACAAAAATTAGCGCAGTGGCATCAGACCACTACGCCTTATCTTTTTTTACATACGCCAGACACTGCCCAGGCACCGGAACTGGTACATACCCTGTGGGAAGACTTACGTAAAACGCTTCCAGAGAt                                                                                                                                                                                                                                 <  2:551878/289‑1 (MQ=255)
                                                gccCGGTTCATGCAGCACGTCCACACAGTGAAGCTATTCGCGACGCTCAACGAAAAAAACCTAAAGTTCCGGTACATGCTGTACTGACGGCGACAAATCCGCTGATCCGTTTTATCGGTAGTGATGATATGACGCAAAACCGGGAATTATTTCAGGTCTGGTTACAAAAATTAGCGCAGTGGCATCAGACCACTACGCCTTATCTTTTTTTACATACGCCAGACACTGCCCAGGCACCGGAACTGGTACATACCCTGTGGGAAGACTTACGTAAAACGCTTCCAGAGAt                                                                                                                                                                                                                                 <  2:551899/289‑1 (MQ=255)
                                                     gTTCATGTGGCACGGCCAGAAAGTGAAGATGTTCCGGAAGCTAAACGAAAAAAACCTAAAGTTCCGGTGCATGCTGTACTGACGGCGACAAATCCGCTGATCCGTTTTATCGGTAGTGATGATATGACGCAAAACCGGGAATTATTTCAGGTCTGGTTACAAAAATTAGCGCAGTGGCATCAGACCACTACGCCTTATCTTTTTTTACATACGCCAGACATTGCCCAGGCACCGGAACTGGTACATACCCTGTGGGAAGACTTACGTAAAACGCTTCCAGAGATCggg                                                                                                                                                                                                                             <  2:345145/288‑2 (MQ=255)
                                                      ttCATGCAGCACGTCCACACAGTGAAGCTATTCGCGACGCTCAACGAAAAAAACCTAAAGTTCCGGTACATGCTGTACTGACGGCGACAAATCCGCTGATCCGTTTTATCGGTAGTGATGATATGACGCAAAACCGGGAATTATTTCAGGTCTGGTTACAAAAATTAGCGCAGTGGCATCAGACCACTACGCCTTATCTTTTTTTACATACGCCAGACATTGCCCAGGCACCGGAACTGGTACATACCCTGTGGGAAGACTTACGTAAAACGCTTCCAGAGATCGGAg                                                                                                                                                                                                                            <  1:304484/288‑1 (MQ=255)
                                                               caCGTCCACACAGTGAAGCTATTCGCGACGCTCAACGAAAAAAACCTAAAGTTCCGGTACATGCTGTACTGACGGCGACAAATCCGCTGATCCGTTTTATCGGTAGTGATGATATGACGCAAAACCGGGAATTATTTCAGGTCTGGTTACAAAAATTAGCGCAGTGGCATCAGACCACTACGCCTTATCTTTTTTTACATACGCCAGACATTGCCCAGGCACCGGAACTGGTACATACCCTGTGGGAAGACTTACGTAAAACGCTTCCAGAGATCGGAGCAGTTCCGGc                                                                                                                                                                                                                  <  2:343892/289‑1 (MQ=255)
                                                               caCGTCCACACAGTGAAGCTATTCGCGACGCTCAACGAAAAAAACCTAAAGTTCCGGTACATGCTGTACTGACGGCGACAAATCCGCTGATCCGCTTTATCGGTAGTGATGATATGACGCAAAACCGGGAATTATTTCAGGTCTGGTTACAAAAATTAGCGCAGTGGCATCAGACCACTACGCCTTATCTTTTTTTACATACGCCAGACATTGCCCAGGCACCGGAACTGGTACATACCCTGTGGGAAGACTTACGTAAAACGCTTCCAGAGATCGGAGCAGTTCCg                                                                                                                                                                                                                    >  1:603278/1‑287 (MQ=255)
                                                                    ccACACAGTGAAGCTATTCGCGACGCTCAACGAAAAAAACCTAAAGTTCCGGTACATGCTGTACTGACGGCGACAAATCCGCTGATCCGTTTTATCGGTAGTGATGATATGACGCAAAACCGGGAATTATTTCAGGTCTGGTTACAAAAATTAGCGCAGTGGCATCAGACCACTACGCCTTATCTTTTTTTACATACGCCAGACATTGCCCAGGCACCGGAACTGGTACATACCCTGTGGGAAGACTTACGTAAAACGCTTCCAGAGATCGGAGCAGTTCCGGCTAtt                                                                                                                                                                                                              >  2:5242/1‑288 (MQ=255)
                                                                              aaGCTATTCGCGACGCTCAACGAAAAAAACCTAAAGTTCCGGTACATGCTGTACTGACGGCGACAAATCCGCTGATCCGTTTTATCGGTAGTGATGATATGACGCAAAACCGGGAATTATTTCAGGTCTGGTTACAAAAATTAGCGCAGTGGCATCAGACCACTACGCCTTATCTTTTTTTACATACGCCAGACATTGCCCAGGCACCGGAACTGGTACATACCCTGTGGGAAGACTTACGTAAAACGCTTCCaga                                                                                                                                                                                                                                    >  1:99944/1‑256 (MQ=255)
                                                                                    ttCGCGACGCTCAACGAAAAAAACCTAAAGTTCCGGTACATGCTGTACTGACGGCGACAAATCCGCTGATCCGTTTTATCGGTAGTGATGATATGACGCAAAACCGGGAATTATTTCAGGTCTGGTTACAAAAATTAGCGCAGTGGCATCAGACCACTACGCCTTATCTTTTTTTACATACGCCAGACATTGCCCAGGCACCGGAACTGGTACATACCCTGTGGGAAGACTTACGTAAAACGCTTTCAGACATAGGAGCAGTTCCGGCTATTCCACAGCAATCctctct                                                                                                                                                                                             >  2:522316/1‑289 (MQ=255)
                                                                                      cgcgACGCTCAACGAAAAAAACCTAAAGTTCCGGTACATGCTGTACTGACGGCGACAAATCCGCTGATCCGTTTTATCGGTAGTGATGATATGACGCAAAACCGGGAATTATTTCAGGTCTGGTTACAAAAATTAGCGCAGTGGCATCAGACCACTACGCCTTATCTTTTTTTACATACGCCAGACATTGCCCAGGCACCGGAACTGGTACATACCCTGTGGGAAGACTTACGTAAAACGCTTCCAGAGATCGGAGCAGtt                                                                                                                                                                                                                       <  2:411017/261‑1 (MQ=255)
                                                                                      cgcgACGCTCAACGAAAAAAACCTAAAGTTCCGGTACATGCTGTACTGACGGCGACAAATCCGCTGATCCGTTTTATCGGTAGTGATGATATGACGCAAAACCGGGAATTATTTCAGGTCTGGTTACAAAAATTAGCGCAGTGGCATCAGACCACTACGCCTTATCTTTTTTTACATACGCCAGACATTGCCCAGGCACCGGAACTGGTACATACCCTGTGGGAAGACTTACGTAAAACGCTTCCAGAGATCGGAGCAGtt                                                                                                                                                                                                                       >  1:411017/1‑261 (MQ=255)
                                                                                                    aaaaaaaCCTAAAGTTCCGGTACATGCTGTACTGACGGCGACAAATCCGCTGATCCGTTTTATCGGTAGTGATGATATGACGCAAAACCGGGAATTATTTCAGGTCTGGTTACAAAAATTAGCGCAGTGGCATCAGACCACTACGCCTTATCTTTTTTTACATACGCCAGACATTGCCCAGGCACCGGAACTGGTACATACCCTGTGGGAAGACTTACGTAAAACGCTTCCAGAGATCGGAGCAGTTCCGGCTATTCCACAGCAATCTTCTCTTTTCTGAATTTGCCAc                                                                                                                                                                             <  1:254516/289‑1 (MQ=255)
                                                                                                    aaaaaaaCCTAAAGTTCCGGTACATGCTGTACTGACGGCGACAAATCCGCTGATCCGTTTTATCGGTAGTGATGATATGACGCAAAACCGGGAATTATTTCAGGTCTGGTTACAAAAATTAGCGCAGTGGCATCAGACCACTACGCCTTATCTTTTTTTACATACGCCAGACATTGCCCAGGCACCGGAACTGGTACATACCCTGTGGGAAGACTTACGTAAAACGCTTCCAGAGATCGGAGCAGTTCCGGCTATTCCACAGCAATCTTCTCTTTTCTGAATTTGCCAc                                                                                                                                                                             <  1:87133/289‑1 (MQ=255)
                                                                                                    aaaaaaaCCTAAAGTTCCGGTACATGCTGTACTGACGGCGACAAATCCGCTGATCCGTGTTATCGGTAGTGATGATATGACGCAAAACCGGGAATTATTTCAGGTCTGGTTACAAAAATTAGCGCAGTGGCATCAGACCACTACGCCTTATCTTTTTTTACATACGCCAGACATTGCCCAGGCACCGGAACTGGTACATACCCTGTGGGAAGACTTACGTAAAACGCTTCCAGAGATCGGAGCAGTTCCGGCTATTCCACAGCAATCTTCTCTTTTCTGAATTTGCCAc                                                                                                                                                                             <  2:408789/289‑1 (MQ=255)
                                                                                                     aaaaaaCCTAAAGTTCCTGTACATGCTGTACTGACGGCGACAAATCCGCTGATCCGTTTTATCGGTAGTGATGCTATGACGCAAAACCGGGAATTATTTCAGGTCTTGTTACAAAAATTAGCGCAGTGGCATCAGACCACTACGCCTTATCTTTTTTTACATACTCCAGACATTGCCCAGGCACCGGAACTGGTACATACCCTGTGGGAAGACTTACGTAAAACGCTTCCAGAGATCGGAGCAGTTCCGGCTATTCCACAGCAATCTTCTCTTTTCTGAATTTGCCAc                                                                                                                                                                             <  1:190780/288‑1 (MQ=255)
                                                                                                                               tGTACTGACGGCGACAAATCCGCTGATCCGTTTTATCGGTAGTGATGATATGACGCAAAACCGGGAATTATTTCAGGTCTGGTTACAAAAATTAGCGCAGTGGCATCAGACCACTACGCCTTATCTTTTTTTACATACGCCAGACATTGCCCAGGCACCGGAACTGGTACATACCCTGTGGGAAGACTTACGTAAAACGCTTCCAGAGATCGGAGCAGTTCCGGCTATTCCACAGCAATCTTCTCTTTTCTGAATTTGCCACCTATCATAGACAGGTGCCATCGGCCAt                                                                                                                                                  >  2:343604/1‑289 (MQ=255)
                                                                                                                                      aCGGCGACAAATCCGCTGATCCGTTTTATCGGTAGTGATGATATGACGCAAAACCGGGAATTATTTCAGGTCTGGTTACAAAAATTAGCGCAGTGGCATCAGACCACTACGCCTTATCTTTTTTTACATACGCCAGACATTGCCCAGGCACCGGAACTGGTACATACCCTGTGGGAAGACTTACGTAAAACGCTTCCAGAGATCGGAGCAGTTCCGGCTATTCCACAGCAAtcttct                                                                                                                                                                                               >  2:571878/1‑237 (MQ=255)
                                                                                                                                      aCGGCGACAAATCCGCTGATCCGTTTTATCGGTAGTGATGATATGACGCAAAACCGGGAATTATTTCAGGTCTGGTTACAAAAATTAGCGCAGTGGCATCAGACCACTACGCCTTATCTTTTTTTACATACGCCAGACATTGCCCAGGCACCGGAACTGGTACATACCCTGTGGGAAGACTTACGTAAAACGCTTCCAGAGATCGGAGCAGTTCCGGCTATTCCACAGCAAtcttct                                                                                                                                                                                               <  1:571878/237‑1 (MQ=255)
                                                                                                                                      aCGGCGACAAATCCGCTGATCCGTTTTATCGGTAGTGATGATATGACGCAAAACCGGGAATTATTTCAGGTCTGGTTACAAAAATTAGCGCAGTGGCATCAGACCACTACGCCTTATCTTTTTTTACATACGCCAGACATTGCCCAGGCACCGGAACTGGTACATACCCTGTGGGAAGACTTACGTAAAACGCTTCCAGAGATCGGAGCAGTTCCGGCTATTCCACAGCAATCTTCTCTTTTCTGAATTTGCCACCTATCATAGACAGGTGCCATCGGCCATTTTAAAg                                                                                                                                           >  1:441583/1‑289 (MQ=255)
                                                                                                                                       cgtcgaacaaTCCGCTTGTCCAGTTGACCGGTTGTGGTGATATGACGCAAAACCGGGAAATATTTGCGGTCTGATTACAAAAATTAGCGCAGTGGCATCAGACCACTACGCCTTATCTGTTTGTACATACGCCAGACATTGCCCAGGCACCGTAACTGGTACATACCCTGGGGGAAGACTTACGTAAAACGCTTCCAGAGATCGGAGCAGTTCCGGCTATTCCACAGCAATCTTCTCTTTTCTGAATTTGCGACCTATCATAGACAGGTGCCATCGGCCATTTTAAAg                                                                                                                                           <  2:186204/280‑1 (MQ=255)
                                                                                                                                         gCGACAAATCCGCTGATCCGTTTTATCGGTAGTGATGATATGACGCAAAACCGGGAATTATTTCAGGTCTGGTTACAAAAATTAGCGCAGTGGCATCAGACCACTACGCCTTATCTTTTTTTACATACGCCAGACATTGCCCAGGCACCGGAACTGGTACATACCCTGTGGGAAGACTTACGTAAAACGCTTCCAGAGATCGGAGCAGTTCCGGCTATTCCACAGCAATCTTCTCTTTTCTGAATTTGCCACCTATCATAGACAGTTGCCATCGGCCATTTTAAAGGGa                                                                                                                                        >  1:584117/1‑289 (MQ=255)
                                                                                                                                           gACAAATCCGCTGATCCGTTTTATCGGTAGTGATGATATGACGCAAAACCGGGAATTATTTCAGGTCTGGTTACAAAAATTAGCGCAGTGGCATCAGACCACTACGCCTTATCTTTTTTTACATACGCCAGACATTGCCCAGGCACCGGAACTGGTACATACCCTGTGGGAAGACTTACGTAAAACGCTTCCAGAGATCGGAGCAGTTCCGGCTATTCCACAGCAATCTTCTCTTTTCTGAATTTGCCACCTATCATAGACAGGTGCCATCGGCCATTTTAAAGGGAg                                                                                                                                       <  2:574903/288‑1 (MQ=255)
                                                                                                                                           gACAAATCCGCTGATCCGTTTTATCGGTAGTGATGATATGACGCAAAACCGGGAATTATTTCAGGTCTGGTTACAAAAATTAGCGCAGTGGCATCAGACCACTACGCCTTATCTTTTTTTACATACGCCAGACATCGCCCAGGCACCGGAACTGGTACATACCCTGTGGGAAGACTTACGTAAAACGCTTCCAGAGATCGGAGCAGTTCCGGCTATTCCACAGCAATCTTCTCTTTTCTGAATTTGCCACCTATCATAGACAGGTGCCATCGGCCATTTTAAAGGGAGt                                                                                                                                      >  2:265163/1‑289 (MQ=255)
                                                                                                                                             cAAATCCGCTGATCCGTTTTATCGGTAGTGATGATATGACGCAAAACCGGGAATTATTTCAGGTCTGGTTACAAAAATTAGCGCAGTGGCATCAGACCACTACGCCTTATCTTTTTTTACATACGCCAGACATTGCCCAGGCACCGGAACTGGTACATACCCTGTGGGAAGACTTACGTAAAACGCTTCCAGAGATCGGAGCAGTTCCGGCTATTCCACAGCACTCTTCTCTTTTCTGAATTTGCCACCTATCATAGACAGGTGCCATCGGCCATTTTAAAGGGAGttt                                                                                                                                    >  1:281407/1‑289 (MQ=255)
                                                                                                                                                      tGATCCGTTTTATCGGTAGTGATGATATGACGCAAAACCGGGAATTATTTCAGGTCTGGTTACAAAAATTAGCGCAGTGGCATCAGACCACTACGCCTTATCTTTTTTTACATACGCCAGACATTGCCCAGGCACCGGAACTGGTACATACCCTGTGGGAAGACTTACGTAAAACGCTTCCAGAGATCGGAGCAGTTCCGGCTATTCCACAGCAATCTTCTCTTTTCTGAATTTGCCACCTATCATAGACAGGTGCCATCGGCCATTTTAAAGGGAGTTTGTATGGTaa                                                                                                                           <  2:441583/289‑1 (MQ=255)
                                                                                                                                                             ttttATCGGTAGTGATGATATGACGCAAAACCGGGAATTATTTCAGGTCTGGTTACAAAAATTAGCGCAGTGGCATCAGACCACTACGCCTTATCTTTTTTTACATACGCCAGACATTGCCCAGGCACCGGAACTGGTACATACCCTGTGGGAAGACTTACGTAAAACGCTTCCAGAGATCGGAg                                                                                                                                                                                                                            <  2:243659/185‑1 (MQ=255)
                                                                                                                                                             ttttATCGGTAGTGATGATATGACGCAAAACCGGGAATTATTTCAGGTCTGGTTACAAAAATTAGCGCAGTGGCATCAGACCACTACGCCTTATCTTTTTTTACATACGCCAGACATTGCCCAGGCACCGGAACTGGTACATACCCTGTGGGAAGACTTACGTAAAACGCTTCCAGAGATCGGAg                                                                                                                                                                                                                            >  1:243659/1‑185 (MQ=255)
                                                                                                                                                               ttATCGGTAGTGATGATATGACGCAAAACCGGGAATTATTTCAGGTCTGGTTACAAAAATTAGCGCAGTGGCATCAGACCACTACGCCTTATCTTTTTTTACATACGCCAGACATTGCCCAGGCACCGGAACTGGTACATACCCTGTGGGAAGACTTACGTAAAACGCTTCCAGAGATCGGAGCAGTTCCGGCTATTCCACAGCAATCTTCTCTTTTCTGAATTTGCCAc                                                                                                                                                                             <  2:600592/230‑1 (MQ=255)
                                                                                                                                                               ttATCGGTAGTGATGATATGACGCAAAACCGGGAATTATTTCAGGTCTGGTTACAAAAATTAGCGCAGTGGCATCAGACCACTACGCCTTATCTTTTTTTACATACGCCAGACATTGCCCAGGCACCGGAACTGGTACATACCCTGTGGGAAGACTTACGTAAAACGCTTCCAGAGATCGGAGCAGTTCCGGCTATTCCACAGCAATCTTCTCTTTTCTGAATTTGCCAc                                                                                                                                                                             >  1:600592/1‑230 (MQ=255)
                                                                                                                                                                   cGGTAGTGATGATATGACGCAAAACCGGGAATTATTTCAGGTCTGGTTACAAAAATTAGCGCAGTGGCATCAGACCACTACGCCTTATCTTTTTTTACATACGCCAGACATTGCCCAGGCACCGGAACTGGTACATACCCTGTGGGAAGACTTACGTAAAACGCTTCCAGAGATCGGAGCAGTTCCGGCTATTCCACAGCAATCTTCTCTTTTCTGAATTTGCCACCTATCATAGACAGGTGCCATCGGCCATTTTAAAGGGAGTTTGTATGGTAAGCGCGCTGTACGc                                                                                                              <  2:621509/289‑1 (MQ=255)
                                                                                                                                                                       aGTGATGATATGACGCAAAACCGGGAATTATTTCAGGTCTGGTTACAAAAATTAGCGCAGTGGCATCAGACCACTACGCCTTATCTTTTTTTACATACGCCAGACATTGCCCAGGCACCGGAACTGGTACATACCCTGTGGGAAGACTTACGTAAAACGCTTCCAGAGATCGGAGCAGTTCCGGCTATTCCACAGCAATCTTCTCTTTTCTGAATTTGCCACCTATCATAGACAGGTGCCATCGGCCATTTTAAAGGGAGTTTGTATGGTAAGCGCGCTGTACGCCGtt                                                                                                          >  1:672083/1‑289 (MQ=255)
                                                                                                                                                                                 tGACGCAAAACCCGGAATTATTTCAGGTCTGGTTACAAAAATTAGCGCAGTTGCATCAGACCACTACGCCTTATCTTTTTTTACATACGCCAGACATTGCCCAGGCACCGGAACTGGTACCTACCCTGTGGGAAGACTTACGTAAAACGCTTCCAGAGATCGGAGCAGTTCCGGCTATTCCACAGCAATCTTCTCTTTTCTGAATTTGCCACCTATCATAGACAGGTGCCATCGGCCATTTTAAAGGGAGTTTGTATGGTAAGCGCGCTGTACGCCGTTTTAAGTGCgt                                                                                                <  1:633784/289‑1 (MQ=255)
                                                                                                                                                                                   aCGCAAAACCGGGAATTATTTCAGGTCTGGTTACAAAAATTAGCGCAGTGGCATCAGACCACTACGCCTTATCTTTTTTTACATACGCCAGACATTGCCCAGGCCCCGGAACTGGTACATACCCTGTGGGCAGACTTACGTAAAACGCTTCCAGGGATCGGAGCAGTTCCGGCTATTCCACAGCAATCTTCCCTTTTCTGAATTTGCCACCTATCATAGACAGCTGCCATCGGCCTTTTTAAAGGGAGTTTGTATTGTAAGCGCGCGGTAGGCGGTTTTAAGTGCgttg                                                                                              >  2:191827/1‑289 (MQ=255)
                                                                                                                                                                                    cGCAAAACCGGGAATTATTTCAGGTCTGGTTACAAAAATTAGCGCAGTGGCATCAGACCACTACGCCTTATCTTTTTTTACATACGCCAGACATTGCCCAGGCACCGGAACGGGTACATACCCTGTGGGAAGACTTACGTAAAACGCTTCCAGAGATCGGAGCAGTTCCGGCTATTCCACAGCAATCTTCTCTTTTCTGAATTTGCCACCTATCATAGACAGGTGACATCGGCCATTTTAAAGGGAGTTTGTATGGTAAGCGAGCTGTACGCCGGTTTAAGTGCgttgt                                                                                             >  1:66308/1‑289 (MQ=255)
                                                                                                                                                                                                                gTTACAAAAATTAGCGCAGTGGCATCAGACCACTACGCCTTATCTTTTTTTACATACGCCAGACATTGCCCAGGCACCGGAACTGGTACATACCCTGTGGGAAGACTTACGTAAAACGCTTCCAGAGATCGGAGCAGTTCCGGC‑‑TTCCACAGCAATCTTCTCTTTTCTGAATTTGCCACCTATCATAGACAGGTGCCATCGGCCATTTTAAAGGGAGTTTGTATGGTAAGCGCGCTGTACGCCGTTTTAAGTGCGTTGTTATTAATGAAGTTCTCTTTTGATGTCGTTc                                                               >  1:118466/1‑289 (MQ=255)
                                                                                                                                                                                                                        aaTTAGCGCAGTGGCATCAGACCACTACGCCTTATCTTTTTTTACATACGCCAGACATTGCCCAGGCACCGGAACTGGTACATACCCTGTGGGAAGACTTACGTAAAACGCTTCCAGAGATCGGAGCAGTTCCGGCTATTCCACAGCAATCTTCTCTTTTCTGAATTTGCCACCTATCATAGACAGGTGCCATCGGCCATTTTAAAGGGAGTTTGTATGGTAAGCGCGCTGTACGCCGTTTTAAGTGCGTTGTTATTAATGAAGTTCTCTTTTGATGTCGTTCGCCTGc                                                         <  1:5242/289‑1 (MQ=255)
                                                                                                                                                                                                                        aaTTAGCACAGTGGCATCAGACCACTACGCCTTATCTTTTTTTACATACGCCAGACATTGCCCAGGCACCGGAACTGGTACATACCCTGTGGGAAGACTTACGTAAAACGCTTCCAGAGATCGGAGCAGTTCCGGCTATTCCACAGCAATCTTCTCTTTTCTGGATTTGCCACCTATCATAGACAGGTGCCATCGGCCATTTTAAAGGGAGTTTGTATGGTAAGCGCGCTGTACGCCGTTTTAAGTGCGTTGTTATTAATGAAGTTCTCTTTTGATGTCGTTCGCCTGc                                                         <  2:603278/289‑1 (MQ=255)
                                                                                                                                                                                                                                   tGGCATCAGACCACTACGCCTTATCTTTTTTTACATACGCCAGACATTGCCCAGGCACCGGAACTGGTACATACCCTGTGGGAAGACTTACGTAAAACGCTTCCAGAGATCGGAGCAGTTCCGGCTATTCCACAGCAATCTTCTCTTTTCTGAATTTGCCACCTATCATAGACAGGTGCCATCGGCCATTTTAAAGGGAGTTTGTATGGTAAGCGCGCTGTACGCCGTTTTAAGTGCGTTGTTATTAATGAAGTTCTCTTTTGATGTCGTTCGCCTGCGAATGCAGTa                                               >  1:536849/1‑288 (MQ=255)
                                                                                                                                                                                                                                   tGGCATCAGACCACTACGCCTTATCTTTTTTTACATACGCCAGACATTGCCCAGGCACCGGAACTGGTACATACCCTGTGGGAAGACTTACGTAAAACGCTTCCAGAGATCGGAGCAGTTCCGGCTATTCCACAGCAATCTTCTCTTTTCTGAATTTGCCACCTATCATAGACAGGTGCCATCGGCCATTTTAAAGGGAGTTTGTATGGTAAGCGCGCTGTACGCCGTTTTAAGTGCGTTGTTATTAATGAAGTTCTCTTTTGATGTCGTTCGCCTGCGAATGCAGTa                                               >  2:347266/1‑288 (MQ=255)
                                                                                                                                                                                                                                   tGGCATCAGACCACTACGCCTTATCTTTTTTTACATACGCCAGACATTGCCCAGGCACCGGAACTGGTACATACCCTGTGGGAAGACTTACGTAAAACGCTTCCAGAGATCGGAGCAGTTCCGGCTATTCCACAGCAATCTTCTCTTTTCTGAATTTGCCACCTATCATAGACAGGTGCCATCGGCCATTTTAAAGGGAGTTTGTATGGTAAGCGCGCTGTACGCCGTTTTAAGTGCGTTGTTATTAATGAAGTTCTCTTTTGATGTCGTTCGCCTGCGAATGCAGTa                                               >  2:650377/1‑288 (MQ=255)
                                                                                                                                                                                                                                    ggCATCAGACCACTACGCCTTATCTTTTTTTACATACGCCAGACATTGCCCAGGCACCGGAACTGGTACATACCCTGTGGGAAGACTTACGTAAAACGCTTCCAGAGATCGGAGCAGTTCCGGCTATTCCACAGCAATCTTCTCTTTTCTGAATTTGCCACCTATCATAGACAGGTGCCATCGGCCATTTTAAAGGGAGTTTGTATGGTAAGCGCGCTGTACGCCGTTTTAAGTGCGTTGTTATTAATGAAGTTCTCTTTTGATGTCGTTCGCCTGCGAATGCAGTa                                               <  1:347266/287‑1 (MQ=255)
                                                                                                                                                                                                                                       aTCAGACCACTACGCCTTATCTTTTTTTACATACGCCAGACATTGCCCAGGCACCGGAACTGGTACATACCCTGTGGGAAGACTTACGTAAAACGCTTCCAGAGATCGGAGCAGTTCCGGCTATTCCACAGCAATCTTCTCTTTTCTGAATTTGCCACCTATCATAGACAGGTGCCATCGGCCATTTTAAAGGGAGTTTGTATGGTAAGCGCGCTGTACGCCGTTTTAAGTGCGTTGTTATTAATGAAGTTCTCTTTTGATGTCGTTCGCCTGCGAATGCAGTACcgca                                          <  1:79779/289‑2 (MQ=255)
                                                                                                                                                                                                                                               aCTACGCCTTATCTTTTTTTACATACGCCAGACATTGCCCAGGCACCGGAACTGGTACATACCCTGTGGGAAGACTTACGTAAAACGCTTCCAGAGATCGGAGCAGTTCCGGCTATTCCACAGCAATCTTCTCTTTTCTGAATTTGCCACCTATCATAGACAGGTGCCATCGGCCATTTTAAAGGGAGTTTGTATGGTAAGCGCGCTGTACGCCGTTTTAAGTGCGTTGTTATTAATGAAGTTCTCTTTTGATGTCGTTCGCCTGCGAATGCAGTACCGCGTTGcc                                     >  2:390233/1‑286 (MQ=255)
                                                                                                                                                                                                                                                cTACGCCTTATCTTTTTTTACATACGCCAGACATTGCCCAGGCACCGGAACTGGTACATACCCTGTGGGAAGACTTACGTAAAACGCTTCCAGAGATCGGAGCAGTTCCGGCTATTCCACAGCAATCTTCTCTTTTCTGAATTTGCCACCTATCATAGACAGGTGCCATCGGCCATTTTAAAGGGAGTTTGTATGGTAAGCGCGCTGTACGCCGTTTTAAGTGCGTTGTTATTAATGAAGTTCTCTTTTGATGTCGTTCGCCTGCGAATGGAGTACGCCATTGCCTATg                                 >  2:79779/1‑289 (MQ=255)
                                                                                                                                                                                                                                                       ttATCTTTTTTTACATACGCCAGACATTGCCCAGGCACCGGAACTGGTACATACCCTGTGGGAAGACTTACGTAAAACGCTTCCAGAGATCGGAGCAGTTCCGGCTATTCCACAGCAATCTTCTCTTTTCTGAATTTGCCACCTATCATAGACAGGTGCCATCGGCCATTTTAAAGGGAGTTTGTATGGTAAGCGCGCTGTACGCCGTTTTAAGTGCGTTGTTATTAATGAAGTTCTCTTTTGATGTCGTTCGCCTGCGAATGCAGTACCGCGTTGCCTATGGCGAcgg                          <  1:522316/289‑1 (MQ=255)
                                                                                                                                                                                                                                                        tATCTTTTTTTACATACGCCAGACATTGCCCAGGCACCGGAACTGGTACATACCCTGTGGGAAGACTTACGTAAAACGCTTCCAGAGATCGGAGCAGTTCCGGCTATTCCACAGCAATCTTCTCTTTTCTGAATTTGCCACCTATCATAGACAGGTGCCATCGGCCATTTTAAAGGGAGTTTGTATGGTAAGCGCGCTGTACGCCGTTTTAAGTGCGTTGTTATTAATGAAGTTCTCTTTTGATGTCGTTCGCCTGCGAATGCAGTACCGCGTTGCCTATGGCGAcggc                         >  1:634876/1‑289 (MQ=255)
                                                                                                                                                                                                                                                         aTCTTTTTTTACATACGCCAGACATTGCCCAGGCACCGGAACTGGTACATACCCTGTGGGAAGACTTACGTAAAACGCTTCCAGAGATCGGAGCAGTTCCGGCTATTCCACAGCAATCTTCTCTTTTCTGAATTTGCCACCTATCATAGACAGGTGCCATCGGCCATTTTAAAGGGAGTTTGTATGGTAAGCGCGCTGTACGCCGTTTTAAGTGCGTTGTTATTAATGAAGTTCTCTTTTGATGTCGTTCGCCTGCGAATGCAGTACAGCGTTGCCTATGGCGAcggcg                        >  2:669034/1‑289 (MQ=255)
                                                                                                                                                                                                                                                               ttttACATACGCCAGACATTGCCCAGGCACCGGAACTGGTACATACCCTGTGGGAAGACTTACGTAAAACGCTTCCAGAGATCGGAGCAGTTCCGGCTATTCCACAGCAATCTTCTCTTTTCTGAATTTGCCAc                                                                                                                                                                             <  1:462488/134‑1 (MQ=255)
                                                                                                                                                                                                                                                               ttttACATACGCCAGACATTGCCCAGGCACCGGAACTGGTACATACCCTGTGGGAAGACTTACGTAAAACGCTTCCAGAGATCGGAGCAGTTCCGGCTATTCCACAGCAATCTTCTCTTTTCTGAATTTGCCAc                                                                                                                                                                             >  2:462488/1‑134 (MQ=255)
                                                                                                                                                                                                                                                                         gCCAGACATTGCCCAGGCACCGGAACTGGTACATACCCTGTGGGAAGACTTACGTAAAACGCTTCCAGAGATCGGAGCAGTTCCGGCTATTCCACAGCAATCTTCTCTTTTCTGAATTTGCCACCTATCATAGACAGGTGCCATCGGCCATTTTAAAGGGAGTTTGTATGGTAAGCGCGCTGTACGCCGTTTTAAGTGCGTTGTTATTAATGAAGTTCTCTTTTGATGTCGTTCGCCTGCGAATGCAGTACCGCGTTGCCTATGGCGACGGCGGTTTTAGCGAACTGCa        <  1:343604/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                               cATTGCCCAGGCACCGGAACTGGTACATACCCTGTGGGAAGACTTACGTAAAACGCTTCCAGAGATCGGAGCAGTTCCGGC‑‑TTCCACAGCAATCTTCTCTTTTCTGAATTTGCCACCTATCATAGACAGGTGCCATCGGCCATTTTAAAGGGAGTTTGTATGGTAAGCGCGCTGTACGCCGTTTTAAGTGCGTTGTTATTAATGAAGTTCTCTTTTGATGTCGTTCGCCTGCGAATGCAGTACCGCGTTGCCTATGGCGACGGCGGTTTTAGCGAACTGCAAAgcgc  <  2:118466/287‑1 (MQ=255)
                                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                                              
ATCGCGGTTTACATCAGCGCGGCGTTAATCGGGTGATTTTAGACAGCCGCCCGGTTCATGCAGCACGTCCACACAGTGAAGCTATTCGCGACGCTCAACGAAAAAAACCTAAAGTTCCGGTACATGCTGTACTGACGGCGACAAATCCGCTGATCCGTTTTATCGGTAGTGATGATATGACGCAAAACCGGGAATTATTTCAGGTCTGGTTACAAAAATTAGCGCAGTGGCATCAGACCACTACGCCTTATCTTTTTTTACATACGCCAGACATTGCCCAGGCACCGGAACTGGTACATACCCTGTGGGAAGACTTACGTAAAACGCTTCCAGAGATCGGAGCAGTTCCGGCTATTCCACAGCAATCTTCTCTTTTCTGAATTTGCCACCTATCATAGACAGGTGCCATCGGCCATTTTAAAGGGAGTTTGTATGGTAAGCGCGCTGTACGCCGTTTTAAGTGCGTTGTTATTAATGAAGTTCTCTTTTGATGTCGTTCGCCTGCGAATGCAGTACCGCGTTGCCTATGGCGACGGCGGTTTTAGCGAACTGCAAAGCGC  >  NC_000913/1951834‑1952393

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑