Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,665,250 G→C 12.1% R206R (CGC→CGG)  PP_1459 ← CBS domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,665,2500GC12.1% 179.9 / 15.8 66R206R (CGC→CGG) PP_1459CBS domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base G (28/30);  new base C (4/4);  total (32/34)
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.07e-01

GGCGCTGTGCAGCGCTTCCAGGGTCAGCTCGGCCTTGGGCAGCAGGTGGCTGATCAGCTTGGTGTTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAGCAGACCCTGGCGAGGCGCTTCGTGCTGTTCATCCTCATCGAAGTCGTC  >  NC_002947/1665103‑1665396
                                                                                                                                                   |                                                                                                                                                  
ggCGCTGTGCAGCGCTTCCAGGGTCAGCTCGGCCTTGGGCAGCAGGTGGCTGATCAGCTTGGTGTTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGc                                                                                                                                                   <  3:3589/149‑1 (MQ=255)
  cgcTGTGCAGCGCTTCCAGGGTCAGCTCGGCCTTGGGCAGCAGGTGGCTGATCAGCTTGGTGTTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGt                                                                                                                                                 <  4:231927/149‑1 (MQ=255)
  cgcTGTGCAGCGCTTCCAGGGTCAGCTCGGCCTTGGGCAGCAGGTGGCTGATCAGCTTGGTGTTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGt                                                                                                                                                 <  2:274374/149‑1 (MQ=255)
          aGCGCTTCCAGGGTCAGCTCGGCCTTGGGCAGCAGGTGGCTGATCAGCTTGGTGTTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGAGCCGATCGAACg                                                                                                                                          >  4:291057/1‑148 (MQ=255)
                 ccAGGGTCAGCTCGGCCTTGGGCAGCAGGTGGCTGATCAGCTTGGTGTTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCtcgtc                                                                                                                                   >  4:104763/1‑148 (MQ=255)
                  cAGGGTCAGCTCGGCCTTGGGCAGCAGGTGGCTGATCAGCTTGGTGTTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGAGCCGATCGAACGCATCGTCGa                                                                                                                                 <  3:291057/149‑1 (MQ=255)
                  cAGGGTCAGCTCGGCCTTGGGCAGCAGGTGGCTGATCAGCTTGGTGTTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGAGCCGATCGAACGCATCGTCGa                                                                                                                                 <  6:146076/149‑1 (MQ=255)
                   aGGGTCAGCTCGGCCTTGGGCAGCAGGTGGCTGATCATCTTGGTGTTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGCCGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAg                                                                                                                                >  1:32161/1‑149 (MQ=255)
                    gggTCAGCTCGGCCTTGGGCAGCAGGTGGCTGATCAGCTTGGTGTTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGAGCCGATCGAACGCATCGTCGAgc                                                                                                                               <  3:157555/149‑2 (MQ=255)
                     ggTCAGCTCGGCCTTGGGCAGCAGGTGGCTGATCAGCTTGGTGTTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGt                                                                                                                              <  5:232181/149‑1 (MQ=255)
                         aGCTCGGCCTTGGGCAGCAGGTGGCTGATCAGCTTGGTGTTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGAt                                                                                                                          <  2:50728/149‑1 (MQ=255)
                         aGCTCGGCCTTGGGCAGCAGGTGGCTGATCAGCTTGGTGTTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGAGCCGATCGAACGCATCGTCGAgcagctg                                                                                                                          <  6:67885/149‑7 (MQ=38)
                           cTCGGCCTTGGGCAGCAGGTGGCTGATCAGCTTGGTGTTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGc                                                                                                                        >  5:282505/1‑149 (MQ=255)
                               gCCTTGGGCAGCAGGTGGCTGATCAGCTTGGTGTTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATc                                                                                                                    <  6:282505/149‑1 (MQ=255)
                                ccTTGGGCAGCAGGTGGCTGATCAGCTTGGTGTTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATcc                                                                                                                   <  6:44849/149‑1 (MQ=255)
                                            ggTGGCTGATCAGCTTGGTGTTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTgcgc                                                                                                       <  3:104763/149‑1 (MQ=255)
                                                 cTGATCAGCTTGGTGTGGAGGATCGCTTCGACCTGGTTGATGGCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCAc                                                                                                  <  7:45490/149‑1 (MQ=255)
                                                   gATCAGCTTGGTGTTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCAc                                                                                                  >  2:54419/1‑147 (MQ=255)
                                                         cTTGGTGTTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGt                                                                                          >  7:276453/1‑149 (MQ=255)
                                                            ggTGTTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGt                                                                                       >  7:252983/1‑149 (MQ=255)
                                                             gtgtTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACTATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGtt                                                                                      <  8:276453/149‑1 (MQ=255)
                                                                  gAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACgg                                                                                 >  6:174856/1‑149 (MQ=255)
                                                                   aGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGt                                                                                >  7:174271/1‑149 (MQ=255)
                                                                         tcTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCtt                                                                          <  7:195385/148‑1 (MQ=255)
                                                                           ttCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCtt                                                                          >  8:195385/1‑147 (MQ=255)
                                                                               aCCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAg                                                                    >  2:230446/1‑149 (MQ=255)
                                                                                ccTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGc                                                                   <  3:284223/149‑1 (MQ=255)
                                                                                 cTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGc                                                                   >  7:21979/1‑148 (MQ=255)
                                                                                 cTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGc                                                                   <  8:21979/148‑1 (MQ=255)
                                                                                 cTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCg                                                                  >  4:229311/1‑149 (MQ=255)
                                                                                    gTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGAt                                                                             >  5:133374/1‑135 (MQ=255)
                                                                                    gTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGAt                                                                             <  6:133374/135‑1 (MQ=255)
                                                                                          gTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTg                                                                               <  5:147389/127‑1 (MQ=255)
                                                                                          gTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTg                                                                               >  6:147389/1‑127 (MQ=255)
                                                                                            cGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCGACTGCAGGcc                                                        <  8:174271/148‑1 (MQ=255)
                                                                                             gTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGt                                                      <  8:252983/149‑1 (MQ=255)
                                                                                                gtgGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGt                                                                        <  3:47644/128‑1 (MQ=255)
                                                                                                gtgGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGt                                                                        >  4:47644/1‑128 (MQ=255)
                                                                                                   gTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAg                                                  >  7:286817/1‑147 (MQ=255)
                                                                                                         cGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCAt                                                                                                                     <  7:9693/74‑1 (MQ=255)
                                                                                                         cGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCAt                                                                                                                     >  8:9693/1‑74 (MQ=255)
                                                                                                           gCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAGCAGACCCTgg                                        >  7:262231/1‑149 (MQ=255)
                                                                                                                cGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAGCAGACCCTTGCGAgg                                   <  5:174856/149‑1 (MQ=255)
                                                                                                                  ggTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAGCAGACCCTGGCGAGgcg                                 <  1:230446/149‑1 (MQ=255)
                                                                                                                   gtgtGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAGCAGACCCTGGCGAGgcgc                                >  6:35232/1‑149 (MQ=255)
                                                                                                                      tgGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAGCAGACCCTGGCGAGGCGCTTc                             <  5:253692/149‑1 (MQ=255)
                                                                                                                      tgGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAGCAGACCCTGGCGAGGCGCTTc                             <  1:54419/149‑1 (MQ=255)
                                                                                                                           gcTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAGCAGACCCTGGc                                       >  6:217222/1‑134 (MQ=255)
                                                                                                                           gcTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAGCAGACCCTGGc                                       <  5:217222/134‑1 (MQ=255)
                                                                                                                            cTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCTTCGAGCAGACCCTGGCGAGGCGCTTCGTGCTg                       <  8:286817/149‑1 (MQ=255)
                                                                                                                              gACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAGCAGACCCTGGCGAGGCGCTTCGTGCTGt                      >  3:107649/1‑148 (MQ=255)
                                                                                                                                   tCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAGCAGACCCTGGCGAGGCGCTTCGTGCTGTTCATcc                >  2:160388/1‑149 (MQ=255)
                                                                                                                                      gctgctCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAGCAGACCCTGGCGAGGCGCTTCGTg                          >  7:26785/1‑136 (MQ=255)
                                                                                                                                      gctgctCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAGCAGACCCTGGCGAGGCGCTTCGTg                          <  8:26785/136‑1 (MQ=255)
                                                                                                                                      gctgctCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAGCAGACCCTGGCGAGGCGCTTCGTGCTGTTCATCCTCa             >  8:298609/1‑149 (MQ=255)
                                                                                                                                          ctCGACGATGAGTTCGATCTGCTCGTCGAGGTTGATGCCATCCTCTTCGTTGCGCCTCACCAGGATGGCGTTGTCGGTTATCTTGTAGAGCGTGTGCAGGCCGTAGAGCATACCCTGGCGAGGCGCTTCGTGCTGTTCATCCTCATc           <  2:32161/147‑1 (MQ=255)
                                                                                                                                          ctCGACGAGCCGATCGAACGCATCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAGCAGACCCTGGCGAGGCGCTTCGTGCTGTTCATCCTCATCg          >  1:12096/1‑148 (MQ=255)
                                                                                                                                           tCGACGAGCCGCTCGAACGCATCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGAACGAGATCTAGTCCAGCGACTCCAGGCCGTCTAGCAGACCCTGGCGAGGAGCTTCGTGCTGTTCATCGTCAACGa         >  2:252867/1‑148 (MQ=255)
                                                                                                                                           tCGACGAGCCGATCGAACGCATCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAGCAGACCCTGGCGAGGCGCTTCGTGCTGTTCCTCCTCATCg          >  8:299491/1‑147 (MQ=255)
                                                                                                                                            cgacgaTGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAGCAGACCCTGGCGAGGCGCTTCGTGCTGTTCATCCTCATCGAAg       >  8:54545/1‑149 (MQ=255)
                                                                                                                                            cgacgaTGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAGCAGACCCTGGCGAGGCGCTTCGTGCTGTTCATCCTCATCGAAg       >  5:66377/1‑149 (MQ=255)
                                                                                                                                            cgacgaTGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAGCAGACCCTGGCGAGGCGCTTCGTGCTGTTCATCCTCATCGAAg       >  6:288778/1‑149 (MQ=255)
                                                                                                                                             gacgaTGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAGCAGACCCTGGCGAGGCGCTTCGTGCTGTTCATCCTCATCGAAg       <  5:288778/148‑1 (MQ=255)
                                                                                                                                                 aTGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAGCAGACCCTGGCGAGGCGCTTCGTGCTGTTCATCCTCATCGAAgtcgtc  <  1:124547/149‑1 (MQ=255)
                                                                                                                                                  tGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAGCAGACCCTGGCGAGgc                                  <  6:305635/116‑1 (MQ=255)
                                                                                                                                                  tGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAGCAGACCCTGGCGAGgc                                  >  5:305635/1‑116 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GGCGCTGTGCAGCGCTTCCAGGGTCAGCTCGGCCTTGGGCAGCAGGTGGCTGATCAGCTTGGTGTTGAGGATCGCTTCGACCTGGTTGATGTCGTTGTGGTAGACCGGCAGGCGGGTGTGGCGGCTGACGATCAGCTGCTCGACGATGCGTTCGATCGGCTCGTCGAGGTTGATGCCATCCACTTCGTTGCGCGGCACCAGGATGTCGTTGACGGTGATCTTGTCCAGCGACTGCAGGCCGTCGAGCAGACCCTGGCGAGGCGCTTCGTGCTGTTCATCCTCATCGAAGTCGTC  >  NC_002947/1665103‑1665396

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