Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,178,886 T→G 7.3% 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 7.3% 204.8 / 6.6 69Q631P (CAA→CCA) PP_1033sulfatase domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base T (32/32);  new base G (2/3);  total (34/35)
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.60e-01

CAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTG‑TCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCTGACGCCAGTTGGTCAAAGAACGTGGACATGAAGAATGTGCGTTTC  >  NC_002947/1178738‑1179029
                                                                                                                                                    |                                                                                                                                                
cAGAAAGCCCGCAATCATGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGtt                                                                                                                                                  >  6:294884/1‑149 (MQ=255)
  gAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTgg                                                                                                                                                <  3:43292/149‑1 (MQ=255)
   aaaGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTggg                                                                                                                                               >  7:157758/1‑149 (MQ=255)
        ccGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGc                                                                                                                                          <  2:148032/149‑1 (MQ=255)
        ccGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCCCAACGCGTGGAGGGTga                                                                                                                                          <  2:77866/149‑2 (MQ=255)
           cATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCttt                                                                                                                                       <  8:157758/149‑1 (MQ=255)
            aggATGGCGTGCTTTCGCGTTCTTGTTCGTTCAGATTGCTTCATCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCCCAACGCGTGGAGGGTgacga                                                                                                                                       <  4:228941/146‑5 (MQ=255)
                   cGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCgg                                                                                                                               <  5:294884/149‑1 (MQ=255)
                    gggCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTgggg                                                                                                                                              >  8:109259/1‑133 (MQ=255)
                    gggCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTgggg                                                                                                                                              <  7:109259/133‑1 (MQ=255)
                    gggCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCggg                                                                                                                              >  6:117487/1‑149 (MQ=255)
                    gggCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCggg                                                                                                                              >  4:133494/1‑149 (MQ=255)
                      gCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGt                                                                                                                             >  5:48049/1‑148 (MQ=255)
                        tttCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCAt                                                                                                                          >  4:114570/1‑149 (MQ=255)
                          tCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATc                                                                                                                         <  4:281648/148‑1 (MQ=255)
                             ggTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATc                                                                                                                         <  6:112165/145‑1 (MQ=255)
                             ggTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATc                                                                                                                         >  5:112165/1‑145 (MQ=255)
                             ggTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGAGGGTCATCCTc                                                                                                                      >  2:46864/1‑148 (MQ=255)
                              gTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTc                                                                                                                      >  7:101617/1‑147 (MQ=255)
                              gTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTc                                                                                                                      >  7:76632/1‑147 (MQ=255)
                              gTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTc                                                                                                                      >  8:88063/1‑147 (MQ=255)
                              gTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTc                                                                                                                      >  2:183907/1‑147 (MQ=255)
                                    gttcgttcAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTggggg                                                                                                                                             >  2:157897/1‑118 (MQ=255)
                                    gttcgttcAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGtt                                                                                                                <  2:298886/147‑1 (MQ=255)
                                    gttcgttcAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGt                                                                                                                 <  6:48049/146‑1 (MQ=255)
                                     ttcgttcAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTa                                                                                                             <  5:69377/149‑1 (MQ=255)
                                         ttcAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTa                                                                                                         >  5:44685/1‑149 (MQ=255)
                                         ttcAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTa                                                                                                         <  8:76632/149‑1 (MQ=255)
                                             gATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCa                                                                                                     <  6:44685/149‑1 (MQ=255)
                                              aTTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGc                                                                                                                                 <  4:257154/120‑1 (MQ=255)
                                              aTTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAg                                                                                                    >  8:223226/1‑149 (MQ=255)
                                              aTTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTACGGGTCATCCTCAAGGTTGTATTTAGTCAg                                                                                                    <  1:183907/149‑1 (MQ=255)
                                               ttGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAgg                                                                                                   <  7:88063/149‑1 (MQ=255)
                                               ttGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAgg                                                                                                   <  8:101617/149‑1 (MQ=255)
                                                tGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAgg                                                                                                   <  7:223226/148‑1 (MQ=255)
                                                   ttATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCATGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGgtg‑t                                                                                              >  6:73532/1‑149 (MQ=255)
                                                   ttATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCT‑GTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTATTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTG‑Tc                                                                                             >  2:238597/1‑149 (MQ=255)
                                                    tATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTG‑Tc                                                                                             >  4:69066/1‑149 (MQ=255)
                                                         ttctGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTG‑TCCGAAc                                                                                        >  2:280109/1‑149 (MQ=255)
                                                           ctGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTG‑TCCGAACt                                                                                       <  5:117487/148‑1 (MQ=255)
                                                            tGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCaa                                                                                                                    >  6:139648/1‑119 (MQ=255)
                                                            tGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCaa                                                                                                                    <  5:139648/119‑1 (MQ=255)
                                                                tttCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCTATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTAAGTCAGGGTG‑TCCGAAATGAATTc                                                                                 >  7:292374/1‑149 (MQ=255)
                                                                             gcgcTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTG‑TCCGAACTGAATTCGTATATGTATTTg                                                                    >  3:172103/1‑149 (MQ=255)
                                                                              cgcTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTG‑TCCGAACTGAATTCGTATATGTATTTg                                                                    >  6:194993/1‑148 (MQ=255)
                                                                                           gCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCCTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTG‑Tcc                                                                                            <  1:201818/111‑1 (MQ=255)
                                                                                           gCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCCTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTG‑Tcc                                                                                            >  2:201818/1‑111 (MQ=255)
                                                                                            cGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTG‑TCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTccacc                                                     <  3:114570/149‑1 (MQ=255)
                                                                                                ataGGAAAGCAGCCTGCTACTCACTGCGATGAACTCGTCAGCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTG‑TCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAg                                                  <  8:292374/148‑1 (MQ=255)
                                                                                                  aGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTG‑TCCGAACTGAATTCGTATATGTAttt                                                                     <  4:74989/127‑1 (MQ=255)
                                                                                                  aGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTG‑TCCGAACTGAATTCGTATATGTAttt                                                                     >  3:74989/1‑127 (MQ=255)
                                                                                                    gAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTa                                                                                                             >  1:5445/1‑86 (MQ=255)
                                                                                                    gAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTa                                                                                                             <  2:5445/86‑1 (MQ=255)
                                                                                                     aaaGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTG‑TCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCTGa                                             <  1:46864/148‑1 (MQ=255)
                                                                                                          agcCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGTTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGCTCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCTGACGCCa                                        <  1:238597/149‑1 (MQ=255)
                                                                                                          agcCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTG‑TCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCTGACGCCAg                                       <  3:89197/149‑1 (MQ=255)
                                                                                                               gCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTG‑TCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCTGACGCCAGTTGGt                                  >  5:136959/1‑149 (MQ=255)
                                                                                                                   ctcACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTG‑TCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCTGACGCCAGTTGGTCaa                               >  8:3654/1‑148 (MQ=255)
                                                                                                                      aCTGCGATGAACTCGTCACCCTCCATGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTG‑TCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCTGACGCCAGTTGGTCAAAGaa                           <  5:73532/149‑1 (MQ=255)
                                                                                                                      aCTGCGATGAACTCGTCACCCTCCACGCGTTGCGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTG‑TCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCTGACGCCAGTTGGTCAAAGaa                           >  3:212162/1‑149 (MQ=255)
                                                                                                                        tGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTG‑TCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCTGACGCCAGTTGGTCAAAGAACg                         >  2:140929/1‑149 (MQ=255)
                                                                                                                          cGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTG‑TCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCTGACGCCAGTTGGTCAAAGAACTTg                       >  3:48759/1‑149 (MQ=255)
                                                                                                                          cGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTG‑TCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCTGACGCCAGTTGGTCAAAGAACGTg                       <  7:3654/149‑1 (MQ=255)
                                                                                                                          cGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTG‑TCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCTGACGCCAGTTGGTCAAAGAACGTg                       <  5:194993/149‑1 (MQ=255)
                                                                                                                                 cTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTG‑TCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCTGACGCCAGTTGGTCAAAGAACGTGGACATg                 >  5:245549/1‑148 (MQ=255)
                                                                                                                                 aaaagcACCCCCAACGCGTGGAGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTG‑TCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCTGACGCCAGTTGGTCAAAGACCGTGGACATg                 >  5:218136/6‑148 (MQ=255)
                                                                                                                                     gcACCCCCAACGCGTGGAGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTG‑TCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCTGACGCCa                                        >  7:11328/2‑121 (MQ=255)
                                                                                                                                     gcACCCCCAACGCGTGGAGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTG‑TCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCTGACGCCa                                        <  8:11328/120‑1 (MQ=255)
                                                                                                                                          cTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTG‑TCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCTGACGCCAGTTGGTCAAAGAACGTGGACATGAAGAATGTGc       <  3:133494/149‑1 (MQ=255)
                                                                                                                                               ttggTTGTGGGTGCTTTTCATTCGGGTCATCATCAAGGTTGTATTTAGTCATGGTG‑TCCGAACTGAATTCGTATATTTATTTGTACGGGTGCTCCACCAAGGCTTACGCCAGTTGGTCAAAGAACGTGGACATGAAGAATGTGCGTTTc  <  4:48759/146‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
CAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTG‑TCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCTGACGCCAGTTGGTCAAAGAACGTGGACATGAAGAATGTGCGTTTC  >  NC_002947/1178738‑1179029

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: ‑