Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,316,022 A→C 10.3% A35A (GCA→GCC)  PP_3788 → non‑ribosomal peptide synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,316,0220AC10.3% 165.4 / 8.7 58A35A (GCA→GCC) PP_3788non‑ribosomal peptide synthetase
Reads supporting (aligned to +/- strand):  ref base A (27/25);  new base C (3/3);  total (30/28)
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.84e-01

ATCAGGGAATCAGTCTCGTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTG‑CGCAGCGGGGATGTGGTTGCGGTGAATGCGACAGTGCAACTGCGCTATCCCGTCG  >  NC_002947/4315875‑4316167
                                                                                                                                                   |                                                                                                                                                  
atcaGGGAATCAGTCTCGTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCa                                                                                                                                                    <  8:80236/148‑1 (MQ=255)
 tcaGGGAATCAGTCTCGTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCCGc                                                                                                                                                  >  1:5907/1‑149 (MQ=255)
   aGGGAATCAGTCTCGTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCgg                                                                                                                                                >  7:258550/1‑149 (MQ=255)
           aGTCTCGTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTcga                                                                                                                                         <  4:176903/148‑1 (MQ=255)
           aGTCTCGTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACACCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTcgac                                                                                                                                        >  7:150170/1‑149 (MQ=255)
            gTCTCGTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCCGCTGCATTcgatg                                                                                                                                       <  7:259969/149‑3 (MQ=255)
            gTCTCGTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCCGCTGCATTcgatg                                                                                                                                       <  2:5907/149‑3 (MQ=255)
              ctcGTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTcgacga                                                                                                                                      <  2:12046/148‑1 (MQ=255)
                 gTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGAc                                                                                                                                  >  2:6890/1‑149 (MQ=255)
                 gTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAAGACATCGAGCGTATCCTCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGAc                                                                                                                                  >  4:276410/1‑149 (MQ=255)
                  tACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTAGCTGGTCGAATACATCGAGCGTATCATCCACCGCTAACCATCGAATGCAGCGGCCTTCGACGAGGAcc                                                                                                                                 >  8:145689/1‑149 (MQ=255)
                    cGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCa                                                                                                                                <  8:150170/148‑1 (MQ=255)
                        acacaGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCt                                                                                                                           >  6:123721/1‑149 (MQ=255)
                          acaGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCtt                                                                                                                          >  1:56092/1‑148 (MQ=255)
                           caGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGa                                                                                                                        >  4:305828/1‑149 (MQ=255)
                           caGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGa                                                                                                                        <  8:258550/149‑1 (MQ=255)
                            aGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGAc                                                                                                                       <  1:6890/149‑1 (MQ=255)
                             gCAAGGAGCGCAGCATGAACCAGGCATTACTGACGCTCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGAc                                                                                                                       >  5:202030/1‑148 (MQ=255)
                                         gcaTGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTAc                                                                                                                   >  4:201812/1‑140 (MQ=255)
                                         gcaTGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTAc                                                                                                                   <  3:201812/140‑1 (MQ=255)
                                                ccAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTTGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATc                                                                                                    >  5:39493/1‑148 (MQ=255)
                                                 cAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCg                                                                                                   >  4:305133/1‑148 (MQ=255)
                                                     cATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTg                                                                                              <  2:93452/149‑1 (MQ=255)
                                                                   aCGGGAGGGCCTGCCTACCCCGTGCACTTCTTGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGc                                                                                <  6:39493/149‑1 (MQ=255)
                                                                         gggCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGaa                                                                                                             <  3:170085/114‑1 (MQ=255)
                                                                         gggCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGaa                                                                                                             >  4:170085/1‑114 (MQ=255)
                                                                          ggCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCt                                                                                                                     <  1:284944/105‑1 (MQ=255)
                                                                          ggCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCt                                                                                                                     >  2:284944/1‑105 (MQ=255)
                                                                           gcctgcctACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGc                                                                        <  5:123721/149‑1 (MQ=255)
                                                                           gcctgcctACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGc                                                                        <  2:56092/149‑1 (MQ=255)
                                                                             ctgcctACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCAt                                                                      >  5:21362/1‑149 (MQ=255)
                                                                                cctACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGaa                                                                   >  6:265470/1‑149 (MQ=255)
                                                                                    ccccGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATGGAATGCAGCGGCATTCGACGAGTACCAGCGGTTGAGCTTCAGCGAACTAAGGATCTAACTGCTCACACTGCACGCTCGGCTGCATTAACTCg                                                               <  1:234092/149‑1 (MQ=255)
                                                                                     cccGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACTCTCCGCTGCATGAACTCgg                                                              >  6:203362/1‑149 (MQ=255)
                                                                                     cccGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATAGAATGCAGCGGCATTAGACGAGGACCAGCGCTTGACCTACAGCGAAATAAGGATCGAACTGATCACACTGCAAGCTACGCTGCATGAACTCg                                                               >  3:256926/1‑148 (MQ=255)
                                                                                              ttCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTG‑CGCAg                                                    >  4:53803/1‑149 (MQ=255)
                                                                                              ttCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTG‑CGCAg                                                    >  2:163340/1‑149 (MQ=255)
                                                                                                ctcGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTg                                                                         >  7:247608/1‑127 (MQ=255)
                                                                                                ctcGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTg                                                                         <  8:247608/127‑1 (MQ=255)
                                                                                                    cTGGTCGAATACATCGACCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGACCTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTG‑CGCAGCGGGGa                                              <  1:129341/149‑1 (MQ=255)
                                                                                                        tCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTG‑CGCAGCGGGGAtgtg                                          <  6:21362/149‑1 (MQ=255)
                                                                                                        tCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTG‑CGCAGCGGGGAtgtg                                          >  1:124003/1‑149 (MQ=255)
                                                                                                          gAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTG‑CGCAGCGGGGATGTgg                                         <  5:265470/148‑1 (MQ=255)
                                                                                                                aTCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTG‑CGCAGCGGGGATGTGGTTGCGGt                                  <  3:305133/149‑1 (MQ=255)
                                                                                                                aTCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTG‑CGCAGCGGGGATGTGGTTGCGGt                                  >  7:299971/1‑149 (MQ=255)
                                                                                                                   gAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTG‑CGCAGCGGGGATGTGGTTGCGGTGaa                               <  3:305828/149‑1 (MQ=255)
                                                                                                                    atcGTATCATCCACCGCTAACCATCGAATGCAGCGGCATTCGAAGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTGTCGCAGCGGGGATGTGGTTGCGGTGaa                               <  1:163340/147‑1 (MQ=255)
                                                                                                                    aGCGTATCATCCACCGCTACCCATCGAATGCAGCAGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTG‑CGCAGCGGGGATGTGGTTGCGGTGAAt                              <  3:143503/149‑1 (MQ=255)
                                                                                                                             tcCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTG‑CGCAGCGGGGATGTGGTTGCg                                    <  3:102375/134‑1 (MQ=255)
                                                                                                                             tcCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTG‑CGCAGCGGGGATGTGGTTGCg                                    >  4:102375/1‑134 (MQ=255)
                                                                                                                               cACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTG‑CGCAGCGGGGATGTGGTTGCGGTGACTGCGACAGTGc                    >  7:291980/1‑148 (MQ=255)
                                                                                                                                          gtCGAATGCCGCTGCATTCGATGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTG‑CGCAGCGGGGATGTGGTTGCGGt                                  <  1:115814/122‑1 (MQ=255)
                                                                                                                                          gtCGAATGCCGCTGCATTCGATGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTG‑CGCAGCGGGGATGTGGTTGCGGt                                  >  2:115814/2‑123 (MQ=255)
                                                                                                                                          aTCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACGTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCTGCGTG‑CGCAGCGGGGATGTGGTTGCGGt                                  <  2:85386/123‑1 (MQ=255)
                                                                                                                                          aTCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTG‑CGCAGCGGGGATGTGGTTGCGGTGAATGCGACAGTGCAACTGCGCTATc        <  4:256926/149‑1 (MQ=255)
                                                                                                                                          aTCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTG‑CGCAGCGGGGATGTGGTTGCGGTGAATGCGACAGTGCAACTGCGCTATc        <  8:299971/149‑1 (MQ=255)
                                                                                                                                          aTCGAATGCAGCGGCATTCGACGAGGAACAGCGCTTGACCTACAGCCAACTAAGGATCGAACTCCTCACACTGCACGCTCCGCTGCATGAACTCGGCCTG‑CGCAGCGGGGATGTGGTTGCGGt                                  >  1:85386/1‑123 (MQ=255)
                                                                                                                                           tCGAATGCCGCTGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTG‑CGCAGCGGGGATGTGGTTGCGGTGAATGCGACAGTGCAACTGCGCTATcc       >  4:34289/1‑149 (MQ=255)
                                                                                                                                             gAATGCAGCGGCATTCGACGAGGACCAGGGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTG‑CGCAGCGGGGATGTGGTTGCGGTGAATGCGACAGTGCAACTGCGCTATCCcg     >  1:177099/1‑149 (MQ=255)
                                                                                                                                                   aGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTG‑CGCAGCGGGGATGTGGTTGCGGTGAATGCGACAGTGCAACTGCGCTATCCcgtcg  <  8:102613/146‑1 (MQ=255)
                                                                                                                                                   aGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTG‑CGCAGCGGGGATGTGGTTGCGGTGAATGCGACAGTGCAACTGCGCTATCCcgtcg  >  7:102613/1‑146 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
ATCAGGGAATCAGTCTCGTACGCCACACAGCAAGGAGCGCAGCATGAACCAGGCATTACTGACGATCACGGGAGGGCCTGCCTACCCCGTGCACTTCTCGCTGGTCGAATACATCGAGCGTATCATCCACCGCTACCCATCGAATGCAGCGGCATTCGACGAGGACCAGCGCTTGACCTACAGCGAACTAAGGATCGAACTGCTCACACTGCACGCTCGGCTGCATGAACTCGGCGTG‑CGCAGCGGGGATGTGGTTGCGGTGAATGCGACAGTGCAACTGCGCTATCCCGTCG  >  NC_002947/4315875‑4316167

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