Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,286,145 C→G 6.7% M403I (ATG→ATC)  ggt ← gamma‑glutamyltranspeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,286,1450CG 6.7% 180.2 / 4.6 60M403I (ATG→ATC) ggtgamma‑glutamyltranspeptidase
Reads supporting (aligned to +/- strand):  ref base C (28/28);  new base G (2/2);  total (30/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 = 6.69e-01

CGATTACCAGCTCGACCTTGCCATCGCGGGTCATCAGGCTGGGGCTCATCGACGACAGCATGCGCTTGCCCGGTTCGATCGCGTTGGCATCACCGCCGACCACGCCAAAGGCGTTCGCCGCCCCCGGCTTGGCGCTGAAGTCATCCATCTCGTCGTTGAGCAGGAAACCGGCACC‑CTTCACCACGACGCCACTGCCGTAGTCGAGGTTGAGGGTGTAGGTATTGCTGACCGCGTTGCCCTGCTTGTCGACGATCGAGAAATGCGTGGTCTGGTGCGGTTCAAGGCCCGGCTTGA  >  NC_002947/5286000‑5286293
                                                                                                                                                 |                                                                                                                                                     
cGATTACCAGCTCGACCTTGCCATCGCGGGTCATCAGGCTGGGGCTCATCGACGACAGCATGCGCTTGCCCGGTTCGATCGCGTTGGCATCACCGCCGACCACGCCAAAGGCGTTCGCCGCCCCCGGCTTGGCGCTGAAGTCATCCAtct                                                                                                                                                   >  1:681336/1‑150 (MQ=255)
    tACCAGCTCGACCTTGCCATCGCGGGTCATCAGGCTGGGGCTCATCGACGACAGCATGCGCTTGCCCGGTTCGATCGCGTTGGCATCACCGCCGACCACGCCAAAGGCGTTTGCCGCCCCCGGCTTGGCGCTGAAGTCATCCATCtcgtc                                                                                                                                               <  2:280551/150‑1 (MQ=255)
         gCTCGACCTTGCCATCGCGGGTCATCAGGCTGGGGCTCATCGACGACAGCATGCGCTTGCCCGGTTCGATCGCGTTGGCATCACCGCCGACCACGCCAAGGGCGTTCGCCGCCCTCGGCTTGGCGCTGAAGTCATCCATCTCGTCGTTGa                                                                                                                                          <  2:1099980/150‑1 (MQ=255)
            cGACCTTGCCATCGCGGGTCACCAGACTGGGGTTCATCGACGACAGCATGCGCTTGCCCGGTTCGTTCGCGTTGGCATCACCGCCGACCACGCCAAAGGCGTTCGCCGCCCCCGGCTTGGCGCTGGAGTCATCCATCTCGTCGTTGAGCa                                                                                                                                       <  2:642588/150‑1 (MQ=255)
                   gCCATCGCGGGTCATCAGGCTGGGGCTCATCGACGACAGCATGCGCTTGCCCGGTTCGATCGCGTTGGCATCACCGCCGACCACGCCAAAGGCGTTCGCCGCCCCCGGCTTGGCGCTGAAGTCATCCATCTCGTCGTTGAGCAGGAAAcc                                                                                                                                <  1:1061308/150‑1 (MQ=255)
                      aTCGCGGGTCATCAGGCTGGGGCTCATCGACGACAGCATGCGCTTGCCCGGTTCGATCGCGTTGGCATCACCGCCGACCACGCCAAAGGCGTTCGCCGCCCCCGGCTTGGCGCTGAAGTCATCCATCTCGTCGTTGAGCAGGAAACCGGc                                                                                                                             <  2:974708/150‑1 (MQ=255)
                       tCGCGGGTCATCAGGCTGGGGCTCATCGACGACAGCATGCGCTTGCCCGGTTCGATCGCGTTGGCATCACCGCCGACCACGCCAAAGGCGTTCGCCGCCCCCGGCTTGGCGCTGAAGTcatccatc                                                                                                                                                    >  2:726069/1‑126 (MQ=255)
                       tCGCGGGTCATCAGGCTGGGGCTCATCGACGACAGCATGCGCTTGCCCGGTTCGATCGCGTTGGCATCACCGCCGACCACGCCAAAGGCGTTCGCCGCCCCCGGCTTGGCGCTGAAGTcatccatc                                                                                                                                                    <  1:726069/126‑1 (MQ=255)
                                atcaGGCTGGGGCTCATCGACGACAGCATGCGCTTGCCCGGTTCGATCGCGTTGGCATCACCGCCGACCACGCCAAAGGCGTTCGCCGCCCCCGGCTTGGCGCTGAAGTCATCCATCTCGTCGTTGAGCAGGAAACCGGCACC‑CTTcacc                                                                                                                  >  1:67719/1‑150 (MQ=255)
                                 tcaGGCTGGGGCTCATCGACGACAGCATGCGCTTGCCCGGTTCGATCGCGTTGGCATCACCGCCGACCACGCCAAAGGCGTTCGCCGCCCCCGGCTTGGCGCTGAAGTCATCCATCTCGTCGTTGAGCAGGAAACCGGCACC‑CTTcacca                                                                                                                 >  2:751857/1‑150 (MQ=255)
                                           gCTCATCGACGACAGCATGCGCTTGCCCGGTTCGATCGCGTTGGCATCACCGCCGACCACGCCAAAGGCGTTCGCCGCCCCCGGCTTGGCGCTGAAGTCATCCATCTCGTCGTTGAGCAGGAAACCGGCACC‑CTTCACCACGACGCCACt                                                                                                       >  2:293505/1‑150 (MQ=255)
                                                     gacAGCATGCGCTTGGCCGGTACGATCGCGTTGGCATCACCGCCGACCACGCCAAAGGCGTTCGCCGCCCCCGGCTTGGCGCTGAAGTCATCCATCTCGTCGTTGAGCAGGAAACCGGCa                                                                                                                            <  2:549325/120‑1 (MQ=255)
                                                     gacAGCATGCGCTTGGCCGGTACGATCGCGTTGGCATCACCGCCGACCACGCCAAAGGCGTTCGCCGCCCCCGGCTTGGCGCTGAAGTCATCCATCTCGTCGTTGAGCAGGAAACCGGCa                                                                                                                            >  1:549325/1‑120 (MQ=255)
                                                     gacAGCATGCGCTTGCCCGGTTCGATCGCGTTGGCATCACCGCCGACCACGCCAAAGGCGTTCGCCGCCCCCGGCTTGGCGCTGAAGTCATCCATCTCGTCGTTGAGCAGGAAACCGGCACC‑CTTCACCACGACGCCACTGCCGTAGTCg                                                                                             >  2:348242/1‑150 (MQ=255)
                                                     gacAGCATGCGCTTGCCCGGTTCGATCGCGTTGGCATCACCGCCGACCACGCCAAAGGCGTTCGCCGCCCCCGGCTTGGCGCTGAAGTCATCCATCTCGTCGTTGAGCAGGAAACCGGCACC‑CTTCACCACGACGCCACTGCCGTAGTCg                                                                                             >  1:344453/1‑150 (MQ=255)
                                                      acAGCATGCGCTTGCCCGGTTCGATCGCGTTGGCATCACCGCCGACCACGCCAAAGGCGTTCGCCGCCCCCGGCTTGGCGCTGAAGTCATCCATCTCGTCGTTGAGCAGGAAACCGGCACC‑CTTCACCACGACGCCACTGCCGTAGTCGa                                                                                            <  1:293505/150‑1 (MQ=255)
                                                      acAGCATGCGCTTGCCCGGTTCGATCGCGTTGGCATCACCGCCGACCACGCCAAAGGCGTTCGCCGCCCCCGGCTTGGCGCTGAAGTCATCCATCTCGTCGTTGAGCAGGAAACCGGCACC‑CTTCACCACGACGCCACTGCCGTAGTCGa                                                                                            <  2:344453/150‑1 (MQ=255)
                                                        tgCATGCGCTTGCCCGGTTCGATCGCGTTGGCATCACCGCCGACCACGCCAAAGGCGTTCGCCGCCCCCGGCTTGGCGCTGAAGTCATCCATCTCGTCGTTGAGCAGGAAACCGGCACC‑CTTCACCACGACGCCACTGCCGTAGTCGAgg                                                                                          <  1:108888/149‑1 (MQ=255)
                                                            tGCGCTTGCCCGGTTCGATCGCGTTGGCATCACCGCCGACCACGCCAAAGGCGTTCGCCGCCCCCGGCTTGGCGCTGAAGTCATCCATCTCGTCGTTGAGCAGGAAACCGGCACC‑CTTCACCACGACGCCACTGCCGTAGTCGAgg                                                                                          >  1:388734/1‑146 (MQ=255)
                                                            tGCGCTTGCCCGGTTCGATCGCGTTGGCATCACCGCCGACCACGCCAAAGGCGTTCGCCGCCCCCGGCTTGGCGCTGAAGTCATCCATCTCGTCGTTGAGCAGGAAACCGGCACC‑CTTCACCACGACGCCACTGCCGTAGTCGAgg                                                                                          <  2:388734/146‑1 (MQ=255)
                                                            tGCGCTTGCCCGGTTCGATCGCGTTGGCATCACCGCCGACCACGCCAAAGGCGTTCGCCGCCCCCGGCTTGGCGCTGAAGTCATCCATCTCGTCGTTGAGCAGGAAACCGGCACC‑CTTCACCACGACGCCACTGCCGTAGTCGAGGTTGa                                                                                      <  1:913637/150‑1 (MQ=255)
                                                                   gCCCGGTTCGATCGCGTTGGCATCACCGCCGACCACGCCAAAGGCGTTCGCCTCCCCCGGCTTGGCGCTGAAGTCATCCATCTCGTCGTTGAGCAGGAAACCGGCACC‑CTTCACCACGACGCCACTGCCGTAGTCGAGGTTGAGGGTGTa                                                                               <  1:1172029/150‑1 (MQ=255)
                                                                   gCCCGGTTCGATCGCGTTGGCATCACCGCCGACCACGCCAAAGGCGTTCGCCGCCCCCGGCTTGGCGCTGAAGTCATCCATCTCGTCGTTGAGCAGGAAACCGGCACC‑CTTCACCACGACGCCACTGCCGTAGTCGAGGTTGAGGGTGTa                                                                               <  2:452303/150‑1 (MQ=255)
                                                                   gCCCGGTTCGATCGCGTTGGCATCACCGCCGACCACGCCAAAGGCGTTCGCCGCCCCCGGCTTGGCGCTGAAGTCATCCATCTCGTCGTTGAGCAGGAAACCGGCACC‑CTTCACCACGACGCCACTGCCGTAGTCGAGGTTGAGCGTGTa                                                                               <  2:1281261/150‑1 (MQ=255)
                                                                    cccGGTTCGATCGCGTTGGCATCACCGCCGACCACGCCAAAGGCGTTCGCCGCCCCCGGCTTGGCGCTGAAGTCATCCATCTCGTCGTTGAGCAGGAAACCGGCACC‑CTTCACCACGACGCCACTGCCGTAGTCGAGGTTGAGGGTGTAg                                                                              <  1:961644/150‑1 (MQ=255)
                                                                       ggTTCGATCGCGTTGGCATCACCGCCGACCACGCCAAAGGCGTTCGCCGCCCCCGGCTTGGCGCTGAAGTCATCCATCTCGTCGTTGAGCAGGAAACCGGCACC‑CTTCACCACGACGCCACTGCCGTAGTCGAGGTTGAGGGTGTAGGTa                                                                           >  1:111210/1‑150 (MQ=255)
                                                                             aTCGCGTTGGCATCACCGCCGACCACGCCAAAGGCGTTCGCCGCCCCCGGCTTGGCGCTGAAGTCATCCATCTCGTCGTTGAGCAGGAAACCGGCACC‑CTTCACCACGACGCCACTGCCGTAGTCGAGGTTGAGGGTGTAGGTATTGCTg                                                                     >  2:1250051/1‑150 (MQ=255)
                                                                              tCGCGTTGGCATCACCGCCGACCACGCCAAAGGCGTTCGCCGCCTCCGGCTTGGCGCTGAAGTCATCCATCTCGTCGTTGAGCAGGACACCGGCACC‑CTTCACCACGACGCCACTGCCGTAGTCGAGGTTGAGGGTGTAGGTATTGCTGa                                                                    <  2:294330/150‑1 (MQ=255)
                                                                              tCGCGTTGGCATCACCGCCGACCACGCCAAAGGCGTTCGCCGCCCCCGGCTTGGCGCTGAAGTCATCCATCTCGTCGTTGAGCAGGAAACCGGCACC‑CTTCACCACGACGCCACTGCCGTAGTCGAGGTTGAGGGTGTAGGTATTGCTGa                                                                    >  1:417781/1‑150 (MQ=255)
                                                                                gcgTTGGCATCACCGCCGACCACGCCAAAGGCGTTCGCCGCCCCCGGCTTGGCGCTGAAGTCATCCATCTCGTCGTTGAGCAGGAAACCGGCACC‑CTTCACCACGACGCCACTGCCGTAGTCGAGGTTGAGGGTGTAGGTATTGCTGAcc                                                                  >  1:609688/1‑150 (MQ=255)
                                                                                gcgTTGGCATCACCGCCGACCACGCCAAAGGCGTTCGCCGCCCCCGGCTTGGCGCTGAAGTCATCCATCTCGTCGTTGAGCAGGAAACCGGCACC‑CTTCACCACGACGCCACTGCCGTAGTCGAGGTTGAGGGTGTAGGTATTGCTGAcc                                                                  >  2:937556/1‑150 (MQ=255)
                                                                                   ttGGCATCACCGCCGACCACGCCAAAGGCGTTCGCCGCCCCCGGCTTGGCGCTGAAGTCATCCATCTCGTCGTTGAGCAGGAAACCGGCACC‑CTTCACCACGACGCCACTGCCGTAGTCGAGGTTGAGGGTGTAGGTATTGCTGACcgcg                                                               >  1:578127/1‑150 (MQ=255)
                                                                                        aTCACCGCCGACCACGCCAAAGGCGTTCGCCGCCCCCGGCTTGGCGCTGAAGTCATCCATCTCGTCGTTGAGCAGGAAACCGGCACC‑CTTCACCACGACGCCACTGCCGTAGTCGAGGTTGAGGGTGTAGGTATTGCTGACCGCGTTGcc                                                          >  1:839719/1‑150 (MQ=255)
                                                                                            ccgccgACCACGCCAAAGGCGTTCGCCGCCCCCGGCTTGGCGCTGAAGTCATCCATCTCGTCGTTGAGCAGGAAACCGGCACC‑CTTCACCACGACGCCACTGCCGTAGTCGAGGTTGAGGGTGTAGGTATTGCTGACCGCGTTGCCCTGc                                                      >  2:1162175/1‑150 (MQ=255)
                                                                                                cgACCACGCCAAAGGCGTTCGCCGCCCCCGGCTTGGCGCTGAAGTCATCCATCTCGTCGTTGAGCAGGAAACCGGCACC‑CTTCCCCACGACGCCACTGCCGTAGTCGAGGTTGAGGGTGTAGGTATTGCTGACCGCGTTGCCCTGCTTGt                                                  <  2:681336/150‑1 (MQ=255)
                                                                                                 gACCACGCCAAAGGCGTTCGCCGCCCCCGGCTTGGCGCTGAAGTCATCCATCTCGTCGTTGAGCAGGAAACCGGCACC‑CTTCACCACGACGCCACTGCCGTAGTCGAGGTTGAGGGTGTAGGTATTGCTGACCGCGTTGCCCTGCTTGTc                                                 <  2:111210/150‑1 (MQ=255)
                                                                                                   ctaCGCCAAAGGCGTTCGCCGCCCCCGGCTTGGCGCTGAAGTCATCCATCTCGTCGTTGAGCAGGAAACCGGCACC‑CTTCACCACGACGCCACTGCCGTAGTCGAGGTTGAGGGTGTAGGTATTGCTGACCGCGTTGCCCTGCTTGTcga                                               <  1:1128371/148‑1 (MQ=255)
                                                                                                    ctcGCCAAAGGCGTTCGCCGCCCCCGGCTTGGCGCTGAAGTCATCCCTCTCGTCGTTTAGCAGGAAACCGGCACC‑CTTCACCACGACGCCACTGCCGTAGTCGAGGTTGAGGGTGTAGGGATTGCTGACCGCGTTGCCCTGCTTGTcgcc                                              >  2:577692/3‑148 (MQ=255)
                                                                                                    cACGCCAAAGGCGTTCGCCGCCCCCGGCTTGGCGCTGAAGTCATCCATCTCGTCGTTGAGCAGGAAACCGGCACC‑CTTCACCACGACGCCACTGCCGTAGTCGAGGTTGAGGGTGTAGGTATTGCTGACCGCGTTGCCCTGCTTGTcgac                                              >  1:134214/1‑150 (MQ=255)
                                                                                                       gCCAAAGGCGTTCGCCGCCCCCGGCTTGGCGCTGAAGTCATCCATCTCGTCGTTGAGCAGGAAACCGGCACC‑CTTCACCACGACGCCACTGCCGTAGTCGAGGTTGAGGGTGCAGGTATTGCTGACCGCGTTGCGGTGTTTGGctacgag                                           >  2:1037547/1‑149 (MQ=255)
                                                                                                         cAAAGGCGTTCGCCGCCCCCGGCTTGGCGCTGAAGTCATCCATCTCGTCGTTGAGCAGGAAACCGGCACC‑CTTCACCACGACGCCACTGCCGTAGTCGAGGTTGAGGGTGTAGGTATTGCTGACCGCGTTGCCCTg                                                       >  1:436455/1‑136 (MQ=255)
                                                                                                         cAAAGGCGTTCGCCGCCCCCGGCTTGGCGCTGAAGTCATCCATCTCGTCGTTGAGCAGGAAACCGGCACC‑CTTCACCACGACGCCACTGCCGTAGTCGAGGTTGAGGGTGTAGGTATTGCTGACCGCGTTGCCCTg                                                       <  2:436455/136‑1 (MQ=255)
                                                                                                         cAAAGGCGTTCGCCGCCCCCGGCTTGGCGCTGAAGTCATCCATCTCGTCGTTGAGCAGGAAACCGGCACC‑CTTCACCACGACGCCACTGCCGTAGTCGAGGTTGAGGGTGTAGGTATTGCTGACCGCGTTGCCCTGCTTGTCGACGATCg                                         >  1:1021448/1‑150 (MQ=255)
                                                                                                                gTTCGCCGCCCCCGGCTTGGCGCTGAAGTCATCCATCTCGTCGTTGAGCAGGAAACCGGCACC‑CTTCACCACGACGCCACTGCCGTAGTCGAGGTTGAGGGTGTAGGTATTGCTGACCGCGTTGCCCTGCTTGTCGACGATCGAGAAATg                                  >  1:271515/1‑150 (MQ=255)
                                                                                                                 ttCGCGGCCCCCGGCTTGGCGCTGAAGTCATCCATCTCGTCGTTGCCTGAGGAACCGGCACC‑CTTCACCACGACGCCACTGCCGTAGTCGAGGTTGAGGGTGTAGGTATTGCTGACCGCGTTGCCCTGCTTGTCGACGATCGTGAAATGc                                 <  1:78909/150‑1 (MQ=255)
                                                                                                                 ttCGCCGCCCCCGGCTTGGCGCTGAAGTCATCCATCTCGTCGTTGAGCAGGAAACCGGCACC‑CTTCACCACGACGCCACTGCCGTAGTCGAGGTTGAg                                                                                     >  2:707016/1‑98 (MQ=255)
                                                                                                                 ttCGCCGCCCCCGGCTTGGCGCTGAAGTCATCCATCTCGTCGTTGAGCAGGAAACCGGCACC‑CTTCACCACGACGCCACTGCCGTAGTCGAGGTTGAg                                                                                     <  1:707016/98‑1 (MQ=255)
                                                                                                                    gccgccCCCGGCTTGGCGCTGAAGTCATCCATCTCGTCGTTGAGCAGGAAACCGGCACC‑CTTCACCACGACGCCACTGCCGTAGTCGAGGTTGAGGGTGTAGGTATTGCTGACCGCGTTGCCCTGCTTGTCGACGATCGAGAAATGCGTg                              >  1:56458/1‑150 (MQ=255)
                                                                                                                           ccGGCTTGGCGCTGAAGTCATCCATCTCGTCGTTGAGCAGGAAACCGGCACC‑CTTCACCACGACGCCACTGCCGTAGTCGAGGTTGAGGGTGTAGGTATTGCTGACCGCGTTGCCCTGCTTGTCGACGATCGAGAAATGCGTGGTCTGGt                       <  2:609688/150‑1 (MQ=255)
                                                                                                                              gCTTGGCGCTGAAGTCATCCATCTCGTCGTTGAGCAGGAAACCGGCACC‑CTTCACCACGACGCCACTGCCGTAGTCGAGGTTGAGGGTGTAGGTATTGCTGACCGCGTTGCCCTGCTTGTCGACGATCGAGAAATGCGTGGTCTGGTGCg                    <  2:134214/150‑1 (MQ=255)
                                                                                                                                  ggCGCTGAAGTCATCCATCTCGTCGTTGAGCAGGAAACCGGCACC‑CTTCACCACGACGCCACTGCCGTAGTCGAg                                                                                           <  1:1008739/75‑1 (MQ=255)
                                                                                                                                  ggCGCTGAAGTCATCCATCTCGTCGTTGAGCAGGAAACCGGCACC‑CTTCACCACGACGCCACTGCCGTAGACGAg                                                                                           >  2:1008739/1‑75 (MQ=255)
                                                                                                                                   gcgcTGAAGTCATCCATCTCGTCGTTGAGCAGGAAACCGGCACC‑CTTCACCACGACGCCACTGCCGTAGTCGAGGTTGAGGGTGTAGGTATTGCTGACCGCGTTGCCCTGCTTGTCGACGATCGAGAAATGCGTGGTCTGGTGCGGTTCa               >  2:691354/1‑150 (MQ=255)
                                                                                                                                    cgcTGAAGTCATCCATCTCGTCGTTGAGCAGGAAACCGGCACC‑CTTCACCACGACGCCACTGCCGTAGTCGAGGTTGAGGGTGTAGGTATTGCTGACCGCGTTGCCCTGCTTGTCGACGTTCGAGAAATGCGTGGTCTGGTGCGGTTCaa              <  1:1037547/150‑1 (MQ=255)
                                                                                                                                      cTGAAGTCATCCATCTCGTCGTTGAGCAGGAAAACGGATCCTCTTC‑CGATCTCGCCACTTGCGTAGTCGAGGTTGAGGGTGTAGGTATTGCTGACCGCGTTGCCCTGCTTGTCGACGATCGAGAAATGCGTGGTCTGGTGCGGTTCAAgg            <  1:1250051/150‑1 (MQ=255)
                                                                                                                                       tGAAGTCATCCATCTCGTCGTTGAGCAGGAAACCGGCACC‑CTTCACCACGACGCCACTGCCGTAGTCGAGGTTGAGGGTGTAGGTATTGCTGACCGCGTTGCCCTGCTTGTCGACGATCGAGAAATGCGTGGTCTGGTGCGGTTCAAGGc           >  1:106108/1‑150 (MQ=255)
                                                                                                                                        gAAGTCATCGATCTCGTCGTTGAGCAGGAAACCGGCACC‑CTTCACCACGACGCCACTGCCGTAGTCGAGGTTGAGGGTGTAGGTATTGCTGACCGCGTTGCCCTGc                                                      >  1:1176568/1‑106 (MQ=255)
                                                                                                                                        gAAGTCATCGATCTCGTCGTTGAGCAGGAAACCGGCACC‑CTTCACCACGACGCCACTGCCGTAGTCGAGGTTGAGGGTGTAGGTATTGCTGACCGCGTTGCCCTGc                                                      <  2:1176308/106‑1 (MQ=255)
                                                                                                                                        gAAGTCATCGATCTCGTCGTTGAGCAGGAAACCGGCACC‑CTTCACCACGACGCCACTGCCGTAGTCGAGGTTGAGGGTGTAGGTATTGCTGACCGCGTTGCCCTGc                                                      >  1:1176308/1‑106 (MQ=255)
                                                                                                                                        gAAGTCATCGATCTCGTCGTTGAGCAGGAAACCGGCACC‑CTTCACCACGACGCCACTGCCGTAGTCGAGGTTGAGGGTGTAGGTATTGCTGACCGCGTTGCCCTGc                                                      <  2:1176568/106‑1 (MQ=255)
                                                                                                                                         aaGTCATCCATCTCGTCGTTGAGCAGGAAACCGGCACC‑CTTCACCACGACGCCACTGCCGTAGTCGAGGTTGAGGGTGTAGGTATTGCTGACCGCGTTGCCCTGCTTGTCGACGATCGAGAAAt                                   <  1:1181917/124‑1 (MQ=255)
                                                                                                                                         aaGTCATCCATCTCGTCGTTGAGCAGGAAACCGGCACC‑CTTCACCACGACGCCACTGCCGTAGTCGAGGTTGAGGGTGTAGGTATTGCTGACCGCGTTGCCCTGCTTGTCGACGATCGAGAAAt                                   >  2:1181917/1‑124 (MQ=255)
                                                                                                                                              atccatcTCGTCGTTGAGCAGGAAACCGGCACC‑CTTCACCACGACGCCACTGCCGTAGTCGAGTATGAGGGTGTAGGTATTTCTGACCGCGTTGCCCTGCTTGTCGACGATCAAGAAATGCGTGGTCTGGTGCGGTTCAAGGCCCGGCtt    >  2:549037/1‑150 (MQ=255)
                                                                                                                                                ccatcTCGTCGTTGAGCAGGAAACCGGCACC‑CTTCACCACGACGCCACTGCCGTAGTCGAGGTTGAGGGTGTAGGTATTGCTGACCGCGTTGCCCTGCTTGTCGACGATCGAGAAATGCGTGGTCTGGTGCGGTTCAAGGCCCGGCTTGa  >  2:704368/1‑150 (MQ=255)
                                                                                                                                                 catcTCGTCGTTGAGCAGGAAACCGGCACC‑CTTCACCACGACGCCACTGCCGTAGTCGAGGTTGAGGGTGTAGGTATTGCTGACCGCGTTGCCCTGCTTGTCGACGATCGAGaaa                                    <  1:864692/115‑1 (MQ=255)
                                                                                                                                                 catcTCGTCGTTGAGCAGGAAACCGGCACC‑CTTCACCACGACGCCACTGCCGTAGTCGAGGTTGAGGGTGTAGGTATTGCTGACCGCGTTGCCCTGCTTGTCGACGATCGAGaaa                                    >  2:864692/1‑115 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
CGATTACCAGCTCGACCTTGCCATCGCGGGTCATCAGGCTGGGGCTCATCGACGACAGCATGCGCTTGCCCGGTTCGATCGCGTTGGCATCACCGCCGACCACGCCAAAGGCGTTCGCCGCCCCCGGCTTGGCGCTGAAGTCATCCATCTCGTCGTTGAGCAGGAAACCGGCACC‑CTTCACCACGACGCCACTGCCGTAGTCGAGGTTGAGGGTGTAGGTATTGCTGACCGCGTTGCCCTGCTTGTCGACGATCGAGAAATGCGTGGTCTGGTGCGGTTCAAGGCCCGGCTTGA  >  NC_002947/5286000‑5286293

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