Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 824,573 T→A 5.1% intergenic (‑76/‑57) ybhP ← / → ybhQ endo/exonuclease/phosphatase family protein/inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913824,5730TA 5.1% 303.0 / 4.3 97intergenic (‑76/‑57)ybhP/ybhQendo/exonuclease/phosphatase family protein/inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (52/40);  new base A (3/2);  total (55/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.13e-01

ACGCGGTAAAGCCTTTGTGAATATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGATTATGTTGCATTTACTGGTAGTGGCGCTGCTGGTGGTGG  >  NC_000913/824428‑824720
                                                                                                                                                 |                                                                                                                                                   
aCGCGGTAAAGCCTTTGTGAATATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTc                                                                                                                                                  >  1:1325067/1‑149 (MQ=255)
  gcgGTAAAGCCTTTGTGAATATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCTTGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCg                                                                                                                                                >  2:994619/1‑149 (MQ=255)
  gcgGTAAAGCCTTTGTGAATATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCg                                                                                                                                                >  2:989264/1‑149 (MQ=255)
   cgGTAAAGCCTTTGTGAATATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAAAATTCCGt                                                                                                                                               <  1:1098867/149‑1 (MQ=255)
          gCCTTTGTGAATATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTc                                                                                                                                        <  2:1227888/149‑1 (MQ=255)
           ccTTTGTGAATATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCg                                                                                                                                       >  2:823958/1‑149 (MQ=255)
                  gAATATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTActct                                                                                                                                <  1:684015/149‑1 (MQ=255)
                   aaTATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTg                                                                                                                               >  2:304117/1‑149 (MQ=255)
                     tatTGATGGTGAGCACCTTGAACGATACTTGTTGGGTTTGATCGGGCATAATTTTCCTGGCTTTGCCTCTTATCGCATTGAAATAGTGTAGTCGGCGGCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGag                                                                                                                             <  1:263045/149‑1 (MQ=255)
                       ttGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCAGTAAAGTTCGGTACTCTGagta                                                                                                                           >  1:1355853/1‑149 (MQ=255)
                        tGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGagtaa                                                                                                                          >  1:756784/1‑149 (MQ=255)
                           tGGGGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGagtaagta                                                                                                                       <  2:1094303/149‑1 (MQ=255)
                                        gAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCATTGCCTATTATCTCATTGAAATAGTGAAGTCGGCGTCACAAAAAGGTGAGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAAttctt                                                                                                          <  1:252172/149‑1 (MQ=255)
                                         aaCGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAAttcttc                                                                                                         <  1:1044350/149‑1 (MQ=255)
                                         aaCGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAAttcttc                                                                                                         >  2:920729/1‑149 (MQ=255)
                                          aCGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCa                                                                                                        <  2:1238960/149‑1 (MQ=255)
                                           cGAATATTGTTGGGTTTGATCTGGCATAATTTTCCTGTGTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAg                                                                                                       <  2:153712/149‑1 (MQ=255)
                                           cGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAg                                                                                                       <  1:1549572/149‑1 (MQ=255)
                                            gAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAgg                                                                                                      >  2:797124/1‑149 (MQ=255)
                                                  tgttgtGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGaa                                                                                                <  1:1546296/149‑1 (MQ=255)
                                                    ttgtGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGc                                                                                              <  1:521372/149‑1 (MQ=255)
                                                    ttgtGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGc                                                                                              <  2:1355853/149‑1 (MQ=255)
                                                      gtgtTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGATATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCa                                                                                            >  1:63171/1‑149 (MQ=255)
                                                      gtgtTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCa                                                                                            >  1:1521518/1‑149 (MQ=255)
                                                       tgtTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCAt                                                                                           <  1:1469784/149‑1 (MQ=255)
                                                       tgtTTGATCGCGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTCAGTAGAGATAAATTCTTCAGGAGAGAAGCCAt                                                                                           <  2:463320/149‑1 (MQ=255)
                                                        gtTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATg                                                                                          >  2:1321953/1‑149 (MQ=255)
                                                         tttGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGa                                                                                         >  2:534837/1‑149 (MQ=255)
                                                         tttGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGa                                                                                         >  2:528787/1‑149 (MQ=255)
                                                            gATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGt                                                                                      <  1:453325/149‑1 (MQ=255)
                                                             aTCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTg                                                                                     <  2:1172215/149‑1 (MQ=255)
                                                              tCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTgg                                                                                    >  1:1351696/1‑149 (MQ=255)
                                                              tCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTgg                                                                                    <  1:304117/149‑1 (MQ=255)
                                                                   cATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCTTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGcaaca                                                                               >  1:1413814/1‑149 (MQ=255)
                                                                    aTAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGcaacaa                                                                              >  2:30653/1‑149 (MQ=255)
                                                                     tAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCaacaac                                                                             <  1:407260/149‑1 (MQ=255)
                                                                     tAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCaacaac                                                                             <  2:1285985/149‑1 (MQ=255)
                                                                     tAATTTTCCTGTCTTTGACTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCaacaac                                                                             <  1:413033/149‑1 (MQ=255)
                                                                      aaTTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGTTACTCTGAGTAAGTAGAGATAAATTCTTCAGTAGAGAAGCCATGCAGATGCAACAACg                                                                            >  2:1100974/1‑149 (MQ=255)
                                                                      aaTTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACg                                                                            >  1:762130/1‑149 (MQ=255)
                                                                        ttttCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAATCCATGAAGTGGCAACAACgtg                                                                          >  2:1218122/1‑149 (MQ=255)
                                                                              tGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCgtg                                                                    >  1:546402/1‑149 (MQ=255)
                                                                               gTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGGACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGCGAAGCCATGAAGTGGCCACAGCGTGTTCgtgt                                                                   >  2:247160/1‑149 (MQ=255)
                                                                                  tttGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGc                                                                >  1:1285116/1‑149 (MQ=255)
                                                                                  tttGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGc                                                                >  1:238810/1‑149 (MQ=255)
                                                                                     gCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAAc                                                             >  2:619878/1‑149 (MQ=255)
                                                                                     gCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAAc                                                             >  1:1391629/1‑149 (MQ=255)
                                                                                      ccTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAGATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAAAGTGTTCGTGTCGCAACg                                                            >  1:255579/1‑149 (MQ=255)
                                                                                      ccTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAAGAACGTGTTCGTGTCGCAACg                                                            >  1:228329/1‑149 (MQ=255)
                                                                                       ctctTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGTAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACgg                                                           <  2:1146286/149‑1 (MQ=255)
                                                                                       ctctTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACgg                                                           <  2:1521518/149‑1 (MQ=255)
                                                                                       ctctTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACgg                                                           <  1:920729/149‑1 (MQ=255)
                                                                                       ctctTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACgg                                                           <  1:888289/149‑1 (MQ=255)
                                                                                       ctctTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACgg                                                           >  1:296530/1‑149 (MQ=255)
                                                                                        tctTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACACCGTGTTCGTGTCGCAACggg                                                          >  2:233272/1‑149 (MQ=255)
                                                                                        tctTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACggg                                                          <  1:774416/149‑1 (MQ=255)
                                                                                        tctTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACggg                                                          >  2:770490/1‑149 (MQ=255)
                                                                                         ctTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGt                                                         >  2:1585955/1‑149 (MQ=255)
                                                                                             tctcATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTaga                                                                                                                     >  1:1239119/1‑85 (MQ=255)
                                                                                             tctcATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTaga                                                                                                                     <  2:1239119/85‑1 (MQ=255)
                                                                                             tctcATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGGTAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTaa                                                     >  2:871892/1‑149 (MQ=255)
                                                                                             tctcATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTaa                                                     >  2:7584/1‑149 (MQ=255)
                                                                                              ctcATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAg                                                    >  2:28358/1‑149 (MQ=255)
                                                                                               tcATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTGTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGt                                                   >  2:1067680/1‑149 (MQ=255)
                                                                                               tcATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGt                                                   >  2:694977/1‑149 (MQ=255)
                                                                                               tcATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGt                                                   >  2:948271/1‑149 (MQ=255)
                                                                                               tcATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCCTGTCGCAACGGGTCTAAGt                                                   >  2:331508/1‑149 (MQ=255)
                                                                                                 aTTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGTTACTCTGAGTAAGTAGAGATAAATTCTTCAGTAGAGAAGACATGAAGTGGCAACAACGTGTTCGTGTCGAAACGGGTCTAAGTTg                                                 >  1:993766/1‑149 (MQ=255)
                                                                                                 aTTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTg                                                 >  2:1457815/1‑149 (MQ=255)
                                                                                                 aTTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTg                                                 >  2:1260647/1‑149 (MQ=255)
                                                                                                    gAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTg                                              <  2:238810/149‑1 (MQ=255)
                                                                                                     aaaTAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTgg                                             >  1:904884/1‑149 (MQ=255)
                                                                                                          gtgtAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGAt                                        <  1:1406650/149‑1 (MQ=255)
                                                                                                           tgtcggcGGCGTCACACAAACGTGCGGTCTTACGGAACTTTCCGTAAAGTTCGGTACTCTGAGTCAGGAGAGAGAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGAtt                                       <  1:247160/143‑1 (MQ=255)
                                                                                                               gTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGATTATGt                                   >  2:961451/1‑149 (MQ=255)
                                                                                                                tCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGATTATGtt                                  <  1:684236/149‑1 (MQ=255)
                                                                                                                  ggCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGATTATGTTGc                                <  1:1093827/149‑1 (MQ=255)
                                                                                                                        acaAAAAGGTGCGGTCTTCCGGATTTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGATTATGTTGCATTTAc                          <  2:328296/149‑1 (MQ=255)
                                                                                                                         caAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGATTATGTTGCATTTACt                         >  2:1467920/1‑149 (MQ=255)
                                                                                                                          aaaaaGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGATTATGTTGCATTTACTg                        <  2:63171/149‑1 (MQ=255)
                                                                                                                           aaaaGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGATTATGTTGCATTTACTgg                       <  1:797124/149‑1 (MQ=255)
                                                                                                                           aaaaGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGATTATGTTGCATTTACTgg                       >  1:596042/1‑149 (MQ=255)
                                                                                                                             caGGGGCGGGCTTACGGAATTTTCCGTAAAGTGCGGTACTCTGAGTAAGTAGAGATAAATGCTTCAGGAGAGAAGCCAGGCAGTGGCAACAACGTGTGCGTGTCGCCACGGGTCTAAGTTGCTGGCAGATTATGTTGCATTTACTGGTa                     <  1:989264/148‑1 (MQ=255)
                                                                                                                             aaGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGGGGCAACAACGTGTTCGTGGCGCAACGGGTCTAAGTTGCTGGCAGATTATGTTGCATTTACTGGTa                     <  1:994619/149‑1 (MQ=255)
                                                                                                                              aGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTATTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGATTATGTTGCATTTACTGGTAg                    <  1:823958/149‑1 (MQ=255)
                                                                                                                                gTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGATTATGTTGCATTTACTGGTAGTg                  >  1:125309/1‑149 (MQ=255)
                                                                                                                                   cGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGATTATGTTGCATTTACTGGTAGTGgcg               <  2:1386148/149‑1 (MQ=255)
                                                                                                                                    aactttACGGAAAATTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGATTATGTTGCATTTACTGGTAGTGgcgc              >  1:1590715/5‑149 (MQ=255)
                                                                                                                                     gTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGATTATGTTGCATTTACTGGTAGTGGCgct             <  2:546402/149‑1 (MQ=255)
                                                                                                                                     actttACGGAAAATTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGATTATGTTGCATTTACTGGTAGTGGCgct             <  2:1590715/146‑1 (MQ=255)
                                                                                                                                       tttACGGAAAATTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGATTATGTTGCATTTACTGGTAGTGGCgctgc           >  1:1212378/2‑149 (MQ=255)
                                                                                                                                       tttACGGAAAATTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGAAACGGGTCTAAGTTGCTGGCAGATTATGTTGCATTTACTGGTAGTGGCgctgc           >  1:57211/2‑149 (MQ=255)
                                                                                                                                        ttACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCGATGAAGTGGCAACAACGTGTTCGTGTCCCAACGGGTCTAAGTTGCTGGCAGATTATATTGCATTTACTGGTAGAGGCgctgat          >  1:334191/1‑147 (MQ=255)
                                                                                                                                           cGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGTCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTTCTGGCAGATTATGTTGCATTTACTGGTAGTGGCGCTGCtggt       >  1:767435/1‑149 (MQ=255)
                                                                                                                                            ggAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGATTATGTTGCATTTACTGGTAGTGGCGCTGCtggtg      >  1:1439613/1‑149 (MQ=255)
                                                                                                                                              aaTTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGATTATGTTGCATTTACTGGTAGTGGCGCTGCtggtggt    <  2:296530/149‑1 (MQ=255)
                                                                                                                                               aTTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGATTATGTTGCATTTACTGGTAGTGGCGCTGCtggtggtg   >  2:597483/1‑149 (MQ=255)
                                                                                                                                                ttttCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGATTATGTTGCATTTACTGGTAGTGGCGCTGCtggtgg     >  2:1260795/1‑146 (MQ=255)
                                                                                                                                                ttttCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGATTATGTTGCATTTACTGGTAGTGGCGCTGCtggtgg     <  1:1260795/146‑1 (MQ=255)
                                                                                                                                                ttttCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGATTATGTTGCATTTACTGGGAGTGGCGCTGCtggtggtgg  <  1:770490/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
ACGCGGTAAAGCCTTTGTGAATATTGATGGTGAGCACCTTGAACGAAAATTGTTGTGTTTGATCGGGCATAATTTTCCTGTCTTTGCCTCTTATCTCATTGAAATAGTGTAGTCGGCGTCACAAAAAGGTGCGGTCTTACGGAATTTTCCGTAAAGTTCGGTACTCTGAGTAAGTAGAGATAAATTCTTCAGGAGAGAAGCCATGAAGTGGCAACAACGTGTTCGTGTCGCAACGGGTCTAAGTTGCTGGCAGATTATGTTGCATTTACTGGTAGTGGCGCTGCTGGTGGTGG  >  NC_000913/824428‑824720

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