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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,866,4880GC 6.5% 862.2 / 41.3 247intergenic (+23/+71)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 = 1.13e-02

Only 100 of 251 total aligned reads displayed.

ACAACCGGCGGGTTACTCTCAGAAGCGGTGATGGGATTGAATCCGATTCAGTGCGTATTGGTGACGCTGGCAGCCTGCTTTGGTGGGCTTGGTGCGTCACATATTAATGACTCAGGGTTCTGGATTGTGACCAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGTACCCCAGCAAACCGAATCGACGTGCCAGCACTTCACGCTGTTTGGCGTTC  >  NC_000913/2866223‑2866754
                                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                          
                          cgTGATGGGATTGAAGCCGTTTCGTGGCGTATTGTGGACGCTGGCAGCCTGCTTTGGTGGGCTGGGTGTGTCACATATTAATGACTCAGGGTTCTGGATTGTGACCAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGc                                                                                                                                                                                                                                         <  4:205892/274‑1 (MQ=255)
                                                                             cTTTGGTGGGCTTGGTGCGTCACATATTAATGACTCAGGGTTCTGGATTGTGACCAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAg                                                                                                                                                                                       >  2:111219/1‑274 (MQ=255)
                                                                                                                     ttCTGGATTGTGACCAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATggg                                                                                                                                                                                                                                                                            >  8:600681/1‑149 (MQ=255)
                                                                                                                     ttCTGGATTGTGACCAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATggg                                                                                                                                                                                                                                                                            >  7:1255766/1‑149 (MQ=255)
                                                                                                                       cTGGATTGTGACCAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCa                                                                                                                                                                                                                                                                          >  7:1216290/1‑149 (MQ=255)
                                                                                                                          gATTGTGACCAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCg                                                                                                                                                                                                                                                                       <  5:914554/149‑1 (MQ=255)
                                                                                                                            ttGTGACCAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGa                                                                                                                                                                                                                                                                     <  8:723270/149‑1 (MQ=255)
                                                                                                                             tgtgACCAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGAc                                                                                                                                                                                                                                                                    <  8:1216290/149‑1 (MQ=255)
                                                                                                                              gtgACCAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGAcc                                                                                                                                                                                                                                                                   <  8:1367738/149‑1 (MQ=255)
                                                                                                                               tgACCAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCt                                                                                                                                                                                                                                                                  >  1:330/1‑149 (MQ=255)
                                                                                                                                gACCAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCtt                                                                                                                                                                                                                                                                 <  5:1374495/149‑1 (MQ=255)
                                                                                                                                gACCAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGCCCATCGGACCtt                                                                                                                                                                                                                                                                 <  7:259878/149‑1 (MQ=255)
                                                                                                                                  ccAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCtttt                                                                                                                                                                                                                                                               >  5:72318/1‑149 (MQ=255)
                                                                                                                                  ccAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCtttt                                                                                                                                                                                                                                                               >  8:329610/1‑149 (MQ=255)
                                                                                                                                  ccAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCtttt                                                                                                                                                                                                                                                               <  5:816563/149‑1 (MQ=255)
                                                                                                                                  ccAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCtttt                                                                                                                                                                                                                                                               >  6:198299/1‑149 (MQ=255)
                                                                                                                                    aaaTATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGTGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGTCTTTCTGACAGATGCTTACTTACCCGCTGAACAGCGCTTCGGTATTCAGCCCCTGCCTTTGCCGGCTTTCGCCCAAACGGCCCAGGCCt                                                                                                                                >  2:14910/1‑274 (MQ=255)
                                                                                                                                    aaaTATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTAt                                                                                                                                                                                                                                                             >  5:1210494/1‑149 (MQ=255)
                                                                                                                                    aaaTATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTAt                                                                                                                                                                                                                                                             >  5:1249701/1‑149 (MQ=255)
                                                                                                                                       tatCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAGTAAAAAGGTCCGATGGGCATCGGACCTTTTATtgt                                                                                                                                                                                                                                                          <  5:1126440/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                                                                                                                                                                                                                                                     >  7:867967/1‑149 (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)
                                                                                                                                                      gTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCa                                                                                                                                                                                                                                           <  8:532627/149‑1 (MQ=255)
                                                                                                                                                           aGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTc                                                                                                                                                                                                                                      >  6:1459404/1‑149 (MQ=255)
                                                                                                                                                            gACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCt                                                                                                                                                                                                                                                                  >  5:512610/1‑120 (MQ=255)
                                                                                                                                                             aCGGTCTGAAAACCTGGACAGTGTTAACTACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGGCAGCCTCGc                                                                                                                                                                                                                                    <  8:1255766/149‑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                                                                                                                                                                                                                              <  7:960622/149‑1 (MQ=255)
                                                                                                                                                                       aaCCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTgg                                                                                                                                                                                                                          >  8:1081225/1‑149 (MQ=255)
                                                                                                                                                                           tGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGc                                                                                                                                                                                                                                             <  5:162518/126‑1 (MQ=255)
                                                                                                                                                                           tGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCAAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCCCTTCGATATTCAGCCCCTTCGGTTGCAGGATTTCGCCCAAACGGCCCAGGCCTTCAACCTGAACCCGGCGAACCCCTTCACGGGTGAGGCCaa                                                                                        >  4:270369/1‑275 (MQ=255)
                                                                                                                                                                            ggACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCttt                                                                                                                                                                                                                     <  7:927308/149‑1 (MQ=255)
                                                                                                                                                                             gACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTc                                                                                                                                                                                                                    <  7:235232/149‑1 (MQ=255)
                                                                                                                                                                             gACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTc                                                                                                                                                                                                                    <  6:498358/149‑1 (MQ=255)
                                                                                                                                                                             gACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTc                                                                                                                                                                                                                    <  8:804940/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)
                                                                                                                                                                                    ttAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGa                                                                                                                                                                                                             >  6:1394898/1‑149 (MQ=255)
                                                                                                                                                                                    ttAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGa                                                                                                                                                                                                             >  6:946660/1‑149 (MQ=255)
                                                                                                                                                                                      aaCGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCtt                                                                                                                                                                                                                                  <  6:1328245/126‑1 (MQ=255)
                                                                                                                                                                                      aaCGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATg                                                                                                                                                                                                           <  8:439332/149‑1 (MQ=255)
                                                                                                                                                                                       aCGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGc                                                                                                                                                                                                          >  6:81742/1‑149 (MQ=255)
                                                                                                                                                                                          aCCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGctta                                                                                                                                                                                                       <  7:329610/149‑1 (MQ=255)
                                                                                                                                                                                            cATTCTCGGTTTTACCGGCTTCTTAATTACCTGGTGCGTATGGGCGGTAATTTTACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAACAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGCCAGATGCTTACTTACTCGCGGAACAgcga                                                                                                                                                                                    >  2:226101/1‑165 (MQ=255)
                                                                                                                                                                                            cATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGcttact                                                                                                                                                                                                     >  8:177717/1‑149 (MQ=255)
                                                                                                                                                                                            cATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGcttact                                                                                                                                                                                                     >  6:892529/1‑149 (MQ=255)
                                                                                                                                                                                            cATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGc                                                                                                                >  4:64127/1‑234 (MQ=255)
                                                                                                                                                                                                 tcGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCTGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCTTTCTTACAGATTCTTACTTACTcg                                                                                                                                                                                               >  8:843163/1‑149 (MQ=255)
                                                                                                                                                                                                 tcGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTc                                                                                                                                                                                                >  6:659395/1‑149 (MQ=255)
                                                                                                                                                                                                    gTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTcgcg                                                                                                                                                                                             <  5:1324262/149‑1 (MQ=255)
                                                                                                                                                                                                               ttCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCtt                                                                                                                                                                                  <  6:104089/149‑1 (MQ=255)
                                                                                                                                                                                                                tCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTc                                                                                                                                                                                 <  7:461291/149‑1 (MQ=255)
                                                                                                                                                                                                                 cTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCg                                                                                                                                                                                <  5:461979/149‑1 (MQ=255)
                                                                                                                                                                                                                       tACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTc                                                                                                                                                                          <  6:422829/149‑1 (MQ=255)
                                                                                                                                                                                                                        aCCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTc                                                                                                                                                                                                                    >  2:362334/1‑106 (MQ=255)
                                                                                                                                                                                                                        aCCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCa                                                                                                                                                                         >  5:386870/1‑149 (MQ=255)
                                                                                                                                                                                                                         ccTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAg                                                                                                                                                                        <  8:359115/149‑1 (MQ=255)
                                                                                                                                                                                                                            ggTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGAc                                                                                                                                                                                                                <  4:362334/106‑1 (MQ=255)
                                                                                                                                                                                                                                 gTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTga                                                                                                                                                                                                                                >  7:561714/1‑85 (MQ=255)
                                                                                                                                                                                                                                 gTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTga                                                                                                                                                                                                                                <  8:561714/85‑1 (MQ=255)
                                                                                                                                                                                                                                   aTGGGCGGTAATTGGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGTACCCCAGCAAACCGAATCGa                                <  2:94785/275‑1 (MQ=255)
                                                                                                                                                                                                                                      ggCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGc                                                                                                                                                           <  7:600681/149‑1 (MQ=255)
                                                                                                                                                                                                                                       gCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCGTCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCa                                                                                                                                                          <  5:261745/149‑1 (MQ=255)
                                                                                                                                                                                                                                       gCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCa                                                                                                                                                          >  7:787349/1‑149 (MQ=255)
                                                                                                                                                                                                                                       gCGGTAATGTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCa                                                                                                                                                          <  7:405828/149‑1 (MQ=255)
                                                                                                                                                                                                                                        cGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAgcgc                                                                                                                                                                                    >  1:20347/1‑122 (MQ=255)
                                                                                                                                                                                                                                               ttGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTcg                                                                                                                                                  <  6:1249701/149‑1 (MQ=255)
                                                                                                                                                                                                                                                   ccTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCa                                                                                                                                              >  7:1195399/1‑149 (MQ=255)
                                                                                                                                                                                                                                                     ttAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCaaa                                                                                                                                            >  7:750854/1‑149 (MQ=255)
                                                                                                                                                                                                                                                     ttAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCaaa                                                                                                                                            >  5:95189/1‑149 (MQ=255)
                                                                                                                                                                                                                                                       aataaAAAGGTCCGCTGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTAGCTTTAGACTGGCCTTTCTGACCGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACTCGTTCACGGGTGAGGCCAATTTCACTCCCTAAATCTTCCAGTGTTGCTGCTTCTCTCCCCACACAACCGAATCCACGTCCACGCACttgcacct             >  4:227745/1‑268 (MQ=255)
                                                                                                                                                                                                                                                         taaAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGgc                                                                                                                                        >  7:1336153/1‑149 (MQ=255)
                                                                                                                                                                                                                                                          aaaaaGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGgcg                                                                                                                                       >  5:950733/1‑149 (MQ=255)
                                                                                                                                                                                                                                                             aaGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAg                                                                                                                                    <  6:674952/149‑1 (MQ=255)
                                                                                                                                                                                                                                                             aaGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAg                                                                                                                                    <  6:16125/149‑1 (MQ=255)
                                                                                                                                                                                                                                                              aGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAg                                                                                                                                                         <  5:1072/127‑1 (MQ=255)
                                                                                                                                                                                                                                                              aGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAg                                                                                                                                                         >  6:1072/1‑127 (MQ=255)
                                                                                                                                                                                                                                                              aGGTCCGATGCCCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAgg                                                                                                                                   >  5:1331110/1‑149 (MQ=255)
                                                                                                                                                                                                                                                              aGGTCCGATGCCCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTATGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAgg                                                                                                                                   >  8:101822/1‑149 (MQ=255)
                                                                                                                                                                                                                                                               ggTCCGATGCCCATCGGACCTTTTATTGTGCACAGAAAAGGCAAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGcttactta                                                                                                                                                                                                   >  6:283693/1‑84 (MQ=38)
                                                                                                                                                                                                                                                                gTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGcc                                                                                                                                 >  5:557367/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                 tCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGCTTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGTACCCCAGCAAACCGAATCGACGTGCCAGCACTTCACGCTGTTTTGAGTTc  >  2:171632/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                   cGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTc                                                                                                                              <  6:559618/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                     aTGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCaa                                                                                                                            <  7:177717/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                       gggCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAAcc                                                                                                                          >  5:642604/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                          
ACAACCGGCGGGTTACTCTCAGAAGCGGTGATGGGATTGAATCCGATTCAGTGCGTATTGGTGACGCTGGCAGCCTGCTTTGGTGGGCTTGGTGCGTCACATATTAATGACTCAGGGTTCTGGATTGTGACCAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGTACCCCAGCAAACCGAATCGACGTGCCAGCACTTCACGCTGTTTGGCGTTC  >  NC_000913/2866223‑2866754

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