Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,111,083 T→G 5.8% E15A (GAA→GCA) ‡ yghD ← putative membrane‑anchored secretion pathway M‑type protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,111,0830TG 5.8% 628.3 / 24.7 172E15A (GAA→GCA) ‡yghDputative membrane‑anchored secretion pathway M‑type protein
Reads supporting (aligned to +/- strand):  ref base T (78/84);  new base G (5/5);  total (83/89)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.29e-01

Only 100 of 172 total aligned reads displayed.

TCCGTTCGCTCAACGGCTGCCAGATGAGGACGTAATATCCGACGCTAAACAGGAACACCACCGCTGCCAGTAACATGCCCTTTTCACGCGGAGAACGCCCCGCCAGGTGTTGTGCCAGCCAGTGTTCGCCACGGCTTAACTGGCGTTCACGCCATTGCTGAAAATAGTGAATAAATTTATCGCGTAACATGTTATTTCCTCCGCAACGTTACGCCGCCGGAAACCGCATCACCCTCTTTCTGTAACGCGTCCTGTTGCACAACATAATCGGTTGCCAGTACACT  >  NC_000913/3110937‑3111220
                                                                                                                                                  |                                                                                                                                         
tCCGTTCGCTCAACGGCTGCCAGATGAGGACGTAATATCCGCCGCTAAACAGGAACACCACCGCTGCCAGTAACATGCCCTTTTCACGCGGAGAACGCCCCGCCAGGTGTTGTGCCAGCCAGTGTTCGCCACGGCTTAACTGGTGTTCa                                                                                                                                         <  2:1224531/149‑1 (MQ=255)
 ccGTTCGCTCAACGGCTGCCAGATGAGGACGTAATATCCGACGCTAAACAGGAACACCACCGCTGCCAGTAACATGCCCTTTTCACGCGGAGAACGCCCCGCCAGGTGTTGTGCCAGCCAGTGTTCGCCACGGCTTAACTGGCGTTCAc                                                                                                                                        >  3:478169/1‑149 (MQ=255)
       gCTCAACGGCTGCCAGATGAGGACGTAATATCCGACGCTAAACAGGAACACCACCGCTGCCAGTAACATGCCCTTTTCACGCGGAGAACGCCCCGCCAGGTGTTGTGCCAGCCAGTGTTCGCCACGGCTTAACTGGCGTTCACGCCAtt                                                                                                                                  >  4:564845/1‑149 (MQ=255)
        cTCAACGGCTGCCAGATGAGGACGTAATATCCGACGCTAAACAGGAACACCACCGCTGCCAGTAACATGCCCTTTTCACGCGGAGAACGCCCCGCCAGGTGTTGTGCCAGCCAGTGTTCGCCACGGCTTAACTGGCGTTCACGCCATTg                                                                                                                                 <  3:692982/149‑1 (MQ=255)
         tCAACGGCTGCCAGATGAGGACGTAATATCCGACGCTAAACAGGAACACCACCGCTGCCAGTAACATGCCCTTTTCACGCGGAGAACGCCCCGCCAGGTGTTGTGCCAGCCAGTGTTCGCCACGGCTTAACTGGCGTTCACGCCATTGc                                                                                                                                >  1:964938/1‑149 (MQ=255)
         tCAACGGCTGCCAGATGAGGACGTAATATCCGACGCTAAACAGGAACACCACCGCTGCCAGTAACATGCCCTTTTCACGCGGAGAACGCCCCGCCAGGTGTTGTGCCAGCCAGTGTTCGCCACGGCTTAACTGGCGTGAACGc                                                                                                                                      <  2:134599/143‑1 (MQ=255)
         tCAACGGCTGCCAGATGAGGACGTAATATCCGACGCTAAACAGGAACACCACCGCTGCCAGTAACATGCCCTTTTCACGCGGAGAACGCCCCGCCAGGTGTTGTGCCAGCCAGTGTTCGCCACGGCTTAACTGGCGTGAACGCCagtt                                                                                                                                 <  1:1319833/148‑4 (MQ=255)
           aaCGGCTGCCAGATGAGGACGTAATATCCGACGCTAAACAGGAACACCACCGCTGCCAGTAACATGCCCTTTTCACGCGGAGAACGCCCCGCCAGGTGTTGTGCCAGCCAGTGTTCGCCACGGCTTAACTGGCGTTCACGCCATTGCTg                                                                                                                              <  2:964938/149‑1 (MQ=255)
              ggCTGCCAGATGAGGACGTAATATCCGACGCTAAACAGGAACACCACCGCTGCCAGTAACATGCCCTTTTCACGCGGAGAACGCCCCGCCAGGTGTTGTGCCAGCCAGTGTTCGCCACGGCTTAACTGGCGTTCa                                                                                                                                         <  3:771992/135‑1 (MQ=255)
               gCTGCCAGATGAGGACGTAATATCCGACGCTAAACAGGAACACCACCGCTGCCAGTAACATGCCCTTTTCACGCGGAGAACGCCCCGCCAGGTGTTGTGCCAGCCAGTGTTCGCCACGGCTTAACTGGCGTTCACGCCATTGCTGaaaa                                                                                                                          <  4:1383581/149‑1 (MQ=255)
                   ccAGATGAGGACGTAATATCCGACGCTAAACAGGAACACCACCGCTGCCAGTAACATGCCCTTTTCACGCGGAGAACGCCCCGCCAGGTGTTGTGCCAGCCAGTGTTCGCCACGGCTTAACTGGCGTTCACGCCATTGCTGAAAATAGt                                                                                                                      <  3:1091365/149‑1 (MQ=255)
                    cAGATGAGGACGTAATATCCGACGCTAAACAGGAACACCACCGCTGCCAGTAACATGCCCTTTTCACGCGGAGAACGCCCCGCCAGGTGTTGTGCCAGCCAGTGTTCGCCACGGCTTAACTGGCGTGAACGCCagt                                                                                                                                  >  4:177927/1‑134 (MQ=255)
                      gATGAGGACGTAATATCCGACGCTAAACAGGAACACCACCGCTGCCAGTAACATGCCCTTTTCACGCGGAGAACGCCCCGCCAGGTGTTGTGCCAGCCAGTGTTCGCCACGGCTTAACTGGCGTTc                                                                                                                                          <  2:690689/126‑1 (MQ=255)
                         gAGGACGTAATATCCGACGCTAAACAGGAACACCACCGCTGCCAGTAACATGCCCTTTTCACGCGGAGAACGCCCCGCCAGGTGTTGTGCCAGCCAGTGTTCGCCACGGCTTAACTGGCGTTCACGCCATTGCTGAAAATAGTGAATa                                                                                                                 <  2:302678/148‑1 (MQ=255)
                          aGGACGTAATATCCGACGCTAAACAGGAACACCACCGCTGCCAGTAACATGCCCTTTTCACGCGGAGAACGCCCCGCCAGGTGTTGTGCCAGCCAGTGTTCGCCACGGCTTAACTGGCGTTCACGCCATTGCTGAAAATAGTGAATaaa                                                                                                               <  3:138491/149‑1 (MQ=255)
                           ggACGTAATATCCGACGCTAAACAGGAACACCACCGCTGCCAGTAACATGCCCTTTTCACGCGGAGAACGCCCCGCCAGGTGTTGTGCCAGCCAGTGTTCGCCACGGCTTAACTGGCGTTCACGCCATTGCTGAAAATAGTGAATAAAt                                                                                                              >  2:1430528/1‑149 (MQ=255)
                              cGTAATATCCGACGCTAAACAGGAACACCACCGCTGCCAGTAACATGCCCTTTTCACGCGGAGAACGCCCCGCCAGGTGTTGTGCCAGCCAGTGTTCGCCACGGCTTAACTGGCGTTCACGCCATTGCTGAAAATAGTGAATAAATTTa                                                                                                           >  1:1008641/1‑149 (MQ=255)
                               gTAATATCCGACGCTAAACAGGAACACCACCGCTGCCAGTAACATGCCCTTTTCACGCGGAGAACGCCCCGCCAGGTGTTGTGCCAGCCAGTGTTCGCCACGGCTTAACTGGCGTTCACGCCATTGCTGAAAATAGTGAATAAATTTAt                                                                                                          <  3:1215973/149‑1 (MQ=255)
                                 aaTATCCGACGCTAAACAGGAACACCACCGCTGCCAGTAACATGCCCTTTTCACGCGGAGAACGCCCCGCCAGGTGTTGTGCCAGCCAGTGTTCGCCACGGCTTAACTGGCGTTCACGCCATTGCTGAAAATAGTGAATAAAt                                                                                                              >  1:147737/1‑143 (MQ=255)
                                 aaTATCCGACGCTAAACAGGAACACCACCGCTGCCAGTAACATGCCCTTTTCACGCGGAGAACGCCCCGCCAGGTGTTGTGCCAGCCAGTGTTCGCCACGGCTTAACTGGCGTTCACGCCATTGCTGAAAATAGTGAATAAAt                                                                                                              <  2:147737/143‑1 (MQ=255)
                                   tatCCGACGCTAAACAGGAACACCACCGCTGCCAGTAACATGCCCTTTTCACGCGGAGAACGCCCCGCCAGGTGTTGTGCCAGCCAGTGTTCGCCACGGCTTAACTGGCGTTCACGCCATTGCTGAAAATAGt                                                                                                                      <  2:529769/133‑1 (MQ=255)
                                   tatCCGACGCTAAACAGGAACACCACCGCTGCCAGTAACATGCCCTTTTCACGCGGAGAACGCCCCGCCAGGTGTTGTGCCAGCCAGTGTTCGCCACGGCTTAACTGGCGTTCACGCCATTGCTGAAAATAGt                                                                                                                      >  1:529769/1‑133 (MQ=255)
                                     tCCGACGCTAAACAGGAACACCACCGCTGCCAGTAACATGCCCTTTTCACGCGGAGAACGCCCCGCCAGGTGTTGTGCCAGCCAGTGTTCGCCACGGCTTAACTGGCGTTCACGCCATTGCTGAAAATAGTGAATAAATTTATCGCGTa                                                                                                    >  4:269791/1‑149 (MQ=255)
                                        gACGCTAAACAGGAACACCACCGCTGCCAGTAACATGCCCTTTTCACGCGGAGAACGCCCCGCCAGGTGTTGTGCCAGCCAGTGTTCGCCACGGCTTAACTGGCGTTCACGCCATTGCTGAAAATAGTGAATAAATTTATCGCGTAAc                                                                                                  <  1:1026277/148‑1 (MQ=255)
                                        gACGCTAAACAGGAACACCACCGCTGCCAGTAACATGCCCTTTTCACGCGGAGAACGCCCCGCCAGGTGTTGTGCCAGCCAGTGTTCGCCACGGCTTAACTGGCGTTCACGCCATTGCTGAAAATAGTGAATAAATTTATCGCGTAAc                                                                                                  <  2:587930/148‑1 (MQ=255)
                                        gACGCTAAACAGGAACACCACCGCTGCCAGTAACATGCCCTTTTCACGCGGAGAACGCCCCGCCAGGTGTTGTGCCAGCCAGTGTTCGCCACGGCTTAACTGGCGTTCACGCCATTGCTGAAAATAGTGAATAAATTTATCGCGTAAc                                                                                                  >  2:1026277/1‑148 (MQ=255)
                                        gACGCTAAACAGGAACACCACCGCTGCCAGTAACATGCCCTTTTCACGCGGAGAACGCCCCGCCAGGTGTTGTGCCAGCCAGTGTTCGCCACGGCTTAACTGGCGTTCACGCCATTGCTGAAAATAGTGAATAAATTTATCGCGTAAc                                                                                                  >  1:587930/1‑148 (MQ=255)
                                        gACGCTAAACAGGAACACCACCGCTGCCAGTAACATGCCCTTTTCACGCGGAGAACGCCCCGCCAGGTGTTGTGCCAGCCAGTGTTCGCCACGGCTTAACTGGCGTTCACGCCATTGCTGAAAATAGTGAATAAATTTATCGCGTAACa                                                                                                 >  1:458016/1‑149 (MQ=255)
                                        gACGCTAAACAGGAACACCACCGCTGCCAGTAACATGCCCTTTTCACGCGGAGAACGCCCCGCCAGGTGTTGTGCCAGCCAGTGTTCGCCACGGCTTAACTGGCGTTCACGCCATTGCTGAAAATAGTGAATAAATTTATCGCGTAACa                                                                                                 >  1:1299912/1‑149 (MQ=255)
                                           gCTAAACAGGAACACCACCGCTGCCAGTAACATGCCCTTTTCACGCGGAGAACGCCCCGCCAGGTGTTGTGCCAGCCAGTGTTCGCCACGGCTTAACTGGCGTTCACGCCATTGCTGAAAATAGTGAATAAATTTATCGCGTAACATGt                                                                                              <  1:556745/149‑1 (MQ=255)
                                           gCTAAACAGGAACACCACCGCTGCCAGTAACATGCCCTTTTCACGCGGAGAACGCCCCGCCAGGTGTTGTGCCAGCCAGTGTTCGCCACGGCTTAACTGGCGTGAACGCCagttaa                                                                                                                               >  1:678846/1‑111 (MQ=255)
                                            cTAAACAGGAACACCACCGCTGCCAGTAACATGCCCTTTTCACGCGGAGAACGCCCCGCCAGGTGTTGTGCCAGCCAGTGTTCGCCACGGCTTAACTGGCGTTCACGCCATTGCTGAAAATAGTGAATAAATTTATCGCGTAACATGtt                                                                                             <  1:178306/149‑1 (MQ=255)
                                               aaCAGGAACACCACCGCTGCCAGTAACATGCCCTTTTCACGCGGAGAACGCCCCGCCAGGTGTTGTGCCAGCCAGTGTTCGCCACGGCTTAACTGGCGTTCACGCCATTGCTGAAAATAg                                                                                                                       <  2:830168/120‑1 (MQ=255)
                                               aaCAGGAACACCACCGCTGCCAGTAACATGCCCTTTTCACGCGGAGAACGCCCCGCCAGGTGTTGTGCCAGCCAGTGTTCGCCACGGCTTAACTGGCGTTCACGCCATTGCTGAAAATAGt                                                                                                                      >  3:942192/1‑121 (MQ=255)
                                               aaCAGGAACACCACCGCTGCCAGTAACATGCCCTTTTCACGCGGAGAACGCCCCGCCAGGTGTTGTGCCAGCCAGTGTTCGCCACGGCTTAACTGGCGTTCACGCCATTGCTGAAAATAGTGAATAAATTTATCGCGTAACATGTTAtt                                                                                          >  4:142158/1‑149 (MQ=255)
                                               aaCAGGAACACCACCGCTGCCAGTAACATGCCCTTTTCACGCGGAGAACGCCCCGCCAGGTGTTGTGCCAGCCAGTGTTCGCCACGGCTTAACTGGCGTTCACGCCATTGCTGAAAATAGTGAATAAATTTATCGCGTAACATGTTAtt                                                                                          <  2:44373/149‑1 (MQ=255)
                                                aCAGGAACACCACCGCTGCCAGTAACATGCCCTTTTCACGCGGAGAACGCCCCGCCAGGTGTTGTGCCAGCCAGTGTTCGCCACGGCTTAACTGGCGTTCACGCCATTGCTGAAAATAGTGAATaa                                                                                                                >  2:921784/1‑126 (MQ=255)
                                                aCAGGAACACCACCGCTGCCAGTAACATGCCCTTTTCACGCGGAGAACGCCCCGCCAGGTGTTGTGCCAGCCAGTGTTCGCCACGGCTTAACTGGCGTTCACGCCATTGCTGAAAATAGTGAATAAATTTATCGCGTAACATGTTAttt                                                                                         <  4:1370253/149‑1 (MQ=255)
                                                     aaCACCACCGCTGCCAGTAACATGCCCTTTTCACGCGGAGAACGCCCCGCCAGGTGTTGTGCCAGCCAGTGTTCGCCACGGCTTAACTGGCGTTCACGCCATTGCTGAAAATAGt                                                                                                                      <  2:281586/115‑1 (MQ=255)
                                                      acacCACCGCTGCCAGTAACATGCCCTTTTCACGCGGAGAACGCCCCGCCAGGTGTTGTGCCAGCCAGTGTTCGCCACGGCTTAACTGGCGTTCACGCCATTGCTGAAAATAGTGAATAAATTTATCGCGTAACATGTTATTTCCTCCg                                                                                   <  2:458016/149‑1 (MQ=255)
                                                        accaccGCTGCCAGTAACATGCCCTTTTCACGCGGAGAACGCCCCGCCAGGTGTTGTGCCAGCCAGTGTTCGCCACGGCTTAACTGGCGTTCACGCCATTGCTGAAAATAGt                                                                                                                      <  2:1134972/112‑1 (MQ=255)
                                                        accaccGCTGCCAGTAACATGCCCTTTTCACGCGGAGAACGCCCCGCCAGGTGTTGTGCCAGCCAGTGTTCGCCACGGCTTAACTGGCGTTCACGCCATTGCTGAAAATAGt                                                                                                                      >  1:1134972/1‑112 (MQ=255)
                                                        accaccGCTGCCAGTAACATGCCCTTTTCACGCGGAGAACGCCCCGCCAGGTGTTGTGCCAGCCAGTGTTCGCCACGGCTTAACTGGCGTTCACGCCATTGCTGAAAATAGTGAATAAATTTa                                                                                                           <  4:384451/123‑1 (MQ=255)
                                                        accaccGCTGCCAGTAACATGCCCTTTTCACGCGGAGAACGCCCCGCCAGGTGTTGTGCCAGCCAGTGTTCGCCACGGCTTAACTGGCGTTCACGCCATTGCTGAAAATAGTGAATAAATTTATCGCGTAACATGTTATTTCCTCCGCa                                                                                 >  1:1248290/1‑149 (MQ=255)
                                                         ccaccGCTGCCAGTAACATGCCCTTTTCACGCGGAGAACGCCCCGCCAGGTGTTGTGCCAGCCAGTGTTCGCCACGGCTTAACTGGCGTTCACGCCATTGCTGAAAATAGTGAATaaa                                                                                                               >  3:1440453/1‑118 (MQ=255)
                                                         ccaccGCTGCCAGTAACATGCCCTTTTCACGCGGAGAACGCCCCGCCAGGTGTTGTGCCAGCCAGTGTTCGCCACGGCTTAACTGGCGTTCACGCCATTGCTGAAAATAGTGAATAAATTTATCGCGTAACAt                                                                                                <  1:1115367/133‑1 (MQ=255)
                                                         ccaccGCTGCCAGTAACATGCCCTTTTCACGCGGAGAACGCCCCGCCAGGTGTTGTGCCAGCCAGTGTTCGCCACGGCTTAACTGGCGTTCACGCCATTGCTGAAAATAGTGAATAAATTTATCGCGTAACATGTTATTTCCTCCGCaa                                                                                >  2:24364/1‑149 (MQ=255)
                                                          caccGCTGCCAGTAACATGCCCTTTTCACGCGGAGAACGCCCCGCCAGGTGTTGTGCCAGCCAGTGTTCGCCACGGCTTAACTGGCGTTCACGCCATTGCTGAAAATAGTGAATAAATTTATCGCGTAACATGTTATTTCCTCCGCAAc                                                                               <  3:269791/149‑1 (MQ=255)
                                                           accGCTGCCAGTAACATGCCCTTTTCACGCGGAGAACGCCCCGCCAGGTGTTGTGCCAGCCAGTGTTCGCCACGGCTTAACTGGCGTTCACGCCATTGCTGAAAATAGTGAATAAATTTATCGCGTAACATGtt                                                                                             <  4:429193/134‑1 (MQ=255)
                                                                  ccAGTAACATGCCCTTTTCACGCGGAGAACGCCCCGCCAGGTGTTGTGCCAGCCAGTGTTCGCCACGGCTTAACTGGCGTTCACGCCATTGCTGAAAATAGTg                                                                                                                     >  1:436383/1‑103 (MQ=255)
                                                                    aGTAACATGCCCTTTTCACGCGGAGAACGCCCCGCCAGGTGTTGTGCCAGCCAGTGTTCGCCACGGCTTAACTGGCGTTCACGCCATTGCTGAAAATAGTGAATAAATTTATCGCGTAACATGTTATTtc                                                                                        <  2:365739/130‑1 (MQ=255)
                                                                    aGTAACATGCCCTTTTCACGCGGAGAACGCCCCGCCAGGTGTTGTGCCAGCCAGTGTTCGCCACGGCTTAACTGGCGTTCACGCCATTGCTGAAAATAGTGAATAAATTTATCGCGTAACATGTTATTtc                                                                                        >  1:365739/1‑130 (MQ=255)
                                                                     gTAACATGCCCTTTTCACGCGGAGAACGCCCCGCCAGGTGTTGTGCCAGCCAGTGTTCGCCACGGCTTAACTGGCGTTCACGCCATTGCTGAAAATAGTGAATAAATTTATCGCGTAACATGTTAt                                                                                           >  3:1361321/1‑126 (MQ=255)
                                                                      tAACATGCCCTTTTCACGCGGAGAACGCCCCGCCAGGTGTTGTGCCAGCCAGTGTTCGCCACGGCTTAACTGGCGTTCACGCCATTGCTGAAAATAGTGAATAAATTTATCGCGTAACATGTTATTTCCTCCGCAACGTTACGccgccg                                                                   >  4:1293621/1‑149 (MQ=255)
                                                                      tAACATGCCCTTTTCACGCGGAGAACGCCCCGCCAGGTGTTGTGCCAGCCAGTGTTCGCCACGGCTTAACTGGCGTTCACGCCATTGCTGAAAATAGTGAATAAATTTATCGCGTAACATGTTATTTCCTCCGCAACGTTACGccgccg                                                                   <  1:1320216/149‑1 (MQ=255)
                                                                          aTGCCCTTTTCACGCGGAGAACGCCCCGCCAGGTGTTGTGCCAGCCAGTGTTCGCCACGGCTTAACTGGCGTTCACGCCATTGCTGAAAATAGTGAATAAATTTATCGCGTAACATGTTATTTCCTCCGCAACGTTACGCCGCCGGaaa                                                               >  1:819074/1‑149 (MQ=255)
                                                                           tGCCCTTTTCACGCGGAGAACGCCCCGCCAGGTGTTGTGCCAGCCAGTGTTCGCCACGGCTTAACTGGCGTTCACGCCATTGCTGAAAATAGTGAATAAATTTATCGCGTAACATGTTATTTCCTCCGCAACGTTACGCCGCCGGAAAc                                                              <  3:1293621/149‑1 (MQ=255)
                                                                            gCCCTTTTCACGCGGAGAACGCCCCGCCAGGTGTTGTGCCAGCCAGTGTTCGCCACGGCTTAACTGGCGTTCa                                                                                                                                         >  1:583426/1‑73 (MQ=255)
                                                                            gCCCTTTTCACGCGGAGAACGCCCCGCCAGGTGTTGTGCCAGCCAGTGTTCGCCACGGCTTAACTGGCGTTCACGCCATTGCTGAAAATAGTGAATAAATTTATCGCGTAACATGTTATTTCCTCCGCAACGTTACGCCGCCGGAAAcc                                                             >  4:1241379/1‑149 (MQ=255)
                                                                             cccTTTTCACGCGGAGAACGCCCCGCCAGGTGTTGTGCCAGCCAGTGTTCGCCACGGCTTAACTGGCGTTCACGCCATTGCTGAAAATAGTGAATAAATTTATCGCGTAACATGTTATTTCCTCCGCAACGTTACGCCGCCGGAAACCg                                                            <  1:1430528/149‑1 (MQ=255)
                                                                              ccTTTTCACGCGGAGAACGCCCCGCCAGGTGTTGTGCCAGCCAGTGTTCGCCACGGCTTAACTGGCGTTCACGCCATTGCTGAAAATAGTGAATAAAt                                                                                                              >  3:941836/1‑98 (MQ=255)
                                                                              ccTTTTCACGCGGAGAACGCCCCGCCAGGTGTTGTGCCAGCCAGTGTTCGCCACGGCTTAACTGGCGTTCACGCCATTGCTGAAAATAGTGAATAAAt                                                                                                              <  4:941836/98‑1 (MQ=255)
                                                                              ccTTTTCACGCGGAGAACGCCCCGCCAGGTGTTGTGCCAGCCAGTGTTCGCCACGGCTTAACTGGCGTTCACGCCATTGCTGAAAATAGTGAATAAATTTATCGCGTAACATGTTATTTCCTCCGCAACGTTACGCCGCCGGAAACCGc                                                           >  4:645589/1‑149 (MQ=255)
                                                                              ccTTTTCACGCGGAGAACGCCCCGCCAGGTGTTGTGCCAGCCAGTGTTCGCCACGGCTTAACTGGCGTTCACGCCATTGCTGAAAATAGTGAATAAATTTATCGCGTAACATGTTATTTCCTCCGCAACGTTACGCCGCCGGAAACCGc                                                           <  2:819074/149‑1 (MQ=255)
                                                                               cTTTTCACGCGGAGAACGCCCCGCCAGGTGTTGTGCCAGCCAGTGTTCGCCACGGCTTAACTGGCGTTCACGCCATTGCTGAAAATAGTGAATAAATTTATCGCGTAACATGTTATTtcctcc                                                                                    <  3:803738/123‑1 (MQ=255)
                                                                                       gcgGAGAACGCCCCGCCAGGTGTTGTGCCAGCCAGTGTTCGCCACGGCTTAACTGGCGTTCACGCCATTGCTGAAAATAGTGAATAAATTTAt                                                                                                          <  1:323669/93‑1 (MQ=255)
                                                                                        cgGAGAACGCCCCGCCAGGTGTTGTGCCAGCCAGTGTTCGCCACGGCTTAACTGGCGTTCACGCCATTGCTGAAAATAGTGAATAAATTTATCGCGTAACATGTTATTTCCTCCGCAACGTTACGCCGCCGGAAACCGCATCACCctct                                                 <  1:921677/149‑1 (MQ=255)
                                                                                         ggAGAACGCCCCGCCAGGTGTTGTGCCAGCCAGTGTTCGCCACGGCTTAACTGGCGTTCACGCCATTGCTGAAAATAGTGAATAAATTTATCGCGTAACATGtt                                                                                             <  4:462383/104‑1 (MQ=255)
                                                                                         ggAGAACGCCCCGCCAGGTGTTGTGCCAGCCAGTGTTCGCCACGGCTTAACTGGCGTTCACGCCATTGCTGAAAATAGTGAATAAATTTATCGCGTAACATGtt                                                                                             >  3:462383/1‑104 (MQ=255)
                                                                                          gagaACGCCCCGCCAGGTGTTGTGCCAGCCAGTGTTCGCCACGGCTTAACTGGCGTTCACGCCATTGCTGAAAATAGTGAATAAATTTATCGCGTAACATGTTATTTCCTCCGCAACGTTACGCCGCCGGAAACCGCATCACCCTCttt                                               <  4:478169/149‑1 (MQ=255)
                                                                                              aCGCCCCGCCAGGTGTTGTGCCAGCCAGTGTTCGCCACGGCTTAACTGGCGTTCACGCCATTGCTGAAAATAGTGAATAAATTTATCGCGTAACATGTTATTTCCTCCGCAACGTTACGCCGCCGGAAACCGCATCACCCTCTTTCTGt                                           <  2:1248290/149‑1 (MQ=255)
                                                                                               cGCCCCGCCAGGTGTTGTGCCAGCCAGTGTTCGCCACGGCTTAACTGGCGTTCACGCCATTGCTGAAAATAGTGAATAAATTTATCGCGTAACATGTTATTTCCTCCGCAACGTTACGCCGCCGGAAACCGCATCACCCTCTTTCTGTa                                          <  2:1099028/149‑1 (MQ=255)
                                                                                                gCCCCGCCAGGTGTTGTGCCAGCCAGTGTTCGCCACGGCTTAACTGGCGTTCACGCCATTGCTGAAAATAGTGAATaaa                                                                                                               >  2:697560/1‑79 (MQ=255)
                                                                                                gCCCCGCCAGGTGTTGTGCCAGCCAGTGTTCGCCACGGCTTAACTGGCGTTCACGCCATTGCTGAAAATAGTGAATAAATTTATCGCGTAACATGTTATTTCCTCCGCAACGTTACGCCGCCGGAAACCGCATCACCCTCTTTCTGTaa                                         >  3:185068/1‑149 (MQ=255)
                                                                                                gCCCCGCCAGGTGTTGTGCCAGCCAGTGTTCGCCACGGCTTAACTGGCGTTCACGCCATTGCTGAAAATAGTGAATAAATTTATCGCGTAACATGTTATTTCCTCCGCAACGTTACGCCGCCGGAAACCGCATCACCCTCTTTCTGTaa                                         <  1:302660/149‑1 (MQ=255)
                                                                                                  cccGCCAGGTGTTGTGCCAGCCAGTGTTCGCCACGGCTTAACTGGCGTTCACGCCATTGCTGAAAATAGTGAATAAATTTa                                                                                                           >  3:347644/1‑81 (MQ=255)
                                                                                                  cccGCCAGGTGTTGTGCCAGCCAGTGTTCGCCACGGCTTAACTGGCGTTCACGCCATTGCTGAAAATAGTGAATAAATTTa                                                                                                           <  4:347644/81‑1 (MQ=255)
                                                                                                   ccGCCAGGTGTTGTGCCAGCCAGTGTTCGCCACGGCTTAACTGGCGTTCACGCCATTGCTGAAAATAGTGAATAAATTTATCGCGTAACATGTTATTTCCTCCGCAACGTTACGCCGCCGGAAACCGCATCACCCTCTTTCTGTAAcgc                                      <  4:185068/149‑1 (MQ=255)
                                                                                                     gCCAGGTGTTGTGCCAGCCAGTGTTCGCCACGGCTTAACTGGCGTTCACGCCATTGCTGAAAATAGTGAATAAATTTATCGCGTAACATGTTATTTCCTCCGCAACGTTACGCCGCCGGAAACCGCATCACCCTCTTTCTGTAACGCGt                                    <  1:888602/149‑1 (MQ=255)
                                                                                                               gtgCCAGCCAGTGTTCGCCACGGCTTAACTGGCGTTCACGCCATTGCTGAAAATAGTGAATAAATTTATCGCGTAACATGTTATTTCCTCCGCAACGTTACGCCGCCGGAAACCGCATCACCCTCTTTCTGTAACGCGTCCTGTTGcac                          <  2:131332/149‑1 (MQ=255)
                                                                                                                     gccagTGTTCGCCACGGCTTAACTGGCGTTCACGCCATTGCTGAAAATAGTGAATAAATTTATCGCGTAACATGTTATTTCCTCCGCAACGTTACGCCGCCGGAAACCGCATCACCCTCttt                                               >  3:411676/1‑122 (MQ=255)
                                                                                                                       cagTGTTCGCCACGGCTTAACTGGCGTTCACGCCATTGCTGAAAATAGTGAATAAATTTATCGCGTAACATGTTATTTCCTCCGCAACGTTACGCCGCCGGAAACCGCATCACCCTCTTTCTGTAACGCGTCCTGTTGCACAACATAAt                  >  3:783740/1‑149 (MQ=255)
                                                                                                                       cagTGTTCGCCACGGCTTAACTGGCGTTCACGCCATTGCTGAAAATAGTGAATAAATTTATCGCGTAACATGTTATTTCCTCCGCAACGTTACGCCGCCGGAAACCGCATCACCCTCTTTCTGTAACGCGTCCTGTTGCACAACATAAt                  >  2:576230/1‑149 (MQ=255)
                                                                                                                        agTGTTCGCCACGGCTTAACTGGCGTTCACGCCATTGCTGAAAATAGTGAATAAATTTATCGCGTAACATGTTATTTCCTCCGCAACGTTACGCCGCCGGAAACCGCATCACCCTCTTTCTGTAACGCGTCCTGTTGCACAACATAATc                 <  1:695318/149‑1 (MQ=255)
                                                                                                                             tCGCCACGGCTTAACTGGCGTTCACGCCATTGCTGAAAATAGTGAATAAATTTATCGCGTAACATGTTATTTCCTCCGCAACGTTACGCCGCCGGAAACCGCATCACCCTCTTTCTGTAACGCGTCCTGTTGCACAACATAATCGGTTg            <  3:1241379/149‑1 (MQ=255)
                                                                                                                               gCCACGGCTTAACTGGCGTTCACGCCATTGCTGAAAATAGTGAATAAATTTATcgc                                                                                                       <  3:381196/56‑1 (MQ=255)
                                                                                                                               gCCACGGCTTAACTGGCGTTCACGCCATTGCTGAAAATAGTGAATAAATTTATcgc                                                                                                       >  4:381196/1‑56 (MQ=255)
                                                                                                                               gCCACGGCTTAACTGGCGTTCACGCCATTGCTGAAAATAGTGAATAAATTTATCGCGTAACATGTTATTTCCTCCGCAACGTTACGCCGCCGGAAACCGCATCACCCTCTTTCTGTAACGCGTCCTGTTGCACAACATaa                   >  2:656837/1‑140 (MQ=255)
                                                                                                                               gCCACGGCTTAACTGGCGTTCACGCCATTGCTGAAAATAGTGAATAAATTTATCGCGTAACATGTTATTTCCTCCGCAACGTTACGCCGCCGGAAACCGCATCACCCTCTTTCTGTAACGCGTCCTGTTGCACAACATAATCGGTTGcc          <  4:783740/149‑1 (MQ=255)
                                                                                                                                ccACGGCTTAACTGGCGTTCACGCCATTGCTGAAAATAGTGAATAAATTTATCGCGTAACATGTTATTTCCTCCGCAACGTTACGCCGCCGGAAACCGCATCACCCTCTTTCTGTAACGCGTCCTGTTGCACAACATAATCGGTTGCCa         <  3:634513/149‑1 (MQ=255)
                                                                                                                                 cACGGCTTAACTGGCGTTCACGCCATTGCTGAAAATAGTGAATAAATTTATCGCGTAACATGTTATTTCCTCCGCAACGTTACGCCGCCGGAAACCGCATCACCCTCTTTCTGTAACGCGTCCTGTTGCACAACATAATCGGTTGCCAg        <  1:25986/149‑1 (MQ=255)
                                                                                                                                       ttAACTGGCGTTCACGCCATTGCTGAAAATAGTGAATAAATTTATCGCGTAACATGTTATTTCCTCCGCAACGTTAcgccgc                                                                     >  1:823747/1‑82 (MQ=255)
                                                                                                                                        gcaatGGCGTGAACGCCATTGCTGAAAATAGTGAATAAATTTATCGCGTAACATGTTATTTCCTCCGCAACGTTACGCCGCCGGAAACCGCATCACCCTCTTTCTGTAACGCGTCCTGTTGCACAACATAATCGGtt             <  4:785109/133‑1 (MQ=255)
                                                                                                                                          aatGGCGTGAACGCCATTGCTGAAAATAGTGAATAAATTTATc                                                                                                         >  4:869609/3‑43 (MQ=255)
                                                                                                                                          aCTGGCGTTCACGCCATTGCTGAAAATAGTGAATAAATTTATCGCGTAACATGTTATTTCCTCCGCAACGTTACGCCGCCGGAAACCGCATCAccc                                                    <  1:911626/96‑1 (MQ=255)
                                                                                                                                          aCTGGCGTTCACGCCATTGCTGAAAATAGTGAATAAATTTATCGCGTAACATGTTATTTCCTCCGCAACGTTACGCCGCCGGAAACCGCATCAccc                                                    >  2:911626/1‑96 (MQ=255)
                                                                                                                                            tGGCGTTCACGCCATTGCTGAAAATAGTGAATAAATTTATCGCGTAACATGTTATTTCCTCCGCAACGTTACGCCGCCGGAAACCGCATCACCCTCTTTCTGTAACGCGTCCTGTTGCACAACATAATCGGTTGCCAGTa      <  2:1075553/140‑1 (MQ=255)
                                                                                                                                            tGGCGTTCACGCCATTGCTGAAAATAGTGAATAAATTTATCGCGTAACATGTTATTTCCTCCGCAACGTTACGCCGCCGGAAACCGCATCACCCTCTTTCTGTAACGCGTCCTGTTGCACAACATAATCGGTTGCCAGTa      >  1:1075553/1‑140 (MQ=255)
                                                                                                                                              gCGTTCACGCCATTGCTGAAAATAGTGAATAAATTTATCGCGTAACATGTTATTTCCTCCGCAACGTTACGCCGCCGGAAACCGCATCACCCTCTTTCTGTAACGCGTCCTGTTGCACAACATAATc                 <  2:324331/127‑1 (MQ=255)
                                                                                                                                              gCGTTCACGCCATTGCTGAAAATAGTGAATAAATTTATCGCGTAACATGTTATTTCCTCCGCAACGTTACGCCGCCGGAAACCGCATCACCCTCTTTCTGTAACGCGTCCTGTTGCACAACATAATc                 >  1:324331/1‑127 (MQ=255)
                                                                                                                                                  |                                                                                                                                         
TCCGTTCGCTCAACGGCTGCCAGATGAGGACGTAATATCCGACGCTAAACAGGAACACCACCGCTGCCAGTAACATGCCCTTTTCACGCGGAGAACGCCCCGCCAGGTGTTGTGCCAGCCAGTGTTCGCCACGGCTTAACTGGCGTTCACGCCATTGCTGAAAATAGTGAATAAATTTATCGCGTAACATGTTATTTCCTCCGCAACGTTACGCCGCCGGAAACCGCATCACCCTCTTTCTGTAACGCGTCCTGTTGCACAACATAATCGGTTGCCAGTACACT  >  NC_000913/3110937‑3111220

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