Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,178,879 C→A 6.8% V633V (GTG→GTT)  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,8790CA 6.8% 178.4 / 3.3 59V633V (GTG→GTT) PP_1033sulfatase domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base C (27/28);  new base A (2/2);  total (29/30)
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.01e-01

AGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCC‑ACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCTGACGCCAGTTGGTCAAAGAACGTGGACATGAAGAATGTG  >  NC_002947/1178739‑1179023
                                                                                                                                            |                                                                                                                                                 
aGACAGCACGCATGATTGCGGGCTTTCTTGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGCACTCGTCACCGTCC‑ACGCGTTg                                                                                                                                          <  2:283341/149‑1 (MQ=255)
     gCCCGCTTGTTTGCGGGCTTTCTGGTTGTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTATGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCC‑ACGCGTTGGGGGt                                                                                                                                     <  6:173756/149‑1 (MQ=255)
        cGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCC‑ACGCGTTGGGGGTGCt                                                                                                                                  >  8:275372/1‑149 (MQ=255)
        cGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCC‑ACGCGTTGGGGGTGCt                                                                                                                                  >  7:192217/1‑149 (MQ=255)
        cGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCCCA‑ACGCGTGGAGGGTgac                                                                                                                                  <  1:140153/149‑3 (MQ=255)
         gCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCC‑ACGCGTTGGGGGTGCtt                                                                                                                                 >  6:204191/1‑149 (MQ=255)
         gCATGATTGCGGGATTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCCCCAACGCGTGGAGGGTgac                                                                                                                                  <  4:255117/149‑3 (MQ=255)
            tGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTCCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCC‑ACGCGTTGGGGGTGCTTTTc                                                                                                                              <  5:329882/149‑1 (MQ=255)
               ttGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCC‑ACGCGTTGGGGGTGCTTTTcc                                                                                                                             >  3:145951/1‑147 (MQ=255)
               ttGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCC‑ACGCGTTGGGGGTGCTTTTcc                                                                                                                             <  4:145951/147‑1 (MQ=255)
                 gCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCC‑ACGCGTTGGGGGTGCTTTTCCTGc                                                                                                                          <  7:233000/148‑1 (MQ=255)
                  cGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCC‑ACGCGTTGGGGGTGCTTTTCCTGCgg                                                                                                                        >  2:147427/1‑149 (MQ=255)
                  cGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCC‑ACGCGTTGGGGGTGCTTTTCCTGCgg                                                                                                                        <  7:79431/149‑1 (MQ=255)
                   gggCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCC‑ACGCGTTGGGGGTGCTTTTCCTGCggg                                                                                                                       <  7:52479/149‑1 (MQ=255)
                             gTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCC‑ACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTc                                                                                                               >  1:306561/1‑147 (MQ=255)
                               tCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCCCA‑ACGCGTGGAGGGTgacgag                                                                                                                               >  2:237898/1‑124 (MQ=255)
                               tCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCCCA‑ACGCGTGGAGGGTgacgag                                                                                                                               <  1:237898/129‑6 (MQ=255)
                                 ttGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGGCTGCTACTCACTGCGATGAAATCGTCACCCTCC‑ACGCGTTGGGGTTGCTTTTCCTGAGGGTCATCCTCAAGttt                                                                                                         >  3:202300/1‑149 (MQ=255)
                                 ttGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCC‑ACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGtt                                                                                                         >  7:51798/1‑149 (MQ=255)
                                    ttcgttcAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCC‑ACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTa                                                                                                      <  4:289389/149‑1 (MQ=255)
                                      cgttcAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCC‑ACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTa                                                                                                      <  1:195242/147‑1 (MQ=255)
                                      cgttcAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCC‑ACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTa                                                                                                      >  2:195242/1‑147 (MQ=255)
                                            gATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCC‑ACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGt                                                                                                >  4:280489/1‑147 (MQ=255)
                                                  ttATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCC‑ACGCGTTGGGGGTGCTTTTCATGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGgtgt                                                                                        >  6:254481/1‑149 (MQ=255)
                                                      tcttctGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCC‑ACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAgg                                                                                            <  2:97064/141‑1 (MQ=255)
                                                      tcttctGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTACGATGAACTCGTCACCCTCC‑ACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAgg                                                                                            >  1:97064/1‑141 (MQ=255)
                                                        ttctGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCC‑ACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAAc                                                                                  <  7:275372/149‑1 (MQ=255)
                                                          ctGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCC‑ACGCGTTgg                                                                                                                                         <  6:255104/92‑1 (MQ=255)
                                                          ctGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCC‑ACGCGTTgg                                                                                                                                         >  5:255104/1‑92 (MQ=255)
                                                            ggCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCC‑ACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGaa                                                                              >  7:64745/1‑149 (MQ=255)
                                                             gCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCC‑ACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAAt                                                                             <  2:306561/149‑1 (MQ=255)
                                                                 tCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCC‑ACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTc                                                                           <  8:51798/147‑1 (MQ=255)
                                                                 tCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCC‑ACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCg                                                                          <  5:204191/148‑1 (MQ=255)
                                                                   tGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCC‑ACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGta                                                                        >  4:59028/1‑148 (MQ=255)
                                                                          tAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCC‑ACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTAtt                                                                <  8:192217/149‑1 (MQ=255)
                                                                            gcgcTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCC‑ACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTAttt                                                               <  1:147427/148‑1 (MQ=255)
                                                                            gcgcTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCC‑ACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTg                                                              >  3:237414/1‑149 (MQ=255)
                                                                                            gTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCC‑ACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGaa                                                                                   >  4:81138/1‑112 (MQ=255)
                                                                                            gTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCC‑ACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGaa                                                                                   <  3:81138/112‑1 (MQ=255)
                                                                                            gTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCC‑ACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTccacc                                               >  5:310750/1‑148 (MQ=255)
                                                                                               ataGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCC‑ACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAg                                            <  5:254481/148‑1 (MQ=255)
                                                                                                  ggAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCC‑ACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGta                                                                        <  4:124157/117‑1 (MQ=255)
                                                                                                  ggAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCC‑ACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGta                                                                        >  3:124157/1‑117 (MQ=255)
                                                                                                             tGCTCCTCACTGCGATGAACTCGTCACCCTCC‑ACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCt                                         >  8:250420/1‑137 (MQ=255)
                                                                                                             tGCTCCTCACTGCGATGAACTCGTCACCCTCC‑ACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCt                                         <  7:250420/137‑1 (MQ=255)
                                                                                                               ctcctcACTGCGATGAACTCGTCACCCTCC‑ACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCTGACGCCAGTTGGt                            <  3:59028/148‑1 (MQ=255)
                                                                                                                tgctcACTGCGATGAACTCGTCACCCTCC‑ACGCGTTGGGGTTGCTTTTCGTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCTGACGc                                    <  8:80694/137‑1 (MQ=255)
                                                                                                                tcctcACTGCGATGAACTCGTCACCCTCC‑ACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTAGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCTGACGc                                    >  7:80694/1‑139 (MQ=255)
                                                                                                                  ctcACTGCGATTAACTCGTCACACTCC‑ACGCGTTGGGGGTGCTTTTCCTGCGTGTCATCCTCCAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACTGGTGCTCCACCAAGGCTGACGCCAGTTGGTCaaa                        <  4:202300/149‑1 (MQ=255)
                                                                                                                  ctcACTGCGATGAACTCGTCACCCTCC‑CCGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGACAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCTGACGCCAGTTGGTCa                          >  6:269962/1‑147 (MQ=255)
                                                                                                                          gATGAACTCGTCACCCTCC‑ACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCTGACGCCAGTTGGTCAAAGAACGTg                 >  5:225015/1‑148 (MQ=255)
                                                                                                                          gATGAACTCGTCACCCTCC‑ACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCTGAAGCCAGTTGGTCAAAGAACGTg                 >  5:119791/1‑148 (MQ=255)
                                                                                                                             gAACTCGTCACCCTCC‑ACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCTGACGCCAGTTGGTCAAAGAACGTGGAc              <  3:100760/148‑1 (MQ=255)
                                                                                                                             gAACTCGTCACCCTCC‑ACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCTGACGCCAGTTGGTCAAAGAACGTGGAc              >  4:100760/1‑148 (MQ=255)
                                                                                                                               gaaaagcACCCCCA‑ACGCGTGGAGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCTGACGCCAGTTGGTCAAAGAACGTGGACAt            >  6:5342/7‑148 (MQ=255)
                                                                                                                               aCTCGTCACCCTCC‑ACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCTGACGCCAGTTGGTCAAAGAACGTGGACATg           <  3:76660/149‑1 (MQ=255)
                                                                                                                                   gTCACCCTCC‑ACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCTGACGCCAGTTGGTCAAAGAATGTGGACATgaa         >  6:296074/1‑147 (MQ=255)
                                                                                                                                   gTCACCCTCC‑ACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCTGACGCCAGTTGGTCAAAGAATGTGGACATgaa         <  5:296074/147‑1 (MQ=255)
                                                                                                                                        ccTCC‑ACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCTGACGCCAGTTGGTCAAAGAACGTGGACATGAAGAAtgtg  >  1:57896/1‑149 (MQ=255)
                                                                                                                                           cc‑ACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCTGACGCCAGTTGGTCAAAGAACGt                  >  1:117071/1‑130 (MQ=255)
                                                                                                                                           cc‑ACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCTGACGCCAGTTGGTCAAAGAACGt                  <  2:117071/130‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                                 
AGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCC‑ACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCTGACGCCAGTTGGTCAAAGAACGTGGACATGAAGAATGTG  >  NC_002947/1178739‑1179023

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