Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,710,161 T→A 9.2% A260A (GCT→GCA ppc → phosphoenolpyruvate carboxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,710,1610TA 9.2% 173.8 / 6.1 62A260A (GCT→GCAppcphosphoenolpyruvate carboxylase
Reads supporting (aligned to +/- strand):  ref base T (30/26);  new base A (2/4);  total (32/30)
Fisher's exact test for biased strand distribution p-value = 4.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

ACTGTGGCATGCCATCCCCAGCCACCTGCGCAAGGTCGACAAGGCCCTGCTGGAGGCCACTGGCCTGCGCCTGCCGCTGGAGGCCGCGCCGATCCGCTTTGCCTCGTGGATGGGTGGGGACCGTGATGGCAACCCCAATGTGACGGCTGCCGTCACCCGCGAAGTACTCCTGCTGGCGCGCTGGATGGCCGCCGACCTGTTCCTGCGCGACATCGATGCACTGGCCGCCGAGTTGTCCATGCAGCAGGCCAACGACACCCTGCGCAAGCAGGTTGGCGACAGTGCCGAACCCTACC  >  NC_002947/1710014‑1710309
                                                                                                                                                   |                                                                                                                                                      
aCTGTGGCATGCCATCCCCAGCCACCTGCGCAAGGTCGACAAGGCCCTGCTGGAGGCCACTGGCCTGCGCCTGCCGCTGGAGGCCGCGCCGATCCGCTTTGCCTCGTGGATGGGTGGGGACCGTGATGGCAACCCCAATGTGACGGCTg                                                                                                                                                       >  3:75389/1‑149 (MQ=255)
    tgGCATGCCATCCCCAGCCACCTGCGCAAGGTCGACAAGGCCCTGCTGGAGGCCACTGGCCTGCGCCTGCCGCTGGAGGCCGCGCCGATCCGCTTTGCCTCGTGGATGGGTGGGGACCGTGATGGCAACCCCAATGTGACGGCTGCCGt                                                                                                                                                   >  8:9032/1‑149 (MQ=255)
        aTGCCATCCCCAGCCACCTGCGCAAGGTCGACAAGGCCCTGCTGGAGGCCACTGGCCTGCGCCTGCCGCTGGAGGCCGCGCCGATCCGCTTTGCCTCGTGGATGGGTGGGGACCGTGATGGCAACCCCAATGTGACGGCAGCCGTCAca                                                                                                                                               >  4:35387/1‑148 (MQ=255)
           ccATCCCCAGCCACCTGCGCAAGGTCGACAAGGCCCTGCTGGAGGCCACTGGCCTGCGCCTGCCGCTGGAGGCCGCGCCGATCCGCTTTGCCTCGTGGATGGGTGGGGACCGTGATGGCAACCCCAATGTGACGGCAGCCGTCAcattg                                                                                                                                            <  3:35387/149‑5 (MQ=255)
           ccATCCCCAGCCACCTGCGCAAGGTCGACAAGGCCCTGCTGGAGGCCACTGGCCTGCGCCTGCCGCTGGAGGCCGCGCCGATCCGCTTTGCCTCGTGGATGGGTGGGGACCGTGATGGCAACCCCAATGTGACGGCAGCCGTCAcatt                                                                                                                                             <  5:132199/148‑4 (MQ=255)
             aTCCCCAGCCACCTGCGCAAGGTCGACAAGGCCCTGCTGGAGGCCACTGGCCTGCGCCTGCCGCTGGAGGCCGCGCCGATCCGCTTTGCCTCGTGGATGGGTGGGGACCGTGATGGCAACCCCAATGTGACGGCAGCCGTCAcattggg                                                                                                                                          <  8:30284/149‑7 (MQ=255)
                 ccAGCCACCTGCGCAAGGTCGACAAGGCCCTGCTGGAGGCCACTGGCCTGCGCCTGCCGCTGGAGGCCGCGCCGATCCGCTTTGCCTCGTGGATGGGTGGGGACCGTGATGGCAACCCCAATGTGACGGCTGCCGTCACCCGCGAAGTa                                                                                                                                      >  8:282485/1‑149 (MQ=255)
                  cAGCCACCTGCGCAAGGTCGACAAGGCCCTGCTGGAGGCCACTGGCCTGCGCCTGCCGCTGGAGGCCGCGCCGATCCGCTTTGCCTCGTGGATGGGTGGGGACCGTGATGGCAACCCCAATGTGACGGCTGCCGTCACCCGCGAAGTa                                                                                                                                      >  8:59210/1‑148 (MQ=255)
                  cAGCCACCTGCGCAAGGTCGACAAGGCCCTGCTGGAGGCCACTGGCCTGCGCCTGCCGCTGGAGGCCGCGCCGATCCGCTTTGCCTCGTGGATGGGTGGGGACCGTGATGGCAACCCCAATGTGACGGCTGCCGTCACCCGCGAAGTAc                                                                                                                                     >  5:16355/1‑149 (MQ=255)
                  cAGCCACCTGCGCAAGGTCGACAAGGCCCTGCTGGAGGCCACTGGCCTGCGCCTGCCGCTGGAGGCCGCGCCGATCCGCTTTGCCTCGTGGATGGGTGGGGACCGTGATGGCAACCCCAATGTGACGGCAGCCGTCAcattg                                                                                                                                            >  4:155098/1‑138 (MQ=255)
                  cAGCCACCTGCGCAAGGTCGACAAGGCCCTGCTGGAGGCCACTGGCCTGCGCCTGCCGCTGGAGGCCGCGCCGATCCGCTTTGCCTCGTGGATGGGTGGGGACCGTGATGGCAACCCCAATGTGACGGCAGCCGTCAcattg                                                                                                                                            <  3:155098/142‑5 (MQ=255)
                      cACCTGCGCAAGGTCGACAAGGCCCTGCTGGAGGCCACTGGCCTGCGCCTGCCGCTGGAGGCCGCGCCGATCCGCTTTGCCTCGTGGATGGGTGGGGACCGTGATGGCAACCCCAATGTGACGGCTGCCGTCACCCGCGAAGTACTCct                                                                                                                                 <  3:70936/149‑1 (MQ=255)
                        ccTGCGCAAGGTCGACAAGGCCCTGCTGGAGGCCACTGGCCTGCGCCTGCCGCTGGAGGCCGCGCCGATCCGCTTTGCCTCGTGGATGGGTGGGGACCGTGATGGCAACCCCAATGTGACGGCTGCCGTCACCCGCGAAGTACTCctgc                                                                                                                               >  1:29901/1‑149 (MQ=255)
                           gcgcAAGGTCGACAAGGCCCTGCTGGAGGCCACTGGCCTGCGCCTGCCGCTGGAGGCCGCGCCGATCCGCTTTGCCTCGTGGATGGGTGGGGACCGTGATGGCAACCCCAATGTGACGGCTGCCGTCACCCGCGAAGTACTCCTGCTgg                                                                                                                            >  8:362931/1‑149 (MQ=255)
                                     gACAAGGCCCTGCTGGAGGCCACTGGCCTGCGCCTGCCGCTGGAGGCCGCGCCGATCCGCTTTGCCTCGTGGATGGGTGGGGACCGTGATGGCAACCCCAATGTGACGGCTGCCGTCACCCGCGAAGTACTCCTGCTGGCGCGCTGGa                                                                                                                   <  5:128624/148‑1 (MQ=255)
                                     gACAAGGCCCTGCTGGAGGCCACTGGCCTGCGCCTGCCGCTGGAGGCCGCGCCGATCCGCTTTGCCTCGTGGATGGGTGGGGACCGTGATGGCAACCCCAATGTGACGGCTGCCGTCACCCGCGAAGTACTCCTGCTGGCGCGCTGGAt                                                                                                                  <  6:249046/149‑1 (MQ=255)
                                      aCAAGGCCCTGCTGGAGGCCACTGGCCTGCGCCTGCCGCTGGAGGCCGCGCCGATCCGCTTTGCCTCGTGGATGGGTGGGGACCGTGATGGCAACCCCAATGTGACGGCTGCCGTCACCCGCGAAGTACTCCTGCTGGCGCGCTGGATg                                                                                                                 <  7:9032/149‑1 (MQ=255)
                                          ggCCCTGCTGGAGGCCACTGGCCTGCGCCTGCCGCTGGAGGCCGCGCCGATCCGCTTTGCCTCGTGGATGGGTGGGGACCGTGATGGCAACCCCAATGTGACGGCTGCCGTCACCCGCGAAGTACTCCTGCTGGCGCGCTGGATGgccg                                                                                                             >  2:39874/1‑149 (MQ=255)
                                                 ctgGAGGCCACTGGCCTGCGCCTGCCGCTGGAGGCCGCGCCGATCCGCTTTGCCTCGTGGATGGGTGGGGACCGTGATGGCAACCCCAATGTGACGGCTGCCGTCACCCGCGAAGTACTCCTGCTGGCGCGCTGGATGGCCGCCGAcc                                                                                                       <  4:96612/148‑1 (MQ=255)
                                                 ctgGAGGCCACTGGCCTGCGCCTGCCGCTGGAGGCCGCGCCGATCCGCTTTGCCTCGTGGATGGGTGGGGACCGTGATGGCAACCCCAATGTGACGGCTGCCGTCACCCGCGAAGTACTCCTGCTGGCGCGCTGGATGGCCGCCGAc                                                                                                        <  7:362931/147‑1 (MQ=255)
                                                 ctgGAGGCCACTGGCCTGCGCCTGCCGCTGGAGGCCGCGCCGATCCGCTTTGCCTCGTGGATGGGTGGGGACCGTGATGGCAACCCCAATGTGACGGCTGCCGTCACCCGCGAAGTACTCCTGCTGGCGCGCTGGATGGCCGCCGAc                                                                                                        <  2:29901/147‑1 (MQ=255)
                                                 ctgGAGGCCACTGGCCTGCGCCTGCCGCTGGAGGCCGCGCCGATCCGCTTTGCCTCGTGGATGGGTGGGGACCGTGATGGCAACCCCAATGTGACGGCTGCCGTCACCCGCGAAGTACTCCTGCTGGCGCGCTGGATGGCCGCCGAc                                                                                                        <  3:327846/147‑1 (MQ=255)
                                                     aGGCCACTGGCCTGCGCCTGCCGCTGGAGGCCGCGCCGATCCGCTTTGCCTCGTGGATGGGTGGGGACCGTGATGGCAACCCCAATGTGACGGCTGCCGTCACCCGCGAAGTACTCCTGCTGGCGCGCTGGATGGCCGCCGACCTGTTc                                                                                                 <  1:165460/149‑1 (MQ=255)
                                                           cTGGCCTGCGCCTGCCGCTGGAGGCCGCGCCGATCCGCTTTGCCTCGTGGATGGGTGGGGACCGTGATGGCAACCCCAATGTGACGGCTGCCGTCACCCGCGAAGTACTCCTGCTGGCGCGCTGGATGGCCGCCGACCTGTTCCTgcgc                                                                                           <  4:75389/149‑1 (MQ=255)
                                                               cctgcgcctgcCGCTGGAGGCCGCGCCGATCCGCTTTGCCTCGTGGATGGGTGGGGACCGTGATGGCAACCCCAATGTGACGGCTGCCGTCACCCGCGAAGTACTCCTGCTGGCGCGCTGGATGGCCGCCGACCTGTTCCTGCGCGACa                                                                                       >  6:47337/1‑149 (MQ=255)
                                                                ctgcgcctgcCGCTGGAGGCCGCGCCGATCCGCTTTGCCTCGTGGATGGGTGGGGACCGTGATGGCAACCCCAATGTGACGGCTGCCGTCACCCGGGAAGTACTCCTGCTGGCGCGCTGGATGGCCGCCGAc                                                                                                        >  1:333784/1‑132 (MQ=255)
                                                                ctgcgcctgcCGCTGGAGGCCGCGCCGATCCGCTTTGCCTCGTGGATGGGTGGGGACCGTGATGGCAACCCCAATGTGACGGCTGCCGTCACCCGCGAAGTACTCCTGCTGGCGCGCTGGATGGCCGCCGAc                                                                                                        <  2:333784/132‑1 (MQ=255)
                                                                     cctgcCGCTGGAGGCCGCGCCGATCCGCTTTGCCTCGTGGATGGGTGGGGACCGTGATGGCAACCCCAATGTGACGGCTGCCGTCACCCGCGAAGTACTCCTGCTGGCGCGCTGGATGGCCGCCGACCTGTTCCTGCGCGACATCGATg                                                                                 >  1:218105/1‑149 (MQ=255)
                                                                     cctgcCGCTGGAGGCCGCGCCGATCCGCTTTGCCTCGTGGATGGGTGGGGACCGTGATGGCAACCCCAATGTGACGGCTGCCGTCACCCGCGAAGTACTCCTGCTGGCGCGCTGGATGGCCGCCGACCTGTTCCTGCGCGACATCGATg                                                                                 <  3:119506/149‑1 (MQ=255)
                                                                      ctgcCGCTGGAGGCCGCGCCGATCCGCTTTGCCTCGTGGATGGGTGGGGACCGTGATGGCAACCCCAATGTGACGGCTGCCGTCACCCGCGAAGTACTCCTGCTGGCGCGCTGGATGGCCGCCGACCTGTTCCTGCGCGACATCGATg                                                                                 >  1:68143/1‑148 (MQ=255)
                                                                      ctgcCGCTGGAGGCCGCGCCGATCCGCTTTGCCTCGTGGATGGGTGGGGACCGTGATGGCAACCCCAATGTGACGGCTGCCGTCACCCGCGAAGTACTCCTGCTGGCGCGCTGGATGGCCGCCGACCTGTTCCTGCGCGACATCGATg                                                                                 >  3:45335/1‑148 (MQ=255)
                                                                                  gCCGCGCCGATCCGCTTTGCCTCGTGGATGGGTGGGGACCGTGATGGCAACCCCAATGTGACGGCTGCCGTCACCCGCGAAGTACTCCTGCTGGCGCGCTGGATGGCCGCCGACCTGTTCCTGCGCGACATCGATGCACTGGCCGCCGa                                                                    <  1:39874/149‑1 (MQ=255)
                                                                                   ccGCGCCGATCCGCTTTGCCTCGTGGATGGGTGGGGACCGTGATGGCAACCCCAATGTGACGGCTGCCGTCACCCGCGAAGTACTCCTGCTGGCGCGCTGGATGGCCGCCGACCTGTTCCTGCGCGACATCGATGCACTGGCCGCCGAg                                                                   >  6:402015/1‑149 (MQ=255)
                                                                                       gcCGATCCGCTTTGCCTCGTGGATGGGTGGGGACCGTGATGGCAACCCCAATGTGACGGCTGCCGTCACCCGCGAAGTACTCCTGCTGGCGCGCTGGATGGCCGCCGACCTGTTCCTGCGCGACATCGATGCACTg                                                                            <  2:140304/136‑1 (MQ=255)
                                                                                       gcCGATCCGCTTTGCCTCGTGGATGGGTGGGGACCGTGATGGCAACCCCAATGTGACGGCTGCCGTCACCCGCGAAGTACTCCTGCTGGCGCGCTGGATGGCCGCCGACCTGTTCCTGCGCGACATCGATGCACTg                                                                            >  1:140304/1‑136 (MQ=255)
                                                                                          gATCCGCTTTGCCTCGTGGATGGGTGGGGACCGTGATGGCAACCCCAATGTGACGGCTGCCGTCACCCGCGAAGTACTCCTGCTGGCGCGCTGGATGGCCGCCGACCTGTTCCTGCGCGACATCGATGCACTGGCCGCCGAGTTGTCCa                                                            <  7:59210/149‑1 (MQ=255)
                                                                                                 tttGCCTCGTGGATGGGTGGGGACCGTGATGGCAACCCCAATGTGACGGCTGCCGTCACCCGCGAAGTACTCCTGCTGGCGCGCTGGATGGCCGCCGACCTGTTCCTGCGCGACATCGATGCACTGGCCGCCGAGTTGTCCATgcagca                                                     <  4:45335/149‑1 (MQ=255)
                                                                                                             aTGGGTGGGGACCGTGATGGCAACCCCAATGTGACGGCTGCCGTCACCCGCGAAGTACTCCTGCTGGCGCGCTGGATGGCCGCCGACCTGTTCCTGCGCGACATCGATGCACTGGCCGCCGAGTTGTCCATGCAGCAGGCCAACGacac                                         >  1:285988/1‑149 (MQ=255)
                                                                                                             aTGGGTGGGGACCGTGATGGCAACCCCAATGTGACGGCTGCCGTCACCCGCGAAGTACTCCTGCTGGCGCGCTGGATGGCCGCCGACCTGTTCCTGCGCGACATCGATGCACTGGCCGCCGAGTTGTCCATGCAGCAGGCCAACGacac                                         >  1:280352/1‑149 (MQ=255)
                                                                                                                 gtgggGACCGTGATGGCAACCCCAATGTGACGGCTGCCGTCACCCGCGAAGTACTCCTGCTGGCGCGCTGGATGGCCGCCGACCTGTTCCTGCGCGACATCGATGCACTGGCCGCCGAGTTGTCCATGCAGCAGGCCAACGACACCCTg                                    >  8:188268/1‑149 (MQ=255)
                                                                                                                  tgggGACCGTGATGGCAACCCCAATGTGACGGCTGCCGTCACCCGCGAAGTACTCCTGCTGGCGCGCTGGATGGCCGCCGACCTGTTCCTGCGCGACATCGATGCACTGGCCGCCGAGTTGTCCATGCAGCAGGCCAACGACACCCTgc                                   >  7:32244/1‑149 (MQ=255)
                                                                                                                   ggggACCGTGATGGCAACCCCAATGTGACGGCTGCCGTCACCCGCGAAGTACTCCTGCTGGCGCGCTGGATGGCCGCCGACCTGTTCCTGCGCGACATCGATGCACTGGCCGCCGAGTTGTCCATGCAGCAGGCCAACGACACCCTgcg                                  >  6:343235/1‑149 (MQ=255)
                                                                                                                      gACCGTGATGGCAACCCCAATGTGACGGCTGCCGTCACCCGCGAAGTACTCCTGCTGGCGCGCTGGATGGCCGCCGACCTGTTCCTGCGCGACATCGATGCACTGGCCGCCGAGTTGTCCATGCAGCAGGCCAACGACACCCTgcgc                                 >  1:201804/1‑147 (MQ=255)
                                                                                                                          gTGATGGCAACCCCAATGTGACGGCTGCCGTCCCCCGCGAAGTACTCCTGCTGGCGCGCTGGATGGCCGCCGACCTCTTCCTGCGCGACATCGATGCACTGGCCGCCGAGATGTCCATGCAGCAGGCCAACGACACCCTGCGCAAGCAg                           >  7:212618/1‑149 (MQ=255)
                                                                                                                           tGATGGCAACCCCAATGTGACGGCTGCCGTCACCCGCGAAGTACTCCTGCTGGCGCGCTGGATGGCCGCCGACCTGTTCCTGAGCGACATCGATGCACTGGCAGCCGAGTTGTCCATGCAGCAGGCCAACGACACCCTGCGCAAGCAgg                          >  7:256013/1‑149 (MQ=255)
                                                                                                                            gATGGCAACCCCAATGTGACGGCTGCCGTCACCCGCGAAGTACTCCTGCTGGCGCGCTGGATGGCCGCCGACCTGTTCCTGCGCGACATCGATGCACTGGCCGCCGAGTTGTCCATGCAGCAGGCCAACGACACCCTgcgc                                 <  8:156966/141‑1 (MQ=255)
                                                                                                                            gATGGCAACCCCAATGTGACGGCTGCCGTCACCCGCGAAGTACTCCTGCTGGCGCGCTGGATGGCCGCCGACCTGTTCCTGCGCGACATCGATGCACTGGCCGCCGAGTTGTCCATGCAGCAGGCCAACGACACCCTgcgc                                 >  7:156966/1‑141 (MQ=255)
                                                                                                                             aTGGCAACCCCAATGTGACGGCTGCCGTCACCCGCGAAGTACTCCTGCTGGCGCGCTGGAt                                                                                                                  >  2:347902/1‑61 (MQ=255)
                                                                                                                             aTGGCAACCCCAATGTGACGGCTGCCGTCACCCGCGAAGTACTCCTGCTGGCGCGCTGGAt                                                                                                                  <  1:347902/61‑1 (MQ=255)
                                                                                                                              tGGCAACCCCAATGTGACGGCTGCCGTCACCCGCGAAGTACTCCTGCTGGCGCGCTGGATGGCCGCCGACCTGTTCCTGCGCGACATCGATGCACTGGCCGCCGAGTTGTCCATGCAGCAGGCCAACGACACCCTGCGCAAGCAGGTTg                       <  6:132156/149‑1 (MQ=255)
                                                                                                                                 cAACCCCAATGTGACGGCTGCCGTCACCCGCGAAGTACTCCTGCTGGCGCGCTGGATGGCCGCCGACCTGTTCCTGCGCGACATCGATGCACTGGCCGCCGAGTTGTCCATGCAGCAGGCCAACGACACCCTGCGCAAGCAGGTTGGc                     >  3:1487/1‑148 (MQ=255)
                                                                                                                                    ccccAATGTGACGGCTGCCGTCACCCGCGAAGTACTCCTGCTGGCGCGCTGGATGGCCGCCGACCTGTTCCTGCGCGACATCGATGCACTGGCCGCCGAGTTGTCCATGCAGCAGGCCAACGACACCCTGCGCAAGCAGGTTGGCGACa                 <  2:52988/149‑1 (MQ=255)
                                                                                                                                     cccAATGTGACGGCTGCCGTCACCCGCGAAGTACTCCTGCTGGCGCGCTGGATGGCCGCCGACCTGTTCCTGCGCGACATCGATGCTCTGGCCGCCGAGTTGTCCATGCAGCAGGCCAACGACACCCTGCGCAAGCAGGTTg                       >  7:325290/1‑142 (MQ=255)
                                                                                                                                     cccAATGTGACGGCTGCCGTCACCCGCGAAGTACTCCTGCTGGCGCGCTGGATGGCCGCCGACCTGTTCCTGCGCGACATCGATGCTCTGGCCGCCGAGTTGTCCATGCAGCAGGCCAACGACACCCTGCGCAAGCAGGTTg                       <  8:325290/142‑1 (MQ=255)
                                                                                                                                     cccAATGTGACGGCTGCCGTCACCCGCGAAGTACTCCTGCTGGCGCGCTGGATGGCCGCCGACCTGTTCCTGCGCGACATCGATGCACTGGCCGCCGAGTTGTCCATGCAGCAGGCCAACGACACCCTGCGCAAGCAGGTTGGCGACAg                >  4:288391/1‑149 (MQ=255)
                                                                                                                                       cAATGTGACGGCTGCCGTCACCCGCGAAGTACTCCTGCTGGCGCGCTGGATGGCCGCCGACCAGGTTACTGCGCGACATCGATGCACTGGCCGCCGAGTTGTACATGCCGCACGCCAAAGCCACCCAGCGCAAGCCGGTTGGCGACact               >  1:320482/1‑147 (MQ=255)
                                                                                                                                       cAATGTGACGGCTGCCGTCACCCGCGAAGTACTCCTGCTGGCGCGCTGGATGGCCGCCGACCTGTTCCTGCGCGACATCGATGCACTGGCCGCCGa                                                                    <  7:143648/96‑1 (MQ=255)
                                                                                                                                       cAATGTGACGGCTGCCGTCACCCGCGAAGTACTCCTGCTGGCGCGCTGGATGGCCGCCGACCTGTTCCTGCGCCACATCGATGCACTGGCCGCCGa                                                                    >  8:143648/1‑96 (MQ=255)
                                                                                                                                               cGGCTGCCGTCACCCGCGAAGTACTCCTGCTGGCGCGCTGGATGGCCGCCGACCTGTTCCTGCGCGACATCGATGCACTGGCCGCCGAGTTGTCCATGCAGCAGGCCAACGACACCCTTGCGCAAGCAGGTTGGCGACAGTGCCGAAcc       <  1:35932/149‑1 (MQ=255)
                                                                                                                                               cGGCTGCCGTCACCCGCGAAGTACTCCTGCTGGCGCGCTGGATGGCCGCCGACCTGTTCCTGCGCGACATCGATGCACTGGCCGCCGAGTTGTCCATGCAGCAGGCCAACGACACCCTGCGCAAGCAGGTTGGCGACAGTGCCGAAccc      <  3:385991/149‑1 (MQ=255)
                                                                                                                                                 gCTGCCGTCACCCGCGAAGTACTCCTGCTGGCGCGCTGGATGGCCGCCGACCTGTTCCTGCGCGACATCGATGCACTGGCCGCCGAGTTGTCCATGCAGCAGGCCAACGACACCCTGCGCAAGCAGGTTGGCGACAGTGCCGAACCCTa    <  2:280352/149‑1 (MQ=255)
                                                                                                                                                 gCTGCCGTCACCCGCGAAGTACTCCTGCTGGCGCGCTGGATGGCCGCCGACCTGTTCCTGCGCGACATCGATGCACTGGCCGCCGAGTTGTCCATGCAGCAGGCCAACGACACCCTGCGCAAGCAGGTTGGCGACAGTGCCGAACCCTa    <  2:285988/149‑1 (MQ=255)
                                                                                                                                                   tGCCGTCACCCGCGAAGTACTCCTGCTGGCGCGCTGGATGGCCGCCGACCTGTTCCTGCGCGACATCGATGCACTGGCCGCCGAGTTGTCCATGCAGCAGGCCAACGACACCCTGCGCAAGCAGGTTGGCGACAGTGCCGAACCCTAcc  >  1:357681/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                      
ACTGTGGCATGCCATCCCCAGCCACCTGCGCAAGGTCGACAAGGCCCTGCTGGAGGCCACTGGCCTGCGCCTGCCGCTGGAGGCCGCGCCGATCCGCTTTGCCTCGTGGATGGGTGGGGACCGTGATGGCAACCCCAATGTGACGGCTGCCGTCACCCGCGAAGTACTCCTGCTGGCGCGCTGGATGGCCGCCGACCTGTTCCTGCGCGACATCGATGCACTGGCCGCCGAGTTGTCCATGCAGCAGGCCAACGACACCCTGCGCAAGCAGGTTGGCGACAGTGCCGAACCCTACC  >  NC_002947/1710014‑1710309

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: