Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 480,237 C→G 7.1% L24L (CTG→CTC hha ← modulator of gene expression, with H‑NS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913480,2370CG 7.1% 297.9 / 12.2 91L24L (CTG→CTChhamodulator of gene expression, with H‑NS
Reads supporting (aligned to +/- strand):  ref base C (43/40);  new base G (3/4);  total (47/44)
Fisher's exact test for biased strand distribution p-value = 7.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.25e-01

CACATAGGGTATGTCTTCATGGCGAAAAAGTATAAAATTCTTAATAAACAGCCGGTTATAGCTCCGAAAGCGAATTATTTATTAGCGAATAAATTTCCATACTGAGGAAGGGATCTTGTCGTACAGTTTATTCATGGTCAATTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCCTGCATATTAATACCATAACTACTGAACAACATAAAGGTGTCATCCAGATATTCGTCGAT  >  NC_000913/480009‑480483
                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                      
cacaTAGGGTATGTCTTCATGGCGAAAAAGTATAAAATTCTTAATAAACAGCCGGTTATAGCTCCGAAAGCGAATTATTTATTAGCGAATAAATTTCCATACTGAGGAAGGGATCTTGTCGTACAGTTTATTCATGGTCAATTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCgg                                                                                                                                                                                                          >  2:404343/1‑275 (MQ=255)
           tGTCTTCATGGCGAAAAAGTATAAAATTCTTAATAAACAGCCGGTTATAGCTCCGAAAGCGAATTATTTATTAGCGAATAAATTTCCATACTGAGGAAGGGGTCTTGTCGTACAGTTTATTCATGGTCAATTCGGCGANGCGGTGATCTGCGGCTGAGGAACATACCGACAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAAACCNCTCACGCGTGTAAATTTTCTTTCCAAGACTNNAACCCATTAATCAACAGGTTTTTGTTtaa                                                                                                                                                                                               >  3:75686/1‑275 (MQ=255)
                                                      gTTATAGCTCCGAAAGCGAATTATTTATTAGCGAATAAATTTCCATACTGAGGAAGGGATCTTGTCGTACAGTTTATTCATGGTCAATTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTaa                                                                                                                                                    <  2:339739/275‑1 (MQ=255)
                                                                                     cGAATAAATTTCCATACTGAGGAAGGGATCTTGTCGTACAGTTTATTCATGGTCAATTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTNTCGATAACACGCTCCATCGTGTCAATTGTTTGTACACGACGTAAACGCATTAAATATCTGGTGTTCTTTNAAGGGTTTTCGGTCATACTTCTACCTATGGTTGTAATAATTCTAANANATTAATATCTCAGGATTCTCTTtctttt                                                                                                                     >  3:40869/1‑271 (MQ=255)
                                                                                        aTAAATTTCCATACTGAGGAAGGGATCTTGTCGTACAGTTTATTCATGGTCAATTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCgtgt                                                                                                                                                                                                                                                 >  5:983742/1‑148 (MQ=255)
                                                                                                      tGAGGAAGGGATCTTGTCGTACAGTTTATTCATGGTCAATTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTGGTAATAATTCAAACAAGATAACCTCGCAGGATTCCCTTTAGTCGCATTAACAAAACAAAgg                                                                                                    >  4:114679/1‑275 (MQ=255)
                                                                                                        aGGAAGGGATCTTGTCGTACAGTTTATTCATGGTCAATTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGaa                                                                                                  >  4:285017/1‑275 (MQ=255)
                                                                                                        aGGAAGGGATCTTGTCGTAAAGTTTAGTCATGGTCAATTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGaa                                                                                                  <  1:75686/275‑1 (MQ=255)
                                                                                                             gggATCTTGTCGTACAGTTTATTCATGGTCAATTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTa                                                                                                                                                                                                                            >  5:341255/1‑148 (MQ=255)
                                                                                                                       cGTACAGTTTATTCATGGTCAATTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTaa                                                                                                                                                                                                                  <  6:341255/148‑1 (MQ=255)
                                                                                                                                tATTCATGGTCAATTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCg                                                                             >  4:299597/1‑272 (MQ=255)
                                                                                                                                   tCATGGTCAATTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCg                                                                                                                                                                                                                                                    <  5:943482/102‑1 (MQ=255)
                                                                                                                                   tCATGGTCAATTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCg                                                                                                                                                                                                                                                    >  6:943482/1‑102 (MQ=255)
                                                                                                                                      tGGTCAATTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATATCACGCTCCAGCGTGTCAATTGTc                                                                                                                                                                                                                                         <  6:633695/110‑1 (MQ=255)
                                                                                                                                      tGGTCAATTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAGCACGCTCCAGCGTGTCAATTGTc                                                                                                                                                                                                                                         >  5:633695/1‑110 (MQ=255)
                                                                                                                                       ggTCATTTCTGCGAGGCGGTGTTCTGCGTCTGAGTAAAATACCGCCAGTTCATTATCTGAAAATTCGTATTTATTTTTCTCGATAACACGCTCCATCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAAGAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGa                                                                   <  4:404342/275‑1 (MQ=255)
                                                                                                                                        gTCAATTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAaataat                                                                                                                                                                                                             >  6:952309/1‑136 (MQ=255)
                                                                                                                                        gTCAATTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAaataat                                                                                                                                                                                                             <  5:952309/136‑1 (MQ=255)
                                                                                                                                        gTCAATTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTa                                                                                                                                                                      >  2:16445/1‑175 (MQ=255)
                                                                                                                                        gTCAATTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGaa                                                                  >  2:5180/1‑275 (MQ=255)
                                                                                                                                            aTTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAacga                                                                                                                                                                                                                                <  6:1025346/113‑1 (MQ=255)
                                                                                                                                            aTTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAacga                                                                                                                                                                                                                                >  5:1025346/1‑113 (MQ=255)
                                                                                                                                            aTTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGt                                                                                                                                                                  <  4:16445/175‑1 (MQ=255)
                                                                                                                                            aTTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGAt                                                               >  2:187568/1‑274 (MQ=255)
                                                                                                                                             ttCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAacg                                                                                                                                                                                                                                 >  6:634993/1‑111 (MQ=255)
                                                                                                                                             ttCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAacg                                                                                                                                                                                                                                 <  5:634993/111‑1 (MQ=255)
                                                                                                                                               cGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCCTg                                                           <  2:285018/275‑1 (MQ=255)
                                                                                                                                                  cGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCCCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCg                                                                                                                                                                                                           >  3:12808/1‑128 (MQ=255)
                                                                                                                                                     ggCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCCTGCATAtt                                                     <  4:187567/275‑1 (MQ=255)
                                                                                                                                                      gCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACACTCCAGCGTGTCAATTGTCTGGCAACGGCGTAAACGCATTAAATAATCGGttt                                                                                                                                                                                                       <  1:12808/128‑1 (MQ=255)
                                                                                                                                                        ggTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCCTGCATATTaa                                                   <  2:299598/274‑1 (MQ=255)
                                                                                                                                                             tCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTc                                                                                                                                                                            <  6:983742/148‑1 (MQ=255)
                                                                                                                                                                       aGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCCTGCATATTAATACCATAactac                                       >  4:233296/1‑271 (MQ=255)
                                                                                                                                                                        gTAAAATACCGCCAGTTCATTATCTTATAATTCGTTGTTATTTTTCTCGATAACACGCTCCAGCGGGTCAATTGGCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCCTGCATATTAATACCATAACTACTGaa                                   <  4:5180/274‑1 (MQ=255)
                                                                                                                                                                        gTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAaataa                                                                                                                                                                                                              <  7:1098937/103‑1 (MQ=255)
                                                                                                                                                                        gTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAaataa                                                                                                                                                                                                              >  8:1098937/1‑103 (MQ=255)
                                                                                                                                                                          aaaaTACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGt                                                                                                                                                               >  8:254449/1‑148 (MQ=255)
                                                                                                                                                                           aaaTACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCCTGCATATTAATACCATAACTACTGaa                                   <  2:233297/271‑1 (MQ=255)
                                                                                                                                                                              tACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTGGAGCgtgt                                                                                                                                                                                                                                                 <  5:340574/62‑1 (MQ=255)
                                                                                                                                                                              tACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTGGAGCgtgt                                                                                                                                                                                                                                                 >  6:340574/1‑62 (MQ=255)
                                                                                                                                                                                ccGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATg                                                                                                                                                                    <  6:52137/137‑1 (MQ=255)
                                                                                                                                                                                ccGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATg                                                                                                                                                                    >  5:52137/1‑137 (MQ=255)
                                                                                                                                                                                 cGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCgtg                                                                                                                                                                                                                                                  >  2:412458/1‑58 (MQ=255)
                                                                                                                                                                                 cGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTtcg                                                                                                                        >  1:137225/1‑180 (MQ=255)
                                                                                                                                                                                 cGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAg                                                                 >  2:379242/1‑235 (MQ=255)
                                                                                                                                                                                 cGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGACAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGaa                                                                  >  2:105179/1‑234 (MQ=255)
                                                                                                                                                                                   ccAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCg                                                                                                                                                                                      >  3:46434/1‑116 (MQ=255)
                                                                                                                                                                                   ccAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGTAATAATTc                                                                                                                                                       >  4:558957/1‑147 (MQ=255)
                                                                                                                                                                                   ccAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTg                                                                    >  4:225439/1‑230 (MQ=255)
                                                                                                                                                                                     aGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCaa                                                                                                                                                                                                                                              <  4:412457/58‑1 (MQ=255)
                                                                                                                                                                                     aGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGc                                                                                                                    <  3:137225/180‑1 (MQ=255)
                                                                                                                                                                                     aGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATcc                                                             <  4:379241/235‑1 (MQ=255)
                                                                                                                                                                                     aGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATc                                                              <  4:105178/234‑1 (MQ=255)
                                                                                                                                                                                       ttCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACa                                                                                                                                                                                  <  1:46434/116‑1 (MQ=255)
                                                                                                                                                                                       ttCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTa                                                                                                                                                                      >  4:327004/1‑128 (MQ=255)
                                                                                                                                                                                       ttCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAAc                                                                                                                                                   <  2:558958/147‑1 (MQ=255)
                                                                                                                                                                                       ttCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGa                                                                <  2:225440/230‑1 (MQ=255)
                                                                                                                                                                                        tCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTgg                                                                                                                                                                                                                                      <  8:159215/63‑1 (MQ=255)
                                                                                                                                                                                        tCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTgg                                                                                                                                                                                                                                      >  7:159215/1‑63 (MQ=255)
                                                                                                                                                                                        tCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCgg                                                                                                                                                                                     <  5:581418/112‑1 (MQ=255)
                                                                                                                                                                                        tCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCgg                                                                                                                                                                                     >  6:581418/1‑112 (MQ=255)
                                                                                                                                                                                         cATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCCTGCa                                                         >  4:428228/1‑235 (MQ=255)
                                                                                                                                                                                           ttATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGt                                                                                                                                                                  <  2:327005/128‑1 (MQ=255)
                                                                                                                                                                                             aTCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCCTGCATAtt                                                     <  2:428229/235‑1 (MQ=255)
                                                                                                                                                                                              tCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCCTGCATATTAATACCATAACTACTGAACAACATAAAGGTGTCATCCAGa             >  1:84813/1‑274 (MQ=255)
                                                                                                                                                                                               cTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAAc                                                                                                                                          >  6:841870/1‑148 (MQ=255)
                                                                                                                                                                                                 gATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCAttt                                                                       >  2:518721/1‑213 (MQ=255)
                                                                                                                                                                                                     aTTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGa                                                                   <  4:518720/213‑1 (MQ=255)
                                                                                                                                                                                                      ttCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTaa                                                                                                                                                                                                <  5:701649/87‑1 (MQ=255)
                                                                                                                                                                                                      ttCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTaa                                                                                                                                                                                                >  6:701649/1‑87 (MQ=255)
                                                                                                                                                                                                      ttCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGa                                                                >  2:204153/1‑215 (MQ=255)
                                                                                                                                                                                                      ttCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACATATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTTGATTCACTGACTTCCGCCTTTTGTGAAGATCCTGCATATTAATACCATAATTACTGAACAAAATAACGGTGTCATCCAGATAttcttc     >  1:158870/1‑274 (MQ=255)
                                                                                                                                                                                                       tCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGa                                                                                                                                  <  7:254449/148‑1 (MQ=255)
                                                                                                                                                                                                        tctatttattTTTCTTGATCACACGGTCGAGAGTGTTAATTGGCTGGAAAAGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTATTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGAGAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCCTGCATATTAATACCATAACTACTGAACAACATAAAGGTGTCATCCAGATATTCGTCGAt  <  4:96744/273‑1 (MQ=255)
                                                                                                                                                                                                        cGTTTTTTTTTTTTTCGATAACACGCTCCAGCGGGTCAATGGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCCTGCATATTAATACCATAACTACTGAACAACATAAAGGTGTCATCCAGATATTCGTCGAt  <  4:174566/275‑1 (MQ=255)
                                                                                                                                                                                                          tatttattTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATtc                                                                                                                               >  8:135500/1‑148 (MQ=255)
                                                                                                                                                                                                          tatttattTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCCt                                                            <  4:204152/215‑1 (MQ=255)
                                                                                                                                                                                                               atttTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACTTAAACGCATTAACTAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCTTCGCATTGACAATACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGATGATCCTGCATATTAATACCATAACTACTTAACAACATAAAGGTGTCATCCAGGTATtcgt      >  2:96745/1‑264 (MQ=255)
                                                                                                                                                                                                                    tctcGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGAt                                                               >  1:96453/1‑202 (MQ=255)
                                                                                                                                                                                                                        gATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCCTg                                                           <  3:96453/202‑1 (MQ=255)
                                                                                                                                                                                                                         aTAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTa                                                                                                                                                     >  8:671387/1‑111 (MQ=255)
                                                                                                                                                                                                                         aTAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTa                                                                                                                                                     <  7:671387/111‑1 (MQ=255)
                                                                                                                                                                                                                          tgacacGCTGGAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTaa                                                                                                                                                                                                >  5:1033114/3‑67 (MQ=255)
                                                                                                                                                                                                                          tgacacGCTGGAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTaa                                                                                                                                                                                                <  6:1033114/65‑1 (MQ=255)
                                                                                                                                                                                                                             cacGCTGGAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGtaata                                                                                                                                                           >  5:380908/1‑101 (MQ=255)
                                                                                                                                                                                                                             cacGCTGGAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGtaata                                                                                                                                                           <  6:380908/101‑1 (MQ=255)
                                                                                                                                                                                                                              acGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGttttt                                                                                                                                                                                        >  6:751641/1‑71 (MQ=255)
                                                                                                                                                                                                                              acGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGttttt                                                                                                                                                                                        <  5:751641/71‑1 (MQ=255)
                                                                                                                                                                                                                              acGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACaaa                                                                                                           <  7:135500/148‑1 (MQ=255)
                                                                                                                                                                                                                               cGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATa                                                                                                                                                                                >  7:290342/1‑78 (MQ=255)
                                                                                                                                                                                                                               cGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATa                                                                                                                                                                                <  8:290342/78‑1 (MQ=255)
                                                                                                                                                                                                                                   ccAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTa                                                                                                                                                                          <  5:347309/80‑1 (MQ=255)
                                                                                                                                                                                                                                   ccAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTa                                                                                                                                                                          >  6:347309/1‑80 (MQ=255)
                                                                                                                                                                                                                                   ccAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAaacaa                                                                                                       >  6:813295/1‑147 (MQ=255)
                                                                                                                                                                                                                                   ccAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAaacaa                                                                                                       <  5:813295/147‑1 (MQ=255)
                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                      
CACATAGGGTATGTCTTCATGGCGAAAAAGTATAAAATTCTTAATAAACAGCCGGTTATAGCTCCGAAAGCGAATTATTTATTAGCGAATAAATTTCCATACTGAGGAAGGGATCTTGTCGTACAGTTTATTCATGGTCAATTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCCTGCATATTAATACCATAACTACTGAACAACATAAAGGTGTCATCCAGATATTCGTCGAT  >  NC_000913/480009‑480483

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