Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,248,398 C→T 7.7% R558H (CGC→CAC)  PP_3724 ← acyl‑CoA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,248,3980CT 7.7% 183.2 / 5.1 65R558H (CGC→CAC) PP_3724acyl‑CoA synthetase
Reads supporting (aligned to +/- strand):  ref base C (32/28);  new base T (3/2);  total (35/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.33e-01

CGCCGGCACCGATCACCAAGGCCAGGGCAAAGTTGTTTTCACGCTTGATGCCCACTATGAGCTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGC  >  NC_002947/4248255‑4248542
                                                                                                                                               |                                                                                                                                                
cgccgGCACCGATCACCAAGGCCAGGGCAAAGTTGTTTTCACGCTTGATGCCCACTATGAGCTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCt                                                                                                                                            >  1:238457/1‑150 (MQ=255)
  ccgGCACCGATCACCAAGGCCAGGGCAAAGTTGTTTTCACGCTTGATGCCCACTATGAGCTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGc                                                                                                                                          <  2:971564/150‑1 (MQ=255)
    ggCACCGATCACCAAGGCCAGGGCAAAGTTGTTTTCACGCTTGATGCCCACTATGAGCTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCAc                                                                                                                                        <  2:565512/150‑1 (MQ=255)
              aCCAAGGCCAGGGCAAAGTTGTTTTCACGCTTGATGCCCACTATGAGCTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCt                                                                                                                              >  2:227974/1‑150 (MQ=255)
              aCCAAGGCCAGGGCAAAGTTGTTTTCACGCTTGATGCCCACTATGAGCTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCt                                                                                                                              >  1:1132054/1‑150 (MQ=255)
                cAAGGCCAGGGCAAAGTTGTTTTCACGCTTGATGCCCACTATGAGCTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTc                                                                                                                            >  2:1074095/1‑150 (MQ=255)
                      cAGGGCAAAGTTGTTTTCACGCTTGATGCCCACTATGAGCTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGgc                                                                                                                      >  1:1021689/1‑150 (MQ=255)
                          gCAAAGTTGTTTTCACGCTTGATGCCCACTATGAGCTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCgg                                                                                                                  <  2:834677/150‑1 (MQ=255)
                            aaaGTTGTTTTCACGCTTGATGCCCACTATGAGCTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGTGCAGATGCGCCAGTGCatc                                                                                                                               >  2:609674/1‑133 (MQ=255)
                            aaaGTTGTTTTCACGCTTGATGCCCACTATGAGCTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGTGCAGATGCGCCAGTGCatc                                                                                                                               <  1:609674/135‑3 (MQ=255)
                            aaaGTTGTTTTCACGCTTGATGCCCACTATGAGCTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGc                                                                                                                >  1:753503/1‑150 (MQ=255)
                                    tttCACGCTTGATGCCCACTATGAGCTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCa                                                                                                        <  1:1074095/150‑1 (MQ=255)
                                         cGCTTGATGCCCACTATGAGCTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGc                                                                                                                                  >  1:1168977/1‑119 (MQ=255)
                                         cGCTTGATGCCCACTATGAGCTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGc                                                                                                                                  <  2:1168977/119‑1 (MQ=255)
                                           cTTGATGCCCACTATGAGCTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCaa                                                                                                 <  2:548826/150‑1 (MQ=255)
                                               aTGCCCACTATGAGCTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTg                                                                                                                                   >  1:1253010/1‑112 (MQ=255)
                                               aTGCCCACTATGAGCTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTg                                                                                                                                   <  2:1253010/112‑1 (MQ=255)
                                                tGCCCACTATGAGCTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCca                                                                                            <  1:159904/150‑1 (MQ=255)
                                                tGCCCACTATGAGCTCGGCCAAGGGCAGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCACCCTGCGGGGCGGAGGCAGCGATGCACTGGCGAAACTGCACCAGTGCAGGTACATGGGCGCGGTCATCGCGAAGGTTGAACGcaa                                                                                            >  2:422025/1‑148 (MQ=255)
                                                  cccACTATGAGCTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCcaca                                                                                          <  1:227974/150‑1 (MQ=255)
                                                   ccACTATGAGCTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAg                                                                                                                         <  1:1267794/118‑1 (MQ=255)
                                                   ccACTATGAGCTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAg                                                                                                                         >  2:1267794/1‑118 (MQ=255)
                                                   ccACTATGAGCTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAg                                                                                         >  2:287984/1‑150 (MQ=255)
                                                       tATGAGCTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCAcc                                                                                     >  2:1142988/1‑150 (MQ=255)
                                                            gCTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAg                                                                                                       <  2:617206/127‑1 (MQ=255)
                                                            gCTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAg                                                                                                       >  1:617206/1‑127 (MQ=255)
                                                             cTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTc                                                                               <  1:290793/150‑1 (MQ=255)
                                                             cTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTc                                                                               <  1:1142988/150‑1 (MQ=255)
                                                              tCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCAc                                                                                                                                        <  1:187197/92‑1 (MQ=255)
                                                              tCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCAc                                                                                                                                        >  2:187197/1‑92 (MQ=255)
                                                              tCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCg                                                                              <  1:166491/150‑1 (MQ=255)
                                                                 gCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTATCGCATCTGCACCACCGCAGCTTTTAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTc                                                                           <  1:158249/150‑1 (MQ=255)
                                                                 gCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAAAATGAAAGGAGTAGTCAACGATGAACTAGCGCATCTTCACGATGTCAGCTTCGAGTGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTc                                                                           >  2:1109601/1‑150 (MQ=255)
                                                                        gcggcTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCg                                                                    >  1:236073/1‑150 (MQ=255)
                                                                                  cATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTgcgc                                                          <  1:122048/150‑1 (MQ=255)
                                                                                     gggCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGTGCAGATGCGCCAGTGc                                                                                                                                  <  1:344916/75‑1 (MQ=37)
                                                                                     gggCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGTGCAGATGCGCCAGTGc                                                                                                                                  >  2:344916/1‑75 (MQ=37)
                                                                                               agcaCCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAg                                                                                                       >  1:439590/1‑92 (MQ=255)
                                                                                               agcaCCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAg                                                                                                       <  2:439590/92‑1 (MQ=255)
                                                                                               agcaCCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTc                                             <  2:236073/150‑1 (MQ=255)
                                                                                                 caCCTGATCAAACGCAACCTGCGGGGCGTAGCCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAg                                           >  1:651524/1‑150 (MQ=255)
                                                                                                    cTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTAc                                        >  2:1060801/1‑150 (MQ=255)
                                                                                                    cTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTAc                                        >  1:975552/1‑150 (MQ=255)
                                                                                                       ataaaaCGCAATCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCaa                                     >  1:1283915/4‑150 (MQ=255)
                                                                                                        tCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAg                                    <  1:1148196/150‑1 (MQ=255)
                                                                                                          aaaCGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAggg                                  >  2:421452/1‑150 (MQ=255)
                                                                                                          aaaCGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAggg                                  <  1:311852/150‑1 (MQ=255)
                                                                                                          aaaCGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAggg                                  >  2:948775/1‑150 (MQ=255)
                                                                                                                   cTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCtt                                              <  2:20836/129‑1 (MQ=255)
                                                                                                                   cTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCtt                                              >  1:20836/1‑129 (MQ=255)
                                                                                                                    tGCGGGGCGTAGTCAGCGATGCCCTGGCTCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGAcg                        >  1:494944/1‑150 (MQ=255)
                                                                                                                    tGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGAcg                        >  1:280054/1‑150 (MQ=255)
                                                                                                                     gCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGAcgc                       <  1:487469/150‑1 (MQ=255)
                                                                                                                             tAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGgcggcg               >  2:99714/1‑150 (MQ=255)
                                                                                                                              aGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGt              >  2:405987/1‑150 (MQ=255)
                                                                                                                                tCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCt            >  2:1268387/1‑150 (MQ=255)
                                                                                                                                tCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCt            <  1:287984/150‑1 (MQ=255)
                                                                                                                                tCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCt            <  1:1060801/150‑1 (MQ=255)
                                                                                                                                    cGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGAGCGCCTTCAGTACGAAGGGGTAACTGACGCGGGCGGCGTCTGCCa        >  2:976601/1‑150 (MQ=255)
                                                                                                                                       tGCACTGGTGCAGATGCGCCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAACGTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACgg     >  1:874172/1‑150 (MQ=255)
                                                                                                                                       tGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACgg     <  2:238457/150‑1 (MQ=255)
                                                                                                                                          aCTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGgcgc  >  2:1095887/1‑150 (MQ=255)
                                                                                                                                          aCTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGgcgc  <  1:1268387/150‑1 (MQ=255)
                                                                                                                                           cTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAggg                                  <  2:297603/117‑1 (MQ=255)
                                                                                                                                           cTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAggg                                  >  1:297603/1‑117 (MQ=255)
                                                                                                                                               |                                                                                                                                                
CGCCGGCACCGATCACCAAGGCCAGGGCAAAGTTGTTTTCACGCTTGATGCCCACTATGAGCTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGC  >  NC_002947/4248255‑4248542

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑