Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,160,521:1 +C 6.9% coding (259/816 nt) pabC → 4‑amino‑4‑deoxychorismate lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,160,5211.C 6.9% 226.5 / 9.6 58coding (259/816 nt)pabC4‑amino‑4‑deoxychorismate lyase
Reads supporting (aligned to +/- strand):  ref base . (23/31);  new base C (2/2);  total (25/33)
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.37e-01

CTTGCTGGACGGCCACCTCGCACGCCTGGCACTGGGCTGCCAGCGCCTGGCGATCAACGCCGATCTGGCGCTGGTGCGCGACGAGCTGCTGCGCTACGCCAGCCAGCTGGGTGAGGGTGTCGTCAAACTCATCCTTACCCGTGGCGATAGCCAGCGTGGCTATGCGCCAGCTGTCGATGCAGCGCCGCGGCGTATCCTGCAAGGCAGCCCATTACCCCGTTATCCTGTCGAGCATGCTGAGCATGGCGTGCGTTTGTTCCCTTGCCATACCCGGCTTGGGGAACAACCACTGCT  >  NC_002947/2160376‑2160669
                                                                                                                                                  |                                                                                                                                                     
cTTGCTGGACGGCCACCTCGCACGCCTGGCACTGGGCTGCCAGCGCCTGGCGATCAACGCCGATCTGGCGCTGGTGCGCGACGAGCTGCTGCGCTACGCCAGCCAGCTGGGTGAGGGTGTCGTCAAACTCATCCTTACCCGTGGCGa                                                                                                                                                      >  3:120198/1‑147 (MQ=255)
cTTGCTGGACGGCAACCTCGCACGCCTGGCACTGGGCTGCCAGCGCCTGGCGATCAACGCCGATCTGGCGCTGGTGCGCGACGAGCTGCTGCGCTACGCCAGCCAGCTGGGTGAGGGTGTCGTCAAACTCATCCTTACCCGTGGCGa                                                                                                                                                      >  8:5229/1‑147 (MQ=255)
aagctgggACGGCCACCTCGCACGCCTGGCACTGGGCTGCCAGCGCCTGGCGATCAACGCCGATCTGGCGCTGGTGCGCGACGAGCTGCTGCGCTACGCCAGCCAGCTGGGTGAGGGTGTCGTCAAACTCATCCTTACCCGTGGCGa                                                                                                                                                      >  8:39288/7‑147 (MQ=255)
aagctgggACGGCCACCTCGCACGCCTGGCACTGGGCTGCCAGCGCCTGGCGATCAACGCCGATCTGGCGCTGGTGCGCGACGAGCTGCTGCGCTACGCCAGCCAGCTGGGTGAGGGTGTCGTCAAACTCATCCTTACCCGTGGCGATa                                                                                                                                                    <  7:39288/143‑1 (MQ=255)
   gCTGGACGGCCACCTCGCACGCCTGGCACTGGGCTGCCAGCGCCTGGCGATCAACGCCGATCTGGCGCTGGTGCGCGACGAGCTGCTGCGCTACGCCAGCCAGCTGGGTGAGGGTGTCGTCAAACTCATCCTTACCCGTGGCGATAg                                                                                                                                                   <  4:220603/147‑1 (MQ=255)
   gCTGGACGGCCACCTCGCACGCCTGGCACTGGGCTGCCAGCGCCTGGCGATCAACGCCGATCTGGCGCTGGTGCGCGACGAGCTGCTGCGCTACGCCAGCCAGCTGGGTGAGGGTGTCGTCAAACTCATCCTTACCCGTGGCGATAg                                                                                                                                                   <  2:208077/147‑1 (MQ=255)
   gCTGGACGGCCACCTCGCACGCCTGGCACTGGGCTGCCAGCGCCTGGCGATCAACGCCGATCTGGCGCTGGTGCGCGACGAGCTGCTGCGCTACGCCAGCCAGCTGGGTGAGGGTGTCGTCAAACTCATCCTTACCCGTGGCGATAGc                                                                                                                                                  <  4:143041/148‑1 (MQ=255)
      ggACGGCCACCTCGCACGCCTGGCACTGGGCTGCCAGCGCCTGGCGATCAACGCCGATCTGGCGCTGGTGCGCGACGAGCTGCTGCGCTACGCCAGCCAGCTGGGTGAGGGTGTCGTCAAACTCATCCTTACCCGTGGCGATAGCCAGc                                                                                                                                              >  2:242284/1‑149 (MQ=255)
          ggCCACCTCGCACGCCTGGCACTGGGCTGCCAGCGCCTGGCGATCAACGCCGATCTGGCGCTGGTGCGCGACGAGCTGCTGCGCTACGCCAGCCAGCTGGGTGAGGGTGTCGTCAAACTCATCCTTACCCGTGGCGATAGCCAGCGTgg                                                                                                                                          >  8:322794/1‑149 (MQ=255)
            ccACCTCGCACGCCTGGCACTGGGCTGCCAGCGCCTGGCGATCAACGCCGATCTGGCGCTGGTGCGCGACGAGCTGCTGCGCTACGCCAGCCAGCTGGGTGAGGGTGTCGTCAAACTCATCCTTACCCGTGGCGATAGCCAGCGTGGc                                                                                                                                         >  1:6275/1‑148 (MQ=255)
             cACCTCGCACGCCTGGCACTGGGCTGCCAGCGCCTGGCGATCAACGCCGATCTGGCGCTGGTGCGCGACGAGCTGCTGCGCTACGCCAGCCAGCTGGGTGAGGGTGTCGTCAAACTCATCCTTACCCGTGGCGCATAGCCacgctggc                                                                                                                                         >  1:162171/1‑141 (MQ=255)
              aCCTCGCACGCCTGGCACTGGGCTGCCAGCGCCTGGCGATCAACGCCGATCTGGCGCTGGTGCGCGACGAGCTGCTGCGCTACGCCAGCCAGCTGGGTGAGGGTGTCGTCAAACTCATCCTTACCCGTGGCGATAGCCAGCGTGGCTAt                                                                                                                                      >  6:170007/1‑149 (MQ=255)
                cTCGCACGCCTGGCACTGGGCTGCCAGTGCCTGGCGATCAACGCCGATCTGGCGCTGGTGCGCGACGAGCTGCTGCGGTACGCCAGCCAGCTGGGTGAGGGTGTCGTCAAACTCATCCTTACCCGTGGCGCATAGCCACGCTGGCTAtc                                                                                                                                     >  1:109206/1‑148 (MQ=255)
                 tCGCACGCCTGGCACTGGGCTGCCAGCGCCTGGCGATCAACGCCGATCTGGCGCTGGTGCGCGACGAGCTGCTGCGCTACGCCAGCCAGCTGGGTGAGGGTGTCGTCAAACTCATCCTTACCCGTGGCGATAGCCAGCGTGGCTATgcg                                                                                                                                   >  3:241516/1‑149 (MQ=255)
                 tCGCACGCCTGGCACTGGGCTGCCAGCCCCTGGCGATCAACCCCGATCTGGCGCTGGTGCGCGACGAGCTGCTGCCCTACGCCAGCCAGCTGGGTGAGGGTGTCGTCAAACTCATCCTTACCCGAGGCGAAAGCCAGCGTGGCTATgcg                                                                                                                                   >  2:300371/1‑149 (MQ=255)
                    cACGCCTGGCACTGGGCTGCCAGTGCCTGGCGATCAACGCCGATCTGGCGCTGGTGCGCGACGAGCTGCTGCGCTACGCCAGCCAGCTGGGTGAGGGTGTCGTCAAACTCATCCTTACCCGTGGCGCATAGCCACGCTGGCTATCGCCa                                                                                                                                <  2:109206/149‑1 (MQ=255)
                    cACGCCTGGCACTGGGCTGCCAGCGCCTGGCGATCAACGCCGATCTGGCGCTGGTGCGCGACGAGCTGCTGCGCTACGCCAGCCAGCTGGGTGAGGGTGTCGTCAAACTCATCCTTACCCGTGGCGCATAGCCACGCTGGCTATCGCCa                                                                                                                                <  2:162171/149‑1 (MQ=255)
                       gCCTGGCACTGGGCTGCCAGCGCCTGGCGATCATCGCCGATGTGGCGCTGGTGCGCGACGAGCTGCTGCGGTACGCCTGCCAGCTGGGTGAGGGTGTCGTCAAACTCATCCTTACCCGTGGCGATAGCCAGCGTGGCTATGCGCCAGCt                                                                                                                             <  5:242079/149‑1 (MQ=255)
                       gCCTGGCACTGGGCTGCCAGCGCCTGGCGATCAACGCCGATCTGGCGCTGGTGCGCGACGAGCTGCTGCGCTACGCCAGCCAGCTGGGTGAGGGTGTCGTCAAACTCATCCTTACCCGTGGCGATAGCCAGCGTGGCTATGCGCCAGCt                                                                                                                             >  8:1176/1‑149 (MQ=255)
                           ggCACTGGGCTGCCAGCGCCTGGCGATCAACGCCGATCTGGCGCTGGTGCGCGACGAGCTGCTGCGCTACGCCAGCCAGCTGGGTGAGGGTGTCGTCAAACTCATCCTTACCCGTGGCGATAGCCAGCGTGGCTATGCGCCAGCTGTc                                                                                                                          >  2:50683/1‑148 (MQ=255)
                            gCACTGGGCTGCCAGCGCCTGGCGATCAACGCCGATCTGGCGCTTGTGCGCGACGAGCTGCTGCGCTACGCCAGCCAGCTGGGTGAGGGTGTCGTCAAACTCATCCTTACCCGTGGCGATAGCCAGCGTGGCTATGCGCCAGCTGTCGa                                                                                                                        <  7:322794/149‑1 (MQ=255)
                            gCACTGGGCTGCCAGCGCCTGGCGATCAACGCCGATCTGGCGCTGGTGCGCGACGAGCTGCTGCGCTACGCCAGCCAGCTGGGTGAGGGTGTCGTCAAACTCATCCTTACCCGTGGCGATAGCCAGCGTGGCTATGCGCCAGCTGTCGa                                                                                                                        <  2:292473/149‑1 (MQ=255)
                                       ccAGCGCCTGGCGATCAACGCCGATCTGGCGCTGGTGCGCGACGAGCTGCTGCGCTACGCCAGCCAGCTGGGTGAGGGTGTCGTCAAACTCATCCTTACCCGTGGCGATAGCCAGCGTGGCTATGCGCCAGCTGTCGATGCAGcgccgc                                                                                                             >  4:190914/1‑149 (MQ=255)
                                        cAGCGCCTGGCGATCAACGCCGATCTGGCGCTGGTGCGCGACGAGCTGCTGCGCTACGCCAGCCAGCTGGGTGAGGGTGTCGTCAAACTCATCCTTACCCGTGGCGATAGCCAGCGTGGCTATGCGCCAGCTGTCGATGCAGCGCcgcg                                                                                                            <  4:120198/149‑1 (MQ=255)
                                         aGCGCCTGGCGATCAACGCCGATCTGGCGCTGGTGCGCGACGAGCTGCTGCGCTACGCCAGCCAGCTGGGTGAGGGTGTCGTCAAACTCATCCTTACCCGTGGCGATAGCCAGCGTGGCTATGCGCCAGTTGTCGATGCAGCGCCgcgg                                                                                                           <  8:50240/149‑1 (MQ=255)
                                          gcgcCTGGCGATCAACGCCGATCTGGCGCTGGTGCGCGACGAGCTGCTGCGCTACGCCAGCCAGCTGGGTGAGGGTGTCGTCAAACTCATCCTTACCCGTGGCGATAGCCAGCGTGGCTATGCGCCAGCTGTCGATGCAGCGCCgcggc                                                                                                          <  6:18509/149‑1 (MQ=255)
                                          gcgcCTGGCGATCAACGCCGATCTGGCGCTGGTGCGCGACGAGCTGCTGCGCTACGCCAGCCAGCTGGGTGAGGGTGTCGTCAAACTCATCCTTACCCGTGGCGATAGCCAGCGTGGCTATGCGCCAGCTGTCGATGCAGCGCCgcggc                                                                                                          <  3:190914/149‑1 (MQ=255)
                                          gcgcCTGGCGATCAACGCCGATCTGGCGCTGGTGCGCGACGAGCTGCTGCGCTACGCCAGCCAGCTGGGTGAGGGTGTCGTCAAACTCATCCTTACCCGTGGCGATAGCCAGCGTGGCTATGCGCCAGCTGTCGATGCAGCGCCgcggc                                                                                                          >  5:249691/1‑149 (MQ=255)
                                            gcCTGGCGATCAACGCCGATCTGGCGCTGGTGCGCGACGAGCTGCTGCGCTACGCCAGCCAGCTGGGTGAGGGTGTCGTCAAACTCATCCTTACCCGTGGCGATAGCCAGCGTGGCTATGCGCCAGCTGTCGATGCAGCGCCGCGGCGt                                                                                                        <  7:111477/149‑1 (MQ=255)
                                                 gCGATCAACGCCGATCTGGCGCTGGTGCGCGACGAGCTGCTGCGCTACGCCAGCCAGCTGGGTGAGGGTGTCGTCAAACTCATCCTTACCCGTGGCGATAGCCAGCGTGGCTATGCGCCAGCTGTCGATGCAGCGCCGCGGCGTATcc                                                                                                    <  1:242284/148‑1 (MQ=255)
                                                  cGATCAACGCCGATCTGGCGCTGGTGCGCGACGAGCTGCTGCGCTACGCCAGCCAGCTGGGTGAGGGTGTCGTCAAACTCATCCTTACCCGTGGCGATAGCCAGCGTGGCTATGCGCCAGCTGTCGATGCAGCGCCGCGGCGTATCCTg                                                                                                  <  1:50683/149‑1 (MQ=255)
                                                   gATCAACGCCGATCTGGCGCTGGTGCGCGACGAGCTGCTGCGCTACGCCAGCCAGCTGGGTGAGGGTGTCGTCAAACTCATCCTTACCCGTGGCGATAGCCAGCGTGGCTATGCGCCAGCTGTCGATGCAGCGCCGCGGCGTATCCTg                                                                                                  >  7:256132/1‑148 (MQ=255)
                                                    aTCAACGCCGATCTGGCGCTGGTGCGCGACGAGCTGCTGCGCTACGCCAGCCAGCTGGGTGAGGGTGTCGTCAAACTCATCCTTACCCGTGGCGATAGCCAGCGTGGCTATGCGCCAGCTGTCGATGCAGCGCCGCGGCGTATCCTGCa                                                                                                <  4:241516/149‑1 (MQ=255)
                                                       aaCGCCGATCTGGCGCTGGTGCGCGACGAGCTGCTGCGCTACGCCAGCCAGCTGGGTGAGGGTGTCGTCAAACTCATCCTTACCCGTGGCGATAGCCAGCGTGGCTATGCGCCAGCTGTCGATGCAGCGCCGCGGCGTATCCTGCAAg                                                                                              <  7:1176/148‑1 (MQ=255)
                                                           ccGATCTGGCGCTGGTGCGCGACGAGCTGCTGCGCTACGCCAGCCAGCTGGGTGAGGGTGTCGTCAAACTCATCCTTACCCGTGGCGATAGCCAGCGTGGCTATGCGCCAGCTGTCGATGCAGCGCCGCGGCGTATCCTGCAAGGCAGc                                                                                         >  5:247649/1‑149 (MQ=255)
                                                                   gcgcTGGTGCGCGACGAGCTGCTGCGCTACGCCAGCCAGCTGGGTGAGGGTGTCGTCAAACTCATCCTTACCCGTGGCGATAGCCAGCGTGGCTATGCGCCAGCTGTCGATGCAGCGCCGCGGCGTATCCTGCAAGGCAGCCCATTAc                                                                                  <  5:97569/148‑1 (MQ=255)
                                                                           aggtgtcgaGCTGCTGCGCTAGGCCAGCCAGCTGGGTGAGGGTGTCGTCAAACTCATCCTTACCCGTGGCGATAGCCAGCGTGGCTATGCGCCAGCTGTCGATGCAGCGCCGCGGCGTATCCTGCAAGGCAGCCCATTACCCCGTTATc                                                                        <  1:295205/143‑1 (MQ=255)
                                                                            cgcgACGAGCTGCTGCGCTACGCCAGCCAGCTGGGTGAGGGTGTCGTCAAACTCATCCTTACCCGTGGCGATAGCCAGCGTGGCTATGCGCCAGCTGTCGATGCAGCGCCGCGGCGTATCCTGCAAGGCAGCCCATTACCCCGTTATc                                                                        <  3:285705/148‑1 (MQ=255)
                                                                               gacgaGCTGCTGCGCTACGCCAGCCAGCTGGGTGAGGGTGTCGTCAAAGTCATCCTTACCCGTGGCGATAGCCAGCGTGGCTATGCGCCAGCTGTCGATGCAGCGCCGCGGCGTATCCTGCAAGGCAGCCCATTACCCCGTTATCCTGt                                                                    <  1:300371/149‑1 (MQ=255)
                                                                                        ctgcGCTACGCCAGCCAGCTGGGTGAGGGTGTCGTCAAACTCATCCTTACCCGTGGCGATAGCCAGCGTGGCTATGCGCCAGCTGTCGATGCAGCGCCGCGGCGTATCCTGCAAGGCAGCCCATTACCCCGTTATCCTGTCGAGCATg                                                            <  7:5229/148‑1 (MQ=255)
                                                                                        ctgcGCTACGCCAGCCAGCTGGGTGAGGGTGTCGTCAAACTCATCCTTACCCGTGGCGATAGCCAGCGTGGCTATGCGCCAGCTGTCGATGCAGCGCCGCGGCGTATCCTGCAAGGCAGCCCATTACCCCGTTATCCTGTCGAGCATGc                                                           >  5:140505/1‑149 (MQ=255)
                                                                                         tgcGCTACGCCAGCCAGCTGGGTGAGGGTGTCGTCAAACTCATCCTTACCCGTGGCGATAGCCAGCGTGGCTATGCGCCAGCTGTCGATGCAGCGCCGCGGCGTATCCTGCAAGGCAGCCCATTACCCCGTTATCCTGTCGa                                                                 <  5:18045/142‑1 (MQ=255)
                                                                                         tgcGCTACGCCAGCCAGCTGGGTGAGGGTGTCGTCAAACTCATCCTTACCCGTGGCGATAGCCAGCGTGGCTATGCGCCAGCTGTCGATGCAGCGCCGCGGCGTATCCTGCAAGGCAGCCCATTACCCCGTTATCCTGTCGa                                                                 >  6:18045/1‑142 (MQ=255)
                                                                                           cgcTACGCCAGCCAGCTGGGTGAGGGTGTCGTCAAACTCATCCTTACCCGTGGCGATAGCCAGCGTGGCTATGCGCCAGCTGTCGATGCAGCGCCGCGGCGTATCCTGCAAGGCAGCCCATTACCCCGTTATCCTGTCGAGCATGCt                                                          <  8:256132/147‑1 (MQ=255)
                                                                                            gcTACGCCAGCCAGCTGGGTGAGGGTGTCGTCAAACTCATCCTTACCCGTGGCGATAGCCAGCGTGGCTATGCGCCAGCTGTCGATGCAGCGCCGCGGCGTATCCTGCAAGGCAGCCCATTACCCCGTTATCCTGTCGAGCATGCTGAg                                                       >  4:68334/1‑149 (MQ=255)
                                                                                              tACGCCAGCCAGCTGGGTGAGGGTGTCGTCAAACTCATCCTTACCCGTGGCGATAGCCAGCGTGGCTATGCGCCAGCTGTCGATGCAGCGCCGCGGCGTATCCTGCAAGGCAGCCCATTACCCCGTTATCCTGTCGAGCATGCTGAGCa                                                     >  7:293718/1‑149 (MQ=255)
                                                                                                cGCCAGCCAGCTGGGTGAGGGTGTCGTCAAACTCATCCTTACCCGTGGCGATAGCCAGCGTGGCTATGCGCCAGCTGTCGATGCAGCGCCGCGGCGTATCCTGCAAGGCAGCCCATTACCCCGTTATCCTGTCGAGCATGCTGAGCATg                                                   <  5:170007/149‑1 (MQ=255)
                                                                                                 gccagccaGCTGGGTGAGGGTGTCGTCAAACTCATCCTTACCCGTGGCGATAGCCAGCGTGGCTATGCGCCAGCTGTCGATGCAGCGCCGCGGCGTATCCTGCAAGGCAGCCCATTACCCCGTTATCCTGTCGa                                                                 >  1:112386/1‑134 (MQ=255)
                                                                                                 gccagccaGCTGGGTGAGGGTGTCGTCAAACTCATCCTTACCCGTGGCGATAGCCAGCGTGGCTATGCGCCAGCTGTCGATGCAGCGCCGCGGCGTATCCTGCAAGGCAGCCCATTACCCCGTTATCCTGTCGa                                                                 <  2:112386/134‑1 (MQ=255)
                                                                                                   cagccagcTGGGTGAGGGTGTCGTCAAACTCATCCTTACCCGTGGCGATAGCCAGCGTGGCTATGCGCCAGCTGTCGATGCAGCGCCGCGGCGTATCCTGCAAGGCAGCCCATTACCCCGTTATCCTGTCGAGCATGCTGAGCATGgcg                                                >  2:235431/1‑149 (MQ=255)
                                                                                                      ccagcTGGGTGAGGGTGTCGTCAAACTCATCCTTACCCGTGGCGATAGCCAGCGTGGCTATGCGCCAGCTGTCGATGCAGCGCCGCGGCGTATCCTGCAAGGCAGcc                                                                                        >  8:237126/1‑107 (MQ=255)
                                                                                                          cTGGGTGAGGGTGTCGTCAAACTCATCCTTACCCGTGGCGATAGCCAGCGTGGCTATGCGCCAGCTGTCGATGCAGCGCCGCGGCGTATCCTGCAAGGCAGCCCATTACCCCGTTATCCTGTCGAGCATGCTGAGCATGGCGTGCGttt                                         <  6:247649/149‑1 (MQ=255)
                                                                                                             ggTGAGGGTGTCGTCAAACTCATCCTTACCCGTGGCGATAGCCAGCGTGGCTATGCGCCAGCTGTCGATGCAGCGCCGCGGCGTATCCTGCAAGGCAGCCCATTACCCCGTTATCCTGTCGAGCATGCTGAGCATGGCGTGCGTTTGtt                                      <  6:140505/149‑1 (MQ=255)
                                                                                                             ggTGAGGGTGTCGTCAAACTCATCCTTACCCGGGGCGAAAGCCAGCGTGGCTATGCGCCAGCTGTCGATGCAGCGGCGCGGGGGAACCTGCAAGGCATCCCCTTACCCCAGTTATCCTGTCGAGCATGCTGAGCATGGCGTGCGttt                                         <  1:235431/147‑1 (MQ=255)
                                                                                                               tGAGGGTGTCGTCAAACTCATCCTTACCCGTGGCGATAGCCAGCGTGGCTATGCGCCAGCTGTCGATGCAGCGCCGCGGCGTATCCTGCAAGGCAGCCCATTACCCCGTTATCCTGTCGAGCATGCTGAGCATGGCGTGCGTTTGTTcc                                    <  3:139237/149‑1 (MQ=255)
                                                                                                                       tcgtcAAACTCATCCTTACCCGTGGCGATAGCCAGCGTGGCTATGCGCCAGCTGGCGATGCAGCGCCGCGGCGTATCCTGCAAGGCAGCCCATTACCCCGTTATCCTGTCGAGCATGCTGAGCATGGCGTGCGTTTGTTCCCTTGCCa                             >  4:135216/1‑148 (MQ=255)
                                                                                                                                  aTCCTTACCCGTGGCGATAGCCAGCGTGGCTATGCGCCAGCTGTCGATGCAGCGCCGCGGCGTATCCTGCAAGGCAGCCCATTACCCCGTTATCCTGTCGAGCATGCTGAGCATGGCGTGCGTTTGTTCCCTTGCCATACCCGGCTTgg                 >  4:312937/1‑149 (MQ=255)
                                                                                                                                         cccGTGGCGATAGCCAGCGTGGCTATGCGCCAGCTGTCGATGCAGCGCCGCGGCGTATCCTGCAAGGCAGCCCATTACCCCGTTATCCTGTCGAGCATGCTGAGCATGGCGTGCGTTTGTTccc                                   >  7:168279/1‑124 (MQ=255)
                                                                                                                                         cccGTGGCGATAGCCAGCGTGGCTATGCGCCAGCTGTCGATGCAGCGCCGCGGCGTATCCTGCAAGGCAGCCCATTACCCCGTTATCCTGTCGAGCATGCTGAGCATGGCGTGCGTTTGTTccc                                   <  8:168279/124‑1 (MQ=255)
                                                                                                                                            gTGGCGATAGCCAGCGTGGCTATGCGCCAGCTGTCGATGCAGCGCCGCGGCGTATCCTGCAAGGCAGCCCATTACCCCGTTATCCTGTCGAGCATGCTGAGCATGGCGTGCGTTTGTTCCCTTGCCATACCCGGCTTgg                 >  2:219876/1‑139 (MQ=255)
                                                                                                                                            gTGGCGATAGCCAGCGTGGCTATGCGCCAGCTGTCGATGCAGCGCCGCGGCGTATCCTGCAAGGCAGCCCATTACCCCGTTATCCTGTCGAGCATGCTGAGCATGGCGTGCGTTTGTTCCCTTGCCATACCCGGCTTgg                 <  1:219876/139‑1 (MQ=255)
                                                                                                                                                 gATAGCCAGCGTGGCTATGCGCCAGCTGTCGATGCAGCGCCGCGGCGTATCCTGCAAGGCAGCCCATTACCCCGTTATCCTGTCGa                                                                 <  7:29013/86‑1 (MQ=255)
                                                                                                                                                 gATAGCCAGCGTGGCTATGCGCCAGCTGTCGATGCAGCGCCGCGGCGTATCCTGCAAGGCAGCCCATTACCCCGTTATCCTGTCGa                                                                 >  8:29013/1‑86 (MQ=255)
                                                                                                                                                 gATAGCCAGCGTGGCTATGCGCCAGCTGTCGATGCAGCGCCGCGGCGTATCCTGCAAGGCAGCCCATTACCCCGTTATCCTGTCGAGCATGCTGAGCATGGCGTGCGTTTGTTCCCTTGCCATACCCGGCTTGGGGAACAACCActgct  <  3:68334/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
CTTGCTGGACGGCCACCTCGCACGCCTGGCACTGGGCTGCCAGCGCCTGGCGATCAACGCCGATCTGGCGCTGGTGCGCGACGAGCTGCTGCGCTACGCCAGCCAGCTGGGTGAGGGTGTCGTCAAACTCATCCTTACCCGTGGCGATAGCCAGCGTGGCTATGCGCCAGCTGTCGATGCAGCGCCGCGGCGTATCCTGCAAGGCAGCCCATTACCCCGTTATCCTGTCGAGCATGCTGAGCATGGCGTGCGTTTGTTCCCTTGCCATACCCGGCTTGGGGAACAACCACTGCT  >  NC_002947/2160376‑2160669

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: