Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,795,166 G→A 6.4% P426L (CCG→CTG)  pvdH ← diaminobutyrate‑2‑oxoglutarate transaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,795,1660GA 6.4% 182.1 / 2.1 62P426L (CCG→CTG) pvdHdiaminobutyrate‑2‑oxoglutarate transaminase
Reads supporting (aligned to +/- strand):  ref base G (29/29);  new base A (2/2);  total (31/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 = 5.26e-01

ACACGCTTGGGTTGGCGGCCAGCCTTGCAGGCCGACCGCGTCAAACAGTGGCATGCAGCGCGGCGTCAAACGCTGCCTACGGCGGCGGCGACTGCAGCGCTGAAGCGCTGGGCAACCTCGTCGATCTGCTCGGCACTGATGATCAACGGCGGCAGGAAACGCACCACCGCCCCATGCCGCCCACCCAACTCCAGGATCAGCCCGCGCTTGAGGCACTCGCGCTGCACCTTCGGCGCCAAGTCGCGATTAGCCGGTGGGTGGCCCAGAGCATCCGGCTGGCCTTGTGGGTCGACCAG  >  NC_002947/4795019‑4795314
                                                                                                                                                   |                                                                                                                                                    
acacGCTTGGGTTGGCGGCCAGCCTTGCAGGCCGACCGCGTCAAACAGTGGCATGCAGCGCGGCGTCAAACGCTGCCTACGGCGGCGGCGACTGCAGCGCTGAAGCGCTGGGCAACCTCGTCGATCTGCTCGGCACTGATGATCAAcgg                                                                                                                                                     >  2:304181/1‑149 (MQ=255)
  acGGTTGGGTTGGCGGCCAGCCTTGCAGGCCGACCGCGTCAAACAGTGGCATGCAGCGCGGCGTCAAACGCTGCCTACGGCGGCGGCGACTGCAGCGCTGAAGCGCTGGGCAACCTCGTCGATCTGCTCGGCACTGATGATCAAcggcg                                                                                                                                                   <  3:190578/149‑1 (MQ=255)
    gCTTGGGTTGGCGGCCAGCCTTGCAGGCCGACCGCGTCAAACAGTGGCATGCAGCGCGGCGTCAAACGCTGCCTACGGCGGCGGCGACTGCAGCGCTGAAGCGCTGGGCAACCTCGTCGATCTGCTCGGCACTGATGATCAAcggcgg                                                                                                                                                  >  7:76884/1‑148 (MQ=255)
    gCTTGGGTTGGCGGCCAGCCTTGCAGGCCGACCGCGTCAAACAGTGGCATGCAGCGCGGCGTCAAACGCTGCCTACGGCGGCGGCGACTGCAGCGCTGAAGCGCTGGGCAACCTCGTCGATCTGCTCGGCACTGATGATCAAcggcgg                                                                                                                                                  >  7:120691/1‑148 (MQ=255)
     cTTGGGTTGGCGGCCAGCCTTGCAGGCCGACCGCGTCAAACAGTGGCATGCAGCGCGGCGTCAAACGCTGCCTACGGCGGCGGCGACTGCAGCGCTGAAGCGCTGGGCAACCTCGTCGATCTGCTCGGCACTGATGATCAACGGCGGCa                                                                                                                                                <  4:96829/149‑1 (MQ=255)
           ttGGCGGCCAGCCTTGCAGGCCGACCGCGTCAAACAGTGGCATGCAGCGCGGCGTCAAACGCTGCCTACGGCGGCGGCGACTGCAGCGCTGAAGCGCTGGGCAACCTCGTCGATCTGCTCGGCACTGATGATCAACGGCGGCAGGAAAc                                                                                                                                          >  2:213808/1‑149 (MQ=255)
               cggcCAGCCTTGCAGGCCGACCGCGTCAAACAGTGGCATGCAGCGCGGCGTCAAACGCTGCCTACGGCGGCGGCGACTGCAGCGCTGAAGCGCTGGGCAACCTCGTCGATCTGCTCGGCACTGATGATCAACGGCGGCAGGAAACGcac                                                                                                                                      <  4:95572/149‑1 (MQ=255)
                  ccAGCCTTGCAGGCCGACCGCGTCAAACAGTGGCATGCAGCGCGGCGTCAAACGCTGCCTACGGCGGCGGCGACTGCAGCGCTGAAGCGCTGGGCAACCTCGTCGATCTGCTCGGCACTGATGATCAACGGCGGCAGGAAACGcaccac                                                                                                                                   >  5:311854/1‑149 (MQ=255)
                            aGGCCGACCGCGTCAAACAGTGGCATGCAGCGCGGCGTCAAACGCTGCCTACGGCGGCGGCGACTGCAGCGCTGAAGCGCTGGGCAACCTCGTCGATCTGCTCGGCACTGATGATCAACGGCGGCAGGAAACGCACCACCGCCCCATgc                                                                                                                         >  5:43645/1‑149 (MQ=255)
                             ggCCGACCGCGTCAAACAGTGGCATGCAGCGCGGCGTCAAACGCTGCCTACGGCGGCGGCGACTGCAGCGCTGAAGCGCTGGGCAACCTCGTCGATCTGCTCGGCACTGATGATCAACGGCGGCAGGAAACGCACCACCGCCCCATgcc                                                                                                                        <  1:304181/149‑1 (MQ=255)
                              gCCGACCGCGTCAAACAGTGGCATGCAGCGCGGCTTCAAACGCTGCCTACGGCGGCGGCGACTGCAGCGCTGAAGCGCTGGGCAACCTCGTCGATCTGCTCGGCACTGATGATCAACGGCGGCAGGAAACGCACCACCGCCCCATgccg                                                                                                                       <  8:120691/149‑1 (MQ=255)
                              gCCGACCGCGTCAAACAGTGGCATGCAGCGCGGCGTCAAACGCTGCCTACGGCGGCGGCGACTGCAGCGCTGAAGCGCTGGGCAACCTCGTCGATCTGCTCGGCACTGATGATCAACGGCGGCAGGAAACGCACCACCGCCCCATgccg                                                                                                                       >  8:110423/1‑149 (MQ=255)
                                     gcgTCAAACAGTGGCATGCAGCGCGGCGTCAAACGCTGCCTACGGCGGCGGCGACTGCAGCGCTGAAGCGCTGGGCAACCTCGTCGATCTGCTCGGCACTGATGATCAACGGCGGCAGGAAACGCACCACCGCCCCATGCCGcccaccc                                                                                                                <  6:311854/149‑1 (MQ=255)
                                      cgTCAAACAGTGGCATGCAGCGCGGCGTCAAACGCTGCCTACGGCGGCGGCGACTGCAGCGCTGAAGCGCTGGGCAACCTCGTCGATCTGCTCGGCACTGATGATCAACGGCGGCAGGAAACGCACCACCGCCCCATGCCGcccaccc                                                                                                                >  2:189438/1‑148 (MQ=255)
                                       gTCAAACAGTGGCATGCAGCGCGGCGTCAAACGCTGCCTACGGCGGCGGCGACTGCAGCGCTGAAGCGCTGGGCAACCTCGTCGATCTGCTCGGCACTGATGATCAACGGCGGCAGGAAACGCACCACCGCCCCATgccgcc                                                                                                                     >  5:153083/1‑142 (MQ=255)
                                       gTCAAACAGTGGCATGCAGCGCGGCGTCAAACGCTGCCTACGGCGGCGGCGACTGCAGCGCTGAAGCGCTGGGCAACCTCGTCGATCTGCTCGGCACTGATGATCAACGGCGGCAGGAAACGCACCACCGCCCCATgccgcc                                                                                                                     <  6:153083/142‑1 (MQ=255)
                                          aaaCAGTGGCATGCAGCGCGGCGTCAAACTCTGCCTACGGCGGCGGCGACTGCAGCGCTGAAGCGCTGGGCAACCTCGTCGATCTGCTCGGCACTGATTATCAACGGCGGCAGGAAACGCACCACCGCCCCATGCCGCCCACCCAACt                                                                                                            <  6:210526/148‑1 (MQ=255)
                                          aaaCAGTGGCATGCAGCGCGGCGTCAAACGCTGCCTACGGCGGCGGCGACTGCAGCGCTGAAGCGCTGGGCAACCTCGTCGATCTGCTCGGCACTGATGATCAACGGCGGCAGGAAACGCACCACCGCCCCATGCCGCCCACCCAACTc                                                                                                           >  7:300398/1‑149 (MQ=255)
                                           aaCAGTGGCATGCAGCGCGGCGTCAAACGCTGCCTACGGCGGCGGCGACTGCAGCGCTGAAGCGCTGGGCAACCTCGTCGATCTGCTCGGCACTGATGATCAACGGCGGCAGGAAACGCACCACCGCCCCATGCCGCCCACCCAACTcc                                                                                                          <  1:213808/149‑1 (MQ=255)
                                           aaCAGTGGCATGCAGCGCGGCGTCAAACGCTGCCTACGGCGGCGGCGACTGCAGCGCTGAAGCGCTGGGCAACCTCGTCGATCTGCTCGGCACTGATGATCAACGGCGGCAGGAAACGCACCACCGCCCCATGCCGCCCACCCAACTc                                                                                                           >  8:136427/1‑148 (MQ=255)
                                           aaCAGTGGCATGCAGCGCGGCGTCAAACGCTGCCTACGGCGGCGGCGACTGCAGCGCTGAAGCGCTGGGCAACCTCGTCGATCTGCTCGGCACTGATGATCAACGGCGGCAGGAAACGCACCACCGCCCCATGCCGCCCACCCAACTc                                                                                                           >  7:295154/1‑148 (MQ=255)
                                            aCAGTGGCATGCAGCGCGGCGTCAAACGCTGCCTACGGCGGCGGCGACTGCAGCGCTGAAGCGCTGGGCAACCTCGTCGATCTGCTCGGCACTGATGATCAACGGCGGCAGGAAACGCACCACCGCCCCATGCCGCCCACCCAACTc                                                                                                           >  3:66588/1‑147 (MQ=255)
                                               gTGGCATGCAGCGCGGCGTCAAACGCTGCCTACGGCGGCGGCGACTGCAGCGCTGAAGCGCTGGGCAACCTCGTCGATCTGCTCGGCACTGATGATCAACGGCGGCAGGAAACGCACCACCGCCCCATGCCGCCCACCCAACTCCAGGa                                                                                                      <  8:24503/149‑1 (MQ=255)
                                                  gCATGCAGCGCGGCGTCAAACGCTGCCTACGGCGGCGGCGACTGCAGCGCTGAAGCGCTGGGCAACCTCGTCGATCTGCTCGGCACTGATGATCAACGGCGGCAGGAAACGCACCACCGCCCCATGCCGCCCACCCAACTCCAGGATCa                                                                                                   <  4:66588/149‑1 (MQ=255)
                                                   cATGCAGCGCGGCGTCAAACGCTGCCTACGNCGGCGGCGACTGCAGCGCTGAAGCGCTGGGCAACCTCGTCGATCTGCTCGGCACTGATGATCAACAGCGGCAGGAAACGCACCACCGCCCCATGCCGCCCACCCAACTCCAGGATCAg                                                                                                  >  1:28493/1‑149 (MQ=255)
                                                   cATGCAGCGCGGCGTCAAACGCTGCCTACGGCGGCGGCGACTGCAGCGCTGAAGCGCTGGGCAACCTCGTCGATCTGCTCGGCACTGATGATCAACAGCGGCAGGAAACGCACCACCGCCCCATGCCGCCCACCCAACTCCAGGATCAg                                                                                                  >  1:23617/1‑149 (MQ=255)
                                                    aTGCAGCGCGGCGTCAAACGCTGCCTACGGCGGCGGCGACTGCAGCGCTGAAGCGCTGGGCAACCTCGTCGATCTGCTCGGCACTGATGATCAACGGCGGCAGGAAACGc                                                                                                                                        <  6:266022/110‑1 (MQ=255)
                                                    aTGCAGCGCGGCGTCAAACGCTGCCTACGGCGGCGGCGACTGCAGCGCTGAAGCGCTGGGCAACCTCGTCGATCTGCTCGGCACTGATGATCAACGGCGGCAGGAAACGc                                                                                                                                        >  5:266022/1‑110 (MQ=255)
                                                    aTGCAGCGCGGCGTCAAACGCTGCCTACGGCGGCGGCGACTGCAGCGCTGAAGCGCTGGGCAACCTCGTCGATCTGCTCGGCACTGATGATCAACGGCGGCAGGAAACGCACCACCGCCCCATGCCGCCCACCCAACTCCAGGATCAGc                                                                                                 >  2:158259/1‑149 (MQ=255)
                                                        aGCGCGGCGTCAAACGCTGCCTACGGCGGCGGCGACTGCAGCGCTGAAGCGCTGGGCAACCTCGTCGATCTGCTCGGCACTGATGATCAACGGCGGCAGGAAACGcaccac                                                                                                                                   <  4:120790/111‑1 (MQ=255)
                                                               cgTCAAACGCTGCCTACGGCGGCGGCGACTGCAGCGCTGAAGCGCTGGGCAACCTCGTCGATCTGCTCGGCACTGATGATCAACGGCGGCAGGAAACGCACCACCGCCCCATGCCGCCCACCCAACTCCAGGATCAGCCCGCGCTTGAg                                                                                      <  6:43645/149‑1 (MQ=255)
                                                                 tCAAACGCTGCCTACGGCGGCGGCGACTGCAGCGCTGAAGCGCTGGGCAACCTCGTCGATCTGCTCGGCACTGATGATCAACGGCGGCAGGAAACGCACCACCGCCCCATGCCGCCCACCCAACTCCAGGATCAGCCCGCGCTTGAgg                                                                                     >  5:228225/1‑148 (MQ=255)
                                                                 tCAAACGCTGCCTACGGCGGCGGCGACTGCAGCGCTGAAGCGCTGGGCAACCTCGTCGATCTGCTCGGCACTGATCATCAAAGGCGGCAGGAAACCCACCACCGCCCCAt                                                                                                                           >  7:2326/1‑110 (MQ=255)
                                                                      cGCTGCCTACGGCGGCGGCGACTGCAGCGCTGAAGCGCTGGGCAACCTCGTCGATCTGCTCGGCACTGATGATCAACGGCGGCAGGAAACGCACCACCGCCCCATGCCGCCCACCCAACTCCAGGATCAGCCCGCGCTTGAGGCACTcg                                                                               <  8:295154/149‑1 (MQ=255)
                                                                       gCTGCCTACGGCGGCGGCGACTGCAGCGCTGAAGCGCTGGGCAACCTCGTCGATCTGCTCGGCACTGATGATCAACGGCGGCAGGAAACGCACCACCGCCCCATGCCGCCCACCCAACTCCAGGATCAGCCCGCGCTTGAGGCACTcgc                                                                              <  4:35422/149‑1 (MQ=255)
                                                                          gCCTACGGCGGCGGCGACTGCAGCGCTGAAGCGCTGGGCAACCTCGTCGATCTGCTCGGCACTGATGATCAACGGCGGCAGGAAACGCACCACCGCCCCATGCCGCCCACCCAACTCCAGGATCAGCCCGCGCTTGAGGCACTCGCCCt                                                                           >  6:332041/1‑149 (MQ=255)
                                                                            cTACGGCGGCGGCTACTGCAGCGCTGAAGCGCTGGGCAACCTCGTCGATCTGCTCGGCAATGATGATCAACGGCGGCAGGAAACGCACCACCGCCCCATGCCGCCCACCCAACTCCAGGATCAGCCCGCGCTTGAGGCACTCGCGCTg                                                                          >  5:283279/1‑148 (MQ=255)
                                                                                  cggcggcgACTGCAGCGCTGAAGCGCTGGGCAACCTCGTCGATCTGCTCGGCACTGATGATCAACGGCGGCAGGAAACGCACCACCGCCCCATGCCGCCCACCCAACTCCAGGATCAGCCCGCGCTTGAGGCACTCGCGCTGCACCTTc                                                                   <  1:189438/149‑1 (MQ=255)
                                                                                     cggcgACTGCAGCGCTGAAGCGCTGGGCAACCTCGTCGATCTGCTCGGCACTGATGATCAACGGCGGCAGGAAACGCACCACCGCCCCATGCCGCCCACCCAACTCCAGGATCAGCCCGCGCTTGAGGCACTCGCGCTGCACCTTCGgc                                                                <  3:785/149‑1 (MQ=255)
                                                                                              cAGCGCTGAAGCGCTGGGCAACCTCGTCGATCTGCTCGGCACTGATGATCAACGGCGGCAGGAAACGCACCACCGCCCCATGCCGCCCACCCAACTCCAGGATCAGCCCGCGCTTGAGGCACTCGCGCTGCACCTTCGGCGCCAAGTcg                                                       >  2:77212/1‑149 (MQ=255)
                                                                                                gcgcTGAAGCGCTGGGCAACCTCGTCGATCTGCTCCGCACTGATGATCAACGGCGGCAGGAAACGCACCACCGCCCCATCCCGCCCCCCCAAATCCAGGAACAGCACGCGCTTGATGCACTCACGGAGACCCTTCGGCGCCAAGTcgcg                                                     >  6:232252/1‑149 (MQ=255)
                                                                                                 cgcTGAAGCGCTGGGCAACCTCGTCGATCTGCTCGGCACTGATGATCAACGGCGGCAGGAAACGCACCACCGCCCCATGCCGCCCACCCAACTCCAGGATCAGCCCGCGCTTGAGGCACTCGCGCTGCACCTTCGGCGCCAAGTCGCGa                                                    <  1:158259/149‑1 (MQ=255)
                                                                                                      aaGCGCTGGGCAACCTCGTCGATCTGCTCGGCACTGATGATCAACGGCGGCAGGAAACGCACCACCGCCCCATGCCGCCCACCCAACTCCAGGATCAGCCCGCGCTTGAGGCACTCGCGCTGCACCTTCGGCGCCAAGTc                                                        <  3:201491/140‑1 (MQ=255)
                                                                                                      aaGCGCTGGGCAACCTCGTCGATCTGCTCGGCACTGATGATCAACGGCGGCAGGAAACGCACCACCGCCCCATGCCGCCCACCCAACTCCAGGATCAGCCCGCGCTTGAGGCACTCGCGCTGCACCTTCGGCGCCAAGTc                                                        >  4:201491/1‑140 (MQ=255)
                                                                                                      aaGCGCTGGGCAACCTCGTCGATCTGCTCGGCACTGATGATCAACGGCGGCAGGAAACGCACCACCGCCCCATGCCGCCCACCCAACTCCAGGATCAGCCCGCGCTTGAGGCACTCGCGCTGCACCTTCGGCGCCAAGTCGCGATTAGc                                               >  8:258191/1‑149 (MQ=255)
                                                                                                        gcgcTGGGCAACCTCGTCGATCTGCTCGGCACTGATGATCAACGGCGGCAGGAAACGCACCACCGCCCCATGCCGCCCACCCAACTCCAGGATCAGCCCGCGCTTGAGGCACTCGCGCTGCACCTTCGGCGCCAAGTCGCGATTAGCCg                                             <  7:110423/149‑1 (MQ=255)
                                                                                                              ggCAACCTCGTCGATCTGCTCGTCACTGATGATCAACGGCGGCAGGAAACGCACCACCGCCCCATGCCGCCCACCCAACTCCAGGATCAGCCCGCGCTTGAGGCACTCGCGCTGCACCTTCGGCGCCAAGTCGCGATTAGCCggtgg                                         <  6:228225/147‑1 (MQ=255)
                                                                                                              ggCAACCTCGTCGATCTGCTCGGCACTGATGATCAACGGCGGCAGGAAACGCACCACCGCCCCATGCCGCCCACCCAACTCCAGGATCAGCCCGCGCTTGAGGCACTCGCGCTGCACCTTCGGCGCCAAGTCGCGATTAGCCggtgggt                                       <  7:136427/149‑1 (MQ=255)
                                                                                                                   ccTCGTCGATCTGCTCGGCACTGATGATCAACGGCGGCAGGAAACGCACCACCGCCCCATGCCGCCCACCCAACTCCAGGATCAGCCCGCGCTTGAGGCACTCGCGCTGCACCTTCGGCGCCAAGTCGCGATTAGCCggt                                           <  2:106475/140‑1 (MQ=255)
                                                                                                                   ccTCGTCGATCTGCTCGGCACTGATGATCAACGGCGGCAGGAAACGCACCACCGCCCCATGCCGCCCACCCAACTCCAGGATCAGCCCGCGCTTGAGGCACTCGCGCTGCACCTTCGGCGCCAAGTCGCGATTAGCCggt                                           >  1:106475/1‑140 (MQ=255)
                                                                                                                     tcgtcgATCTGCTCGGCACTGATGATCAACGGCGGCAGGAAACGCACCACCGCCCCATGCCGCCCACCCAACTCCAGGATCAGCCCGCGCTTGAGGCACTCGCGCTGCACCTTCGGCGCCAAGTCGCGATTAGCCGGTGGGTGGcc                                   >  7:81630/1‑146 (MQ=255)
                                                                                                                     tcgtcgATATGCTCGGCACTGATGATCAACGGCGGCAGGAAACGCACCAACGCCCCATGCCGCCCACCCAACTCCAGGATCAGCCAGCGCTTGAGGCACTCGCGCTCCAACTTAGGCGCCAAGTCGCCATTAGCCGgtgggtggac                                   >  6:177985/1‑144 (MQ=255)
                                                                                                                          gATCTGCTCGGCACTGATGATCAACGGCGGCAGGAAACGCACCACCGCCCCATGCCGCCCACCCAACTCCAGGATCAGCCCGCGCTTGAGGCACTCGCGCTGCACCTTCGGCGCCAAGTCGCGATTAGCCGGTGGGTGGCCCAGAGCa                            <  8:76884/148‑1 (MQ=255)
                                                                                                                                   ggCACTGATGATCAACGGCGGCAGGAAACGCACCACCGCCCCATGCCGCCCACCCAACTCCAGGATCAGCCCGCGCTTGAGGCACTCGCGCTGCACCTTCGGCGCCAAGTCGCGATTAGCCGGTGGGTGGCCCAGAGCATCCGGCTGGc                  <  8:81630/149‑1 (MQ=255)
                                                                                                                                      aCTGATGATCAACGGCGGCAGGAAACGCACCACCGCCCCATGCCGCCCACCCAACTCCAGGATCAGCCCGCGCTTGAGGGACTCGCGCTGCACCTTCGGCGCCAAGTCGCGATTAGCCGGTGGGTGGCCCAGAGCATCCGGCTGGCCtt               <  6:283279/149‑1 (MQ=255)
                                                                                                                                        tgatgaTCAACGGCGGCAGGAAACGCACCACCGCCCCATGCCGCCCACCCAACTCCAGGATCAGCCCGAGCATGAGGCACTCGCGATGCACCTTCGGCGCCAAGTCGCGATTAGCCGGTGGGTGGCCCAGAGCATCCGGCTCGCCTtgt             >  1:299852/1‑149 (MQ=255)
                                                                                                                                         gatgatCAACGGCGGCAGGAAACGCACCACCGCCCCATGCCGCCCACCCAACTCCAGGATCAGCCCGCGCTTGAGGCACTCGCGCTGCACCTTCGGCGCCAAGTCGCGATTAGCCGGTGGGTGGCCCAGAGCATCCGGCTGGCCTtgtg            >  3:72632/1‑149 (MQ=255)
                                                                                                                                            gatCAACGGCGGCAGGAAACGCACCACCGCCCCATGCCGCCCACCCAACTCCAGGATCAGCCCGCGCTTGAGGCACTcgc                                                                              >  4:58332/1‑80 (MQ=255)
                                                                                                                                            gatCAACGGCGGCAGGAAACGCACCACCGCCCCATGCCGCCCACCCAACTCCAGGATCAGCCCGCGCTTGAGGCACTcgc                                                                              <  3:58332/80‑1 (MQ=255)
                                                                                                                                             atCAACGGCGGCAGGAAACGCACCACCGCCCCATGCCGCCCACCCAACTCCAGGATCAGCCCGCGCTTGAGGCACTCGCGCTGCACCTTCGGCGCCAAGTCGCGATTAGCCGGTGGGTGGCCCAGAGCATCCGGCTGGCCTTGTGGGt         <  8:300398/148‑1 (MQ=255)
                                                                                                                                             atCAACAGCGGCAGGAAACGCACCACCGCCCCATGCCGCCCACCCAACTCCAGGATCAGCCCGCGCTTGAGGCACTCGCGCTGCACCTTCGGCGCCAAGTCGCGATTAGCCGGTGGGTGGCCCAGAGCATCCGGCTGGCCTTGTGGGTc        <  2:28493/149‑1 (MQ=255)
                                                                                                                                             atCAACAGCGGCAGGAAACGCACCACCGCCCCATGCCGCCCACCCAACTCCAGGATCAGCCCGCGCTTGAGGCACTCGCGCTGCACCTTCGGCGCCAAGTCGCGATTAGCCGGTGGGTGGCCCAGAGCATCCGGCTGGCCTTGTGGGTc        <  2:23617/149‑1 (MQ=255)
                                                                                                                                               ccaCGGCGGCAGGAAACGCACCAACGCCCCATGCCGCCCCCCCAACTCCAGGATCAGCCCGCGCTAGAGGCACTCGCGCTGCACCTTCGCCGACAAGTAGCGATTAGCCGGTGGGTGGCCCTGAGCATCCGGCTCGGCTTGTGGGTCGa      >  4:226655/3‑149 (MQ=255)
                                                                                                                                                 aCGGCGGCAGGAAACGCACCACCGCCCCATGCCGCCCACCCAACTCCAGGATCAGCCCGCGCTTGAGGCACTCGCGCTGCACCTTCGGCGCCAAGTCGCGATTAGCCGGTGGGTGGCCCAGAGc                             <  3:168872/124‑1 (MQ=255)
                                                                                                                                                 aCGGCGGCAGGAAACGCACCACCGCCCCATGCCGCCCACCCAACTCCAGGATCAGCCCGCGCTTGAGGCACTCGCGCTGCACCTTCGGCGCCAAGTCGCGATTAGCCGGTGGGTGGCCCAGAGCATCCGGCTGGCCTTGTGGGTCGAc     >  6:327287/1‑148 (MQ=255)
                                                                                                                                                 aCGGCGGCAGGAAACGCACCACCGCCCCATGCCGCCCACCCAACTCCAGGATCAGCCCGCGCTTGAGGCACACGCGCTGCACCATCGGCGCCAAGTCGCGATTAGCAGGAGGGTGGCCCAGAGc                             >  4:168872/1‑124 (MQ=255)
                                                                                                                                                   ggcggcAGGAAACGCACCACCGCCCCATGCCGCCCACCCAACTCCAGGATCAGCCCGCGCTTGAGGCACTCGCGCTGCACCTTCGGCGCCAAGTCGCGATTAGCCGGTGGGTGGCCCAGAGCATCCGGCTGGCCTTGTGGGTCGACCAg  <  5:203016/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
ACACGCTTGGGTTGGCGGCCAGCCTTGCAGGCCGACCGCGTCAAACAGTGGCATGCAGCGCGGCGTCAAACGCTGCCTACGGCGGCGGCGACTGCAGCGCTGAAGCGCTGGGCAACCTCGTCGATCTGCTCGGCACTGATGATCAACGGCGGCAGGAAACGCACCACCGCCCCATGCCGCCCACCCAACTCCAGGATCAGCCCGCGCTTGAGGCACTCGCGCTGCACCTTCGGCGCCAAGTCGCGATTAGCCGGTGGGTGGCCCAGAGCATCCGGCTGGCCTTGTGGGTCGACCAG  >  NC_002947/4795019‑4795314

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