Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 298,117 C→G 5.7% intergenic (‑114/+40) gshA ← / ← PP_0244 glutamate‑‑cysteine ligase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947298,1170CG 5.7% 226.8 / 3.5 70intergenic (‑114/+40)gshA/PP_0244glutamate‑‑cysteine ligase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (35/31);  new base G (2/2);  total (37/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 = 8.41e-01

GCTCAGGCGGCGGTTGAGGAGTTCGCTCAAGATGTATTCCTTCACGCGTCAGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTCGGCGTCGACCACCAGGGTGCGCCGCCCGGCGTGCAGCACACGAGCGGTG  >  NC_002947/297974‑298265
                                                                                                                                               |                                                                                                                                                    
gCTCAGGCGGCGGTTGAGGAGTTCGCTCAAGATGTATTCCTTCACGCGTCAGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAgcg                                                                                                                                                  >  5:60748/1‑148 (MQ=255)
  tCAGGCGGCGGTTGAGGAGTTCGCTCAAGATGTATTCCTTCACGCGTCAGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCAc                                                                                                                                               >  8:104740/1‑149 (MQ=255)
    aGGCGGCGGTTGAGGAGTTCGCTCAAGATGTATTCCTTCACGCGTCAGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCActg                                                                                                                                             <  4:210379/149‑1 (MQ=255)
      gcggcgGTTGAGGAGTTCGCTCAAGATGTATTCCTTCACGCGTCAGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCActgct                                                                                                                                           <  1:261556/149‑1 (MQ=255)
         gcggTTGAGGAGTTCGCTCAAGATGTATTCCTTCACGCGTCAGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGaa                                                                                                                                        <  5:17011/149‑1 (MQ=255)
           ggTTGAGGAGTTCGCTCAAGATGTATTCCTTCACGCGTCAGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGaaaa                                                                                                                                      <  4:133198/149‑1 (MQ=255)
                     ttCGCTCAAGATGTATTCCTTCACGCGTCAGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCGCTGGCActgct                                                                                                                                           <  7:20000/134‑1 (MQ=255)
                     ttCGCTCAAGATGTATTCCTTCACGCGTCAGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCGCTGGCActgct                                                                                                                                           >  8:20000/1‑134 (MQ=255)
                     ttCGCTCAAGATGTATTCCTTCACGCGTCAGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGa                                                                                                                             >  4:31051/1‑148 (MQ=255)
                     ttCGCTCAAGATGTATTCCTTCACGCGTCAGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGa                                                                                                                             >  7:272963/1‑148 (MQ=255)
                     ttCGCTCAAGATGTATTCCTTCACGCGTCAGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGa                                                                                                                             >  8:166149/1‑148 (MQ=255)
                       cGCTCAAGATGTATTCCTTCACGCGTCAGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATAc                                                                                                                          >  2:108711/1‑149 (MQ=255)
                            aaGATGTATTCCTTCACGCGTCAGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCtt                                                                                                                     >  7:194091/1‑149 (MQ=255)
                             aGATGTATTCCTTCACGCGTCAGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGCGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTg                                                                                                                    <  6:60748/149‑1 (MQ=255)
                              gATGTATTCCTTCACGCGTCAGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTgg                                                                                                                   >  2:295498/1‑149 (MQ=255)
                                    ttCCTTCATGCGTCAGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCt                                                                                                             <  5:269342/149‑1 (MQ=255)
                                       cTTCACGCGTCAGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTa                                                                                                            >  2:286715/1‑147 (MQ=255)
                                       cTTCACGCGTCAGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTa                                                                                                            >  6:63737/1‑147 (MQ=255)
                                        ttCACGCGTCAGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCa                                                                                                                               <  3:257396/127‑1 (MQ=255)
                                        ttCACGCGTCAGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCa                                                                                                                               >  4:257396/1‑127 (MQ=255)
                                        ttCACGCGTCAGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTa                                                                                                            >  8:207284/1‑146 (MQ=255)
                                            cgcgTCAGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAAcc                                                                                                     <  7:258959/149‑1 (MQ=255)
                                                tCAGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGa                                                                                                 <  7:104740/149‑1 (MQ=255)
                                                  aGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAg                                                                                               <  5:63737/149‑1 (MQ=255)
                                                  aGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAg                                                                                               <  1:286715/149‑1 (MQ=255)
                                                       ccccAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGaa                                                                                                <  8:204332/143‑1 (MQ=255)
                                                       ccccAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGaa                                                                                                >  7:204332/1‑143 (MQ=255)
                                                        cccAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCt                                                                                         >  7:140616/1‑149 (MQ=255)
                                                             tatGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTACCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTgcgc                                                                                    <  1:108711/149‑1 (MQ=255)
                                                             tatGGGGGTAGAGATAACGGGCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTgcgc                                                                                    <  4:810/149‑1 (MQ=255)
                                                                gggggTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAgg                                                                                              >  6:311004/1‑136 (MQ=255)
                                                                gggggTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAgg                                                                                              <  5:311004/136‑1 (MQ=255)
                                                                  gggTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTg                                                                                <  1:295498/148‑1 (MQ=255)
                                                                       gagaTAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGccacc                                                                          <  2:243819/149‑1 (MQ=255)
                                                                        agaTAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGccacca                                                                         <  8:272963/149‑1 (MQ=255)
                                                                        agaTAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGccac                                                                           >  1:310915/1‑147 (MQ=255)
                                                                             aCGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTagag                                                                                                         <  5:238551/112‑1 (MQ=255)
                                                                             aCGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTagag                                                                                                         >  6:238551/1‑112 (MQ=255)
                                                                                gTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCg                                                                 <  7:207284/149‑1 (MQ=255)
                                                                                 tCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCAGCAGTGCGCTGGCACTGCTGgaa                                                                                                                                       <  2:121958/77‑1 (MQ=25)
                                                                                 tCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCAGCAGTGCGCTGGCACTGCTGgaa                                                                                                                                       >  1:121958/1‑77 (MQ=25)
                                                                                    aCAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCt                                                             >  5:298541/1‑149 (MQ=255)
                                                                                     cAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGccc                                                              <  8:140616/147‑1 (MQ=255)
                                                                                     cAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCt                                                             <  3:31051/148‑1 (MQ=255)
                                                                                           tAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTAc                                                      >  2:345311/1‑149 (MQ=255)
                                                                                           tAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACAGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTAc                                                      >  7:51689/1‑149 (MQ=255)
                                                                                             ggAAGAAAGGAACTAGAGTTGTCTCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACAGCAGATACTGCTTGGTGAGCTAGAGAAACGCGAAGGTTCCTAGCGCTTTGGACACCAGCTTGTCGCCCTGCACTActa                                                    >  2:27681/1‑147 (MQ=255)
                                                                                                    aGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTCGGCGTc                                             >  2:139407/1‑149 (MQ=255)
                                                                                                    aGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTCGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGCCgggcgg                                             >  1:110523/1‑145 (MQ=255)
                                                                                                      gAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTCGGCGt                                              <  6:235223/146‑1 (MQ=255)
                                                                                                      gAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTCGGCGt                                              >  5:235223/1‑146 (MQ=255)
                                                                                                         cTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTCGGCGTCGacca                                        <  8:194091/149‑1 (MQ=255)
                                                                                                            gCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTCGGCGTCGaccac                                       >  3:126750/1‑147 (MQ=255)
                                                                                                              gTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTCGGCGTCGACCACCAgg                                   >  5:245662/1‑149 (MQ=255)
                                                                                                                tGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTCGGCGTCGACCACCAGGGt                                 >  6:75407/1‑149 (MQ=255)
                                                                                                                   cGCCTGGTCTACGCGCTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTCGGCGTCGACCACCAGGGTgcg                              <  8:51689/149‑1 (MQ=255)
                                                                                                                    gCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTc                                                                                           <  8:89467/87‑1 (MQ=255)
                                                                                                                    gCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTc                                                                                           >  7:89467/1‑87 (MQ=255)
                                                                                                                    gCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGACCTGCACTACGTCGGCGTCGACCACCAGGGTgcgc                             >  8:94219/1‑149 (MQ=255)
                                                                                                                    gCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACAAGCTTGTCGCCCTGAACTACGTCGGCGTCGACCACCAGGGTgcgc                             >  8:94221/1‑149 (MQ=255)
                                                                                                                      cTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTCGGCGTCGACCACCAGGGTGcgccg                           >  7:247502/1‑149 (MQ=255)
                                                                                                                        ggTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTCGGCGTCGACCACCAGGGTGCgccgcc                         >  5:80525/1‑149 (MQ=255)
                                                                                                                         gTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTCGGCGTCGACCACCAGGGTGCGCCGccc                        >  3:161625/1‑149 (MQ=255)
                                                                                                                               gcgGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTCGGCGTCGACCACCAGGGTGCGCCGCCCGGCGTg                  >  3:149822/1‑149 (MQ=255)
                                                                                                                                cgGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTCGGCGTCGACCACCAGGGTGCGCCGCCCGGCGTgc                 <  2:310915/149‑1 (MQ=255)
                                                                                                                                 ggTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTCGGCGTCGACCACCAGGGTGCGCCGCCCGGCGTg                  >  7:166067/1‑147 (MQ=255)
                                                                                                                                  gTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTCGGCGTCGACCACCAGGGTGCGCCGCCCGGCGTgcag               <  6:298541/149‑1 (MQ=255)
                                                                                                                                   ttCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTCGGCGTCGACCACCAGGGTGCGCCGCCCGGCGTgcagc              <  6:245662/149‑1 (MQ=255)
                                                                                                                                     ccAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTCGGCGTCGACCACCAGGGTGCGCCGCCCGGCGTgcagca             <  1:345311/148‑1 (MQ=255)
                                                                                                                                          agTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTCGGCGTCGACCACCAGGGTGCGCCGCCCGGCGTGCAGCACACGAg       <  4:126750/149‑1 (MQ=255)
                                                                                                                                              ccAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTCGGCGTCGACCACCAGGGTGCGCCGCCCGGCGTGCAGCACACGAGCGGt   <  8:247502/149‑1 (MQ=255)
                                                                                                                                               cAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTCGGCGTCGACCACCAGGGTGCGCCGCCCGGCGTGCAGCACACGAGCGGTg  >  4:77951/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
GCTCAGGCGGCGGTTGAGGAGTTCGCTCAAGATGTATTCCTTCACGCGTCAGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTCGGCGTCGACCACCAGGGTGCGCCGCCCGGCGTGCAGCACACGAGCGGTG  >  NC_002947/297974‑298265

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: