Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,866,487 G→C 6.5% 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 6.5% 865.1 / 40.4 247intergenic (+22/+72)ygbN/rpoSputative transporter/RNA polymerase, sigma S (sigma 38) factor
Reads supporting (aligned to +/- strand):  ref base G (122/108);  new base C (11/5);  total (133/114)
Fisher's exact test for biased strand distribution p-value = 3.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.02e-01

Only 100 of 254 total aligned reads displayed.

ACAACCGGCGGGTTACTCTCAGAAGCGGTGATGGGATTGAATCCGATTCAGTGCGTATTGGTGACGCTGGCAGCCTGCTTTGGTGGGCTTGGTGCGTCACATATTAATGACTCAGGGTTCTGGATTGTGACCAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGTACCCCAGCAAACCGAATCGACGTGCCAGCACTTCACGCTGTTTGGCGTTC  >  NC_000913/2866223‑2866754
                                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                            
              aCTCTCAGAAGCGGTGATGGGATTGAATCCGATTCAGTGCGTATTGGTGACGCTGGCAGCCTGCTTTGGTGGGCTTGGTGCGTCACATATTAATGACTCAGGGTTCTGGATTGTGACCAAATAGCTGTGTTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGCTAACGACCATTCCAGGGTGTAACGGATTCTTAATTTCCGGGTGCGTGTGGGCGGTGATTTTAACTTAATAAAAAGGTCCGCTGGGCATCTGATCTGTTaacacattag                                                                                                                                                                                                                                                     >  4:2756/1‑266 (MQ=255)
                          cgTGATGGGATTGAAGCCGTTTCGTGGCGTATTGTGGACGCTGGCAGCCTGCTTTGGTGGGCTGGGTGTGTCACATATTAATGACTCAGGGTTCTGGATTGTGACCAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGc                                                                                                                                                                                                                                         <  4:205892/274‑1 (MQ=255)
                                                                             cTTTGGTGGGCTTGGTGCGTCACATATTAATGACTCAGGGTTCTGGATTGTGACCAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAg                                                                                                                                                                                       >  2:111219/1‑274 (MQ=255)
                                                                                                                    gTTCTGGATTGTGACCAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATgg                                                                                                                                                                                                                                                                              <  5:52149/149‑1 (MQ=255)
                                                                                                                     ttCTGGATTGTGACCAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATggg                                                                                                                                                                                                                                                                             >  8:600681/1‑149 (MQ=255)
                                                                                                                       cTGGATTGTGACCAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCa                                                                                                                                                                                                                                                                           >  7:1216290/1‑149 (MQ=255)
                                                                                                                         ggATTGTGACCAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATc                                                                                                                                                                                                                                                                         <  7:293935/149‑1 (MQ=255)
                                                                                                                          gATTGTGACCAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCg                                                                                                                                                                                                                                                                        <  5:914554/149‑1 (MQ=255)
                                                                                                                            ttGTGACCAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGa                                                                                                                                                                                                                                                                      <  8:723270/149‑1 (MQ=255)
                                                                                                                             tgtgACCAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGAc                                                                                                                                                                                                                                                                     <  8:1216290/149‑1 (MQ=255)
                                                                                                                               tgACCAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCt                                                                                                                                                                                                                                                                  >  1:330/1‑149 (MQ=255)
                                                                                                                                gACCAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCtt                                                                                                                                                                                                                                                                 <  5:1374495/149‑1 (MQ=255)
                                                                                                                                  ccAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCtttt                                                                                                                                                                                                                                                               >  5:72318/1‑149 (MQ=255)
                                                                                                                                  ccAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCtttt                                                                                                                                                                                                                                                               <  5:816563/149‑1 (MQ=255)
                                                                                                                                  ccAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCtttt                                                                                                                                                                                                                                                               >  6:198299/1‑149 (MQ=255)
                                                                                                                                  ccAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCtttt                                                                                                                                                                                                                                                               <  5:1274573/149‑1 (MQ=255)
                                                                                                                                  ccAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCtttt                                                                                                                                                                                                                                                               >  5:1335057/1‑149 (MQ=255)
                                                                                                                                    aaaTATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGTGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGTCTTTCTGACAGATGCTTACTTACCCGCTGAACAGCGCTTCGGTATTCAGCCCCTGCCTTTGCCGGCTTTCGCCCAAACGGCCCAGGCCt                                                                                                                                >  2:14910/1‑274 (MQ=255)
                                                                                                                                    aaaTATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTAt                                                                                                                                                                                                                                                             >  5:1210494/1‑149 (MQ=255)
                                                                                                                                     aaTATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTa                                                                                                                                                                                                                                                              <  7:330097/147‑1 (MQ=255)
                                                                                                                                     aaTATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTAtt                                                                                                                                                                                                                                                            <  6:143320/149‑1 (MQ=255)
                                                                                                                                       tatCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATtgt                                                                                                                                                                                                                                                          <  6:120767/149‑1 (MQ=255)
                                                                                                                                            ggggTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTAtt                                                                                                                                                                                                                                                            >  5:792006/1‑142 (MQ=255)
                                                                                                                                            ggggTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTAtt                                                                                                                                                                                                                                                            <  6:792006/142‑1 (MQ=255)
                                                                                                                                            ggggTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGcaca                                                                                                                                                                                                                                                     <  6:1283902/149‑1 (MQ=255)
                                                                                                                                                ttGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGCGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAACAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCCGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCATTCAACCTGAAtt                                                                                                                    >  2:276596/1‑273 (MQ=255)
                                                                                                                                                  gTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAg                                                                                                                                                                                                                                               >  8:407776/1‑149 (MQ=255)
                                                                                                                                                  gTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAg                                                                                                                                                                                                                                               <  5:665140/149‑1 (MQ=255)
                                                                                                                                                  gTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTNTATTGTNCACAGAAAAGGCCAGCCTCGCTTtag                                                                                                                                                                                                                               >  3:21544/1‑165 (MQ=255)
                                                                                                                                                   tCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAgg                                                                                                                                                                                                                                              <  8:867967/149‑1 (MQ=255)
                                                                                                                                                    cGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGc                                                                                                                                                                                                                                             <  6:56062/149‑1 (MQ=255)
                                                                                                                                                      gTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCa                                                                                                                                                                                                                                           <  8:532627/149‑1 (MQ=255)
                                                                                                                                                           aGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTc                                                                                                                                                                                                                                      >  6:1459404/1‑149 (MQ=255)
                                                                                                                                                            gACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCt                                                                                                                                                                                                                                                                  <  6:512610/120‑1 (MQ=255)
                                                                                                                                                              cGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCt                                                                                                                                                                                                                                   <  7:753076/149‑1 (MQ=255)
                                                                                                                                                               ggTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCtt                                                                                                                                                                                                                                  >  7:439332/1‑149 (MQ=255)
                                                                                                                                                                 tCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGaaaa                                                                                                                                                                                                                                                >  8:643628/1‑133 (MQ=255)
                                                                                                                                                                 tCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTga                                                                                                                                                                                                                                <  5:222695/149‑1 (MQ=255)
                                                                                                                                                                  cTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTgag                                                                                                                                                                                                                               >  7:386094/1‑149 (MQ=255)
                                                                                                                                                                   tGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTgaga                                                                                                                                                                                                                              >  8:235232/1‑149 (MQ=255)
                                                                                                                                                                   tGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTgaga                                                                                                                                                                                                                              >  8:1266358/1‑149 (MQ=255)
                                                                                                                                                                       aaCCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTgg                                                                                                                                                                                                                          >  6:324095/1‑149 (MQ=255)
                                                                                                                                                                           tGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGc                                                                                                                                                                                                                                             >  6:162518/1‑126 (MQ=255)
                                                                                                                                                                            ggACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCttt                                                                                                                                                                                                                     <  7:927308/149‑1 (MQ=255)
                                                                                                                                                                             gACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTc                                                                                                                                                                                                                    <  8:804940/149‑1 (MQ=255)
                                                                                                                                                                              aCAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCt                                                                                                                                                                                                                   <  5:466635/149‑1 (MQ=255)
                                                                                                                                                                                aGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGTTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGa                                                                                                                                                                                                                 >  8:785147/1‑149 (MQ=255)
                                                                                                                                                                                aGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGa                                                                                                                                                                                                                 <  8:1272127/149‑1 (MQ=255)
                                                                                                                                                                                  tgtTAACGACCATTCTCGGTTTTACCGGATTCTTACTTACCTGGTGCGTATGGGCGGTCATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTCCCCGCGGACCCGCGCTTAGATATTCATCCCCTTCGTTTGCAGGCTTTCGCGCAACCGGCCCCGGCCTTCAACCTGAATCTGGCGCACACTTTCCCGGCTCAGGCCAATTTCAc                                                                                  >  4:163621/1‑274 (MQ=255)
                                                                                                                                                                                  tgtTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCt                                                                                                                                                                                                                                   <  7:107685/129‑1 (MQ=255)
                                                                                                                                                                                  tgtTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTc                                                                                                                                                                                                                    <  7:1272536/144‑1 (MQ=255)
                                                                                                                                                                                     tAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGAt                                                                                                                                                                                                            <  7:1202012/149‑1 (MQ=255)
                                                                                                                                                                                      aaCGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCtt                                                                                                                                                                                                                                  <  6:1328245/126‑1 (MQ=255)
                                                                                                                                                                                      aaCGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATg                                                                                                                                                                                                           <  8:439332/149‑1 (MQ=255)
                                                                                                                                                                                      aaCGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATg                                                                                                                                                                                                           >  7:1399357/1‑149 (MQ=255)
                                                                                                                                                                                            cATTCTCGGTTTTACCGGCTTCTTAATTACCTGGTGCGTATGGGCGGTAATTTTACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAACAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGCCAGATGCTTACTTACTCGCGGAACAgcga                                                                                                                                                                                    >  2:226101/1‑165 (MQ=255)
                                                                                                                                                                                            cATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGcttact                                                                                                                                                                                                     >  8:177717/1‑149 (MQ=255)
                                                                                                                                                                                            cATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGcttact                                                                                                                                                                                                     >  6:892529/1‑149 (MQ=255)
                                                                                                                                                                                                 tcGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCTGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCTTTCTTACAGATTCTTACTTACTcg                                                                                                                                                                                               >  8:843163/1‑149 (MQ=255)
                                                                                                                                                                                                 tcGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTc                                                                                                                                                                                                >  6:659395/1‑149 (MQ=255)
                                                                                                                                                                                                    gTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTcgcg                                                                                                                                                                                             <  5:1324262/149‑1 (MQ=255)
                                                                                                                                                                                                       ttACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGaa                                                                                                                                                                                          <  5:1172614/149‑1 (MQ=255)
                                                                                                                                                                                                               ttCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCtt                                                                                                                                                                                  <  6:104089/149‑1 (MQ=255)
                                                                                                                                                                                                                tCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTc                                                                                                                                                                                 <  7:461291/149‑1 (MQ=255)
                                                                                                                                                                                                                  ttAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCCGGTTTTCCGCTTCGttc                                                  >  4:94785/1‑272 (MQ=255)
                                                                                                                                                                                                                   tAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGTAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGat                                                                                                                                                                              <  6:1150968/149‑1 (MQ=255)
                                                                                                                                                                                                                    aaTTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGata                                                                                                                                                                             <  7:868477/149‑1 (MQ=255)
                                                                                                                                                                                                                       tACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTc                                                                                                                                                                          <  6:422829/149‑1 (MQ=255)
                                                                                                                                                                                                                        aCCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTc                                                                                                                                                                                                                    >  2:362334/1‑106 (MQ=255)
                                                                                                                                                                                                                         ccTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAg                                                                                                                                                                        <  8:359115/149‑1 (MQ=255)
                                                                                                                                                                                                                            ggTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGAc                                                                                                                                                                                                                <  4:362334/106‑1 (MQ=255)
                                                                                                                                                                                                                                 gTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTga                                                                                                                                                                                                                                <  8:561714/85‑1 (MQ=255)
                                                                                                                                                                                                                                 gTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTga                                                                                                                                                                                                                                >  7:561714/1‑85 (MQ=255)
                                                                                                                                                                                                                                   aTGGGCGGTAATTGGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGTACCCCAGCAAACCGAATCGa                                <  2:94785/275‑1 (MQ=255)
                                                                                                                                                                                                                                      ggCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGc                                                                                                                                                           <  7:600681/149‑1 (MQ=255)
                                                                                                                                                                                                                                       gCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCa                                                                                                                                                          >  7:787349/1‑149 (MQ=255)
                                                                                                                                                                                                                                       gCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCa                                                                                                                                                          <  5:397257/149‑1 (MQ=255)
                                                                                                                                                                                                                                       gCGGTAATGTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCa                                                                                                                                                          <  7:405828/149‑1 (MQ=255)
                                                                                                                                                                                                                                        cGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAgcgc                                                                                                                                                                                    >  1:20347/1‑122 (MQ=255)
                                                                                                                                                                                                                                          gTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGa                                                                                                                                                       >  8:1000772/1‑149 (MQ=255)
                                                                                                                                                                                                                                          gTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGa                                                                                                                                                       >  5:175378/1‑149 (MQ=255)
                                                                                                                                                                                                                                                 gACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTcgcg                                                                                                                                                >  8:259647/1‑149 (MQ=255)
                                                                                                                                                                                                                                                     ttAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGc                                                                                                                                                                       <  6:832738/122‑1 (MQ=255)
                                                                                                                                                                                                                                                     ttAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCaaa                                                                                                                                            >  8:119148/1‑149 (MQ=255)
                                                                                                                                                                                                                                                     ttAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCaaa                                                                                                                                            >  7:750854/1‑149 (MQ=255)
                                                                                                                                                                                                                                                     ttAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCaaa                                                                                                                                            >  5:95189/1‑149 (MQ=255)
                                                                                                                                                                                                                                                     cacaataaaaGGTCCGATGCCCATCGGACCTTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAg                                                                                                                                                         <  6:986435/131‑1 (MQ=38)
                                                                                                                                                                                                                                                       aataaAAAGGTCCGCTGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTAGCTTTAGACTGGCCTTTCTGACCGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACTCGTTCACGGGTGAGGCCAATTTCACTCCCTAAATCTTCCAGTGTTGCTGCTTCTCTCCCCACACAACCGAATCCACGTCCACGCACttgcacct             >  4:227745/1‑268 (MQ=255)
                                                                                                                                                                                                                                                         taaAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGgc                                                                                                                                        >  6:759229/1‑149 (MQ=255)
                                                                                                                                                                                                                                                          aaaaaGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGgcg                                                                                                                                       >  5:950733/1‑149 (MQ=255)
                                                                                                                                                                                                                                                              aGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAg                                                                                                                                                         <  5:1072/127‑1 (MQ=255)
                                                                                                                                                                                                                                                              aGGTCCGATGCCCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTATGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAgg                                                                                                                                   >  8:101822/1‑149 (MQ=255)
                                                                                                                                                                                                                                                               ggTCCGATGCCCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGc                                                                                                                                  >  7:209723/1‑149 (MQ=255)
                                                                                                                                                                                                                                                               ggTCCGATGCCCATCGGACCTTTTATTGTGCACAGAAAAGGCAAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGcttactta                                                                                                                                                                                                   <  5:283693/84‑1 (MQ=38)
                                                                                                                                                                                                                                                                gTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGcc                                                                                                                                 >  5:557367/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                 tCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGCTTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGTACCCCAGCAAACCGAATCGACGTGCCAGCACTTCACGCTGTTTTGAGTTc  >  2:171632/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                   cGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTcg                                                                                                                                                                                               >  8:1418274/1‑84 (MQ=255)
                                                                                                                                                                                                                                                                      tGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAAc                                                                                                                           <  7:1266358/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                       gggCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAAcc                                                                                                                          >  5:642604/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                        ggCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCt                                                                                                                         <  7:240698/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                            
ACAACCGGCGGGTTACTCTCAGAAGCGGTGATGGGATTGAATCCGATTCAGTGCGTATTGGTGACGCTGGCAGCCTGCTTTGGTGGGCTTGGTGCGTCACATATTAATGACTCAGGGTTCTGGATTGTGACCAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGTACCCCAGCAAACCGAATCGACGTGCCAGCACTTCACGCTGTTTGGCGTTC  >  NC_000913/2866223‑2866754

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