Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,178,888 G→A 8.0% P630P (CCC→CCT)  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,8880GA 8.0% 172.8 / 4.3 64P630P (CCC→CCT) PP_1033sulfatase domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base G (28/29);  new base A (2/3);  total (31/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.48e-01

GAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTG‑TCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCTGACGCCAGTTGGTCAAAGAACGTGGACATGAAGAATGTGCGTTTCCCCTGTT  >  NC_002947/1178740‑1179036
                                                                                                                                                    |                                                                                                                                                     
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:76632/1‑147 (MQ=255)
                            gTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTc                                                                                                                             >  2:183907/1‑147 (MQ=255)
                            gTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTc                                                                                                                             >  8:88063/1‑147 (MQ=255)
                            gTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTc                                                                                                                             >  7:101617/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                                                                                                          <  8:101617/149‑1 (MQ=255)
                                             ttGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAgg                                                                                                          <  7:88063/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                                                                                                                           <  5:139648/119‑1 (MQ=255)
                                                          tGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCaa                                                                                                                           >  6:139648/1‑119 (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                                               <  8:11328/120‑1 (MQ=255)
                                                                                                                                   gcACCCCCAACGCGTGGAGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTG‑TCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCTGACGCCa                                               >  7:11328/2‑121 (MQ=255)
                                                                                                                                        cTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTG‑TCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCTGACGCCAGTTGGTCAAAGAACGTGGACATGAAGAATGTGc              <  3:133494/149‑1 (MQ=255)
                                                                                                                                             ttggTTGTGGGTGCTTTTCATTCGGGTCATCATCAAGGTTGTATTTAGTCATGGTG‑TCCGAACTGAATTCGTATATTTATTTGTACGGGTGCTCCACCAAGGCTTACGCCAGTTGGTCAAAGAACGTGGACATGAAGAATGTGCGTTTc         <  4:48759/146‑1 (MQ=255)
                                                                                                                                                    ggggTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTG‑TCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCTGACGCCAGTTGGTCAAAGAACGTGGACATGAAGAATGTGCGTTTCCCCTGtt  >  2:114715/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTG‑TCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCTGACGCCAGTTGGTCAAAGAACGTGGACATGAAGAATGTGCGTTTCCCCTGTT  >  NC_002947/1178740‑1179036

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