Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,178,888 G→A 8.6% 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.6% 159.9 / 4.9 58P630P (CCC→CCT) PP_1033sulfatase domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base G (27/26);  new base A (2/3);  total (29/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 = 4.62e-01

GCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCC‑ACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCTGACGCCAGTTGGTCAAAGAACGTGGACATGAAGAATGTG  >  NC_002947/1178744‑1179023
                                                                                                                                                 |                                                                                                                                       
gCCCGCTTGTTTGCGGGCTTTCTGGTTGTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTATGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCC‑ACGCGTTGGGGGt                                                                                                                                     <  6:173756/149‑1 (MQ=255)
   cGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCC‑ACGCGTTGGGGGTGCt                                                                                                                                  >  7:192217/1‑149 (MQ=255)
   cGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCC‑ACGCGTTGGGGGTGCt                                                                                                                                  >  8:275372/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                                                                                                                             <  4:145951/147‑1 (MQ=255)
          ttGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCC‑ACGCGTTGGGGGTGCTTTTcc                                                                                                                             >  3:145951/1‑147 (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                                                                                                                               <  1:237898/129‑6 (MQ=255)
                          tCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCCCA‑ACGCGTGGAGGGTgacgag                                                                                                                               >  2:237898/1‑124 (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                                                                        >  3:124157/1‑117 (MQ=255)
                                                                                             ggAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCC‑ACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGta                                                                        <  4:124157/117‑1 (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              >  4:100760/1‑148 (MQ=255)
                                                                                                                        gAACTCGTCACCCTCC‑ACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCTGACGCCAGTTGGTCAAAGAACGTGGAc              <  3:100760/148‑1 (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)
                                                                                                                                                 |                                                                                                                                       
GCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCC‑ACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCTGACGCCAGTTGGTCAAAGAACGTGGACATGAAGAATGTG  >  NC_002947/1178744‑1179023

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