Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,535,876 T→A 5.2% intergenic (‑18/+61) yddE ← / ← narV PhzC‑PhzF family protein/nitrate reductase 2 (NRZ), gamma subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,535,8760TA 5.2% 470.3 / 12.2 134intergenic (‑18/+61)yddE/narVPhzC‑PhzF family protein/nitrate reductase 2 (NRZ), gamma subunit
Reads supporting (aligned to +/- strand):  ref base T (67/60);  new base A (3/4);  total (70/64)
Fisher's exact test for biased strand distribution p-value = 7.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.34e-01

Only 100 of 135 total aligned reads displayed.

GCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGATATCGACGAGTAAAGTACTCAAACGGCGCGCTCCACACATGCACCAGTCGGGTGAACGGGAAG  >  NC_000913/1535730‑1536024
                                                                                                                                                  |                                                                                                                                                    
gCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTgg                                                                                                                                                    >  3:755186/1‑149 (MQ=255)
  cTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGgag                                                                                                                                                  >  4:346849/1‑149 (MQ=255)
  cTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGgag                                                                                                                                                  >  1:436537/1‑149 (MQ=255)
   tAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCAAAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGgaga                                                                                                                                                 >  2:314015/1‑149 (MQ=255)
      cTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACt                                                                                                                                              <  2:721593/149‑1 (MQ=255)
       tCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTa                                                                                                                                             >  2:808998/1‑149 (MQ=255)
       tCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTa                                                                                                                                             >  1:86754/1‑149 (MQ=255)
         gcgGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATc                                                                                                                                           <  1:1128819/149‑1 (MQ=255)
          cgGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCAGTCACAGGAGAAcgttaa                                                                                                                                          <  4:994692/149‑6 (MQ=255)
          cgGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCAGTCACAGGAGAAcgttaa                                                                                                                                          <  1:991150/149‑6 (MQ=255)
            ggCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATa                                                                                                                                        <  4:1366616/149‑1 (MQ=255)
              cGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGc                                                                                                                                      <  3:457475/149‑1 (MQ=255)
               gATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGcc                                                                                                                                     >  3:203424/1‑149 (MQ=255)
               gATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGcc                                                                                                                                     >  2:1391220/1‑149 (MQ=255)
                     cTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTg                                                                                                                                                     <  3:1423597/127‑1 (MQ=255)
                     cTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGGTGTGAGGTAAAAGCATCGACGTGGTAAACCTGCGGTTTCATTAACGTTCTCCTGTGACTg                                                                                                                                                     >  4:1423597/1‑127 (MQ=255)
                      tGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGt                                                                                                                              <  2:522959/149‑1 (MQ=255)
                       gCATTTGCGCTTCGCTAAGACTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTg                                                                                                                             <  2:12486/149‑1 (MQ=255)
                         aTTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGgag                                                                                                                                                  >  2:179624/1‑126 (MQ=255)
                          tttGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCACCGACGTGGTACACCTGCGGTTTCCTTAACGTTCTCCTGGGACTGGGGGACTATTATAGCCTGCTAGTGGcc                                                                                                                          >  4:74082/1‑149 (MQ=255)
                               gcTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCt                                                                                                                                    <  2:44521/134‑1 (MQ=255)
                               gcTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCt                                                                                                                                    >  1:44521/1‑134 (MQ=255)
                               gcTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGgaga                                                                                                                     <  4:22288/149‑1 (MQ=255)
                               gcTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGgaga                                                                                                                     <  2:1303386/149‑1 (MQ=255)
                                 ttCGCTAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGc                                                                                                                   >  4:1194238/1‑149 (MQ=255)
                                    gCTAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGaa                                                                                                                <  1:322518/149‑1 (MQ=255)
                                      tAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCGTTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAgg                                                                                                              >  4:565010/1‑149 (MQ=255)
                                      tAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAgg                                                                                                              <  1:584262/149‑1 (MQ=255)
                                      tAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAgg                                                                                                              <  4:1063623/149‑1 (MQ=255)
                                      tAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAgg                                                                                                              >  1:829938/1‑149 (MQ=255)
                                                 gcggGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGcc                                                                                                   >  2:808957/1‑149 (MQ=255)
                                                   ggggAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAg                                                                                                 >  3:620685/1‑149 (MQ=255)
                                                    gggAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAgg                                                                                                >  1:1340374/1‑149 (MQ=255)
                                                           caacaCCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAAt                                                                                         >  4:342287/1‑149 (MQ=255)
                                                             acacCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTa                                                                                       >  2:1614540/1‑149 (MQ=255)
                                                                ccGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATcg                                                                                    <  3:346849/149‑1 (MQ=255)
                                                                  ggCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATcgcc                                                                                  >  2:660019/1‑149 (MQ=255)
                                                                        aTTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAAc                                                                            >  4:1170993/1‑149 (MQ=255)
                                                                           gCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGca                                                                         <  4:620685/149‑1 (MQ=255)
                                                                               cgAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAAt                                                                     >  2:1397880/1‑149 (MQ=255)
                                                                                gAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAAtt                                                                    >  2:1243069/1‑149 (MQ=255)
                                                                                 aaaCGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAAttt                                                                   >  4:702714/1‑149 (MQ=255)
                                                                                  aaCGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTg                                                                  >  3:100608/1‑149 (MQ=255)
                                                                                   aCGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGa                                                                 <  1:1668630/149‑1 (MQ=255)
                                                                                   aCGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGa                                                                 <  2:829938/149‑1 (MQ=255)
                                                                                    cGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGat                                                                >  1:540654/1‑149 (MQ=255)
                                                                                            aGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGATATcgacga                                                        <  1:940028/149‑1 (MQ=255)
                                                                                              gTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGATATCGACGAGt                                                      >  2:1637089/1‑149 (MQ=255)
                                                                                              gTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGATATCGACGAGt                                                      >  3:1335842/1‑149 (MQ=255)
                                                                                                 aaaGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGATATCGACGAGTaaa                                                   <  3:237568/149‑1 (MQ=255)
                                                                                                  aaGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGATATCGACGAGTAAAg                                                  >  4:1282011/1‑149 (MQ=255)
                                                                                                  aaGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGATATCGACGAGTAAAg                                                  >  3:1204668/1‑149 (MQ=255)
                                                                                                  aaGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGATATCGACGAGTAAAg                                                  <  1:609747/149‑1 (MQ=255)
                                                                                                   aGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATcgccgc                                                                                <  2:505877/118‑1 (MQ=255)
                                                                                                   aGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGCTGAAATTATcgccgc                                                                                >  1:505877/1‑118 (MQ=255)
                                                                                                    gCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGATATCGACGAGTAAAGTa                                                <  1:1614540/149‑1 (MQ=255)
                                                                                                      aTCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGATATCGACGAGTAAAGTACt                                              <  3:1055022/149‑1 (MQ=255)
                                                                                                          aCGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGATATCGACGAGTAAAGTACTCaaa                                          >  2:600126/1‑149 (MQ=255)
                                                                                                              ggTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGATATCGACGAGTAAAGTACTCAAACGgc                                      <  4:1681921/149‑1 (MQ=255)
                                                                                                              ggTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGATATCGACGAGTAAAGTACTCAAACGgc                                      <  3:65422/149‑1 (MQ=255)
                                                                                                                     cTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGATATCGACGAGTAAAGTACTCAAACGGCGCGCTcc                               <  4:674819/149‑1 (MQ=255)
                                                                                                                     cTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGATATCGACGAGTAAAGTACTCAAACGGCGCGCTcc                               >  3:1246704/1‑149 (MQ=255)
                                                                                                                      tGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGATATCGACGAGTAAAGTACTCAAACGGCGCGCTCca                              >  2:1624856/1‑149 (MQ=255)
                                                                                                                       gCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGATATCGACGAGTaaa                                                   >  4:473061/1‑127 (MQ=255)
                                                                                                                       gCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGATATCGACGAGTaaa                                                   <  3:473061/127‑1 (MQ=255)
                                                                                                                        cGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGATATCGACGAGTAAAGTACTCAAACGGCGCGCTCcaca                            <  1:1397880/149‑1 (MQ=255)
                                                                                                                         ggTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGATATCGACGAGTAAAGTACTCAAACGGCGCGCTCcacac                           >  1:826290/1‑149 (MQ=255)
                                                                                                                           tttCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGATATCGACGAGTAAAGTACTCAAACGGCGCGCTCCACACAt                         <  3:690286/149‑1 (MQ=255)
                                                                                                                             tCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGATATCGACGAGTAAAGTACTCAAACGGCGCGCTCCACACATGc                       <  1:1615927/149‑1 (MQ=255)
                                                                                                                              cATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGATATc                                                             <  1:1015527/110‑1 (MQ=255)
                                                                                                                              cATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGATATc                                                             >  2:1015527/1‑110 (MQ=255)
                                                                                                                               aTTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAg                                                                                                 >  1:356967/1‑73 (MQ=255)
                                                                                                                               aTTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGa                                                                 >  1:75235/1‑105 (MQ=255)
                                                                                                                               aTTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGATATCGACGAGTAAAGTACTCAAACGGCGCGCTCCACACATGCAc                     >  3:1317694/1‑149 (MQ=255)
                                                                                                                                ttAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGATATCGACGAGTAAAGTAc                                               <  4:368418/122‑1 (MQ=255)
                                                                                                                                ttAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGATATCGACGAGTAAAGTAc                                               >  3:368418/1‑122 (MQ=255)
                                                                                                                                ttAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGATATCGACGAGTAAAGTACTCAAACGGCGCGCTCCACACATGCAcc                    >  3:814690/1‑149 (MQ=255)
                                                                                                                                ttAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGATATCGACGAGTAAAGTACTCAAACGGCGCGCTCCACACATGCAcc                    >  1:760917/1‑149 (MQ=255)
                                                                                                                                ttAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGATATCGACGAGTAAAGTACTCAAACGGCGCGCTCCACACATGCAcc                    >  2:1359984/1‑149 (MQ=255)
                                                                                                                                  aaCGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGATATCGACGAGTAAAGTACTCAAACGGCGCGCTCCACACATGCACCAg                  >  2:400011/1‑149 (MQ=255)
                                                                                                                                   aCGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGa                                                                 >  4:495330/1‑101 (MQ=255)
                                                                                                                                   aCGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGa                                                                 <  3:495330/101‑1 (MQ=255)
                                                                                                                                     gTTCTCCAGTCACAGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGATATCGACGAGTAAAGTACTCAAACGGCGCGCTCCACACATGCACCAGTCg               >  1:1023162/1‑149 (MQ=255)
                                                                                                                                      ttCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGATATCGACGAGTAAAGTACTCAAACGGCGCGCTCCACACATGCACCAGTCgg              <  4:100608/149‑1 (MQ=255)
                                                                                                                                      ttCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGATATCGACGAGTAAAGTACTCAAACGGCGCGCTCCACACATGCACCAGTCgg              <  1:1419180/149‑1 (MQ=255)
                                                                                                                                      ttCTCCTGTGACAGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGATATCGACGAGTAAAGTACTCAAACGGCGCGCTCCACACATGCACCAGTCgg              <  4:950781/149‑1 (MQ=255)
                                                                                                                                      ttCTCCAGTCACAGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAggg                                                                                                             <  4:705823/54‑1 (MQ=37)
                                                                                                                                         tcCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGATATCGACGAGTAAAGTACTCAAACGGCGCGCTCCACACATGCACCAGTCGGGTg           <  4:1318077/149‑1 (MQ=255)
                                                                                                                                         tcCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGATATCGACGAGTAAAGTACTCAAACGGCGCGCTCCACACATGCACCAGTCGGGTg           >  2:433090/1‑149 (MQ=255)
                                                                                                                                           cTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGATATCGACGAGTAAAGTACTCAAACGgc                                      <  2:271043/120‑1 (MQ=255)
                                                                                                                                           cTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGATATCGACGAGTAAAGTACTCAAACGgc                                      >  1:271043/1‑120 (MQ=255)
                                                                                                                                           cTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGATATCGACGAGTAAAGTACTCAAACGGCGCGCTCCACACATGCACCAGTCGGGTGaa         <  4:208624/149‑1 (MQ=255)
                                                                                                                                           cTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGATATCGACGAGTAAAGTACTCAAACGGCGCGCTCCACACATGCACCAGTCGGGTGaa         <  3:74082/149‑1 (MQ=255)
                                                                                                                                            tgtgACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTCTCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGATATCGACGAGTAAAGTACTCAAACGGCGCGCTCCACACATGCACCAGTCGGGTGAAc        >  4:1134070/1‑149 (MQ=255)
                                                                                                                                            tgtgACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGATATCGACGAGTAAAGTACTCAAACGGCGCGCTCCACACATGCACCAGTCGGGTGAAc        >  4:460974/1‑149 (MQ=255)
                                                                                                                                             gtgACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGATATCGACGAGTAAAGTACTCAAACGGCGCGCTCCACACATGCACCAGTCGGGTGAACg       <  3:1283103/149‑1 (MQ=255)
                                                                                                                                               gACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGATATCGACGAGTAAAGTACTCAAACGGCGCGCTCCACACATGCACCAGTCGGGTGAACggg     <  4:1204668/149‑1 (MQ=255)
                                                                                                                                                aCTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGATATCGACGAGTAAAGTACTCAAACGGCGCGCTCCACACATGCa                      >  3:995140/1‑131 (MQ=255)
                                                                                                                                                aCTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGATATCGACGAGTAAAGTACTCAAACGGCGCGCTCCACACATGCACCAGTCGGGTGAACGGGa    <  3:1525101/149‑1 (MQ=255)
                                                                                                                                                  tGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGATATCGACGAGTAAAGTACTCAAACGGCGCGCTCCACACATGCACCAGTCGGGTGAACGGGAag  >  3:1588574/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGATATCGACGAGTAAAGTACTCAAACGGCGCGCTCCACACATGCACCAGTCGGGTGAACGGGAAG  >  NC_000913/1535730‑1536024

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