Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,178,886 T→G 6.1% Q631P (CAA→CCA)  PP_1033 ← sulfatase domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,178,8860TG 6.1% 250.0 / 6.2 82Q631P (CAA→CCA) PP_1033sulfatase domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base T (39/38);  new base G (2/3);  total (41/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.08e-01

AGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGC‑TTTCTGGCAAA‑TAGCGCT‑TTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCTGACGCCAGTTGGTCAAAGAACGTGGACATGAAGAATGTGCGTTTCCCCT  >  NC_002947/1178739‑1179033
                                                                                                                                                      |                                                                                                                                                   
aGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGC‑TTTCTGGCAAA‑TAGCGCT‑TTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTg                                                                                                                                                    >  7:34899/1‑149 (MQ=255)
 gAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGC‑TTTCTGGCAAA‑TAGCGCT‑TTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTgg                                                                                                                                                   >  4:250040/1‑149 (MQ=255)
  aaaGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGC‑TTTCTGGCAAA‑TAGCGCT‑TTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTggg                                                                                                                                                  >  5:141024/1‑149 (MQ=255)
   aaGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGC‑TTTCTGGCAAA‑TAGCGCT‑TTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTgggg                                                                                                                                                 <  1:246324/149‑1 (MQ=255)
     gCCCGCTTGTTTGGGGGTTTTTTGGTTTTTTTTGGTTACTATTCCTTATTTTTTTGGT‑TTTTTGGAAAATTCGCTCTCTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTgggg                                                                                                                                                 <  1:101359/149‑1 (MQ=255)
      cccGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGC‑TTTCTGGCAAA‑TAGCGCT‑TTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTg                                                                                                                                              <  5:305672/149‑1 (MQ=255)
      cccGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGC‑TTTCTGGCAAA‑TAGCGCT‑TTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTg                                                                                                                                              >  6:154906/1‑149 (MQ=255)
        cGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGC‑TTTCTGGCAAA‑TAGCGCT‑TTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCCCAACGCGTGGAGGGTgac                                                                                                                                            <  5:294953/149‑3 (MQ=255)
          cATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGC‑TTTCTGGCAAA‑TAGCGCT‑TTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCCCAACGCGTGGAGGGTgacga                                                                                                                                          <  3:317666/149‑5 (MQ=255)
                 gCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCGTTTCTGGCAAA‑TAGCGCT‑TTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGc                                                                                                                                    >  1:71990/1‑149 (MQ=255)
                   gggCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGC‑TTTCTGGCAAA‑TAGCGCT‑TTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCggg                                                                                                                                 <  7:162640/149‑1 (MQ=255)
                       tttCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGC‑TTTCTGGCAAA‑TAGCGCT‑TTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCAt                                                                                                                             >  2:180279/1‑149 (MQ=255)
                         tCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGC‑TTTCTGGCAAA‑TAGCGCT‑TTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTTCGGGTCATc                                                                                                                            <  2:46163/148‑1 (MQ=255)
                             gTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGC‑TTTCTGGCAAA‑TAGCGCT‑TTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCaa                                                                                                                       <  2:250649/149‑1 (MQ=255)
                                    ttcgttcAGATTGCTTATTCTTCTGGC‑TTTCTGGCAAA‑TAGCGCT‑TTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGt                                                                                                                 <  5:154906/148‑1 (MQ=255)
                                      cgttcAGATTGCTTATTCTTCTGGC‑TTTCTGGCAAA‑TAGCGCT‑TTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTAt                                                                                                               <  3:270911/148‑1 (MQ=255)
                                        ttcAGATTGCTTATTCTTCTGGC‑TTTCTGGCAAA‑TAGCGCT‑TTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCATGCGGGTCATCCTCAAGGTTGTATTTa                                                                                                            <  1:180279/149‑1 (MQ=255)
                                            gATTGCTTATTCTTCTGGC‑TTTCTGGCAAA‑TAGCGCT‑TTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCa                                                                                                        <  6:186477/149‑1 (MQ=255)
                                              ttGCTTATTCTTCTGGC‑TTTCTGGCAAA‑TAGCGCT‑TTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAgg                                                                                                      <  8:34899/149‑1 (MQ=255)
                                                  ttATTCTTCTGGC‑TTTCTGGCAAA‑TAGCGCT‑TTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGgtgt                                                                                                  >  7:1810/1‑149 (MQ=255)
                                                  gtaatcttctGGA‑TTTCTGGCAAA‑TAGCGCT‑TTCTGGTAGGCGTTATAGTAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTACACGCGTTGGGGGTGCTTTTCATGCGTGTCATCCTCAACTTTGTATTTATTCAgggag                                                                                                   >  4:224203/5‑146 (MQ=255)
                                                    aTTCTTCTGGC‑TTTCTGGCAAA‑TAGCGCT‑TTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTcc                                                                                                >  4:138673/1‑149 (MQ=255)
                                                    aTTCTTCTGGC‑TTTCTGGCAAA‑TAGCGCT‑TTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGCGTCATCCTCAAGGTTGAATTTAGTCAGGGTGTcc                                                                                                >  4:238163/1‑149 (MQ=255)
                                                        ttctGGC‑TTTCTGGCAAA‑TAGCGCT‑TTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCGAGGTTGTATTTAGTCAGGGTGTCCGAAc                                                                                            >  7:2627/1‑149 (MQ=255)
                                                        ttctGGC‑TTTCTGGCAAA‑TAGCGCT‑TTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAAc                                                                                            >  7:165799/1‑149 (MQ=255)
                                                        ttctGGC‑TTTCTGGCAAA‑TAGCGCT‑TTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAAc                                                                                            >  1:131761/1‑149 (MQ=255)
                                                        ttctGGC‑TTTCTGGCAAA‑TAGCGCT‑TTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAAc                                                                                            >  5:326136/1‑149 (MQ=255)
                                                        ttctGGC‑TTTCTGGCAAA‑TAGCGCT‑TTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAAc                                                                                            >  4:149621/1‑149 (MQ=255)
                                                        ttctGGC‑TTTCTGGCAAA‑TAGCGCT‑TTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAAc                                                                                            >  1:288620/1‑149 (MQ=255)
                                                        ttctGGC‑TTTCTGGCAAA‑TAGCGCT‑TTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTAAGGGTGTCCGAAc                                                                                            >  1:274197/1‑149 (MQ=255)
                                                           tGGC‑TTTCTGGCAAA‑TAGCGCT‑TTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTg                                                                                          >  5:287776/1‑148 (MQ=255)
                                                             gC‑TTTCTGGCAAA‑TAGCGCT‑TTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCTAGGTTGTATTTAGTCAGGGTGTCCGAACTGa                                                                                         <  8:2627/147‑1 (MQ=255)
                                                             gC‑TTTCTGGCAAA‑TAGCGCT‑TTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAAt                                                                                       >  5:107421/1‑149 (MQ=255)
                                                                    tGGCAAA‑TAGCGCT‑TTCT‑GTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGtata                                                                                >  2:290284/1‑149 (MQ=255)
                                                                    tGGCAAA‑TAGCGAT‑TTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGta                                                                                  >  2:323294/1‑148 (MQ=255)
                                                                     ggCAAA‑TAGCGCT‑TTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCgg                                                                                                                                  >  2:63671/1‑99 (MQ=255)
                                                                     ggCAAA‑TAGCGCT‑TTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCgg                                                                                                                                  <  1:63671/99‑1 (MQ=255)
                                                                     ggCAAA‑TAGCGCT‑TTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAgg                                                                                                      >  8:140793/1‑127 (MQ=255)
                                                                     ggCAAA‑TAGCGCT‑TTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAgg                                                                                                      <  7:140793/127‑1 (MQ=255)
                                                                      gCAAA‑TAGCGCT‑TTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGta                                                                                  <  2:71990/146‑1 (MQ=255)
                                                                      gCAAA‑TAGCGCT‑TTCTCGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGtat                                                                                 <  8:143202/147‑1 (MQ=255)
                                                                          a‑TAGCGCT‑TTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTAt                                                                           <  2:240473/149‑1 (MQ=255)
                                                                             aGCGCT‑TTCTGGTCGGCGTTATTGTAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCTCGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCATGGTTGTATTTAGTCTGTGTTTCCGAACTGTATTCGTATATGTAttt                                                                         <  2:288620/149‑1 (MQ=255)
                                                                               cgcT‑TTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTg                                                                        >  1:233037/1‑148 (MQ=255)
                                                                                      cTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACggg                                                                  <  2:233037/148‑1 (MQ=255)
                                                                                      cTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACGGGt                                                                 >  3:83701/1‑149 (MQ=255)
                                                                                             gCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTccac                                                          <  8:1810/149‑1 (MQ=255)
                                                                                               gTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTccacc                                                         >  2:122217/1‑148 (MQ=255)
                                                                                                ttATAGGAAAGGTGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGATGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCaa                                                       <  3:238163/149‑1 (MQ=255)
                                                                                                   taGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAgg                                                     <  3:250040/148‑1 (MQ=255)
                                                                                                     ggAAAGCAGCCTGCTCCTCACTGCGATGAACGCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCt                                                   >  1:29618/1‑148 (MQ=255)
                                                                                                      gAAAGCAGCCTGCTCCTGACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGTGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCTGa                                                 <  3:224203/149‑1 (MQ=255)
                                                                                                      gAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCTGa                                                 >  2:256067/1‑149 (MQ=255)
                                                                                                      gAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCTGa                                                 <  1:323294/149‑1 (MQ=255)
                                                                                                             gcCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGGTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCTGACGCCAGt                                          <  8:38157/149‑1 (MQ=255)
                                                                                                              ccTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCTGACGCCAGtt                                         <  6:326136/149‑1 (MQ=255)
                                                                                                               cTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCTGACGCCAGTTg                                        >  5:179345/1‑149 (MQ=255)
                                                                                                                  ctcctcACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCTGACGCCAGTTgg                                       <  4:91646/147‑1 (MQ=255)
                                                                                                                  ctcctcACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCTGACGCCAGTTGGt                                      <  6:141024/148‑1 (MQ=255)
                                                                                                                     ctcACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCTGACGCCAGTTg                                        >  5:153640/1‑143 (MQ=255)
                                                                                                                     ctcACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCTGACGCCAGTTg                                        <  6:153640/143‑1 (MQ=255)
                                                                                                                      tcACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCTGACGCCAGTTGGTCa                                    >  3:160308/1‑146 (MQ=255)
                                                                                                                      tcACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCTGACGCCAGTTGGTCa                                    >  3:94103/1‑146 (MQ=255)
                                                                                                                      tcACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCTGACGCCAGTTGGTCa                                    >  1:180038/1‑146 (MQ=255)
                                                                                                                         cTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCTGACGCCAGTTGGTCAAAGAAc                              >  2:105234/1‑149 (MQ=255)
                                                                                                                            cGATGACCTCGTCACCCTCCACGCGTTGGGGGTGCGTTTCCTGCGGGTCATCCTCAAGGGTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCTGACGCCAGTTGGTCAAAGAACGTg                           <  6:109279/149‑1 (MQ=255)
                                                                                                                            cGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCTGACGCCAGTTGGTCAAAGAACGt                            <  4:83701/148‑1 (MQ=255)
                                                                                                                             gATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCTGACGCCAGTTGGTCAAAGAACGTgg                          <  3:138673/149‑1 (MQ=255)
                                                                                                                                 accTCGTCACCCTGCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTTGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACgg                                                                   <  2:115223/102‑1 (MQ=255)
                                                                                                                                 aaCTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCTGACGCCAGTTGGTCAAAGAACGTGGACa                       <  3:149621/148‑1 (MQ=255)
                                                                                                                                 aaCTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCTGACGCCAGTTGGTCAAAGAACGTGGACAt                      >  1:147692/1‑149 (MQ=255)
                                                                                                                                 aaCTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCTGACGCCAGTTGGTCAAAGAACGTGGACAt                      >  5:109308/1‑149 (MQ=255)
                                                                                                                                     aagcACCCCCAACGCGTGGAGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCTGACGCCAGTTGGTCAAAGAACGTGGACATg                     >  2:207090/4‑146 (MQ=255)
                                                                                                                                          cccTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCTGACGCCAGTTGGTCAAAGAACGTGGACATGAAGAAtgt             <  6:107421/149‑1 (MQ=255)
                                                                                                                                           ccccaaCGCGTGGAGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCTGACGCCAGTTGGTCAAAGAACGTGGACATGAAGATtgtg            >  6:45003/6‑149 (MQ=255)
                                                                                                                                           ccTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCTGACGCCAGTTGGTCAAAGAACGTGGACATGAAGAAtgtg            >  4:250009/1‑149 (MQ=255)
                                                                                                                                           ccTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCTGACGCCAGTTGGTCAAAGAACGTGGACATGAAGAAtgtg            <  2:274197/149‑1 (MQ=255)
                                                                                                                                            cTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCTGACGCCAGTTGGTCAAAGAACGTGGACATGAAGAATGTGc           >  7:218867/1‑149 (MQ=255)
                                                                                                                                             tCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCTGACGCCAGTTGGTCAAAGAACGTGGACATGAAGAATGTGCg          >  3:278787/1‑149 (MQ=255)
                                                                                                                                               cACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCTGACGCCAGTTGGTCAAAGAACGTGGACATGAAGAATGTGCGt         <  1:290284/148‑1 (MQ=255)
                                                                                                                                                  gcgTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCTGACGCCAGTTGGTCAAAGAACGTGGACATGAAGAATGTGCGTTTcc     <  6:179345/149‑1 (MQ=255)
                                                                                                                                                  gcgTGGAGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCTGACGCCAGTTGGTCAAAGAACGTGGACATGAAGAATGTGCGTTTcc     <  1:207090/149‑1 (MQ=255)
                                                                                                                                                     ttGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCTGACGCCAGTTGGTCAAAGAACGTGGACATGAAGAATGTGCGTTTCCCCt  <  4:94103/149‑1 (MQ=255)
                                                                                                                                                     ttGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCTGACGCCAGTTGGTCAAAGAACGTGGACATGAAGAATGTGCGTTTCCCCt  <  4:278787/149‑1 (MQ=255)
                                                                                                                                                     ttGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCTGACGCCAGTTGGTCAAAGAACGTGGACATGAAGAATGTGCGTTTCCCCt  <  6:109308/149‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                   
AGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGC‑TTTCTGGCAAA‑TAGCGCT‑TTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCTGACGCCAGTTGGTCAAAGAACGTGGACATGAAGAATGTGCGTTTCCCCT  >  NC_002947/1178739‑1179033

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