Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 504,094 C→G 6.1% G49G (GGC→GGG rpe → ribulose‑5‑phosphate 3‑epimerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947504,0940CG 6.1% 201.8 / 4.4 66G49G (GGC→GGGrperibulose‑5‑phosphate 3‑epimerase
Reads supporting (aligned to +/- strand):  ref base C (32/30);  new base G (2/2);  total (34/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 = 7.19e-01

CCTACGCTATTGCCCCCTCCATTCTCTCTGCCGATTTCGCCCGCCTGGGTGAGGACGTCGACAAGGTATTGGCTGCGGGTGCAGACATCGTCCACTTCGATGTCATGGACAATCACTATGTCCCCAACCTGACCATCGGCCCGATGGTTTGCACCGCTCTGCGCAAGTACGGCGTTACTGCGCCGATCGACGTGCACCTGATGGTCAGCCCGGTTGACCGCATCATTGGCGACTTTATCGAAGCGGGTGCCACCTACATCACCTTCCACCCTGAAGCCACGCAGCACA  >  NC_002947/503955‑504242
                                                                                                                                           |                                                                                                                                                    
ccTACGCTATTGCCCCCTCCATTCTCTCTGCCGATTTCGCCCGCCTGGGTGAGGACGTCGACAAGGTATTGGCTGCGGGTGCAGACATCGTCCACTTCGATGTCATGGACAATCACTATGTCCCCAACCTGACCATCGGGCCGATGGtc                                                                                                                                             <  3:149659/149‑2 (MQ=255)
ccTACGCTATTGCCCCCTCCATTCTCTCTGCCGATTTCGCCCGCCTGGGTGAGGACGTCGACAAGGTATTGGCTGCGGGTGCAGACATCGTCCACTTCGATGTCATGGACAATCACTATGTCCCCAACCTGACCATCGGCCCGATgg                                                                                                                                               >  8:139564/1‑147 (MQ=255)
ccTACGCTATTGCCCCCTCCATTCTCTCTGCCGATTTCGCCCGCCTGGGTGAGGACGTCGACAAGGTATTGGCTGCGGGTGCAGACATCGTCCACTTCGATGTCATGGACAATCACTATGTCCCCAACCTGACCATCGGCCCGATgg                                                                                                                                               <  7:139564/147‑1 (MQ=255)
  tACGCTATTGCCCCCTCCATTCTCTCTGCCGATTTCGCCCGCCTGGGTGAGGACGTCGACAAGGTATTGGCTGCGGGTGCAGACATCGTCCACTTCGATGTCATGGACAATCACTATGTCCCCAACCTGACCATCGGGCCGATGGt                                                                                                                                              >  4:157594/1‑146 (MQ=255)
    cGCTATTGCCCCCTCCATTCTCTCTGCCGATTTCGCCCGCCTGGGTGAGGACGTCGACAAGGTATTGGCTGCGGGTGCAGACATCGTCCACTTCGATGTCATGGACAATCACTATGTCCCCAACCTGACCATCGGGCCGATGGtca                                                                                                                                            <  3:157594/146‑3 (MQ=255)
       tATTGCCCCCTCCATTCTCTCTGCAGATTTCGCCAGACTGGGTGAGGACGTCGAAAAGGTATTGGCTGCGGGTGCAGACATCGTCCACTTCGATGTCATGGACAATCACTATGTCCACAACCTGACCATCGGCCCGATGGTTTGAACCg                                                                                                                                      >  4:252915/1‑149 (MQ=255)
                   cATTCTCTCTGCCGATTTCGCCCGCCTGGGTGAGGACGTCGACAAGGTATTGGCTGCGGGTGCAGACATCGTCCACTTCGATGTCATGGACAATCACTATGTCCCCAACCTGACCATCGGCCCGATGGTTTGCACCGCTCTGCGCAAGt                                                                                                                          <  7:147127/149‑1 (MQ=255)
                             gCCGATTTCGCCCGCCTGGGTGAGGACGTCGACAAGGTATTGGCTGCGGGTGCAGACATCGTCCACTTCGATGTCATGGACAATCACTATGTCCCCAACCTGACCATCGGCCCGATGGTTTGCACCGCTCTGCGCAAGTACGGCGTTAc                                                                                                                >  2:94050/1‑149 (MQ=255)
                                   ttCGCCCGCCTGGGTGAGGACGTCGACAAGGTATTCGCTGCCGGTGCAGACATCGTCCACTTCGAAGTCATGGACAATCACTATCTCCCCAACCTCACCATCGGCCCGATGGTTTGCCCCGCTCTGCGCAAGCACGCCGTTACTGCGcc                                                                                                          <  4:45746/149‑1 (MQ=255)
                                       cccgccTGGGTGAGGACGTCGACAAGGTATTGGCTGCGGGTGCAGACATCGTCCACTTCGATGTCATGGACAATCACTATGTCCCCAACCTGACCATCGGCCCGATGGTTTGCACCGCTCTGCGCAAGTACGGCGTTACTGCGCCGATc                                                                                                      >  8:271647/1‑149 (MQ=255)
                                       cccgccTGGGTGAGGACGTCGACAAGGTATTGGCTGCGGGTGCAGACATCGTCCACTTCGATGTCATGGACAATCACTATGTCCCCAACCTGACCATCGGCCCGATGGTTTGCACCGCTCTGCGCAAGTACGGCGTTACTGCGCCGATc                                                                                                      >  8:264090/1‑149 (MQ=255)
                                         cgccTGGGTGAGGACGTCGACAAGGTATTGGCTGCGGGTGCAGACATCGTCCACTTCGATGTCATGGACAATCACTATGCCCCCAACCTGACCATCGGCCCGATGGTTTGCACCGCTCTGCGCAAGTACGGCGTTACTGCGCCGATCGa                                                                                                    <  7:66275/149‑1 (MQ=255)
                                          gccTGGGTGAGGACGTCGACAAGGTATTGGCTGCGGGTGCAGACATCGTCCACTTCGATGTCATGGACAATCACTATGTCCCCAACCTGACCATCGGCCCGATGGTTTGCACCGCTCTGCGCAAGTACGGCGTTACTGCGCCGATCGAc                                                                                                   <  6:292760/149‑1 (MQ=255)
                                           ccTGGGTGAGGACGTCGACAAGGTATTGGCTGCGGGTGCAGACATCGTCCACTTCGATGTCATGGACAATCACTATGTCCCCAACCTGACCATCGGCCCGATGGTTTGCACCGCTCTGCGCAAGTACGGCGTTACTGCGCCGATCGACg                                                                                                  >  5:274831/1‑149 (MQ=255)
                                              gggTGAGGACGTCGACAAGGTATTGGCTGCGGGTGCAGACATCGTCCACTTCGATGTCATGGACAATCACTATGTCCCCAACCTGACCATCGGCCCGATGGTTTGCACCGCTCTGCGCAAGTACGGCGTTACTGCGCCGATCGACGTGc                                                                                               >  2:61089/1‑149 (MQ=255)
                                               ggTGAGGACGTCGACAAGGTATTGGCTGCGGGTGCAGACATCGTCCACTTCGATGTCATGGACAATCACTATGTCCCCAACCTGACCATCGGCCCGATGGTTTGCACCGCTCTGCGCAAGTACGGCGTTACTGCGCCGATCGACGTGCa                                                                                              >  4:53417/1‑149 (MQ=255)
                                                    ggACGTCGACAAGGTATTGGCTGCGGGTGCAGACATCGTCCACTTCGATGTCATGGACAATCACTATGTCCCCAACCTGACCATCGGCCCGATGGTTTGCACCGCTCTgc                                                                                                                                >  5:115010/1‑110 (MQ=255)
                                                    ggACGTCGACAAGGTATTGGCTGCGGGTGCAGACATCGTCCACTTCGATGTCATGGACAATCACTATGTCCCCAACCTGACCATCGGCCCGATGGTTTGCACCGCTCTgc                                                                                                                                <  6:115010/110‑1 (MQ=255)
                                                    ggACGTCGACAAGGTATTGGCTGCGGGTGCAGACATCGTCCACTTCGATGTCATGGACAATCACTATGTCCCCAACCTGACCATCGGCCCGATGGTTTGCACCGCTCTGCGCAAGTACGGCGTTACTGCGCCGATCGACGTGCACCTGa                                                                                         <  2:89555/149‑1 (MQ=255)
                                                    ggACGTCGACAAGGTATTGGCTGCGGGTGCAGACATCGTCCACTTCGATGTCATGGACAATCACTATGTCCCCAACCTGACCATCGGCCCGATGGTTTGCACCGCTCTGCGCAAGTACGGCGTTACTGCGCCGATCGACGTGCACCTGa                                                                                         <  1:94050/149‑1 (MQ=255)
                                                      aCGTCGACAAGGTATTGGCTGCGGGTGCAGACATCGTCCACTTCGATGTCATGGACAATCACTATGTCCCCAACCTGACCATCGGCCCGATGGTTTGCACCGCTCTGCGCAAGTACGGCGTTACTGCGCCGATCGACGTGCACCTGATg                                                                                       <  2:247699/149‑1 (MQ=255)
                                                             cAAGGTATTGGCTGCGGGTGCAGACATCGTCCACTTCGATGTCATGGACAATCACTATGTCCCCAACCTGACCATCGGCCCGATGGTTTGCACCGCTCTGCGCAAGTACGGCGTTACTGCGCCGATCGACGTGCACCTGATGGTCAGcc                                                                                >  8:133088/1‑149 (MQ=255)
                                                             cAAGGTATTGGCTGCGGGTGCAGACATCGTCCACTTCGATGTCATGGACAATCACTATGTCCCCAACCTGACCATCGGCCCGATGGTTTGCACCGCTCTGCGCAAGTACGGCGTTACTGCGCCGATCGACGTGCACCTGATGGTCAGcc                                                                                <  3:53417/149‑1 (MQ=255)
                                                                ggTATTGGCTGCGGGTGCAGACATCGTCCACTTCGATGTCATGGACAATCACTATGTCCCCAACCTGACCATCGGCCCGATGGTTTGCACCGCTCTGCGCAAGTACGGCGTTACTGCGCCGATCGACGt                                                                                                 <  3:31474/129‑1 (MQ=255)
                                                                ggTATTGGCTGCGGGTGCAGACATCGTCCACTTCGATGTCATGGACAATCACTATGTCCCCAACCTGACCATCGGCCCGATGGTTTGCACCGCTCTGCGCAAGTACGGCGTTACTGCGCCGATCGACGt                                                                                                 >  4:31474/1‑129 (MQ=255)
                                                                ggTATTGGCTGCGGGTGCAGACATCGTCCACTTCGATGTCATGGACAATCACTATGTCCCCAACCTGACCATCGGCCCGATGGTTTGCACCGCTCTGCGCAAGTACGGCGTTACTGCGCCGATCGACGTGCACCTGATGGTCAGCCCgg                                                                             >  6:39925/1‑149 (MQ=255)
                                                                 gTATTGGCTGCGGGTGCAGACATCGTCCACTTCGATGTCATGGACAATCACTATGTCCCCAACCTGACCATCGGCCCGATGGTTTGCACCGCTCTGCGCAAGTACGGCGTTACTGCGCCGATCGACGTGCACCTGATGGTCAGCCCGGt                                                                            >  6:280907/1‑149 (MQ=255)
                                                                  tATTGGCTGCGGGTGCAGACATCGTCCACTTCGATGTCATGGACAATCACTATGTCCCCAACCTGACCATCGGCCCGATGGTTTGCACCGCTCTGCGCAAGTACGGCGTTACTGCGCCGATCGACGTGCACCTGATGGTCAGCCCGGtt                                                                           <  7:271647/149‑1 (MQ=255)
                                                                      ggCTGCGGGTGCAGACATCGTCCACTTCGATGTCATGGACAATCACTATGTCCCCAACCTGACCATCGGCCCGATGGTTTGCACCGCTCTGCGCAAGTACGGCGTTACTGCGCCGATCGACGTGCACCTGATGGTCAGCCCGGTTGAc                                                                        <  8:83204/148‑1 (MQ=255)
                                                                      ggCAGCGGGTGCAGACATCGTCCACTTCGATGTCATGGACAATCACTATGTCCCCAACCTGACCATCGGCCCGATGGTTTGCACCGCTCTGCGCAAGTACGGCGTTACTGCGCCGATCGACGTGCACCTGATGGTCAGCCCGGTTGAcc                                                                       >  7:53451/1‑149 (MQ=255)
                                                                       gCTGCGGGTGCAGACATCGTCCACTTCGATGTCATGGACAATCACTATGTCCCCAACCTGACCATCGGCCCGATGGTTTGCACCGCTCTGCGCAAGTACGGCGTTACTGCGCCGATCGACGTGCACCTGATGGTCAGCCCGGTTGAcc                                                                       <  8:53451/148‑1 (MQ=255)
                                                                         tGCGGGTGCAGACATCGTCCACTTCGATGTCATGGACAATCACTATGTCCCCAACCTGACCATCGGCCCGATGGTTTTCACCGCTCTGCGCAAGTACGGAGTTACTGCGCCGATCGACGTGCACCTGATGGTCAGCCCGGTTGACCGca                                                                    >  3:231528/1‑149 (MQ=255)
                                                                          gCGGGTGCAGACATCGTCCACTTCGATGTCATGGACAATCACTATGTCCCCAACCTGACCATCGGCCCGATGGTTTGCACCGCTCTGCGCAAGTACGGCGTTACTGCGCCGATCGACGTGCACCTGATGGTCAGCCCGGTTGACCGcat                                                                   >  2:95230/1‑149 (MQ=255)
                                                                              gTGCAGACATCGTCCACTTCGATGTCATGGACAATCACTATGTCCCCAACCTGACCATCGGCCCGATGGTTTGCACCGCTCTGCGCAAGTACGGCGTTACTGCGCCGATCGACGTGCACCTGATGGTCAGCCCGGTTGACCGCATCAtt                                                               <  6:30100/149‑1 (MQ=255)
                                                                                 cAGACATCGTCCACTTCGATGTCATGGACAATCACTATGTCCCCAACCTGACCATCGGCCCGATGGTTTGCACCGCTCTGCGCAAGTACGGCGTTACTGCGCCGATCGACGTGCACCTGATGGTCAGCCCGGTTGACCGCATCATTgg                                                             <  7:264090/148‑1 (MQ=255)
                                                                                 cAGACATCGTCCACTTCGATGTCATGGACAATCACTATGTCCCCAACCTGACCATCGGCCCGATGGTTTGCACCGCTCTGCGCAAGTACGGCGTTACTGCGCCGATCGACGTGCACCTGATGGTCAGCCCGGTTGACCGCATCATTGGc                                                            >  6:157421/1‑149 (MQ=255)
                                                                                  aGACATCGTCCACTTCGATGTCATGGACAATCACTATGTCCCCAACCTGACCATCGGCCCGATGGTTTGCACCGCTCTGCGCAAGTACGGCGTTACTGCGCCGATCGACGTGCACCTGATGGTCAGCCCGGTTGACCGCATCATTGGc                                                            >  4:48605/1‑148 (MQ=255)
                                                                                      aTCGTCCACTTCGATGTCATGGACAATCACTATGTCCCCAACCTGACCATCGGCCCGATGGTTTGCACCGCTCTGCGCAAGTACGGCGTTACTGCGCCGATCGACGTGCACCTGATGGTCAACCGGGCTGAccatcagg                                                                 <  1:8229/139‑7 (MQ=255)
                                                                                      aTCGTCCACTTCGATGTCATGGACAATCACTATGTCCCCAACCTGACCATCGGCCCGATGGTTTGCACCGCTCTCCGCAAGTACGGCGTTACTGCGCCGATCGACGTGCACCTGATGGTCAACCGGGCTGAccatcagg                                                                 >  2:8229/1‑133 (MQ=255)
                                                                                         gTCCACTTCGATGTCATGGACAATCACTATGTCCCCAACCTGACCATCGGCCCGATGGTTTGCACCGCTCTGCGCAAGTACGGCGTTACTGCGCCGATCGACGTGCACCTGATGGTCAGCCCGGTTGACCGCATCATTGGCGACTTTa                                                     <  7:133088/148‑1 (MQ=255)
                                                                                          tCCACTTCGATGTCATGGACAATCACTATGTCCCCAACCTGACCATCGGCCCGATGGTTTGCACCGCTCTGCGCAAGTACGGCGTTACTGCGCCGATCGACGTGCACCTGATGGTCAGCCCGGTTGACCGCATCATTGGCGACTTTATc                                                   >  1:187844/1‑149 (MQ=255)
                                                                                             aCTTCGATGTCATGGACAATCACTATGTCCCCAACCTGACCATCGGCCCGATGGTTTGCACCGCTCTGCGCAAGTACGGCGTTACTGCGCCGATCGACGTGCACCTGATGGTCAGCCCGGTTGACCGCATCATTGGCGACTTTATCg                                                  >  7:31886/1‑147 (MQ=255)
                                                                                               ttCGATGTCATGGACAATCACTATGTCCCCAACCTGACCATCGGCCCGATGGTTTGCACCGCTCTGCGCAAGTACGGCGTTACTGCGCCGATCGACGTGCACCTGATGGTCAGCCCGGTTGACCGCATCATTGGCGACTTTATCGAAGc                                              >  4:239801/1‑149 (MQ=255)
                                                                                                 cGATGTCATGGACAATCACTATGTCCCCAACCTGACCATCGGCCCGATGGTTTGCACCGCTCTGCGCAAGTACGGCGTTACTGCGCCGATCGACGTGCACCTGATGGTCAGCCCGGTTGACCGCATCATTGGCGACTTTATCGAAGc                                              <  1:61089/147‑1 (MQ=255)
                                                                                                        attgACAATCACTATGTCCCCAACCTGACCATCGGCCCGTTGGTTTGCAGCGGTCTGCGCAAGTACGGCGTTACTGCTCCGATAGACGTGCACCTGATGGTCAGCCCGGTTGACCGCATCATTg                                                              <  4:143669/121‑1 (MQ=255)
                                                                                                        aTGGACAATCACTATGTCCCCAACCTGACCATCGGCCCGATGGTTTGCACCGCTCTGCGCAAGTACGGCGTTACTGCGCCGATCGACGTGCACCTGATGGTCAGCCCGGTTGACCGCATCATTg                                                              >  3:143669/1‑124 (MQ=255)
                                                                                                         tGGACAATCACTATGTCCCCTTCCTGACCATCGGCCCGATGGTTTGCACCGCTCTGCGCAAGTACGGCGTTACTGCGCCGATCGACGTTCACCTGATGGTCAGCGCGGTTGACCGCATCATTGGCGACTTTATCGAAGCGGGTGCCTcc                                    <  4:231528/149‑1 (MQ=255)
                                                                                                         tGGACAATCACTATGTCCCCAACCTGACCATCGGCCCGATGGTTTGCACCGCTCTGCGCAAGTACGGCGTTACTGCGCCGATCGACGTGCACCTGATGGTCAGCCCGGTTGACCGCATCATTGGCGACTTTATCGAAGCGGGt                                          >  4:152272/1‑143 (MQ=255)
                                                                                                         tGGACAATCACTATGTCCCCAACCTGACCATCGGCCCGATGGTTTGCACCGCTCTGCGCAAGTACGGCGTTACTGCGCCGATCGACGTGCACCTGATGGTCAGCCCGGTTGACCGCATCATTGGCGACTTTATCGAAGCGGGt                                          <  3:152272/143‑1 (MQ=255)
                                                                                                          ggACAATCACTATGTCCCCAACCTGACCATCGGCCCGATGGTTTGCACCGCTCTGCGCAAGTACGGCGTTACTGCGCCGATCGACGTGCACCTGATGGTCAGCCCGGTTGACCGCATCATTGGCGACTTTATCGaa                                                <  3:269619/136‑1 (MQ=255)
                                                                                                          ggACAATCACTATGTCCCCAACCTGACCATCGGCCCGATGGTTTGCACCGCTCTGCGCAAGTACGGCGTTACTGCGCCGATCGACGTGCACCTGATGGTCAGCCCGGTTGACCGCATCATTGGCGACTTTATCGaa                                                >  4:269619/1‑136 (MQ=255)
                                                                                                            aCAATCACTATGTCCCCAACCTGACCATCGGCCCGATGGTTTGCACCGCTCTGCGCAAGTACGGCGTTACTGCGCCGATCGACGTGCACCTGATGGTCAGCCCGGTTGACCGCATCATTgg                                                             <  2:112748/121‑1 (MQ=255)
                                                                                                            aCAATCACTATGTCCCCAACCTGACCATCGGCCCGATGGTTTGCACCGCTCTGCGCAAGTACGGCGTTACTGCGCCGATCGACGTGCACCTGATGGTCAGCCCGGTTGACCGCATCATTgg                                                             >  1:112748/1‑121 (MQ=255)
                                                                                                             cAATCACTATGTCCCCAACCTGACCATCGGCCCGATGGTTTGCACCGCTCTGCGCAAGTACGGCGTTACTGCGCCGATCGACGTGCACCTGATGGTCAGCCCGGTTGACCGCATCATTGGCGACTTTATCGAAGCGGGTGCCACCTACa                                <  3:48605/149‑1 (MQ=255)
                                                                                                               aTCACTATGTCCCCAACCTGACCATCGGCCCGATGGTTTGCACCGCTCTGCGCAAGTACGGCGTTACTGCGCCGATCGACGTGCACCTGATGGTCAGCCCGGTTGACCGCATCATTGGCGACTTTATCGAAGCGGGTGCCACCTACATc                              <  5:39925/149‑1 (MQ=255)
                                                                                                                  aCTATGTCCCCAACCTGACCATCGGCCCGATGGTTTGCACCGCTCTGCGCAAGTACGGCGTTACTGCGCCGATCGACGTGCACCTGATGGTCAGCCCGGTTGACCGCATCATTGGCGACTTTATCGAAGCGGGTGCCACCTACATCAcc                           <  5:280907/149‑1 (MQ=255)
                                                                                                                          cccAACCTGACCATCGGCCCGATGGTTTGCACCGCTCTGCGCAAGTACGGCGTTACTGCGCCGATCGACGTGCACCTGATGGTCAGCCCGGTTGACCGCATCATTGGCGACTTTATCGAAGCGGGTGCCACCTACATCACCTTGCAccc                   >  4:144177/1‑149 (MQ=255)
                                                                                                                          cccAACCTGACCATCGGCCCGATGGTTTGCACCGCTCTGCGCAAGTACGGCGTTACTGCGCCGATCGACGTGCACCTGATGGTCAGCCCGGTTGACCGCATCATTGGCGACTTTATCGAAGCGGGTGCCACCTACATCACCTTCCAccc                   >  4:178511/1‑149 (MQ=255)
                                                                                                                          cccAACCTGACCATCGGCCCGATGGTTTGCACCGCTCTGCGCAAGTACGGCGTTACTGCGCCAATCGACGTGCACCTGATGGTCAGCCCGGTTGACCGCATCATTGGCGACTTTATCGAAGCGGGTGCCACCTACATCACCTTCCAccc                   >  1:127107/1‑149 (MQ=255)
                                                                                                                              aCCTGACCATCGGCCCGATGGTTTGCACCGCTCTGCGCAAGTACGGCGTTACTGCGCCGATCGACGTGCACCTGATGGTCAGCCCGGTTGACCGCATCATTGTGGACTTTATCGAAGCGGGTGCCACCTACATCACCTTCCACCCTGaa               <  3:252915/149‑1 (MQ=255)
                                                                                                                              aCCTGACCATCGGCCCGATGGTTTGCACCGCTCTGCGCAAGTACGGCGTTACTGCGCCGATCGACGTGCACCTGATGGTCAGCCCGGTTGACCGCATCATTGGCGACTTTATCGAAGCGGGTGCCACCTACATCACCTTCCACCCTGaa               <  5:157421/149‑1 (MQ=255)
                                                                                                                                  gACCATCGGCCCGATGGTTTGCACCGCTCTGCGCAAGTACGGCGTTACTGCGCCGATCGACGTGCACCTGATGGTCAGCCCGGTTGACCGCATCATTGGCGACTTTATCGAAGCGGGTGCCACCTACATCACCTTCCACCCTGAAGCCa           >  8:36815/1‑149 (MQ=255)
                                                                                                                                  aaCCATCGGGCCGATGGTTTGCACCGCTCTGCGCAAGTACGGCGTTACTGCGCCGATCGACGTGCACCTGATGGTCAGCCCGGTTGACCGCATCATTGGCGACTTTATCGAAGCGGGTGCCACCTACATCACCTTCCACCCTGAAGCCa           >  3:29780/2‑149 (MQ=255)
                                                                                                                                   aCCATCGGCCCGATGGTTTGCACCGCTCTGCGCAAGTACGGCGTTACTGCGCCGATCGACGTGCACCTGATGGTCAGCCCGGTTGACCGCATCATTGGCGACTTTATCGAAGCGGGTGCCACCTACATCACCTTCCACCCTGAAGCCa           <  2:187844/148‑1 (MQ=255)
                                                                                                                                        cGGCCCGATGGTTTGCACCGCTCTGCGCAAGTACGGCGTTACTGCGCCAATCGACGTGCACCTGATGGTCAGCCCGGTTGACCGCATCATTGGCGACTTTATCGAAGCGGGTGCCACCTACATCACCTTCCACCCTGAAGCCACgcag      >  7:170202/1‑148 (MQ=255)
                                                                                                                                         ggCCCGATGGTTTGCACCGCTCTGCGCAAGTACGGCGTTACTGCGCCGATCGACGTGCACCTGATGGTCAGCCCGGTTGACCGCATCATTGGCGACTTTATCGAAGCGGGTGCCACCTACATCACCTTCCACCCTGAAGCCACgcagca    >  6:70870/1‑149 (MQ=255)
                                                                                                                                           cccGATGGTTTGCACCGCTCTGCGCAAGTACGGCGTTACTGCGCCGATCGACGTGCACCTGATGGTCAGCCCGGTTGACCGCATCATTGGCGACTTTATCGAAGCGGGTGCCACCTACATCACCTTCCACCCTGAAGCCACGCAGcaca  <  3:239801/149‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                                    
CCTACGCTATTGCCCCCTCCATTCTCTCTGCCGATTTCGCCCGCCTGGGTGAGGACGTCGACAAGGTATTGGCTGCGGGTGCAGACATCGTCCACTTCGATGTCATGGACAATCACTATGTCCCCAACCTGACCATCGGCCCGATGGTTTGCACCGCTCTGCGCAAGTACGGCGTTACTGCGCCGATCGACGTGCACCTGATGGTCAGCCCGGTTGACCGCATCATTGGCGACTTTATCGAAGCGGGTGCCACCTACATCACCTTCCACCCTGAAGCCACGCAGCACA  >  NC_002947/503955‑504242

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: