Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,460,087 T→C 6.3% L371S (TTG→TCG)  trkH‑I → potassium transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,460,0870TC 6.3% 181.4 / 2.6 64L371S (TTG→TCG) trkH‑Ipotassium transporter
Reads supporting (aligned to +/- strand):  ref base T (32/27);  new base C (2/2);  total (34/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 = 9.91e-01

GGTGATTACACACTCTGGGGCAGCTTCGCCATCATGCTGTTCTTTTACCTGACGTTCGTCGGCGGGTGCTCGGGCTCCACCGCCGGCGGTTTGAAAATATTCCGTTTCCAGGTGGCCGCGGCCTTGCTTGTCAGCAGCCTGAAGCAGTTGATC‑CATCCACGGGCAACGATCTCGAAGAAATACAATGGGCACCCCATCGATGAAGAGATCGTACGGTCACTGCTGACGTTCTCGTTCTTCTTCACCGTCACCATTGCGGTCATCGCCCTGGGCCTGTCTTTGATAGGACTCGACTG  >  NC_002947/4459939‑4460234
                                                                                                                                                    |                                                                                                                                                    
ggTGATTACACACTCTGGGGCAGCTTCGCCATCATGCTGTTCTTTTACCTGACGTTCGTCGGCGGGTGCTCGGGCTCCACCGCCGGCGGTTTGAAAATATTCCGTTTCCAGGTGGCCGCGGCCTTGCTTGTCAGCAGCCTGAAGCAGtt                                                                                                                                                      <  6:239106/149‑1 (MQ=255)
       acacacTCTGGGGCAGCTTCGCCATCATGCTGTTCTTTTACCTGACGTTCGTCGGCGGGTGCTCGGGCTCCACCGCCGGCGGTTTGAAAATATTCCGTTTCCAGGTGGCCGCGGCCTTGCTTGTCAGCAGCCTGAAGCAGTTGatc‑cat                                                                                                                                              >  5:134990/1‑149 (MQ=255)
         acacTCTGGGGCAGCTTCGCCATCATGCTGTTCTTTTACCTGACGTTCGTCGGCGGGTGCTCGGGCTCCACCGCCGGCGGTTTGAAAATATTCCGTTTCCAGGTGGCCGCGGCCTTGCTTGTCAGCAGCCTGAAGCAGTTGatc‑catc                                                                                                                                             >  4:91438/1‑148 (MQ=255)
               tGGGGCAGCTTCGCCATCATGCTGTTCTTTTACCTGACGTTCGTCGGCGGGTGCTCGGGCTCCACCGCCGGCGGTTTGAAAATATTCCGTTTCCAGGTGGCCGCGGCCTTGCTTGTCAGCAGCCTGAAGCAGTTGATC‑CATCCACgg                                                                                                                                        <  3:86739/147‑1 (MQ=255)
                 gggCAGCTTCGCCATCATGCTGTTCTTTTACCTGACGTTCGTCGGCGGGTGCTCGGGCTCCACCGCCGGCGGTTTGAAAATATTCCGTTTCCAGGTGGCCGCGGCCTTGCTTGTCAGCAGCCTGAAGCAGTTGatc‑cat                                                                                                                                              <  4:210219/139‑1 (MQ=255)
                 gggCAGCTTCGCCATCATGCTGTTCTTTTACCTGACGTTCGTCGGCGGGTGCTCGGGCTCCACCGCCGGCGGTTTGAAAATATTCCGTTTCCAGGTGGCCGCGGCCTTGCTTGTCAGCAGCCTGAAGCAGTTGatc‑cat                                                                                                                                              >  3:210219/1‑139 (MQ=255)
                 gggCAGCTTCGCCATCATGCTGTTCTTTTACCTGACGTTCGTCGGCGGGTGCTCGGGCTCCACCGCCGGCGGTTTGAAAATATTCCGTTTCCAGGTGGCCGCGGCCTTGCTTGTCAGCAGCCTGAAGCAGTTGATC‑CATCCACggg                                                                                                                                       <  5:52777/146‑1 (MQ=255)
                   gCAGCTTCGCCATCATGCTGTTCTTTTACCTGACGTTCGTCGGCGGGTGCTCGGGCTCCACCGCCGGCGGTTTGAAAATATTCCGTTTCCAGGTGGCCGCGGCCTTGCTTGTCAGCAGCCTGAAGCAGTTGATC‑CATCCACGGGCAAc                                                                                                                                   >  7:140392/1‑148 (MQ=255)
                     aGCTTCGCCATCATGCTGTTCTTTTACCTGACGTTCGTCGGCGGGTGCTCGGGCTCCACCGCCGGCGGTTTGAAAATATTCCGTTTCCAGGTGGCCGCGGCCTTGCTTGTCAGCAGCCTGAAGCAGTTGATC‑CATCCACGGGCAACGAt                                                                                                                                <  8:140392/149‑1 (MQ=255)
                        ttCGCCATCATGCTGTTCTTTTACCTGACGTTCGTCGGCGGGTGCTCGGGCTCCACCGCCGGCGGTTTGAAAATATTCCGTTTCCAGGTGGCCGCGGCCTTGCTTGTCAGCAGCCTGAAGCAGTTGATC‑CATCCACGGGCAACGAtc                                                                                                                               <  3:338373/147‑1 (MQ=255)
                        ttCGCCATCATGCTGTTCTTTTACCTGACGTTCGTCGGCGGGTGCTCGGGCTCCACCGCCGGCGGTTTGAAAATATTCCGTTTCCAGGTGGCCGCGGCCTTGCTTGTCAGCAGCCTGAAGCAGTTGATC‑CATCCACGGGCAACGAtc                                                                                                                               >  4:338373/1‑147 (MQ=255)
                           gCCATCATGCTGTTCTTTTACCTGACGTTCGTCGGCGGGTGCTCGGGCTCCACCGCCGGCGGTTTGAAAATATTCCGTTTCCAGGTGGCCGCGGCCTTGCTTGTCAGCAGCCTGAAGCAGTTGATC‑CATCCACGGGCAACGATCTCgaa                                                                                                                          <  7:255734/149‑1 (MQ=255)
                           gCCATCATGCTGTTCTTTTACCTGACGTTCGTCGGCGGGTGCTCGGGCTCCACCGCCGGCGGTTTGAAAATATTCCGTTTCCAGGTGGCCGCGGCCTTGCTTGTCAGCAGCCTGAAGCAGTTGATC‑CATCCACGGGCAACGATCTCg                                                                                                                            >  5:290065/1‑147 (MQ=255)
                             catcatGCTGTTCTTTTACCTGACGTTCGTCGGCGGGTGCTCGGGCTCCACCGCCGGCGGTTTGAAAATATTCCGTTTCCAGGTGGCCGCGGCCTTGCTTGTCAGCAGCCTGAAGCAGTTGATC‑CATCCACGGGCAACGATCTCgaaga                                                                                                                        <  3:292521/149‑1 (MQ=255)
                                   gCTGTTCTTTTACCTGACGTTCGTCGGCGGGTGCTCGGGCTCCACCGCCGGCGGTTTGAAAATATTCCGTTTCCAGGTGGCCGCGGCCTTGCTTGTCAGCAGCCTGAAGCAGTTGATC‑CATCCACGGGCAACGATCTCGAAGAAATAc                                                                                                                   >  1:337134/1‑148 (MQ=255)
                                     tGTTCTTTTACCTGACGTTCGTCGGCGGGTGCTCGGGCTCCACCGCCGGCGGTTTGAAAATATTCCGTTTCCAGGTGGCCGCGGCCTTGCTTGTCAGCAGCCTGAAGCAGTTGATC‑CATCCACGGGCAACGATCTCGAAGAAATACAAt                                                                                                                >  1:36764/1‑149 (MQ=255)
                                     tGTTCTTTTACCTGACGTTCGTCGGCGGCTGCTCGGGCTCCACCGCCGGCGGTTTGAAAATATTCCGTTTCCAGGTGGCCGCGGCCTTGCTTTTCAGAAGCCTGAAGAAGTTGATC‑CCTCCACGGGCAACGATCTCGAAGAAATACAAt                                                                                                                >  7:77424/1‑149 (MQ=255)
                                       ttCTTTTACCTGACGTTCGTCGGCGGGTGCTCGGGCTCCACCGCCGGCGGTTTGAAAATATTCCGTTTCCAGGTGGCCGCGGCCTTGCTTGTCAGCAGCCTGAAGCAGTTGATC‑CATCCACgg                                                                                                                                        <  2:326051/123‑1 (MQ=255)
                                       ttCTTTTACCTGACGTTCGTCGGCGGGTGCTCGGGCTCCACCGCCGGCGGTTTGAAAATATTCCGTTTCCAGGTGGCCGCGGCCTTGCTTGTCAGCAGCCTGAAGCAGTTGATC‑CATCCACgg                                                                                                                                        >  1:326051/1‑123 (MQ=255)
                                       ttCTTTTACCTGACGTTCGTCGGCGGGTGCTCGGGCTCCACCGCCGGCGGTTTGAAAATATTCCGTTTCCAGGTGGCCGCGGCCTTGCTTGTCAGCAGCCTGAAGCAGTTGATC‑CATCCACGGGCAACGATCTCGAAGAAATACAATgg                                                                                                              <  7:310265/149‑1 (MQ=255)
                                       ttCTTTTACCTGACGTTCGTCGGCGGGTGCTCGGGCTCCACCGCCGGCGGTTTGAAAATATTCCGTTTCCAGGTGGCCGCGGCCTTGCTTGTCAGCAGCCTGAAGCAGTTGATC‑CATCCACGGGCAACGATCTCGAAGAAATACAATg                                                                                                               <  6:134990/148‑1 (MQ=255)
                                                          tcgGCGGGTGCTCGGGCTCCACCGCCGGCGGTTTGAAAATATTCCGTTTCCAGGTGGCCGCGGCCTTGCTTGTCAGCAGCCTGAAGCAGTTGATC‑CATCCACGGGCAACGAtctc                                                                                                                             >  4:170638/1‑115 (MQ=255)
                                                          tcgGCGGGTGCTCGGGCTCCACCGCCGGCGGTTTGAAAATATTCCGTTTCCAGGTGGCCGCGGCCTTGCTTGTCAGCAGCCTGAAGCAGTTGATC‑CATCCACGGGCAACGAtctc                                                                                                                             <  3:170638/115‑1 (MQ=255)
                                                                  tGCTCGGGCTCCACCGCCGGCGGTTTGAAAATATTCCGTTTCCAGGTGGCCGCGGCCTTGCTTGTCAGCAGCCTGAAGCAGTTGATC‑CATCCACGGGCAACGATCTCGAAGAAATACAATGGGCACCCCATCGATGAAGAGATCGTACg                                                                                   <  6:290065/149‑1 (MQ=255)
                                                                   gCTCGGGCTCCACCGCCGGCGGTTTGAAAATATTCCGTTTCCAGGTGGCCGCGGCCTTGCTTGTCAGCAGCCTGAAGCAGTTGATC‑CATCCACGGGCAACGATCTCGAAGAAATACAATGGGCACCCCATCGATGAAGAGATCGTACgg                                                                                  >  5:188666/1‑149 (MQ=255)
                                                                       gggCTCCACCGCCGGCGGTTTGAAAATATTCCGTTTCCAGGTGGCCGCGGCCTTGCTTGTCAGCAGCCTGAAGCAGTTGATC‑CATCCACGGGCAACGATCTCGAAGAAATACAATGGGCACCCCATCGATGAAGAGCTCGTACGGTCAc                                                                              >  8:310006/1‑149 (MQ=255)
                                                                                cgccgGCGGTTTGAAAATATTCCGTTTCCAGGTGGCCGCGGCCTTGCTTGTCAGCAGCCTGAAGCAGTTGATC‑CATCCACGGGCAACGATCTCGAAGAAATACAATGGGCAc                                                                                                          >  5:19481/1‑112 (MQ=255)
                                                                                cgccgGCGGTTTGAAAATATTCCGTTTCCAGGTGGCCGCGGCCTTGCTTGTCAGCAGCCTGAAGCAGTTGATC‑CATCCACGGGCAACGATCTCGAAGAAATACAATGGGCAc                                                                                                          <  6:19481/112‑1 (MQ=255)
                                                                                 gccgGCGGTTTGAAAATATTCCGTTTCCAGGTGGCCGCGGCCTTGCTTGTCAGCAGCCTGAAGCAGTTGATC‑CATCCACGGGCAACGATCTCGAAGAAATACAATGGGCACCCCATCGATGAAGAGATCGTACGGTCACTGCTGAcgtt                                                                    <  2:337134/149‑1 (MQ=255)
                                                                                    ggcggTTTGAAAATATTCCGTTTCCAGGTGGCCGCGGCCTTGCTTGTCAGCAGCCTGAAGCAGTTGATC‑CATCCACGGGCAACGATCTCGAAGAAATACAATGGGCACCCCATCGATGAAGAGATCGTACGGTCACTGCTGAcgttctc                                                                 >  5:283014/1‑149 (MQ=255)
                                                                                    ggcggTTTGAAAATATTCCGTTTCCAGGAGGCACCGGCCTTGCTTGTCATCAGCCTGAAGCAGTGGATC‑CATCCACTGGCAACGATCTCGAAGAAATACAATCGGCACCCCATCTATGAATAGATCGTACGGTCACTGCTGAcgttctc                                                                 <  8:77424/149‑1 (MQ=255)
                                                                                        gTTTGAAAATTTTCCGTTTCCAGGTGGCCGCGGCCTTGCTTGTCAGCAGCCTGAAGCAGTTGATC‑CATCCACGGGCAACGATCTCGAAGAAATACAATGGGCACCCCATCGATGAAGAGATCGTACGGTCACTGCTGAcgttctcgttc                                                             <  6:211355/149‑1 (MQ=255)
                                                                                        gTTTGAAAATATTCCGTTTCCAGGTGGCCGCGGCCTTGCTTGTCAGCAGCCTGAAGCAGTTGATC‑CATCCACGGGCAACGATCTCGAAGAAATACAATGGGCACCCCATCGATGAAGAGATCGTACGGTCACTGCTGAcgttctcgttc                                                             <  3:91438/149‑1 (MQ=255)
                                                                                        gTTTGAAAATATTCCGTTTCCAGGTGGCCGCGGCCTTGCTTGTCAGCAGCCTGAAGCAGTTGATC‑CATCCACGGGCAACGATCTCGAAGAAATACAATGGGCACCCCATCGATGAAGAGATCGTACGGTCACTGCTGAcgttctcgttc                                                             <  3:333086/149‑1 (MQ=255)
                                                                                            gAAAATATTCCGTTTCCAGGTGGCCGCGGCCTTGCTTGTCAGCAGCCTGAAGCAGTTGATC‑CATCCACGGGCAACGATCTCGAAGAAATACAATGGGCACCCCATCGATGAAGAGATCGTACGGTCACTGCTGACGTTCTCGttcttc                                                          <  6:188666/148‑1 (MQ=255)
                                                                                               aaTATTCCGTTTCCAGGTGGCCGCGGCCTTGCTTGTCAGCAGCCTGAAGCAGTTGATC‑CATCCACGGGCAACGATCTCGAAGAAATACAATGGGCACCCCATCGATGAAGAGATCGTACGGTCACTGCTGACGTTCTCGTTCTTCTTCa                                                      >  5:64895/1‑149 (MQ=255)
                                                                                                atatTCCGTTTCCAGGTGGCCGCGGCCTTGCTTGTCAGCAGCCTGAAGCAGTTGATC‑CATCCACGGGCAACGATCTCGAAGAAATACAATGGGCACCCCATCGATGAAGAGATCGTACGGTCACTGCTGACGTTCTCGTTCTTCTTCAc                                                     <  2:36764/149‑1 (MQ=255)
                                                                                                      cGTTTCCAGGTGGCCGCGGCCTTGCTTGTCAGCAGCCTGAAGCAGTTGATC‑CATCCACGGGCAACGATCTCGAAGAAATACAATGGGCACCCCATCGATGAAGAGATCGTACGGTCACTGCTGACGTTCTCGTTCTTCTTCACCGTCa                                                <  6:64895/148‑1 (MQ=255)
                                                                                                      cGTTTCCAGGTGGCCGCGGCCTTGCTTGTCAGCAGCCTGAAGCAGTTGATC‑CATCCACGGGCAACGATCTCGAAGAAATACAATGGGCACCCCATCGATGAAGAGATCGTACGGTCACTGCTGACGTTCTCGTTCTTCTTCACCGTCAc                                               >  3:38092/1‑149 (MQ=255)
                                                                                                       gTTTCCAGGTGGCCGCGGCCTTGCTTGTCAGCAGCCTGAAGCAGTTGATC‑CATCCACGGGCAACGATCTCGAAGAAATACAATGGGCAGCCCATCGATGAAGAGCTCGTACGGTCACTGCTGACGTTCTCGTGCTTCTTCACCGTCAc                                               >  7:185453/1‑148 (MQ=255)
                                                                                                       gTTTCCAGGTGGCCGCGGCCTTGCTTGTCAGCAGCCTGAAGCAGTCGATC‑CATCCACGGGCAACGATCTCGAAGAAATACAATGGGCACCCCATCGATGAAGAGATCGTACGGTCAct                                                                             <  8:32487/118‑1 (MQ=255)
                                                                                                       gTTTCCAGGTGGCCGCGGCCTTGCTTGTCAGCAGCCTGAAGCAGTCGATC‑CATCCACGGGCAACGATCTCGAAGAAATACAATGGGCACCCCATCGATGAAGAGATCGTACGGTCAct                                                                             <  8:38076/118‑1 (MQ=255)
                                                                                                       gTTTCCAGGTGGCCGCGGCCTTGCTTGTCAGCAGCCTGAAGCAGTCGATC‑CATCCACGGGCAACGATCTCGAAGAAATACAATGGGCACCCCATCGATGAAGAGATCGTACGGTCAct                                                                             >  7:38076/1‑118 (MQ=255)
                                                                                                       gTTTCCAGGTGGCCGCGGCCTTGCTTGTCAGCAGCCTGAAGCAGTCGATC‑CATCCACGGGCAACGATCTCGAAGAAATACAATGGGCACCCCATCGATGAAGAGATCGTACGGTCAct                                                                             >  7:32487/1‑118 (MQ=255)
                                                                                                        tttCCAGGTGGCCGCGGCCTTGCTTGTCAGCAGCCTGAAGCAGTTGATCGCATCCACGGGCAACGATCTCGAAGAAATACAATGGGCACCCCATCGATGAAGAGATCGTACGGTCACTGCTGACGTTCTCGTTCTTCTTCACCGTCAc                                               >  1:84034/1‑148 (MQ=255)
                                                                                                           ccAGGTGGCCGCGGCCTTGCTTGTCAGCAGCCTGAAGCAGTTGATC‑CATCCACGGGCAACGATCTCGAAGAAATACAATGGGCACCCCATCGATGAAGAGATCGTACGGTCACTGCTGAc                                                                       <  8:311752/120‑1 (MQ=255)
                                                                                                           ccAGGTGGCCGCGGCCTTGCTTGTCAGCAGCCTGAAGCAGTTGATC‑CATCCACGGGCAACGATCTCGAAGAAATACAATGGGCACCCCATCGATGAAGAGATCGTACGGTCACTGCTGAc                                                                       >  7:311752/1‑120 (MQ=255)
                                                                                                             aGGTGGCCGCGGCCTTGCTTGTCAGCAGCCTGAAGCAGTTGATC‑CATCCACGGGCAACGATCTCGAAGAAATACAATGGGCACCCCATCGATGAAGAGATCGTACGGTCACTGCTGAcgttctcgttc                                                             >  7:16989/1‑128 (MQ=255)
                                                                                                             aGGTGGCCGCGGCCTTGCTTGTCAGCAGCCTGAAGCAGTTGATC‑CATCCACGGGCAACGATCTCGAAGAAATACAATGGGCACCCCATCGATGAAGAGATCGTACGGTCACTGCTGAcgttctcgttc                                                             <  8:16989/128‑1 (MQ=255)
                                                                                                             aGGTGGCCGCGGCCTTGCTTGTCAGCAGCCTGAAGCAGTTGATC‑CATCCACGGGCAACGATCTCGAAGAAATACAATGGGCACCCCATCGATGAAGAGATCGTACGGTCACTGCTGACGTTCTCGTTCTTCTTCACCGTCACCATTGCg                                        >  6:295457/1‑149 (MQ=255)
                                                                                                                tGGCCGCGGCCTTGCTTGTCAGCAGCCTGAAGCAGTTGATC‑CATCCACGGGCAACGATCTCGAAGAAATACAATGGGCACCCCATCGATGAAGAGATCGTACGGTCACTGCTGAcgttc                                                                   >  2:358303/1‑119 (MQ=255)
                                                                                                                tGGCCGCGGCCTTGCTTGTCAGCAGCCTGAAGCAGTTGATC‑CATCCACGGGCAACGATCTCGAAGAAATACAATGGGCACCCCATCGATGAAGAGATCGTACGGTCACTGCTGAcgttc                                                                   <  1:358303/119‑1 (MQ=255)
                                                                                                                    cgcgGCCTTGCTTGTCAGCAGCCTGAAGCAGTTGATC‑CATCCACGGGCAACGATCTCGAAGAAATACAATGGGCACCCCATCGATGAAGAGATCGTACGGTCACTGCTGACGTTCTCGTTCTTCTTCACCGTCACCATTGCGGTCATCg                                 >  2:184762/1‑149 (MQ=255)
                                                                                                                     gcgGCCTTGCTTGTCAGCAGCCTGAAGCAGTTGATC‑CATCCACGGGCAACGATCTCGAAGAAATACAATGGGCACCCCATCGATGATGAGATCGTACGGTCACTGCTGACGTTCTCGTTCTTCTTCACCGTCACCATTGCGGTCATCGc                                >  4:219610/1‑149 (MQ=255)
                                                                                                                     gcgGCCTTGCTTGTCAGCAGCCTGAAGCAGTTGATC‑CATCCACGGGCAACGATCTCGAAGAAATACAATGGGCACCCCATCGATGAAGAGATCGTACGGTCAc                                                                              >  7:255318/1‑103 (MQ=255)
                                                                                                                     gcgGCCTTGCTTGTCAGCAGCCTGAAGCAGTTGATC‑CATCCACGGGCAACGATCTCGAAGAAATACAATGGGCACCCCATCGATGAAGAGATCGTACGGTCAc                                                                              <  8:255318/103‑1 (MQ=255)
                                                                                                                         ccTTGCTTGTCAGCAGCCTGAAGCAGTTGATC‑CATCCACGGGCAACGATCTCGAAGAAATACAATGGGCACCCCATCGATGAAGAGATCGTACGGTCACTGCTGACGTTCTCGTTCTTCTTCACCGTCACCATTGCGGTCATCGCCCTg                            >  2:336833/1‑149 (MQ=255)
                                                                                                                             gcttgTCAGCAGCCTGAAGCAGTTGATC‑CATCCACGGGCAACGATCTCGAAGAAATACAATGGGCACCCCATCGATGAAGAGATCGTACGGTCACTGCTGACGTTCTCGTTCTTCTTCACCGTCACCATTGCGGTCATCGCCCTGGGc                         <  5:295457/148‑1 (MQ=255)
                                                                                                                             gcttgTCAGCAGCCTGAAGCAGTTGATC‑CATCCACGGGCAACGATCTCGAAGAAATACAATGGGCACCCCATCGATGAAGAGATCGTACGGTCACTGCTGACGTTCTCGTTCTTCTTCACCGTCACCATTGCGGTCATCGCCCTGGGc                         <  1:184762/148‑1 (MQ=255)
                                                                                                                               ttgTCAGCAGCCTGAAGCAGTTGATC‑CATCCACGGGCAACGATCTCGAAGAAATACAATGGGCACCCCATCGATGAAGAGATCGTACGGTCACTGCTGCCGTTCTCGTTCTTCTTCTCCGTCACCATTGCGGTCATCGCACTGGGCCTg                      >  6:278609/1‑149 (MQ=255)
                                                                                                                                     gcagcCTGAAGCAGTTGATC‑CATCCACGGGCAACGATCTCGAAGAAATACAATGGGCACCCCATCGATGAAGAGATCGTACGGTCACTGCTGACGTTCTCGTTCTTCTTCACCGTCACCATTGCGGTCATCGCCCTGGGCCTGTCttt                 >  8:163571/1‑148 (MQ=255)
                                                                                                                                       agcCTGAAGCAGTTGATC‑CATCCACGGGCAACGATCTCGAAGAAATACAATGGGCACCCCATCGATGAAGAGATCGTACGGTCACTGCTGACGTTCTCGTTCTTCTTCACCGTCACCATTGCGGTCATCGCCCTGGGCCTGTCTTTGAt              <  7:163571/149‑1 (MQ=255)
                                                                                                                                            gAAGCAGTTGATC‑CATCCACGGGCAACGATCTCGAAGAAATACAATGGGCACCCCATCGATGAAGAGATCGTACGGTCACTGCTGACGTTCTCGTTCTTCTTCACCGTCACCATTGCGGTCATCGCCCTGGGCCTGTCTTTGATAGGAc         >  7:317716/1‑149 (MQ=255)
                                                                                                                                               gCAGTTGATC‑CATCCACGGGCAACGATCTCGAAGAAATACAATGGGCACCCCATCGATGAAGAGATCGTACGGTCACTGCTGACGTTCTCGTTCTTCTTCACCGTCACCATTGCGGTCATCGCCCTGGGCCTGTCTTTGATAGGACTc       >  7:253242/1‑148 (MQ=255)
                                                                                                                                                  gTTGATC‑CATCCACGGGCAACGATCTCGAAGAAATACAATGGGCACCCCATCGATGAAGAGATCGTACGGTCACTGCTGACGTTCTCGTTCTTCTTCACCGTCACCATTGCGGTCATCGCCCTGGGCCTGTCTTTGATAGGACTCGACt   >  8:194965/1‑149 (MQ=255)
                                                                                                                                                   ttGATC‑CATCCACGGGCAACGATCTCGAAGAAATACAATGGGCACCCCATCGATGAAGAGATCGTACGGTCACTGCTGACGTTCTCGTTCTTCTTCACCGTCACCATTGCGGGCATCGCCCTGGGCCTGTCTTTGATAGGACTCGACTg  <  6:283014/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GGTGATTACACACTCTGGGGCAGCTTCGCCATCATGCTGTTCTTTTACCTGACGTTCGTCGGCGGGTGCTCGGGCTCCACCGCCGGCGGTTTGAAAATATTCCGTTTCCAGGTGGCCGCGGCCTTGCTTGTCAGCAGCCTGAAGCAGTTGATC‑CATCCACGGGCAACGATCTCGAAGAAATACAATGGGCACCCCATCGATGAAGAGATCGTACGGTCACTGCTGACGTTCTCGTTCTTCTTCACCGTCACCATTGCGGTCATCGCCCTGGGCCTGTCTTTGATAGGACTCGACTG  >  NC_002947/4459939‑4460234

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