Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,060,884 C→T 6.5% S30S (AGC→AGT)  PP_0918 → 3‑beta hydroxysteroid dehydrogenase/isomerase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,060,8840CT 6.5% 181.0 / 2.2 61S30S (AGC→AGT) PP_09183‑beta hydroxysteroid dehydrogenase/isomerase family protein
Reads supporting (aligned to +/- strand):  ref base C (30/27);  new base T (2/2);  total (32/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

GCCCTTGTAGCCATGGCTGGCGTCAGTTGCATCTGAAATCTGGCAA‑GACGGAGTGTTACATGCGAATTCTGGTCACCGGCGCGAGTGGCTTCAT‑TGGCGGGCGTTTTGCGCGTTTCGCCCTGGAGCAAGGCCTGGACGTGCGCGTC‑AGCGGCCGCCGTGCCGAAGGGGTCGAGCATCTGGTCAAGCGTGGCGCACAGTTCATCCCGGGTGATCTGGGTGACGCCGAGCTGGCCCGGCGTTTGTGCCAGGGTGTCGAAGCCGTGGTGCATTGCGCTGGCGCGGTGGGTAACTGGGGG  >  NC_002947/1060736‑1061031
                                                                                                                                                       |                                                                                                                                                   
gCCCTTGTAGCCATGGCTGGCGTCAGTTGCATCTGAAATCTGGCAA‑GACGGAGTGTTACATGCGAATTCTGGTCACCGGCGCGAGTGGCTTCAT‑TGGCGGGCGTTTTGCGCGTTTCGCCCTGGAGCAAGGCCTGGACGTGCGCGTC‑AGc                                                                                                                                                     >  8:43703/1‑149 (MQ=255)
 cccTTGTAGCCATGGCTGGCGTCAGTTGCATCTGAAATCTGGCAA‑GACGGAGTGTTACATGCGAATTCTGGTCACCGGCGCGAGTGGCTTCAT‑TGGCGGGCGTTTTGCGCGTTTCGCCCTGGAGCAAGGCCTGGACGTGCGCGTC‑AGCg                                                                                                                                                    >  7:81795/1‑149 (MQ=255)
      gTAGCCATGGCTGGCGTCAGTTGCATCTGAAATCTGGCAA‑GACGGAGTGTTACATGCGAATTCTGGTCACCGGCGCGAGTGGCTTCAT‑TGGCGGGCGTTTTGCGCGTTTCGCCCTGGAGCAAGGCCTGGACGTGCGCGTC‑AGCGgcc                                                                                                                                                 <  3:186960/147‑1 (MQ=255)
          ccATGGCTGGCGTCAGTTGCATCTGAAATCTGGCAATGACGGAGTGTTACATGCGAATTCTGGTCACCGGCGCTAGTGGCTTCAT‑TGGCGGGCGTTTTGCGCGTTTCGCCCTGGAGCAAGGCCTGGACGTGCGCGTC‑AGCGGCCGCCGt                                                                                                                                            >  1:53585/1‑149 (MQ=255)
              ggctggcGTCAGTTGCATCTGAAATCTGGCAA‑GACGGAGTGTTACATGCGAATTCTGGTCACCGGCGCGAGTGGCTTCAT‑TGGCGGGCGTTTTGCGCGTTTCGCCCTGGAGCAAGGCCTGGACGTGCGCGTC‑AGCGGCCGCCGTGCCg                                                                                                                                        >  2:291967/1‑148 (MQ=255)
                ctggcGTCAGTTGCATCTGAAATCTGGCAA‑GACGGAGTGTTACATGCGAATTCTGGTCACCGGCGCGAGTGGCTTCAT‑TGGCGGGCGTTTTGCGCGTTTCGCCCTGGAGCAAGGCCTGGACGTGCGCGTC‑AGCGGCCGCCGTGCCGAAg                                                                                                                                     <  8:299525/149‑1 (MQ=255)
                ctggcGTCAGTTGCATCTGAAATCTGGCAA‑GACGGAGTGTTACATGCGAATTCTGGTCACCGGCGCGAGTGGCTTCAT‑TGGCGGGCGTTTTGCGCGTTTCGCCCTGGAGCAAGGCCTGGACGTGCGCGTC‑AGCGGCCGCCGTGCCGAAg                                                                                                                                     >  1:6271/1‑149 (MQ=255)
                     gTCAGTTGCATCTGAAATCTGGCAA‑GACGGAGTGTTACATGCGAATTCTGGTCACCGGCGCGAGTGGCTTCAT‑TGGCGGGCGTTTTGCGCGTTTCGCCCTGGAGCAAGGCCTGGACGTGCGCGTC‑AGCGgccg                                                                                                                                                <  1:294210/133‑1 (MQ=255)
                     gTCAGTTGCATCTGAAATCTGGCAA‑GACGGAGTGTTACATGCGAATTCTGGTCACCGGCGCGAGTGGCTTCAT‑TGGCGGGCGTTTTGCGCGTTTCGCCCTGGAGCAAGGCCTGGACGTGCGCGTC‑AGCGgccg                                                                                                                                                >  2:294210/1‑133 (MQ=255)
                     gTCAGTTGCATCTGAAATCTGGCAA‑GACGGAGTGTTACATGCGAATTCTGGTCACCGGCGCGAGTGGCTTCAT‑TGGCGGGCGTTTTGCGCGTTTCGCCCTGGAGCAAGGCCTGGACGTGCGCGTC‑AGCGGCCGCCGTGCCGAAgg                                                                                                                                    >  1:295224/1‑145 (MQ=255)
                     gTCAGTTGCATCTGAAATCTGGCAA‑GACGGAGTGTTACATGCGAATTCTGGTCACCGGCGCGAGTGGCTTCAT‑TGGCGGGCGTTTTGCGCGTTTCGCCCTGGAGCAAGGCCTGGACGTGCGCGTC‑AGCGGCCGCCGTGCCGAAgg                                                                                                                                    <  2:295224/145‑1 (MQ=255)
                                 tGAAATCTGGCAA‑GACGGAGTGTTACATGCGAATTCTGGTCACCGGCGCGAGTGGCTTCAT‑TGGCGGGCGTTTTGCGCGTTTCGCCCTGGAGCAAGGCCTGGACGTGCGCGTC‑AGCGGCCGCCGTGCCGAAGGGGTCGAGCATCTGGt                                                                                                                     >  2:155252/1‑148 (MQ=255)
                                  gAAATCTGGCAA‑GACGGAGTGTTACATGCGAATTCTGGTCACCGGCGCGAGTGGCTTCAT‑TGGCGGGCGTTTTGCGCGTTTCGCCCTGGAGCAAGGCCTGGACGTGCGCGTC‑AGCGGCCGCCGTGCCGAAGGGGTCGAGCATCTGGTc                                                                                                                    >  8:104650/1‑148 (MQ=255)
                                        tGGCAA‑GACGGAGTGTTACATGCGAATTCTGATCACCGGCGCGAGTGGCTTCAT‑TGGCGGGCGTTTTGCGCGTTTCGCCCTGGAGCAAGGCCTGGACGTGCGCGTC‑AGCGGCCGCCGTGCCGAAGGGGTCGAGCATCTGGTCAAGCGTg                                                                                                             >  7:193100/1‑149 (MQ=255)
                                          gCAA‑GACGGAGTGTTACATGCGAATTCTGATCACCGGCGCGAGTGGCTTCAT‑TGGCGGGCGTTTTGCGCGTTTCGCCCTGGAGCAAGGCCTGGACGTGCGCGTC‑AGCGGCCGCCGTGCCGAAGGGGTCGAGCATCTGGTCAAGCGTg                                                                                                             <  8:193100/147‑1 (MQ=255)
                                          gCAA‑GACGGAGGTTTACATGCGAATGCTGGTCGCCGGCGCGAGTGGCTTCAT‑TGGCGGGCGTTTTGCGCGTTTCGCCCTGGAGCAAGGCCTGGACGTGCGCGTC‑AGCGGCCGCCGTGCCGAAGGGGTCGAGCATCTGGTCAAGCGt                                                                                                              <  8:134259/146‑1 (MQ=255)
                                            aa‑GACGGAGTGTTACATGCGAATTCTGGTCACCGGCGCGAGTGGCTTCATATGGCGGGCGTTTTGCGCGTTTCGCCCTGGAGCAAGGCCTGGACGTGCGCGTC‑AGCGGCCGCCGTGCCGAAGGGGTCGAGCATCTGGTCAAGCGTGgcg                                                                                                          >  1:79630/1‑149 (MQ=255)
                                            aa‑GACGGAGTGTTACATGCGAATTCTGGTCACCGGCGCGAGTGGCTTCAT‑TGGCGGGCGTTTTGCGCGTTTCGCCCTGGAGCAAGGCCTGGACGTGCGCGTC‑AGCGGCCGCCGTGCCGAAGGGGTCGAGCATCTGGTCAAGCGTGgcg                                                                                                          >  6:307141/1‑148 (MQ=255)
                                                              gCGAATTCTGGTCACCGGCGCGAGTGGCTTCAT‑TGGCGGGCGTTTTGCGCGTTTCGCCCTGGAGCAAGGCCTGGACGTGCGCGTCTAGCGGCCGCCGTGCCGAAGGGGTCGa                                                                                                                              <  1:167245/112‑1 (MQ=255)
                                                              gCGAATTCTGGTCACCGGCGCGAGTGGCTTCAT‑TGGCGGGCGTTTTGCGCGTTTCGCCCTGGAGCAAGGCCTGGACGTGCGCGTC‑AGCGGCCGCCGTGCCGAAGGGGTCGa                                                                                                                              >  2:167245/1‑111 (MQ=255)
                                                              gCGAATTCTGGTCACCGGCGCGAGTGGCTTCAT‑TGGCGGGCGTTTTGCGCGTTTCGCCCTGGAGCAAGGCCTGGACGTGCGCGTC‑AGCGGCCGCCGTGCCGAAGGGGTCGAGCATc                                                                                                                         >  4:179138/1‑116 (MQ=255)
                                                               cGAATTCTGGTCACCGGCGCGAGTGGCTTCAT‑TGGCGGGCGTTTTGCGCGTTTCGCCCTGGAGCAAGGCCTGGACGTGCGCGTC‑AGCGGCCGCCGTGCCGAAGGGGTCGAGCATCTGGTCAAGCGTGGCGCACAGTTCATCCCGGGTg                                                                                        <  7:7609/148‑1 (MQ=255)
                                                               cGAATTCTGGTCACCGGCGCGAGTGGCTTCAT‑TGGCGGGCGTTTTGCGCGTTTCGCCCTGGAGCAAGGCCTGGACGTGCGCGTC‑AGCGGCCGCCGTGCCGAAGGGGTCGAGCATCTGGTCAAGCGTGGCGCACAGTTCATCCCGGGTGa                                                                                       <  3:41234/149‑1 (MQ=255)
                                                                  aTTCTGGTCACCGGCGCGAGTGGCTTCAT‑TGGCGGGCGTTTTGCGCGTTTCGCCCTGGAGCAAGGCCTGGACGTGCGCGTC‑AGCGGCCGCCGTGCCGAAGGGGTCGAGCATCTGGTCAAGCGTGGCGCACAGTTCATCCCGGGt                                                                                         <  4:121373/144‑1 (MQ=255)
                                                                  aTTCTGGTCACCGGCGCGAGTGGCTTCAT‑TGGCGGGCGTTTTGCGCGTTTCGCCCTGGAGCAAGGCCTGGACGTGCGCGTC‑AGCGGCCGCCGTGCCGAAGGGGTCGAGCATCTGGTCAAGCGTGGCGCACAGTTCATCCCGGGTGagat                                                                                    >  3:121373/1‑146 (MQ=255)
                                                                   ttCTGGTCACCGGCGCGAGTGGCTTCAT‑TGGTGGGCGTTTTGCGCGTTTCGCCCTGGAGCAAGGCCTGGACGTGCGCGTC‑AGCGGCCGCCGTGCCGAAGGGGTCGAGCATCTGGTCAAGCGTGGCGCACAGTTCATCCCGGGTGATCTg                                                                                   >  2:98983/1‑149 (MQ=255)
                                                                        gTCACCGGCGCGAGTGGCTTCAT‑TGGCGGGCGTTTTGCGCGTTTCGCCCTGGAGCAAGGCCTGGACGTGCGCGTC‑AGCGGCCGCCGTGCCGAAGGGGTCGAGCATCTGGTCAAGCGTGGCGCACAGTTCATCCCGGGTGATCTGGGTGa                                                                              >  6:144/1‑149 (MQ=255)
                                                                             cGGCGCGAGTGGCTTCAT‑TGGCGGGCGTTTTGCGCGTTTCGCCCTGGAGCAAGGCCTGGACGTGCGCGTC‑AGCGGCCGCCGTGCCGAAGGGGTCGAGCATCTGGTCAAGCGTGGCGCACAGTTCATCCCGGGTGATCTGGGTGACGcc                                                                          >  5:58620/1‑148 (MQ=255)
                                                                              ggCGCGAGTGGCTTCAT‑TGGCGGGCGTTTTGCGCGTTTCGCCCTGGAGCAAGGCCTGGACGTGCGCGTC‑AGCGGCCGCCGTGCCGAAGGGGTCGAGCATCt                                                                                                                        <  1:109524/101‑1 (MQ=255)
                                                                              ggCGCGAGTGGCTTCAT‑TGGCGGGCGTTTTGCGCGTTTCGCCCTGGAGCAAGGCCTGGACGTGCGCGTC‑AGCGGCCGCCGTGCCGAAGGGGTCGAGCATCt                                                                                                                        >  2:109524/1‑101 (MQ=255)
                                                                                            cAT‑TGGCGGGCGTTTTGCGCGTTTCGCCCTGGAGCAAGGCCTGGACGTGCGCGTC‑AGCGGCCGCCGTGCCGAAGGGGTCGAGCATCTGGTCAAGCGTGGCGCACAGTTCATCCCggg                                                                                          >  5:114474/1‑117 (MQ=255)
                                                                                            cAT‑TGGCGGGCGTTTTGCGCGTTTCGCCCTGGAGCAAGGCCTGGACGTGCGCGTC‑AGCGGCCGCCGTGCCGAAGGGGTCGAGCATCTGGTCAAGCGTGGCGCACAGTTCATCCCggg                                                                                          <  6:114474/117‑1 (MQ=255)
                                                                                             aT‑TGGCGGGCGTTTTGCGCGTTTCGCCCTGGAGCAAGGCCTGGACGTGCGCGTC‑AGCGGCCGCCGTGCCGAAGGGGTCGAGCATCTGGTCAAGCGTGGCGCACAGTTCATCCCGGGTGATCTGGGTGACGCCGGGCCAGCTCGGCGtca                                                         <  7:104650/149‑3 (MQ=255)
                                                                                             aT‑TGGCGGGCGTTTTGCGCGTTTCGCCCTGGAGCAAGGCCTGGACGTGCGCGTC‑AGCGGCCGCCGTGCCGAAGGGGTCGAGCATCTGGTCAAGCGTGGCGCACAGTTCATCCCGGGTGATCTGGGTGACGCCGAGCTGGCCCGGCGttt                                                         >  2:151068/1‑149 (MQ=255)
                                                                                                 ggcgggcgTTTTGCGCGTTTCGCCCTGGAGCAAGGCCTGGACGTGCGCGTC‑AGCGGCCGCCGTGCCGAAGGGGTCGAGCATCTGGTCAAGCGTGGCGCACAGTTCATCCCGGGTGATCTGGGTGACGCCGGGCCAGCTCGGCGtcacc                                                       <  1:155252/148‑5 (MQ=255)
                                                                                                 ggcgggcgTTTTGCGCGTTTCGCCCTGGAGCAAGGCCTGGACGTGCGCGTC‑AGCGGCCGCCGTGCCGAAGGGGTCGAGCATCTGGTCAAGCGTGGCGCACAGTTCATCCCGGGTGATCTGGGTGACGCCGAGCTGGCCCGGCGtcacc                                                       <  6:58620/148‑5 (MQ=255)
                                                                                                  gcgggcgTTTTGCGCGTTTCGCCCTGGAGCAAGGCCTGGACGTGCGCGTC‑AGCGGCCGCCGTGCCGAAGGGGTCGAGCATCTGGTCAAGCGTGGCGCACAGTTCATCCCGGGTGATCTGGGTGACGCCGAGCTGGCCCGGCGTTTGTGc                                                     >  5:29747/1‑149 (MQ=255)
                                                                                                      gcgTTTTGCGCGTTTCGCCCTGGAGCAAGGCCTGGACGTGCGCGTC‑AGCGGCCGCCGTGCCGAAGGGGTCGAGCATCTGGTCAAGCGTGGCGCACAGTTCATCCCGGGTGATCTGGGTGACGCCGAGCTGGCCCGGCGTTTGTGCCAgg                                                 <  1:291967/149‑1 (MQ=255)
                                                                                                              cgcgTTTCGCCCTGGAGCAAGGCCTGGACGTGCGCGTC‑AGCGGCCGCCGTGCCGAAGGGGTCGAGCATCTGGTCAAGCGTGGCGCACAGTTCATCCCGGGTGATCTGGGTGACGCCGAGCTGGCCCGGCGTTTGTGCCAGGGTGTCGa                                          <  6:29747/148‑1 (MQ=255)
                                                                                                                     cGCCCTGGAGCAAGGCCTGGACGTGCGCGTC‑AGCGGCCGCCGTGCCGAAGGGGTCGAGCATCTGGTCAAGCGTGGCGCACAGTTCATCCCGGGTGATCTGGGTGACGCCGAGCTGGCCCGGCGTTTGTGCCAGGGTGTCGAAGCCgtg                                   <  5:307141/148‑1 (MQ=255)
                                                                                                                     cGCCCTGGAGCAAGGCCTGGACGTGCGCGTC‑AGCGGCCGCCGTGCCGAAGGGGTCGAGCATCTGGTCAAGCGTGGCGCACAGTTCATCCCGGGTGATCTGGGTGACGCCGAGCTGGCCCGGCGTTTGTGCCAGGGTGTCGAAGCCgtg                                   <  5:144/148‑1 (MQ=255)
                                                                                                                      gCCCTGGAGCAAGGCCTGGACGTGCGCGTC‑AGCGGCCGCCGTGCCGAAGGGGTCGAGCATCTGGTCAAGCGTGGCGCACAGTTCATCCCGGGTGATCTGGGTGACGCCGAGCTGGCCCGGCGTTTGTGCCAGGGTGTCGAAGCCgtggt                                 <  1:98983/149‑1 (MQ=255)
                                                                                                                        ccTGGAGCAAGGCCTGGACGTGCGCGTC‑AGCGGCCGCCGTGCCGAAGGGGTCGAGCATCTGGTCAAGCGTGGCGCACAGTTCATCCCGGGTGATCTGGGTGACGCCGAGCTGGCCCGGCGTTTGTGCCAGGGTGTCGAAGCCgtggt                                 <  6:184176/147‑1 (MQ=255)
                                                                                                                        ccTGGAGCAAGGCCTGGACGTGCGCGTC‑AGCGGCCGCCGTGCCGAAGGGGTCGAGCATCTGGTCAAGCGTGGCGCACAGTTCATCCCGGGTGATCTGGGTGACGCCGAGCTGGCCCGGCGTTTGTGCCAGGGTGTCGAAGCCgtggt                                 >  5:184176/1‑147 (MQ=255)
                                                                                                                         cTGGAGCAAGGCCTGGACGTGCGCGTC‑AGCGGCCGCCGTGCCGAAGGGGTCGAGCATCTGGTCAAGCGTGGCGCACAGTTCATCCCGGGTGATCTGGGTGACGCCGAGCTGGCCCGGCGTTTGTGCCAGGGTGTCGAAGCCGTGGTGCa                              >  7:109857/1‑149 (MQ=255)
                                                                                                                         cTGGAGCAAGGCCTGGACGTGCGCGTC‑AGCGGCCGCCGTGCCGAAGGGGTCGAGCATCTGGTCAAGCGTGGCGCACAGTTCATCCCGGGTGATCTGGGTGACGCCGAGCTGGCCCGGCGTTTGTGCCAGGGTGTCGAAGCCGTGGTGCa                              >  4:301940/1‑149 (MQ=255)
                                                                                                                               cAAGGCCTGGACGTGCGCGTC‑AGCGGCCGCCGTGCCGAAGGGGTCGAGCATCTGGTCAAGCGTGGCGCACAGTTCATCCCGGGTGATCTGGGTGACGCCGAGCTGGCCCGGCGTTTGTGCCAGGGTGTCGAAGCCGTGGTGCATTgcgc                        >  5:270194/1‑149 (MQ=255)
                                                                                                                                 aGGCCTGGACGTGCGCGTC‑AGCGGCCGCCGTGCCGAAGGGGTCGAGCATCTGGTCAAGCGTGGCGCACAGTTCATCCCGGGTGATCTGGGTGACGCCGAGCTGGCCCGGCGTTTGTGCCAGGGTGTCGAAGCCGTGGTGCATTGCGCTg                      >  1:304760/1‑149 (MQ=255)
                                                                                                                                 aGGCCTGGACGTGCGCGTC‑AGCGGCCGCCGTGCCGAAGGGGTCGAGCATCTGGTCAAGCGTGGCGCACAGTTCATCCCGGGTGATCTGGGTGACGCCGAGCTGGCCCGGCGTTTGTGCCAGGGTGTCGAAGCCGTGGTGCATTGCGCTg                      >  6:117020/1‑149 (MQ=255)
                                                                                                                                   gCCTGGACGTGCGCGTC‑AGCGGCCGCCGTGCCGAAGGGGTCGAGCATCTGGTCAAGCGTGGCGCACAGTTCATCCCGGGTGATCTGGGTGACGCCGAGCTGGCCCGGCGTTTGTGCCAGGGTGTCGAAGCCGTGGTGCATTGCGCTGgc                    <  8:109857/149‑1 (MQ=255)
                                                                                                                                     cTGGACGTGCGCGTC‑AGCGGCCGCCGTGCCGAAGGGGTCGAGCATCTGGTCAAGCGTGGCGCACAGTTCATCCCGGGTGATCTGGGTGACGCCGAGCTGGCCCGGCGTTTGTGCCAGGGTGTCGAAGCCGTGGTGCAtt                            >  7:238002/1‑139 (MQ=255)
                                                                                                                                     cTGGACGTGCGCGTC‑AGCGGCCGCCGTGCCGAAGGGGTCGAGCATCTGGTCAAGCGTGGCGCACAGTTCATCCCGGGTGATCTGGGTGACGCCGAGCTGGCCCGGCGTTTGTGCCAGGGTGTCGAAGCCGTGGTGCAtt                            <  8:238002/139‑1 (MQ=255)
                                                                                                                                        gACGTGCGCGTC‑AGTGGCCGCCGTGCCGAAGGGGTCGAGCATCTGGTCAAGCGTGGCGCACAGTTCATCCCGGGTGATCTGGGTGACGCCGAGCTGGCCCGGCGTTTGTGCCAGGGTGTCGAAGCCGTGGTGCATTgcg                         <  7:61441/139‑1 (MQ=255)
                                                                                                                                        gACGTGCGCGTC‑AGTGGCCGCCGTGCCGAAGGGGTCGAGCATCTGGTCAAGCGTGGCGCACAGTTCATCCCGGGTGATCTGGGTGACGCCGAGCTGGCCCGGCGTTTGTGCCAGGGTGTCGAAGCCGTGGTGCATTgcg                         >  8:57107/1‑139 (MQ=255)
                                                                                                                                        gACGTGCGCGTC‑AGTGGCCGCCGTGCCGAAGGGGTCGAGCATCTGGTCAAGCGTGGCGCACAGTTCATCCCGGGTGATCTGGGTGACGCCGAGCTGGCCCGGCGTTTGTGCCAGGGTGTCGAAGCCGTGGTGCATTgcg                         >  8:61441/1‑139 (MQ=255)
                                                                                                                                        gACGTGCGCGTC‑AGTGGCCGCCGTGCCGAAGGGGTCGAGCATCTGGTCAAGCGTGGCGCACAGTTCATCCCGGGTGATCTGGGTGACGCCGAGCTGGCCCGGCGTTTGTGCCAGGGTGTCGAAGCCGTGGTGCATTgcg                         <  7:57107/139‑1 (MQ=255)
                                                                                                                                        gACGTGCGCGTC‑AGCGGCCGCCGTGCCGAAGGGGTCGAGCATCTGGTCAAGCGTGGCGCACAGTTCATCCCGGGTGATCTGGGTGACGCCGAGCTGGCCCGGCGTTTGTGCCAGGGTGTCGAAGCCGTGGTGCATTGCGCTGGCGCGGt               <  1:151068/149‑1 (MQ=255)
                                                                                                                                         aCGTGCTCGTC‑AGCGGCCGCCGTGCCGAAGGGGTCGAGCATCTGGTCAAGCGTGGCGCACAGTTCATCCCGGGTGATCTGGGTGACGCCGAGCTGGCCCGGCGTTTGTGCCAGGGTGTCGAAGCCGTGGTGCATTGCGCTGGCGCGGTg              <  2:79630/149‑1 (MQ=255)
                                                                                                                                            tGCGCGTC‑AGCGGCCGCCGTGCCGAAGGGGTCGAGCATCTGGTCAAGCGTGGCGCACAGTTCATCCCGGGTGATCTGGGTGACGCCGAGCTGG‑CCGGCGGTTGGTCCCGGGTGGCCCAGCCCTGGTTGCTTGCGCTGGCGCGGGGGGt           >  4:78619/1‑148 (MQ=255)
                                                                                                                                              cgcgTC‑AGCGGCCGCCGTGCCGAAGGGGTCGAGCATCTGGTCAAGCGTGGCGCACAGTTCATCCCGGGTGATCTGGGTGACGCCGAGCTGGCCCGGCGTTTGTGCCAgg                                                 <  4:32839/109‑1 (MQ=255)
                                                                                                                                              cgcgTC‑AGCGGCCGCCGTGCCGAAGGGGTCGAGCATCTGGTCAAGCGTGGCGCACAGTTCATCCCGGGTGATCTGGGTGACGCCGAGCTGGCCCGGCGTTTGTGCCAgg                                                 >  3:32839/1‑109 (MQ=255)
                                                                                                                                                      gCGGCCGCCGTGCCGAAGGGGTCGAGCATCTGGTCAAGCGTGGCGCACAGTTCATCCCGGGTGATCTGGGTGACGCCGAGCTGGCCCGGCGTTTGTGCCAGGGTGTCGAAGCCGTGGTGCATTGCGCTGGCGCGGTGGGTAACTggggg  <  5:117020/149‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                   
GCCCTTGTAGCCATGGCTGGCGTCAGTTGCATCTGAAATCTGGCAA‑GACGGAGTGTTACATGCGAATTCTGGTCACCGGCGCGAGTGGCTTCAT‑TGGCGGGCGTTTTGCGCGTTTCGCCCTGGAGCAAGGCCTGGACGTGCGCGTC‑AGCGGCCGCCGTGCCGAAGGGGTCGAGCATCTGGTCAAGCGTGGCGCACAGTTCATCCCGGGTGATCTGGGTGACGCCGAGCTGGCCCGGCGTTTGTGCCAGGGTGTCGAAGCCGTGGTGCATTGCGCTGGCGCGGTGGGTAACTGGGGG  >  NC_002947/1060736‑1061031

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