Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,006,162 G→A 6.3% Q572* (CAG→TAG)  PP_3535 ← gamma‑glutamyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,006,1620GA 6.3% 188.0 / 2.4 63Q572* (CAG→TAG) PP_3535gamma‑glutamyltransferase
Reads supporting (aligned to +/- strand):  ref base G (30/29);  new base A (2/2);  total (32/31)
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.25e-01

TTCGATTTCACGCTGATAGACATCTAGATATTCAACTACTGAATGGCAAGTCTGCGTTCTATCACGGCTAGTCACCGAGGGCCACTCCTTCCCGCCGGGGATCGGCACCGCCTGACCAACCCTCGCCGCTGCGCTGGATGATCTGCATGCCGCTGGTCATGGTCATCGGTGTCACCTCATGGCCCCAGGCTGCCAACTGGCGGATCAGCGCCGGGCTGCTCAGCCCTGCTTCCACTTCGGTGCCGGCATTGCGGCTGCCGAAGTTGGGCAGGTCGGCGGCCTGCTGCGGGT  >  NC_002947/4006018‑4006308
                                                                                                                                                |                                                                                                                                                  
ttCGATTTCACGCTGATAGACATCTAGATATTCAACTACTGAATGGCAAGTCTGCGTTCTATCACGGCTAGTCACCGAGGGCCACTCCTTCCCGCCGGGGATCGGCACCGCCTGACCAACCCTCGCCGCTGCGCTGGATGATCTGCATg                                                                                                                                                >  7:23438/1‑149 (MQ=255)
  cGATTTCACGCTGATAGACATCTAGATATTCAACTACTGAATGGCAAGTCTGCGTTCTATCACGGCTAGTCACCGAGGGCCACTCCTTCCCGCCGGGGATCGGCACCGCCTGACCAACCCTCGCCGCTGCGCTGGATGATCTGCATGc                                                                                                                                               <  8:40427/148‑1 (MQ=255)
   gATTTCACGCTGATAGACATCTAGATATTCAACTACTGAATGGCAAGTCTGCGTTCTATCACGGCTAGTCACCGAGGGCCACTCCTTCCCGCCGGGGATCGGCACCGCCTGACCAACCCTCGCCGCTGCGCTGGATGATCTGCATGCCg                                                                                                                                             >  8:132008/1‑149 (MQ=255)
       tCACGCTGATAGACATCTAGATATTCAACTACTGAATGGCAAGTCTGCGTTCTATCACGCCAAGTCACCGAGGGCCACTCCTTCCCGCCGGGGATCGGCACCGCCTGACCAACCCTCGCCGCTGCGCTGGATGATATGCATGCCGCtgg                                                                                                                                         >  7:119726/1‑149 (MQ=255)
                tAGACATCTAGATATTCAACTACTGAATGGCAAGTCTGCGTTCTATCACGGCTAGTCACCGAGGGCCACTCCTTCCCGCCGGGGATCGGCACCGCCTGACCAACCCTCGCCGCTGCGCTGGATGATCTGCATGCCGCtggtcatggtca                                                                                                                                <  1:204357/149‑1 (MQ=255)
                       cTAGATATTCAACTACTGAATGGCAAGTCTGCGTTCTATCACGGCTAGTCACCGAGGGCCACTCCTTCCCGCCGGGGATCGGCACCGCCTGACCAACCCTCGCCGCTGCGCTGGATGATCTGCATGCCGCTGGTCATGGTCATCGgtg                                                                                                                          >  8:181869/1‑148 (MQ=255)
                        tAGATATTCAACTACTGAATGGCAAGTCTGCGTTCTATCACGGCTAGTCACCGAGGGCCACTCCTTCCCGCCGGGGATCGGCACCGCCTGACCAACCCTCGCCGCTGCGCTGGATGATCTGCATGCCGCTGGTCATGGTCATCGGTGTc                                                                                                                        >  7:161175/1‑149 (MQ=255)
                                     actGAATGGCAAGTCTGCGTTCTATCACGGCTAGTCACCGAGGGCCACTCCTTCCCGCCGGGGATCGGCACCGCCTGACCAACCCTCGCCGCTGCGCTGGATGATCTGCATGCCGCTGGTCATGGTCATCGGTGTCACCTCATGGcccc                                                                                                           <  6:180351/149‑1 (MQ=255)
                                      ctGAATGGCAAGTCTGCGTTCTATCACGGCTAGTCACCGAGGGCCACTCCTTCCCGCCGGGGATCGGCACCGCCTGACCAACCCTCGCCGCTGCGCTGGATGATCTGCATGCCGCTGGTCATGGTCATCGGTGTCACCTCATGGCCCCa                                                                                                          <  8:161175/149‑1 (MQ=255)
                                         aaTGGCAAGTCTGCGTTCTATCACGGCTAGTCACCGAGGGCCACTCCTTCCCGCCGGGGATCGGCACCGCCTGACCAACCCTCGCCGCTGCGCTGGATGATCTGCATGCCGCTGGTCATGGTCATCGGTGTCACCTCATGGCCCCAGGc                                                                                                       <  2:181132/149‑1 (MQ=255)
                                          aTGGCAAGTCTGCGTTCTATCACGGCTAGTCACCGAGGGCCACTCCTTCCCGCCGGGGATCGGCACCGCCTGACCAACCCTCGCCGCTGCGCTGGATGATCTGCATGCCGCTGGTCATGGTCATCGGTGTCACCTCATGGCCCCAGGCt                                                                                                      <  5:290004/149‑1 (MQ=255)
                                            ggCAAGTCTGCGTTCTATCACGGCTAGTCACCGAGGGCCACTCCTTCCCGCCGGGGATCGGCACCGCCTGACCAACCCTCGCCGCTGCGCTGGATGATCTGCATGCCGCTGGTCATGGTCATCGGTGTCACCTCATGGCCCCAGGCTGc                                                                                                    >  3:218460/1‑149 (MQ=255)
                                             gCAAGTCTGCGTTCTATCACGGCTAGTCACCGAGGGCCACTCCTTCCCGCCGGGGATCGGCACCGCCTGACCAACCCTCGCCGCTGCGCTGGATGATCTGCATGCCGCTGGTCATGGTCATCGGTGTCACCTCATGGCCCCAGGCTGc                                                                                                    >  3:250148/1‑148 (MQ=255)
                                               aaGTCTGCGTTCTATCACGGCTAGTCACCGAGGGCCACTCCTTCCCGCCGGGGATCGGCACCGCCTGACCAACCCTCGCCGCTGCGCTGGATGATCTGCATGCCGCTGGTCATGGTCATCGGTGTCACCTCATGGCCCCAGGCTGCCaa                                                                                                 <  1:235876/149‑1 (MQ=255)
                                                  tCTGGGTTCTTTCACGGCTAGTCACCGAGGGCCACTCCTTACCTCCGGGGATCGGCACCGCCTGACCAACCCTCGCCGCTGCGCTGGATGATCTGCATGCCGCTGGTGATGGTCATCGGTGTCACCTCATGGCCCCAGGCTGCGAACt                                                                                               <  8:119726/148‑1 (MQ=255)
                                                    tGCGTTCTATCACGGCTAGTCACCGAGGGCCACTCCTTCCCGCCGGGGATCGGCACCGCCTGACCAACCCTCGCCGCTGCGCTGGATGATCTGCATGCCGCTGGTCATGGTCATCGGTGTCACCTCATGGCCCCAGGCTGCCAACTGGc                                                                                            >  5:50853/1‑149 (MQ=255)
                                                      cGTTCTATCACGGCTAGTCACCGAGGGCCACTCCTTCCCGCCGGGGATCGGCACCGCCTGACCAACCCTCGCCGCTGCGCTGGATGATCTGCATGCCGCTGGTCATGGTCATCGGTGTCACCTCATGGCCCCAGGCTGCCAACTGGCgg                                                                                          <  8:271311/149‑1 (MQ=255)
                                                        ttCTATCACGGCTAGTCACCGAGGGCCACTCCTTCCCGCCGGGGATCGGCACCGCCTGACCAACCCTCGCCGCTGCGCTGGATGATCTGCATGCCGCTGGTCATGGTCATCGGTGTCACCTCATGGCCCCAGGCTGCCAACTGGCGGAt                                                                                        <  8:130685/149‑1 (MQ=255)
                                                            aTCACGGCTAGTCACCGAGGGCCACTCCTTCCCGCCGGGGATCGGCACCGCCTGACCAACCCTCGCCGCTGCGCTGGATGATCTGCATGCCGCTGGTCATGGTCATCGGTGTCACCTCATGGCCCCAGGCTGCCAACTGGCGGATCAgc                                                                                    >  1:150960/1‑149 (MQ=255)
                                                             tCACGGCTAGTCACCGAGGGCCACTCCTTCCCGCCGGGGATCGGCACCGCCTGACCAACCCTCGCCGCTGCGCTGGATGATCTGCATGCCGCTGGTCATGGTCATCGGTGTCACCTCATGGCCCCAGGCTGc                                                                                                    <  1:151756/132‑1 (MQ=255)
                                                             tCACGGCTAGTCACCGAGGGCCACTCCTTCCCGCCGGGGATCGGCACCGCCTGACCAACCCTCGCCGCTGCGCTGGATGATCTGCATGCCGCTGGTCATGGTCATCGGTGTCACCTCATGGCCCCAGGCTGCCAACTGGCGGATCa                                                                                      <  4:187218/146‑1 (MQ=255)
                                                             tCACGGCTAGTCACCGAGGGCCACTCCTTCCCGCAGGGGATCGGCACCGCCTGACCAACCCTCGCCGCTGCGCTGGATGATCTGCATGCCGCTGGTCATGGTCATCGGTGTCACCTCATGGCCCCAGGCTGc                                                                                                    >  2:151756/1‑132 (MQ=255)
                                                                 ggCTAGTCACCGAGGGCCACTCCTTCCCGCCGGGGATCGGCACCGCCTGACCAACCCTCGCCGCTGCGCTGGATGATCTGCATGCCGCTGGTCATGGTCATCGGTGTCACCTCATGGCCCCAGGCTGCCAACTGGCGGATCAGCGCCgg                                                                               <  3:147129/149‑1 (MQ=255)
                                                                   cTAGTCACCGAGGGCCACTCCTTCCCGCCGGGGATCGGCACCGCCTGACCAACCCTCGCCGCTGCGCTGGATGATCTGCATGCCGCTGGTCATGGTCATCGGTGTCACCTCATGGCCCCAGGCTGCCAACTGGCGGATCAGCGCCGGgc                                                                             <  2:150960/149‑1 (MQ=255)
                                                                          ccGAGGGCCACTCCTTCCCGCCGGGGATCGGCACCGCCTGACCAACCCTCGCCGCTGCGCTGGATGATCTGCATGCCGCTGGTCATGGTCATCGGTGTCACCTCATGGCCCCAGGCTGCCAACTGGCGGATCAGCGCCGGGCTGCTCAg                                                                      >  3:189171/1‑149 (MQ=255)
                                                                             aGGGCCACTCCTTCCCGCCGGGGATCGGCACCGCCTGACCAACCCTCGCCGCTGCGCTGGATGATCTGCATGCCGCTGGTCATGGTCATCGGTGTCACCTCATGGCCCCAGGCTGCCAACTGGCg                                                                                           <  4:156767/125‑1 (MQ=255)
                                                                               ggCCACTCCTTCCCGCCGGGGATCGGCACCGCCTGACCAACCCTCGCCGCTGCGCTGGATGATCTGCATGCCGCTGGTCATGGTCATCGGTGTCACCTCATGGCCCCAGGCTGCCAACTGGCGGATCAGCGCCGGGCTGCTCAGCCCTg                                                                 >  4:122272/1‑149 (MQ=255)
                                                                                   aCTCCTTCCCGCCGGGGATCGGCACCGCCTGACCAACCCTCGCCGCTGCGCTGGATGATCTGCATGCCGCTGGTCATGGTCATCGGTGTCACCTCATGGCCCCAGGCTGCCAACTGGCGGATCAGCGCCGGGCTGCTCAGCCCTGCTTc                                                             <  7:132008/149‑1 (MQ=255)
                                                                                       cTTCCCGCCGGGGATCGGCACCGCCTGACCAACCCTCGCCGCTGCGCTGGATGATCTGCATGCCGCTGGTCATGGTCATCGGTGTCAc                                                                                                                      >  2:453/1‑88 (MQ=255)
                                                                                       cTTCCCGCCGGGGATCGGCACCGCCTGACCAACCCTCGCCGCTGCGCTGGATGATCTGCATGCCGCTGGTCATGGTCATCGGTGTCAc                                                                                                                      <  1:453/88‑1 (MQ=255)
                                                                                          cccGCCGGGGATCGGCACCGCCTGACCAACCCTCGCCGCTGCGCTGGATGATCTGCATGCCGCTGGTCATGGTCATCGGTGTCACCTCATGGCCCCAGGCTGCCAACTGGCGGATCAGCGCCGGGCTGCTCAGCCCTGCTTCCACTTc                                                       <  1:241022/148‑1 (MQ=255)
                                                                                            cgccgGGGATCGGCACCGCCTGACCAACCCTCGCCGCTGCGCTGGATGATCTGCATGCCGCTGGTCATGGTCATCGGTGTCACCTCATGGCCCCAGGCTGCCAACTGGCGGATCAGAGCCGGGCTCCTCAGCCCTGCTTCAACTTCGGt                                                    >  4:268931/1‑149 (MQ=255)
                                                                                               cgGGGATCGGCACCGCCTGACCAACCCTCGCCGCTGCGCTGGATGATCTGCATGCCGCTGGTCATGGTCATCGGTGTCACCTCATGGCCCCAGGCTGCCAACTGGCGGATCAGCGCCGGGCTGCTCAGCCCTGCt                                                               <  7:2968/135‑1 (MQ=255)
                                                                                               cgGGGATCGGCACCGCCTGACCAACCCTCGCCGCTGCGCTGGATGATCTGCATGCCGCTGGTCATGGTCATCGGTGTCACCTCATGGCCCCAGGCTGCCAACTGGCGGATCAGCGCCGGGCTGCTCAGCCCTGCt                                                               >  8:2968/1‑135 (MQ=255)
                                                                                               cgGGGATCGGCACCGCCTGACCAACCCTCGCCGCTGCGCTGGATGATCTGCATGCCGCTGGTCATGGTCATCGGTGTCACCTCATGGCCCCAGGCTGCCAACTGGCGGATCAGCGCCGGGCTGCTCAGCCCTGCTTCCACTTCGGTGcc                                                 <  4:218460/149‑1 (MQ=255)
                                                                                                     tCGGCACCGCCTGACCAACCCTCGCCGCTGCGCTGGATGATCTGCATGCCGCTGGTCATGGTCATCGGTGTCACCTCATGGCCCCAGGCTGCCAACTGGCGGATCAGCGCCGGGCTGCTCAGCCCTGCTTCCACTTCGGTGCCGGCAt                                            <  8:23438/148‑1 (MQ=255)
                                                                                                        gCACCGCCTGACCAACCCTCGCCGCTGCGCTGGATGATCTGCATGCCGCTGGTCTTGGTCATCGGTGTCACCTCATGGCCCCAGGCTGCCAACTGGCGGATCAGCGCCGGGCTGCTCAGCCCTGCTTCCACTTCGGTGCCGGCATTGCg                                        >  2:30641/1‑149 (MQ=255)
                                                                                                        gCACCGCCTGACCAACCCTCGCCGCTGCGCTGGATGATCTGCATGCCGCTGGTCATGGTCATCGGTGTCACCTCATGGCCCCAGGCTGCCAACTGGCGGATCAGCGCCGGGCTGCTCAGCCCTGCTTCCACTTCGGTGCCGGCATTGCg                                        >  2:73520/1‑149 (MQ=255)
                                                                                                        gCAACGCCTGACCAACCCTCGCCGCTGCGCTGGATGATCTGCATGCCGCTGGTCATGGTCATCGGTGTCACCTCATGGCCCCAGGCTGCCAACTGGCGGATCAGCGCCGGGCTGCTCAGCCCTGCTTCCACTTCGGTGCCGGCATTGCg                                        >  6:216542/1‑149 (MQ=255)
                                                                                                              ccTGACCAACCCTCGCCGCTGCGCTGGATGATCTGCATGCCGCTGGTCATGGTCATCGGTGTCACCTCATGGCCCCAGGCTGCCAACTGGCGGATCAGCGCCGGGCTGCTCAGCCCTGCTTCCACTTCGGTGCCGGCATTGCGCCTGcc                                  >  1:119133/1‑149 (MQ=255)
                                                                                                                 gACCAACCCTCGCCGCTGCGCTGGATGATCTACATGCCGCTGGTCATGGTCATCGGTGTCACCTCATGGCCCCAGGCTGCCAACTGGCGGATCAGCGCCGGGCTGCTCAGCCCTGCTTCCACTTCGGTGCCGGCATTGCGGCTGCCg                                 >  1:114074/1‑147 (MQ=255)
                                                                                                                  aCCAACCCTCGCCGCTGCGCTGGATGATCTGCATGCCGCTGGTCATGGTCATCGGTGTCACCTCATGGCCCCAGGCTGCCAACTGGCGGAt                                                                                        >  7:142738/1‑91 (MQ=255)
                                                                                                                  aCCAACCCTCGCCGCTGCGCTGGATGATCTGCATGCCGCTGGTCATGGTCATCGGTGTCACCTCATGGCCCCAGGCTGCCAACTGGCGGAt                                                                                        <  8:142738/91‑1 (MQ=255)
                                                                                                                  aCCAACCCTCGCCGCTGCGCTGGATGATCTGCATGCCGCTGGTCATGGTCATCGGTGTCACCTCATGGCCCCAGGCTGCCAACTGGCGGATCAGCGCCGGGCTGCTCAGCCCTGCTTCCACTTc                                                       <  8:256414/124‑1 (MQ=255)
                                                                                                                  aCCAACCCTCGCCGCTGCGCTGGATGATCTGCATGCCGCTGGTCATGGTCATCGGTGTCACCTCATGGCCCCAGGCTGCCAACTGGCGGATCAGCGCCGGGCTGCTCAGCCCTGCTTCCACTTc                                                       >  7:256414/1‑124 (MQ=255)
                                                                                                                  aCCAACCCTCGCCGCTGCGCTGGATGATCTGCATGCCGCTGGTCATGGTCATCGGTGTCACCTCATGGCCCCAGGCTGCCAACTGGCGGATCAGCGCCGGGCTGCTCAGCCCTGCTTCCACTTCGGTGCCGGCATTGCGGCTGCCGAAg                              >  7:166481/1‑149 (MQ=255)
                                                                                                                     aaCCCTCGCCGCTGCGCTGGATGATCTGCATGCCGCTGGTCATGGTCATCGGTGTCACCTCATGGCCCCAGGCTGCCAACTGGCGGATCAGCGCCGGGCTGCTCAGCCCTGCTTCCACTTCGGTGCCGGCATTGCGGCTGCCGAAGTTg                           >  1:211816/1‑149 (MQ=255)
                                                                                                                     aaCCCTCGCCGCTGCGCTGGATGATCTGCATGCCGCTGGTCATGGTCATCGGTGTCACCTCATGGCCCCAGGCTGCCAACTGGCGGATCAGCGCCGGGCTGCTCAGCCCTGCTTCCACTTCGGTGCCGGCATTGCGGCTGCCGAAGTTg                           <  4:250148/149‑1 (MQ=255)
                                                                                                                     aaCACTCGCCGCTTCGCTGGATGATCTGCTTGCCGCTGGTCATGGTCATCTGGGTCACCTCATGGCCCCAGGCTGCCAACTGGCGGATCAGCGCCGGGCTGCTCAGCCCTGCTTCCACTTCGGTGCCGGCATTGCGGCTGCCGAAGTTg                           <  4:189171/149‑1 (MQ=255)
                                                                                                                       cccTCGCCGCTGCGCTGGATGATCTACATGCCGCTGGTCATGGTCATCGGTGTCACCTCATGGCCCCAGGCTGCCAACTGGCGGATCAGCGCCGGGCTGCTCAGCCCTGCTTCCACTTCGGTGCCGGCATTGCGGCTGCCGAAGTTggg                         <  2:114074/149‑1 (MQ=255)
                                                                                                                         cTCGCCGCTGCGCTGGATGATCTGCATGCCGCTGGTCATGGTCATCGGTGTCACCTCATGGCCCCAGGCTGCCAACTGGCGGATCAGCGCCGGGCTGCTCAGCCCTGCTTcc                                                            >  7:248007/1‑112 (MQ=255)
                                                                                                                         cTCGCCGCTGCGCTGGATGATCTGCATGCCGCTGGTCATGGTCATCGGTGTCACCTCATGGCCCCAGGCTGCCAACTGGCGGATCAGCGCCGGGCTGCTCAGCCCTGCTTcc                                                            <  8:248007/112‑1 (MQ=255)
                                                                                                                          tCGCCGCTGCGCTGGATGATCTACATGCCGCTGGTCATGGTCATCGGTGTCACCTCATGGCCCCAGGCTGCCAACTGGCGGATCAGCGCCGGGCTGCTCAGCCCTGCTTCCACTTCGGTGCCGGCATTGCGGCTGCCGAAGt                             >  8:20674/1‑142 (MQ=255)
                                                                                                                          tCGCCGCTGCGCTGGATGATCTACATGCCGCTGGTCATGGTCATCGGTGTCACCTCATGGCCCCAGGCTGCCAACTGGCGGATCAGCGCCGGGCTGCTCAGCCCTGCTTCCACTTCGGTGCCGGCATTGCGGCTGCCGAAGt                             <  7:20674/142‑1 (MQ=255)
                                                                                                                            gccgcTGCGCTGGATGATCTGCATGCCGCTGGTCATGGTCATCGGTGTCACCTCATGGCCCCAGGCTGCCAACTGGCGGATCAGCGCCGGGCTGCTCAGCCCTGCTTCCACTTCGGTGCCGGCATTGCGGCTGCCGAAGTTGGGCAGGt                    <  3:122272/149‑1 (MQ=255)
                                                                                                                                    gctgGATGATCTGCATGCCGCTGGTCATGGTCATCGGTGTCACCTCATGGCCCCAGGCTGCCAACTGGCGGATCAGCGCCGGGCTGCTCAGCCCTGCTTCCACTTCGGTGCCGGCATTGCGGCTGCCGAAGTTGGGCAGGTCGGCGGcc            >  4:193152/1‑149 (MQ=255)
                                                                                                                                     ctgGATGATCTGCATGCCGCTGGTCATGGTCATCGGTGTCACCTCATGGCCCCAGGCTGCCAACTGGCGGATCAGCGCCGGGCTGCTCAGCCCTGCTTCCACTTCGGTGCCGGCATTGCGGCTGCCGAAGTTGGGCAGGTCGGCGGCct           >  2:69643/1‑149 (MQ=255)
                                                                                                                                        gatgatCTGCATGCCGCTGGTCATGGTCATCGGTGTCACCTCATGGCCCCAGGCTGCCAACTGGCGGATCAGCGCCGGGCTGCTCAGCCCTGCTTCCACTTCGGTGCCGGCATTGCGGCTGCCGAAGTTGGGCAGGTCGGCGGCctg          <  1:69643/147‑1 (MQ=255)
                                                                                                                                           gatCTGCATGCCGCTGGTCATGGTCATCGGTGTCACCTCATGGCCCCAGGCTGCCAACTGGCGGATCAGCGCCGGGCTGCTCAGCCCTGCTTCCACTTCGGTGCCGGCATTGCGGCTGCCGAAGTTGGGCAGGTCGGCGGcc            >  3:93396/1‑142 (MQ=255)
                                                                                                                                           gatCTGCATGCCGCTGGTCATGGTCATCGGTGTCACCTCATGGCCCCAGGCTGCCAACTGGCGGATCAGCGCCGGGCTGCTCAGCCCTGCTTCCACTTCGGTGCCGGCATTGCGGCTGCCGAAGTTGGGCAGGTCGGCGGcc            <  4:93396/142‑1 (MQ=255)
                                                                                                                                              cTGCATGCCGCTGGTCATGGTCATCGGTGTCACCTCATGGCCCCAGGCTGCCAACTGGCGGATCAGCGCCGGGCTGCTCAGCCCTGCTTCCACTTCGGTGCCGGCATTGCGGCTGCCGAAGtt                            >  4:215357/1‑123 (MQ=255)
                                                                                                                                              cTGCATGCCGCTGGTCATGGTCATCGGTGTCACCTCATGGCCCCAGGCTGCCAACTGGCGGATCAGCGCCGGGCTGCTCAGCCCTGCTTCCACTTCGGTGCCGGCATTGCGGCTGCCGAAGTTGGGCAGGTCGGCGGCCTGCTGCGGGt  >  5:64792/1‑149 (MQ=255)
                                                                                                                                              cTGCATGCCGCTGGTCATGGTCATCGGTGTCACCTCATGGCCCCAGGCTGCCAACTGGCGGATCAGCGCCGGGCTGCTCAGCCCTGCTTCCACTTCGGTGCCGGCATTGCGGCTGCCGAAGTTGGGCAGGTCGGCGGCCTGCTGCGGGt  >  5:43270/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
TTCGATTTCACGCTGATAGACATCTAGATATTCAACTACTGAATGGCAAGTCTGCGTTCTATCACGGCTAGTCACCGAGGGCCACTCCTTCCCGCCGGGGATCGGCACCGCCTGACCAACCCTCGCCGCTGCGCTGGATGATCTGCATGCCGCTGGTCATGGTCATCGGTGTCACCTCATGGCCCCAGGCTGCCAACTGGCGGATCAGCGCCGGGCTGCTCAGCCCTGCTTCCACTTCGGTGCCGGCATTGCGGCTGCCGAAGTTGGGCAGGTCGGCGGCCTGCTGCGGGT  >  NC_002947/4006018‑4006308

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