Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,866,488 G→C 6.9% intergenic (+23/+71) ygbN → / ← rpoS putative transporter/RNA polymerase, sigma S (sigma 38) factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,866,4880GC 6.9% 502.6 / 23.9 145intergenic (+23/+71)ygbN/rpoSputative transporter/RNA polymerase, sigma S (sigma 38) factor
Reads supporting (aligned to +/- strand):  ref base G (72/63);  new base C (5/5);  total (77/68)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.47e-01

Only 100 of 149 total aligned reads displayed.

TTCTGGATTGTGACCAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCT  >  NC_000913/2866340‑2866635
                                                                                                                                                    |                                                                                                                                                    
ttCTGGATTGTGACCAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATggg                                                                                                                                                      >  1:1581618/1‑149 (MQ=255)
    ggATTGTGACCAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATc                                                                                                                                                  >  3:234082/1‑149 (MQ=255)
    ggATTGTGACCAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATc                                                                                                                                                  >  3:221754/1‑149 (MQ=255)
     gATTGTGACCAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCg                                                                                                                                                 >  2:1302327/1‑149 (MQ=255)
     gATTGTGACCAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGCCCATCg                                                                                                                                                 >  2:647205/1‑149 (MQ=255)
      aTTGTGACCAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCgg                                                                                                                                                <  1:430065/149‑1 (MQ=255)
       ttGTGACCAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGa                                                                                                                                               <  2:1034269/149‑1 (MQ=255)
        tgtgACCAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGCCCATCGGAc                                                                                                                                              <  2:1329646/149‑1 (MQ=255)
        tgtgACCAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGCCCATCGGAc                                                                                                                                              <  1:647205/149‑1 (MQ=255)
         gtgACCAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGAcc                                                                                                                                             <  3:102306/149‑1 (MQ=255)
          tgACCAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑t                                                                                                                                           <  2:1192443/149‑1 (MQ=255)
            aCCAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGATAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑ttt                                                                                                                                         <  2:27175/149‑1 (MQ=255)
            aCCAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCTGTTTTACCGGATTCTTAATTACCCGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGCTGGGCATCGGACC‑ttt                                                                                                                                         <  2:411898/149‑1 (MQ=255)
            aCCAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑ttt                                                                                                                                         <  2:560805/149‑1 (MQ=255)
            aCCAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑ttt                                                                                                                                         >  1:606322/1‑149 (MQ=255)
             ccAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑tttt                                                                                                                                        >  4:204412/1‑149 (MQ=255)
              cAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATCAAAAGGTCCGATGGGCATCGGACC‑TTTTa                                                                                                                                       <  4:42408/149‑1 (MQ=255)
              cAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTa                                                                                                                                       <  2:705324/149‑1 (MQ=255)
              cAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTa                                                                                                                                       >  3:730363/1‑149 (MQ=255)
               aaaTATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTAt                                                                                                                                      >  3:749821/1‑149 (MQ=255)
               aaaTATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTAt                                                                                                                                      <  3:1547315/149‑1 (MQ=255)
               aaaTATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTAt                                                                                                                                      <  4:794453/149‑1 (MQ=255)
                aaTATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTAtt                                                                                                                                     <  2:748139/149‑1 (MQ=255)
                  tatCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATtgt                                                                                                                                   <  4:315415/149‑1 (MQ=255)
                      tGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATtgt                                                                                                                                   <  3:1537035/145‑1 (MQ=255)
                      tGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATtgt                                                                                                                                   >  4:1537035/1‑145 (MQ=255)
                       ggggTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGcaca                                                                                                                              <  1:508645/149‑1 (MQ=255)
                       ggggTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGcaca                                                                                                                              >  2:638512/1‑149 (MQ=255)
                          gTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGaa                                                                                                                           >  1:860998/1‑149 (MQ=255)
                           ttGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGaaa                                                                                                                          >  4:1303215/1‑149 (MQ=255)
                                ggTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGcc                                                                                                                     >  2:269050/1‑149 (MQ=255)
                                ggTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGcc                                                                                                                     >  1:1470507/1‑149 (MQ=255)
                                   aGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGc                                                                                                                  <  1:1581018/149‑1 (MQ=255)
                                         cGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGCCCATCGGACC‑t                                                                                                                                           <  3:367151/118‑1 (MQ=255)
                                          ggTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCtt                                                                                                           >  4:549215/1‑149 (MQ=255)
                                            tCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTga                                                                                                         <  4:234082/149‑1 (MQ=255)
                                               gAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGAc                                                                                                      >  1:658192/1‑149 (MQ=255)
                                                aaaaCCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTa                                                                                                                                       <  3:1487674/115‑1 (MQ=255)
                                                aaaaCCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTa                                                                                                                                       >  4:1487674/1‑115 (MQ=255)
                                                 aaaCCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTg                                                                                                    >  3:636038/1‑149 (MQ=255)
                                                     cTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGc                                                                                                                      >  3:795096/1‑127 (MQ=255)
                                                     cTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCt                                                                                                >  2:1495618/1‑149 (MQ=255)
                                                        gACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTc                                                                                             <  3:780259/149‑1 (MQ=255)
                                                         aCAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCt                                                                                            <  4:1260455/149‑1 (MQ=255)
                                                         aCAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCt                                                                                            <  4:336284/149‑1 (MQ=255)
                                                          cAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGa                                                                                                                            >  4:1197398/1‑116 (MQ=255)
                                                           aGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTAtt                                                                                                                                     >  4:1282191/1‑106 (MQ=255)
                                                             tgtTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACa                                                                                        <  4:891759/149‑1 (MQ=255)
                                                             tgtTAACGACCATTCTCGGTTTTACCGGATGCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCGAGCCTCGCTTGAGACTGGCCTTTCTGACa                                                                                        <  3:1276295/149‑1 (MQ=255)
                                                              gtTAACGCCCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAg                                                                                       <  4:636038/149‑1 (MQ=255)
                                                                    gACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGAc                                                                                                                                              <  3:1303591/89‑1 (MQ=255)
                                                                    gACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGc                                                                                                             >  1:118046/1‑121 (MQ=255)
                                                                    gACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGctt                                                                                 <  4:221754/149‑1 (MQ=255)
                                                                    gACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGctt                                                                                 <  2:367828/149‑1 (MQ=255)
                                                                     aCCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGa                                                                                          >  4:1558955/1‑139 (MQ=255)
                                                                      ccATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGcttac                                                                               <  4:803199/149‑1 (MQ=255)
                                                                       cATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGcttact                                                                              >  4:1172835/1‑149 (MQ=255)
                                                                       cATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGcttact                                                                              >  1:297744/1‑149 (MQ=255)
                                                                         ttCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGc                                                                                                                                 <  4:776486/96‑1 (MQ=255)
                                                                            tcGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGTCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTc                                                                         >  1:1229166/1‑149 (MQ=255)
                                                                            tcGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTc                                                                         >  2:1106438/1‑149 (MQ=255)
                                                                            tcGGTTTTACCGGATTCTTAATTACCAGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTc                                                                         >  2:1528222/1‑149 (MQ=255)
                                                                                      cGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAgc                                                               >  3:944275/1‑149 (MQ=255)
                                                                                       ggATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAgcg                                                              >  1:662680/1‑149 (MQ=255)
                                                                                       ggATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAgcg                                                              <  2:337354/149‑1 (MQ=255)
                                                                                           tCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTc                                                          <  3:1268010/149‑1 (MQ=255)
                                                                                             ttAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTAtt                                                                                                                                     >  3:1157606/1‑72 (MQ=255)
                                                                                             ttAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGc                                                                                                             >  3:812996/1‑96 (MQ=255)
                                                                                             ttAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGc                                                                                                             <  4:812996/96‑1 (MQ=255)
                                                                                             ttAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGa                                                                                          <  1:1420703/115‑1 (MQ=255)
                                                                                             ttAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGa                                                                                          >  2:1420703/1‑115 (MQ=255)
                                                                                             ttAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGa                                                        >  1:760827/1‑149 (MQ=255)
                                                                                             ttAATTACCTGGTGCGTATGGGCAGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTAtt                                                                                                                                     <  4:1157606/72‑1 (MQ=255)
                                                                                              tAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGat                                                       <  2:125445/149‑1 (MQ=255)
                                                                                                 ttACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATAtt                                                    >  4:1354896/1‑149 (MQ=255)
                                                                                                 ttACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATACTTACTTACTCGCGGAACAGCGCTTCGATAtt                                                    <  2:1470507/149‑1 (MQ=255)
                                                                                                     cTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCGTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGc                                                >  3:1292455/1‑149 (MQ=255)
                                                                                                      tGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGcc                                               >  2:397439/1‑149 (MQ=255)
                                                                                                         tGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCt                                            <  3:1172835/149‑1 (MQ=255)
                                                                                                          gCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTg                                           >  4:69436/1‑149 (MQ=255)
                                                                                                             tATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGt                                        >  3:495919/1‑149 (MQ=255)
                                                                                                                  gCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCa                                   <  4:944275/149‑1 (MQ=255)
                                                                                                                  gCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCa                                   <  1:540861/149‑1 (MQ=255)
                                                                                                                            gACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTcgcg                         <  3:662260/149‑1 (MQ=255)
                                                                                                                               cTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGcttac                                                                               >  3:1052781/1‑92 (MQ=255)
                                                                                                                               cTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGcttac                                                                               <  4:1052781/92‑1 (MQ=255)
                                                                                                                               cTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCaa                      >  1:1096173/1‑149 (MQ=255)
                                                                                                                               cTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCaa                      >  4:1440726/1‑149 (MQ=255)
                                                                                                                               cTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCaa                      >  2:454296/1‑149 (MQ=255)
                                                                                                                               cTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCaa                      >  3:375886/1‑149 (MQ=255)
                                                                                                                               cTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCaa                      >  1:1419927/1‑149 (MQ=255)
                                                                                                                                ttAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGTCTCAAGCGAGGCTGGCCttttct                                                                                           <  4:1371973/79‑4 (MQ=18)
                                                                                                                                ttAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCaaa                     >  1:884427/1‑149 (MQ=255)
                                                                                                                                       aaaGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCtt                                                                                                           <  3:473450/56‑1 (MQ=255)
                                                                                                                                         aGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAgg            >  4:697524/1‑149 (MQ=255)
                                                                                                                                         aGGTCCGATGCCCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAgg            >  1:1465017/1‑149 (MQ=255)
                                                                                                                                           gTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGcc          >  4:203139/1‑149 (MQ=255)
                                                                                                                                            tCCGATGCCCATCGGACCTTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGcc          >  1:751071/1‑149 (MQ=255)
                                                                                                                                                 tGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAAc    <  4:379747/149‑1 (MQ=255)
                                                                                                                                                  gggCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAAcc   >  3:904568/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TTCTGGATTGTGACCAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCT  >  NC_000913/2866340‑2866635

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