Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 252,594 G→C 5.5% W98C (TGG→TGC)  yafO → mRNA interferase toxin of the YafO‑YafN toxin‑antitoxin system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913252,5940GC 5.5% 644.7 / 22.9 181W98C (TGG→TGC) yafOmRNA interferase toxin of the YafO‑YafN toxin‑antitoxin system
Reads supporting (aligned to +/- strand):  ref base G (85/86);  new base C (6/4);  total (91/90)
Fisher's exact test for biased strand distribution p-value = 7.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.30e-02

Only 100 of 184 total aligned reads displayed.

ACTTATTCGCCTGCAACTTACAGCAGAGGAACTTGATGCGTTAACGGCGGATTTTATTTCCTATAAGCGTGACGGTGTTTTGCCAGATATATTTGGTCGCGATGCACTCTACGACGACTCCTTTACCTGGCCATTAATCAAATTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATACAAATAAGAAACTATCAGCCTGGCGATTTTCAGCAACTATGCGCTATTTTCATTAGAGCGGTTACGATGACCGCCAGTCAGCATTATTCACCACAACAAATTTCCGCCTGGGCGCAGATTGACGAATCTC  >  NC_000913/252321‑252850
                                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                
    attttccTGCAACGTACAGGAGAGTAACTNNTTTCGTTAACGGCGGATTTTATTTCCTATANGCGGTATGGTTTTTTACCAGTTATATTTGGTCGCGATGCAGTCTACGTCGACTCCTTTACCTGGCCATTAATCANATTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCt                                                                                                                                                                                                                                                             <  3:88455/270‑1 (MQ=255)
                                                       aTTTCCTATAAGCGTGACGGTGTTTTGCCAGATATATTTGGTCGCGATGCACTCTACGACGACTCCTTTACCTGGCCATTAATCAAATTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGCACGAATGACTAAGCGCATTTGGTATATTGTCAGGGGGCGTTTTATGCGCAAGCATGGTTTCTCATTGCAATTCTGAAACCTGAACCTCATAAACTGGCTCGCGCTACCACCCa                                                                                                                                                                                                          >  1:79526/1‑275 (MQ=255)
                                                        tttCCTATAAGCGTGACGTTGTTTTGCCAGATATATTTGGTCTCGATGCACTCTACGACGACTCCTTTACCTGGCCATTAATCAAATTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCaa                                                                                                                                                                                                         <  1:136248/275‑1 (MQ=255)
                                                                      gACGGTGTTTTGCCAGATATATTTGGTCGCGATGCACTCTACGACGACTCCTTTACCTGGCCATTAATCAAATTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTgg                                                                                                                                                                                           >  2:463567/1‑275 (MQ=255)
                                                                                    agagATAGTTGGTCGCGATGCACTCTACGACGACTCCTTTACCTTGCCATTAATCAAATTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAg                                                                                                                                                                             <  2:181497/275‑1 (MQ=255)
                                                                                            ttGGTCGCGATGCACTCTACGACGACTCCTTTACCTGGCCATTAATCAAATTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGc                                                                                                                                                                     >  2:393826/1‑275 (MQ=255)
                                                                                                   cgATGCACTCTACGACGACTCCTTTACCTGGCCATTAATCAAATTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCACTTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTTCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTg                                                                                                                                                                                            >  2:630479/1‑245 (MQ=255)
                                                                                                                tatcaaaaaaagaCCTGTCCATTAATCAAATTTGTGCGCGTTTCTCATATTCAACTGGCAAATAAGAATAATCCATNTCCGCCACAATTGAGCCAATTCAGCATAACGAATGACGAATCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTTGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTtata                                                                                                                                                  <  3:160990/262‑1 (MQ=255)
                                                                                                                     cTCCTTTACCTGGCCATTAATCAAATTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCCTAAAATTGGAAAAATGGCAGAAGCGTTTCGAATGCGTTTTTGAATTta                                                                                                                                                    >  4:99477/1‑267 (MQ=255)
                                                                                                                              cTGGCCATTAATCAAATTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGt                                                                                                                                                                                                                                                                 <  8:673041/149‑1 (MQ=255)
                                                                                                                                  ccATTAATCAAATTTGAGCGCGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTCTCAGGGGGCGTTTGATGAGAAANCATGGTTGCTCATTGCCATTTTGAACCCTGAACCTCATAAACTGGCTCGAGAGAAAACACAACTGAAAAAAATTGGGAAANAGTCAAAACCGTTACGAAAAGGTTGNTGAATTTTTATTATAAAAAAcaaacaaaa                                                                                                                               >  3:118673/1‑274 (MQ=255)
                                                                                                                                    aTTAATCAAATTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCa                                                                                                                                                                                                                                                           <  5:379802/149‑1 (MQ=255)
                                                                                                                                    aTTAATCAAATTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCa                                                                                                                                                                                                                                                           >  5:275728/1‑149 (MQ=255)
                                                                                                                                    aTTAATCAAATTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAAttt                                                                                                                                                     <  4:68350/251‑1 (MQ=255)
                                                                                                                                     ttAATCAAATTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCAt                                                                                                                                                                                                                                                          <  4:105682/149‑1 (MQ=255)
                                                                                                                                      tAATCAAATTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCAtt                                                                                                                                                                                                                                                         <  7:383468/149‑1 (MQ=255)
                                                                                                                                      tAATCAAATTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCAtt                                                                                                                                                                                                                                                         <  8:632851/149‑1 (MQ=255)
                                                                                                                                       aaTCAAATTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAAtt                                                                                                                                                      >  1:160990/1‑247 (MQ=255)
                                                                                                                                         tCAAATTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGcc                                                                                                                                                                                                                                                      >  7:1049840/1‑149 (MQ=255)
                                                                                                                                         tCAAATTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGcc                                                                                                                                                                                                                                                      >  8:607941/1‑149 (MQ=255)
                                                                                                                                         tCAAATTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGcc                                                                                                                                                                                                                                                      >  8:343290/1‑149 (MQ=255)
                                                                                                                                          cAAATTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCa                                                                                                                                                                                                                                                     <  5:557162/149‑1 (MQ=255)
                                                                                                                                             aTTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTc                                                                                                                                                                                                                                                  >  7:471073/1‑149 (MQ=255)
                                                                                                                                              tttGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCTCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCt                                                                                                                                                                                                                                                 >  6:692375/1‑149 (MQ=255)
                                                                                                                                              tttGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCGTGGTTGCTCATTGCCATTCt                                                                                                                                                                                                                                                 >  7:887669/1‑149 (MQ=255)
                                                                                                                                              tttGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCt                                                                                                                                                                                                                                                 >  6:896999/1‑149 (MQ=255)
                                                                                                                                                  aGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGaaa                                                                                                                                                                                                                                             <  8:692962/149‑1 (MQ=255)
                                                                                                                                                  aGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGaaa                                                                                                                                                                                                                                             <  8:210303/149‑1 (MQ=255)
                                                                                                                                                        ttGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGaa                                                                                                                                                                                                                                       <  7:316909/149‑1 (MQ=255)
                                                                                                                                                          gCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAAcc                                                                                                                                                                                                                                     >  5:106242/1‑149 (MQ=255)
                                                                                                                                                           cTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCt                                                                                                                                                                                                                                    >  8:215754/1‑149 (MQ=255)
                                                                                                                                                           cTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCt                                                                                                                                                                                                                                    >  5:906133/1‑149 (MQ=255)
                                                                                                                                                            tCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAACCATGCTTGCTCa                                                                                                                                                                                                                                                           >  7:748281/1‑125 (MQ=255)
                                                                                                                                                             cATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCa                                                                                                                                                                                                                                  >  6:234954/1‑149 (MQ=255)
                                                                                                                                                                   cATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACt                                                                                                                                                                                                                            <  5:502068/149‑1 (MQ=255)
                                                                                                                                                                    aTCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTg                                                                                                                                                                                                                           <  5:1044474/149‑1 (MQ=255)
                                                                                                                                                                     tCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTgg                                                                                                                                                                                                                          <  5:896999/149‑1 (MQ=255)
                                                                                                                                                                        ggCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTc                                                                                                                                                                                                                       <  8:155763/149‑1 (MQ=255)
                                                                                                                                                                         gCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTAtattat                                                                                                                                              >  4:43756/1‑221 (MQ=255)
                                                                                                                                                                            aaTGTGAATATACCATTTCCGCCACAGTTTCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATACAAATAAGAAACTATCAGCCTGGCGATTTTCAGCAACTATGCGCTAt                                                                                     <  4:463567/275‑1 (MQ=255)
                                                                                                                                                                            aaTGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCgaga                                                                                                                                                                                                                   <  5:365901/149‑1 (MQ=255)
                                                                                                                                                                             aTGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGAt                                                                                                                                                                                                                  >  6:1024331/1‑149 (MQ=255)
                                                                                                                                                                                 gaATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAACCATGCTTGCTCa                                                                                                                                                                                                                                                           <  5:635717/104‑1 (MQ=255)
                                                                                                                                                                                  aataatCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATaacaa                                                                                                                                                                                                             >  7:562685/1‑149 (MQ=255)
                                                                                                                                                                                      atCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTc                                                                                                                                                                                                                                                  <  7:726564/108‑1 (MQ=255)
                                                                                                                                                                                      atCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCaa                                                                                                                                                                                                         >  8:872554/1‑149 (MQ=255)
                                                                                                                                                                                       tCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTg                                                                                                                                                                                                                                                               >  8:118114/1‑94 (MQ=255)
                                                                                                                                                                                       tCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTg                                                                                                                                                                                                                                                               <  7:118114/94‑1 (MQ=255)
                                                                                                                                                                                       tCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCaaa                                                                                                                                                                                                        >  6:107688/1‑149 (MQ=255)
                                                                                                                                                                                                gCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATaaaa                                                                                                                                                                                               >  6:855581/1‑149 (MQ=255)
                                                                                                                                                                                                 ccACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATCACAACCAAATGCATAAAAt                                                                                                                                                                                              >  8:736432/1‑149 (MQ=255)
                                                                                                                                                                                                  cacaGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAAtt                                                                                                                                                                                             <  5:692375/149‑1 (MQ=255)
                                                                                                                                                                                                   acaGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTg                                                                                                                                                                                            >  6:616635/1‑149 (MQ=255)
                                                                                                                                                                                                   acaGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTg                                                                                                                                                                                            <  8:558502/149‑1 (MQ=255)
                                                                                                                                                                                                         gcgcCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGaaaa                                                                                                                                                                                      <  7:215754/149‑1 (MQ=255)
                                                                                                                                                                                                          cgcCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAAt                                                                                                                                                                                     <  7:447468/149‑1 (MQ=255)
                                                                                                                                                                                                          cgcCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAAt                                                                                                                                                                                     <  8:565236/149‑1 (MQ=255)
                                                                                                                                                                                                            ccAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATgg                                                                                                                                                                                   <  8:887669/149‑1 (MQ=255)
                                                                                                                                                                                                            ccAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATgg                                                                                                                                                                                   <  7:599651/149‑1 (MQ=255)
                                                                                                                                                                                                             cAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGc                                                                                                                                                                                  <  8:761164/149‑1 (MQ=255)
                                                                                                                                                                                                               aTTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTAtattat                                                                                                                                              >  1:100026/1‑183 (MQ=255)
                                                                                                                                                                                                                   agcagAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAAt                                                                                                                                          <  3:100026/183‑1 (MQ=255)
                                                                                                                                                                                                                    gcagAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTc                                                                                                                                                                                                                                   <  7:930396/93‑1 (MQ=255)
                                                                                                                                                                                                                          cGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGc                                                                                                                                                                                                                         >  4:653316/1‑97 (MQ=255)
                                                                                                                                                                                                                            aaTGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCAt                                                                                                                                                                   >  7:261271/1‑149 (MQ=255)
                                                                                                                                                                                                                              tGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCga                                                                                                                                                                                                                     <  2:653316/97‑1 (MQ=255)
                                                                                                                                                                                                                                  gAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGttt                                                                                                                                                             >  8:944100/1‑149 (MQ=255)
                                                                                                                                                                                                                                   aaGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGtttt                                                                                                                                                            <  7:607941/149‑1 (MQ=255)
                                                                                                                                                                                                                                     gcgcATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGc                                                                                                                                                                                                                         <  6:764460/86‑1 (MQ=255)
                                                                                                                                                                                                                                        cATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAAt                                                                                                                                                       <  6:224281/149‑1 (MQ=255)
                                                                                                                                                                                                                                          tttGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTGATCTCCATTCTGCAACCTGACCCCCATATCCTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAACATGGCAGAAGCGTCTCGCACGCGTTTTTCAAttt                                                                                                                                                     <  8:377934/149‑1 (MQ=255)
                                                                                                                                                                                                                                           ttGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTta                                                                                                                                                    <  8:1035493/149‑1 (MQ=255)
                                                                                                                                                                                                                                              gTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTatat                                                                                                                                                 >  6:937358/1‑149 (MQ=255)
                                                                                                                                                                                                                                               tataTTGTCAGGGGACGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTAtatt                                                                                                                                                >  8:237314/1‑149 (MQ=255)
                                                                                                                                                                                                                                                  atTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATg                                                                                                                                             >  6:974697/1‑149 (MQ=255)
                                                                                                                                                                                                                                                   ttGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAAc                                                                                                                                       >  3:171307/1‑154 (MQ=255)
                                                                                                                                                                                                                                                     gTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAAt                                                                                                                                          <  8:43692/149‑1 (MQ=255)
                                                                                                                                                                                                                                                      tCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTggg                                                                                                                                                                                          >  8:1042777/1‑100 (MQ=255)
                                                                                                                                                                                                                                                      tCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATa                                                                                                                                         >  6:642340/1‑149 (MQ=255)
                                                                                                                                                                                                                                                       cAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATAc                                                                                                                                   <  1:171307/154‑1 (MQ=255)
                                                                                                                                                                                                                                                        aGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCg                                                                                                                                                                      >  4:752736/1‑118 (MQ=255)
                                                                                                                                                                                                                                                         gggggCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGaa                                                                                                                                                        >  6:611129/1‑131 (MQ=255)
                                                                                                                                                                                                                                                         gggggCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGaa                                                                                                                                                        <  5:611129/131‑1 (MQ=255)
                                                                                                                                                                                                                                                          ggggCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACAt                                                                                                                                     >  7:491307/1‑149 (MQ=255)
                                                                                                                                                                                                                                                          ggggCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACAt                                                                                                                                     >  7:335963/1‑149 (MQ=255)
                                                                                                                                                                                                                                                           gggCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATa                                                                                                                                    <  8:976507/149‑1 (MQ=255)
                                                                                                                                                                                                                                                           gggCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATa                                                                                                                                    <  7:872554/149‑1 (MQ=255)
                                                                                                                                                                                                                                                            ggCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATg                                                                                                                                                                  <  2:752736/118‑1 (MQ=255)
                                                                                                                                                                                                                                                            ggCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAAttt                                                                                                                                                     <  2:715415/131‑1 (MQ=255)
                                                                                                                                                                                                                                                            ggCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATAc                                                                                                                                   >  5:972633/1‑149 (MQ=255)
                                                                                                                                                                                                                                                               tggcaatgaGCAACCATGCTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATACaaa                                                                                                                                >  5:76981/6‑149 (MQ=255)
                                                                                                                                                                                                                                                               gTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATACaaa                                                                                                                                >  6:523800/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                   gatgaGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATACAAATAAg                                                                                                                            >  8:577037/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                   gatgaGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATACAAATAAg                                                                                                                            >  5:642301/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                   aatgaGCAACCATGCTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATaaaa                                                                                                                                                                                               <  7:418715/81‑1 (MQ=255)
                                                                                                                                                                                                                                                                    atgaGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTatat                                                                                                                                                 <  1:111536/127‑1 (MQ=255)
                                                                                                                                                                                                                                                                    atgaGCAACCATGCTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATACAAATAAGa                                                                                                                           >  6:40312/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                      gaGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAAtt                                                                                                                                                                                             <  8:286377/81‑1 (MQ=255)
                                                                                                                                                                                                                                                                           agcaTGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATACAAATAAGAAACTATc                                                                                                                    <  8:731561/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                                 gTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATACAAATAAGAAACTATCAGCCTGGCGATTTTCAGCAACTa                                                                                             >  2:191923/1‑166 (MQ=255)
                                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                
ACTTATTCGCCTGCAACTTACAGCAGAGGAACTTGATGCGTTAACGGCGGATTTTATTTCCTATAAGCGTGACGGTGTTTTGCCAGATATATTTGGTCGCGATGCACTCTACGACGACTCCTTTACCTGGCCATTAATCAAATTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATACAAATAAGAAACTATCAGCCTGGCGATTTTCAGCAACTATGCGCTATTTTCATTAGAGCGGTTACGATGACCGCCAGTCAGCATTATTCACCACAACAAATTTCCGCCTGGGCGCAGATTGACGAATCTC  >  NC_000913/252321‑252850

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