Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 480,236 C→G 6.6% E25Q (GAG→CAG)  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,2360CG 6.6% 305.4 / 11.8 91E25Q (GAG→CAG) hhamodulator of gene expression, with H‑NS
Reads supporting (aligned to +/- strand):  ref base C (44/41);  new base G (3/3);  total (47/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.62e-02

AAAAGTCAGTGAAGAGAACGTTCCCATTACATCTGTGTGATAAAGATCACATAGGGTATGTCTTCATGGCGAAAAAGTATAAAATTCTTAATAAACAGCCGGTTATAGCTCCGAAAGCGAATTATTTATTAGCGAATAAATTTCCATACTGAGGAAGGGATCTTGTCGTACAGTTTATTCATGGTCAATTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCCTGCATATTAATACCATAACTACTGAACAACATAAAGGTGTCATCCAGATATTCGTCGAT  >  NC_000913/479962‑480483
                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                       
aaaaGTCAGTGAAGAGAACGTTCCCATTACATCTGTGTGATAAAGATCACATAGGGTATGTCTTCATGGCGAAAAAGTATAAAATTCTTAATAAACAGCCGGTTATAGCTCCGAAAGCGAATTATTTATTAGCGAATAAATTTCCATACTGAGGAAGGGATCTTGTCGTACAGTTTATTCATGGTCAATTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTc                                                                                                                                                                                                                                                         <  4:182959/275‑1 (MQ=255)
                                               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                                                                                                                                                                                                                                 <  5:634993/111‑1 (MQ=255)
                                                                                                                                                                                            ttCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAacg                                                                                                                                                                                                                                 >  6:634993/1‑111 (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)
                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                       
AAAAGTCAGTGAAGAGAACGTTCCCATTACATCTGTGTGATAAAGATCACATAGGGTATGTCTTCATGGCGAAAAAGTATAAAATTCTTAATAAACAGCCGGTTATAGCTCCGAAAGCGAATTATTTATTAGCGAATAAATTTCCATACTGAGGAAGGGATCTTGTCGTACAGTTTATTCATGGTCAATTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGTAATAATTCTAACAAGATAAACTCGCAGGATTCTCTTTCGTCGCATTGACAAAACAACGGAATAGTCGATTACCTGACTTCCGCCATTTCTGAAGATCCTGCATATTAATACCATAACTACTGAACAACATAAAGGTGTCATCCAGATATTCGTCGAT  >  NC_000913/479962‑480483

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