Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 298,120 C→G 6.2% intergenic (‑117/+37) 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,1200CG 6.2% 209.6 / 4.3 65intergenic (‑117/+37)gshA/PP_0244glutamate‑‑cysteine ligase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (31/30);  new base G (2/2);  total (33/32)
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.11e-01

TCAGGCGGCGGTTGAGGAGTTCGCTCAAGATGTATTCCTTCACGCGTCAGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTCGGCGTCGACCACCAGGGTGCGCCGCCCGGCGTGCAGCACACGAGCGGTGCA  >  NC_002947/297976‑298267
                                                                                                                                                |                                                                                                                                                   
tCAGGCGGCGGTTGAGGAGTTCGCTCAAGATGTATTCCTTCACGCGTCAGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCAc                                                                                                                                                 >  4:202477/1‑149 (MQ=255)
             gaggagTTCGCTCAAGATGTATTCCTTCACGCGTCAGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTTCAGCAGTGCGCTGGCACTGCTGgaaccgc                                                                                                                                     <  2:88800/148‑5 (MQ=255)
               ggagTTCGCTCAAGATGTATTCCTTCACGCGTCAGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCg                                                                                                                                   >  5:278256/1‑148 (MQ=255)
                  gTTCGCTCAAGATGTATTCCTTCACGCGTCAGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGa                                                                                                                               <  3:202477/149‑1 (MQ=255)
                   ttCGCTCAAGATGTATTCCTTCACGCGTCAGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGa                                                                                                                               >  7:4689/1‑148 (MQ=255)
                   ttCGCTCAAGATGTATTCCTTCACGCGTCAGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGa                                                                                                                               >  7:39392/1‑148 (MQ=255)
                   ttCGCTCAAGATGTATTCCTTCACGCGTCAGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGa                                                                                                                               >  3:41414/1‑148 (MQ=255)
                      gCTCAAGATGTATTCCTTCACGCGTCAGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGTTACt                                                                                                                           >  8:302144/1‑149 (MQ=255)
                        tCAAGATGTATTCCTTCACGCGTCAGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGc                                                                                                                         >  8:218358/1‑149 (MQ=255)
                          aaGATGTATTCCTTCACGCGTCAGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTTCAGCAGTGCGCTGGCACTGCTGgaaccgc                                                                                                                                     >  4:149304/1‑131 (MQ=255)
                          aaGATGTATTCCTTCACGCGTCAGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTTCAGCAGTGCGCTGGCACTGCTGgaaccgc                                                                                                                                     <  3:149304/135‑5 (MQ=255)
                          aaGATGTATTCCTTCACGCGTCAGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCtt                                                                                                                       <  7:18770/149‑1 (MQ=255)
                           aGATGTATTCCTTCACGCGTCAGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTg                                                                                                                      <  5:132230/149‑1 (MQ=255)
                           aGATGTATTCCTTCACGCGTCAGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTg                                                                                                                      <  2:204152/149‑1 (MQ=255)
                            gATGTATTCCTTCACGCGTCAGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTgg                                                                                                                     >  1:239018/1‑149 (MQ=255)
                               gTATTCCTTCACGCGTCAGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGGGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGt                                                                                                                    >  5:311568/1‑147 (MQ=255)
                                 aTTCCTTCACGCGTCAGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGc                                                                                                                >  2:184264/1‑149 (MQ=255)
                                    ccTTCACGCGTCAGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTa                                                                                                              <  5:270069/148‑1 (MQ=255)
                                    ccTTCACGCGTCAGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTa                                                                                                              <  6:278256/148‑1 (MQ=255)
                                    ccTTCACGCGTCAGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTa                                                                                                              <  4:51383/148‑1 (MQ=255)
                                     cTTGACGCGTCAGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTagt                                                                                                            <  4:241091/149‑2 (MQ=255)
                                      ttCACGCGTCAGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGa                                                                                                                               >  3:326528/1‑129 (MQ=255)
                                      ttCACGCGTCAGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGa                                                                                                                               <  4:326528/129‑1 (MQ=255)
                                      ttCACGCGTCAGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTg                                                                                                                      <  6:189718/138‑1 (MQ=255)
                                      ttCACGCGTCAGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTg                                                                                                                      >  5:189718/1‑138 (MQ=255)
                                      ttCACGCGTCAGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGt                                                                                                                    >  5:336612/1‑140 (MQ=255)
                                      ttCACGCGTCAGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGt                                                                                                                    <  6:336612/140‑1 (MQ=255)
                                                aGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGaa                                                                                                  <  7:218358/148‑1 (MQ=255)
                                                 gTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAgg                                                                                                >  5:230180/1‑149 (MQ=255)
                                                     ccccAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTcc                                                                                            >  2:175664/1‑149 (MQ=255)
                                                      cccAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTg                                                                                                                   >  1:127839/1‑125 (MQ=255)
                                                      cccAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTg                                                                                                                   <  2:127839/125‑1 (MQ=255)
                                                      cccAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCt                                                                                           >  7:99034/1‑149 (MQ=255)
                                                       ccAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCtt                                                                                          <  8:99034/149‑1 (MQ=255)
                                                         aaTATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTgc                                                                                        <  4:41414/149‑1 (MQ=255)
                                                              gggggTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCtt                                                                                          >  6:7193/1‑142 (MQ=255)
                                                              gggggTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCtt                                                                                          <  5:7193/142‑1 (MQ=255)
                                                               ggggTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTg                                                                                  <  8:4689/149‑1 (MQ=255)
                                                               ggggTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTg                                                                                  <  7:144400/149‑1 (MQ=255)
                                                               ggggTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTg                                                                                  <  4:85079/149‑1 (MQ=255)
                                                                     gagaTAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGccac                                                                             <  1:184264/148‑1 (MQ=255)
                                                                        aTAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGc                                                                         >  1:102035/1‑149 (MQ=255)
                                                                           aCGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTg                                                                      <  6:311568/149‑1 (MQ=255)
                                                                                   cAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCt                                                               <  8:39392/148‑1 (MQ=255)
                                                                                   cAAGCGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCt                                                               >  3:50905/1‑148 (MQ=255)
                                                                                      gggTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAAc                                                           <  3:123324/149‑1 (MQ=255)
                                                                                              agaaAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTCgg                                                   >  3:212751/1‑149 (MQ=255)
                                                                                                  aGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTCGGCGTc                                               >  3:61468/1‑149 (MQ=255)
                                                                                                  aGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCACTTTGGCCACCAGCTTGTCGCCCTGAACTACGTCGGCGTc                                               >  1:302521/1‑149 (MQ=255)
                                                                                                      aCTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTCGGCGTCGac                                            >  5:60415/1‑148 (MQ=255)
                                                                                                          gCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTCGGCGTCGaccac                                         >  5:260007/1‑147 (MQ=255)
                                                                                                            gTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTc                                                     <  2:67905/133‑1 (MQ=255)
                                                                                                            gTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTc                                                     >  1:67905/1‑133 (MQ=255)
                                                                                                             ttGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGTTACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTCGGCGTCGACCACCAggg                                    <  7:302144/149‑1 (MQ=255)
                                                                                                               gTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTCGGCGTCGACCACCAGGGt                                   <  3:246544/148‑1 (MQ=255)
                                                                                                                  gCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTCGGCGTCGACCACCAGGGTgcgc                               <  2:239018/149‑1 (MQ=255)
                                                                                                                   ccTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTgg                                                                                 >  2:131676/1‑98 (MQ=255)
                                                                                                                   ccTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTgg                                                                                 <  1:131676/98‑1 (MQ=255)
                                                                                                                    cTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTCGGCGTCGACCACCAGGGTGcgccg                             >  2:237228/1‑149 (MQ=255)
                                                                                                                     tGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCGGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTCGGCGTCGACCACCAGGGTGcgccgc                            >  2:50951/1‑149 (MQ=255)
                                                                                                                              cgGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTCGGCGTCGACCACCAGGGTGCGCCGCCCGGCGTgc                   <  4:212751/149‑1 (MQ=255)
                                                                                                                               ggTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTCGGCGTCGACCACCAGGGTGCGCCGCCCGGCGTg                    >  4:90566/1‑147 (MQ=255)
                                                                                                                                ttttcagcagTGCGCTGGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTCGGCGTCGACCACCAGGGTGCGCCGCCCGGCGTgcag                 >  8:338472/5‑149 (MQ=255)
                                                                                                                                             cAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTCGGCGTCGACCACCAGGGTGCGCCGCCCGGCGTGCAGCACACGAGCGGTg    <  2:102035/149‑1 (MQ=255)
                                                                                                                                             cAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTCGGCGTCGACCACCAGGGTGCGCCGCCCGGCGTGCAGCACACGAGCGGTg    >  5:366900/1‑149 (MQ=255)
                                                                                                                                               gcgcACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTCGGCGTCGACCACCAGGGTGCGCCGCCCGGCGTGCAGCACACGAGCGGTGca  >  6:118117/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
TCAGGCGGCGGTTGAGGAGTTCGCTCAAGATGTATTCCTTCACGCGTCAGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTCGGCGTCGACCACCAGGGTGCGCCGCCCGGCGTGCAGCACACGAGCGGTGCA  >  NC_002947/297976‑298267

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