Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,866,487 G→C 5.1% intergenic (+22/+72) 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,4870GC 5.1% 894.0 / 28.6 255intergenic (+22/+72)ygbN/rpoSputative transporter/RNA polymerase, sigma S (sigma 38) factor
Reads supporting (aligned to +/- strand):  ref base G (120/120);  new base C (7/6);  total (128/127)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.44e-01

Only 100 of 261 total aligned reads displayed.

TGACAACCGGCGGGTTACTCTCAGAAGCGGTGATGGGATTGAATCCGATTCAGTGCGTATTGGTGACGCTGGCAGCCTGCTTTGGTGGGCTTGGTGCGTCACATATTAATGACTCAGGGTTCTGGATTGTGACCAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGTACCCCAGCAAACCGAATCGACGTGCCAGCACTTCA  >  NC_000913/2866221‑2866739
                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                             
                                                                 aCGCTGGCAGCCTGCTTTGGTGGGCTTGGTGCGTCACATATTAATGACTCAGGGTTCTGGATTGTGACCACATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTTGTTCGTATGGGCGGTAATTTGACCTTAATAAAATGGTCAGATGGGCATCTGAACTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTGGGACACTGCCTTTCGGACAAATTcttaat                                                                                                                                                                                      >  4:516320/1‑272 (MQ=255)
                                                                    cgggCAGCCTGGTTGGGTGGGCTTGGTGCGTGACAGTTTAATGACTCAGGGTTCTGGATTGTGACCAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAAACTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTcgcg                                                                                                                                                                              <  1:122709/277‑1 (MQ=255)
                                                                                tttGGTGGGCTTGGTGCGTCACATATTAATGACTCAGGGTTCTGGATTGTGACCAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAgcg                                                                                                                                                                      <  2:277817/275‑1 (MQ=255)
                                                                                     tggGCTTGGTGCGTCACATATTAATGACTCAGGGTTCTGGATTGTGACCAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTa                                                                                                                                                                                                                                                <  2:654577/197‑1 (MQ=255)
                                                                                                                cTCAGGGTTCTGGATTGTGACCAAATATCTGGGGTTGCCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATCAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGTAAAGGCCAGCCTCGCTTTAGACTGGCATTTCTGCCAGCTTCTTACTTCCCCCCGGACCACGGCTTCGATATTCACCCCCTCCGTTTGCAGGAt                                                                                                                                       >  4:131326/1‑274 (MQ=255)
                                                                                                                  cAGGGTTCTGGATTGTGACCAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTc                                                                                                                                    >  2:581502/1‑275 (MQ=255)
                                                                                                                   aGGGTTCTGGATTGTGACCAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTGTGCAGGAGTTCGCGCa                                                                                                                               >  1:187338/1‑279 (MQ=255)
                                                                                                                       ttCTGGATTGTGACCAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATggg                                                                                                                                                                                                                                                              >  7:1507401/1‑149 (MQ=255)
                                                                                                                       ttCTGGATTGTGACCAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCgg                                                                                                                                                                                                                                                        <  4:572155/155‑1 (MQ=255)
                                                                                                                           ggATTGTGACCAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATc                                                                                                                                                                                                                                                          <  7:1378496/149‑1 (MQ=255)
                                                                                                                              ttGTGACCAAATGTCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTGTACCGGATTCTTAATTACCTGGTGCGTAGGGGCGGGAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGa                                                                                                                                                                                                                                                       <  6:204669/149‑1 (MQ=255)
                                                                                                                              ttGTGACCAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGa                                                                                                                                                                                                                                                       <  8:328231/149‑1 (MQ=255)
                                                                                                                               tgtgACCAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGCCCATCGGAc                                                                                                                                                                                                                                                      <  6:560831/149‑1 (MQ=255)
                                                                                                                                gtgACCAAATATCTGGGGTTTTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCa                                                                                                                      <  2:435197/275‑1 (MQ=255)
                                                                                                                                gtgACCAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGAcc                                                                                                                                                                                                                                                     >  6:1354438/1‑149 (MQ=255)
                                                                                                                                  gACCAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCtt                                                                                                                                                                                                                                                   >  7:1138450/1‑149 (MQ=255)
                                                                                                                                  gACCAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCtt                                                                                                                                                                                                                                                   >  7:518136/1‑149 (MQ=255)
                                                                                                                                   aCCAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCttt                                                                                                                                                                                                                                                  <  5:1633639/149‑1 (MQ=255)
                                                                                                                                   aCCAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCttt                                                                                                                                                                                                                                                  >  5:125538/1‑149 (MQ=255)
                                                                                                                                   aCCAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCttt                                                                                                                                                                                                                                                  >  5:421554/1‑149 (MQ=255)
                                                                                                                                   aCCAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCttt                                                                                                                                                                                                                                                  >  8:627850/1‑149 (MQ=255)
                                                                                                                                     cAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTa                                                                                                                                                                                                                                                <  5:1103867/149‑1 (MQ=255)
                                                                                                                                      aaaTATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATc                                                                                                                                                                                                                                                          <  6:113494/138‑1 (MQ=255)
                                                                                                                                             tGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGcac                                                                                                                                                                                                                                        <  6:280937/149‑1 (MQ=255)
                                                                                                                                                   tGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCt                                                                                                                                                                                                                         >  2:505009/1‑158 (MQ=255)
                                                                                                                                                    gTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAg                                                                                                                                                                                                                                 <  6:640374/149‑1 (MQ=255)
                                                                                                                                                    gTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGACCGATGGGCATCGGACCTTTTATTGTGCCCAGAAAAGGCCAGCCTCGCTGCAGACTGGCCTTTc                                                                                                                                                                                                     >  4:441885/1‑176 (MQ=255)
                                                                                                                                                      cGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGc                                                                                                                                                                                                                               <  7:237048/149‑1 (MQ=255)
                                                                                                                                                        gTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCa                                                                                                                                                                                                                             <  5:719148/149‑1 (MQ=255)
                                                                                                                                                              gACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATgg                                                                                                                                                                                                                                                               <  6:445520/109‑1 (MQ=255)
                                                                                                                                                              gACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCg                                                                                                                                                                                                                       >  5:951550/1‑149 (MQ=255)
                                                                                                                                                              gACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGACAGCCTCg                                                                                                                                                                                                                       >  7:17563/1‑149 (MQ=255)
                                                                                                                                                               aCGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCa                                                                                                                                                                                                                                                            >  4:468277/1‑111 (MQ=255)
                                                                                                                                                                 ggTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCtt                                                                                                                                                                                                                   <  5:1304401/149‑1 (MQ=255)
                                                                                                                                                                   tCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTga                                                                                                                                                                                                                 <  5:313861/149‑1 (MQ=255)
                                                                                                                                                                   tCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACTTTCACGGGTGa                                                                               >  1:80895/1‑279 (MQ=255)
                                                                                                                                                                     tGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTgaga                                                                                                                                                                                                               >  5:18808/1‑149 (MQ=255)
                                                                                                                                                                     tGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTgaga                                                                                                                                                                                                               >  7:1276412/1‑149 (MQ=255)
                                                                                                                                                                      gAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGACGGGCATCGGACCTTTTATtgtg                                                                                                                                                                                                                                           <  8:1362961/121‑1 (MQ=255)
                                                                                                                                                                       aaaaCCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAg                                                                                                                                                                                                                                 <  6:700172/130‑1 (MQ=255)
                                                                                                                                                                       aaaaCCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAg                                                                                                                                                                                                                                 >  5:700172/1‑130 (MQ=255)
                                                                                                                                                                         aaCCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCttt                                                                                                                                                                                                                                                  >  8:1098471/1‑111 (MQ=255)
                                                                                                                                                                          aCCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGca                                                                                                                                                                                                                                         >  5:973035/1‑119 (MQ=255)
                                                                                                                                                                          aCCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGc                                                                                                                                                                                                          <  6:159247/149‑1 (MQ=255)
                                                                                                                                                                            cTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCt                                                                                                                                                                                                        >  7:1236578/1‑149 (MQ=255)
                                                                                                                                                                             tGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCtt                                                                                                                                                                                                       >  7:412597/1‑149 (MQ=255)
                                                                                                                                                                             tGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCtt                                                                                                                                                                                                       >  6:744695/1‑149 (MQ=255)
                                                                                                                                                                             tGGACAGTGTTAACGACCATTCTCGGTTTTACAGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGAGGCTTACTTACTCTCGAAACAGCGCTTCGa                                                                                                                                                                >  4:313716/1‑188 (MQ=255)
                                                                                                                                                                               gACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTc                                                                                                                                                                                                     <  6:147830/149‑1 (MQ=255)
                                                                                                                                                                               gACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTc                                                                                                                                                                                                     <  7:1359335/149‑1 (MQ=255)
                                                                                                                                                                                 cAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATAtt                                                                                                                                                            <  2:313720/188‑1 (MQ=255)
                                                                                                                                                                                   gtgtTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGAc                                                                                                                                                                                                 <  6:1077373/149‑1 (MQ=255)
                                                                                                                                                                                    tgtTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACCGATGCTTACTTACTCGCTGAACAGCGCTTCGATATTCCACCCCTGCGTTTCCAGGATTTCGCGCAAACCGCCCAGGCCTTCAACCTGAATCTGGCCAACCACTTACCGGGTGAGGCCAATTTCAca                                                                  >  4:500332/1‑274 (MQ=255)
                                                                                                                                                                                       tAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGAt                                                                                                                                                                                             <  5:1160352/149‑1 (MQ=255)
                                                                                                                                                                                        aaCGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGc                                                                                                                                                                                           <  4:741863/150‑1 (MQ=255)
                                                                                                                                                                                          cGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGct                                                                                                                                                                                          <  7:746051/149‑1 (MQ=255)
                                                                                                                                                                                          cGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAgcg                                                                                                                                                                      <  2:522088/169‑1 (MQ=255)
                                                                                                                                                                                           gACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGctt                                                                                                                                                                                         >  6:123471/1‑149 (MQ=255)
                                                                                                                                                                                           gACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGctt                                                                                                                                                                                         >  7:1003812/1‑149 (MQ=255)
                                                                                                                                                                                           gACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGCCCATCGGACCttt                                                                                                                                                                                                                                                  >  5:593294/1‑93 (MQ=255)
                                                                                                                                                                                             ccATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGcttac                                                                                                                                                                                       <  8:1626331/149‑1 (MQ=255)
                                                                                                                                                                                              cATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTTACNGATGCTTACTGTCTCGCGGAAAAGCGCTTCAATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAACCGGCCCCNGCCCTCAACCTGACTCTGGCGACAACGTTCNNCGGTGAGGCCAACTTCACCAACAACAGCTTc                                                      >  3:106077/1‑277 (MQ=255)
                                                                                                                                                                                              cATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGc                                                                                                                                                                                           >  4:703759/1‑144 (MQ=255)
                                                                                                                                                                                                  ctcGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGcttac                                                                                                                                                                                       <  2:703764/144‑1 (MQ=255)
                                                                                                                                                                                                   tcGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTc                                                                                                                                                                                 >  8:67752/1‑149 (MQ=255)
                                                                                                                                                                                                             cGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAgc                                                                                                                                                                       >  6:1698822/1‑149 (MQ=255)
                                                                                                                                                                                                              ggATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAgcg                                                                                                                                                                      >  8:1162934/1‑149 (MQ=255)
                                                                                                                                                                                                                aTTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCt                                                                                                                                                                    <  5:1412672/149‑1 (MQ=255)
                                                                                                                                                                                                                    ttAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGa                                                                                                                                                                <  5:1543505/149‑1 (MQ=255)
                                                                                                                                                                                                                    ttAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGa                                                                                                                                                                >  8:563006/1‑149 (MQ=255)
                                                                                                                                                                                                                          gccTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTATGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGTGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGTACCCCAGCa                            <  2:131328/274‑1 (MQ=255)
                                                                                                                                                                                                                            cTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGc                                                                                                                                                        >  8:1500804/1‑149 (MQ=255)
                                                                                                                                                                                                                            cTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGc                                                                                                                                                        <  6:1436308/149‑1 (MQ=255)
                                                                                                                                                                                                                             tGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGcc                                                                                                                                                       >  5:802660/1‑149 (MQ=255)
                                                                                                                                                                                                                                tGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGc                                                                                                                                                                                                          >  2:182476/1‑95 (MQ=255)
                                                                                                                                                                                                                                    tATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCttt                                                                                                                                                                                                      <  4:182473/95‑1 (MQ=255)
                                                                                                                                                                                                                                    tATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGt                                                                                                                                                >  5:346140/1‑149 (MQ=255)
                                                                                                                                                                                                                                    tATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACg                                                                                          >  4:417636/1‑203 (MQ=255)
                                                                                                                                                                                                                                      tGGGAGTTAATTTGACCTTAATAAAAAGGTCCGATAGGCAGCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGTACCCCAGCAAACCGAATCGACGTGc            <  1:264747/279‑1 (MQ=255)
                                                                                                                                                                                                                                                tttGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTc                                                                                                                                    <  5:835493/149‑1 (MQ=255)
                                                                                                                                                                                                                                                tttGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGt                                                                                 <  2:450369/200‑1 (MQ=255)
                                                                                                                                                                                                                                                      cTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGa                                                                                                                                        <  6:392906/139‑1 (MQ=255)
                                                                                                                                                                                                                                                      cTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCaa                                                                                                                              >  5:1018022/1‑149 (MQ=255)
                                                                                                                                                                                                                                                       ttAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCaaa                                                                                                                             >  7:375066/1‑149 (MQ=255)
                                                                                                                                                                                                                                                         aataaAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGCAAAGGCCAGCCTCGCTTGAGACTCGCTTTTCTGACAGATTCTTACTTACTCCCGTAACCACGCTTCGACATTCTACCCCCGCCTTTGCCGGATTNCCCGCACACGGCGCAGCCGTTCAACCTGAATTTGGCCAac                                                                                             >  3:237039/1‑179 (MQ=255)
                                                                                                                                                                                                                                                         aataaAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAac                                                                                             <  1:237039/179‑1 (MQ=255)
                                                                                                                                                                                                                                                         aataaAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGACCAGCCTCGCTCGAGACTGGCCTTTc                                                                                                                                                                                                     >  7:975527/1‑75 (MQ=255)
                                                                                                                                                                                                                                                           taaAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGgc                                                                                                                         <  8:1236578/149‑1 (MQ=255)
                                                                                                                                                                                                                                                            aaaaaGGTCCGATGGGCATCGGACCTTTTATTGTGCGCAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTcgc                                                                                                                                                                               <  6:1401834/94‑1 (MQ=255)
                                                                                                                                                                                                                                                             aaaaGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTACAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTAc                                                            <  2:391300/208‑1 (MQ=255)
                                                                                                                                                                                                                                                               aaGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAg                                                                                                                     >  8:277239/1‑149 (MQ=255)
                                                                                                                                                                                                                                                               aaGGTCCGATGCCCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAg                                                                                                                     >  8:815430/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                 ggTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGc                                                                                                                   <  5:1354438/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                 ggTCCGATGCCCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTc                                                                                                                                                           <  6:881230/109‑1 (MQ=255)
                                                                                                                                                                                                                                                                  gTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGcc                                                                                                                  >  5:928335/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                   tCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCt                                                                                                                 <  5:123471/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                     cGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAg                                                                                                                     >  6:526539/1‑143 (MQ=255)
                                                                                                                                                                                                                                                                     cGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAg                                                                                                                     <  5:526539/143‑1 (MQ=255)
                                                                                                                                                                                                                                                                      gATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCa                                                                                                              >  7:161704/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                          ggCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCt                                                                                                          <  7:563006/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                             
TGACAACCGGCGGGTTACTCTCAGAAGCGGTGATGGGATTGAATCCGATTCAGTGCGTATTGGTGACGCTGGCAGCCTGCTTTGGTGGGCTTGGTGCGTCACATATTAATGACTCAGGGTTCTGGATTGTGACCAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGTACCCCAGCAAACCGAATCGACGTGCCAGCACTTCA  >  NC_000913/2866221‑2866739

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