Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,098,681 G→A 6.8% H295H (CAC→CAT)  kdsD ← D‑arabinose 5‑phosphate isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,098,6810GA 6.8% 172.9 / 2.4 59H295H (CAC→CAT) kdsDD‑arabinose 5‑phosphate isomerase
Reads supporting (aligned to +/- strand):  ref base G (28/27);  new base A (2/2);  total (30/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.08e-01

CCAGCTTGATGCCTTTGCCGCGTTGCATGAGGTCCTGGTTCATTCCATCGCTCCTTTACATTACGCCGGCGCGC‑AGCAGGTCATGCATGTTCAGGGCGCCGGTCGGGCGGTCTTCACGGTCCACCACCACCAGGGCGTTGATCTTGTGGT‑CTTCCATTATCTTCAGGGCTTCGGCGGCGAGCATGTCGGCACGCGCAGTCTTGCCATGAACGGTCATGACCTGGTCGATCAGGGTGGTGTGCACATCGATATTGCGGTCCAGGCTGCGGCGCAGGTCACCGTCGGTGAAGATCC  >  NC_002947/1098536‑1098829
                                                                                                                                                  |                                                                                                                                                     
ccAGCTTGATGCCTTTGCCGCGTTGCATGAGGTCCTGGTTCATTCCATCGCTCCTTTACATTACGCCGGCGCGC‑AGCAGGTCATGCATGTTCAGGGCGCCGGTCGGGCGGTCTTCACGGTCCACCACCACCAGGGCGTTGATCTTGTgg                                                                                                                                                    >  6:353997/1‑149 (MQ=255)
       gATGCCTTTGCCGCGTTGCATGAGGTCCTGGTTCATTCCATCGCTCCTTTACATTACGCCGGCGCGC‑AGCAGGTCATGCATGTTCAGGGCGCCGGTCGGGCGGTCTTCACGGTCCACCACCACCAGGGCGTTGATCTTGTGGT‑CTTcc                                                                                                                                             <  8:5140/148‑1 (MQ=255)
         tGCCTTTGCCGCGTTGCATGAGGTCCTGGTTCATTCCATCGCTCCTTTACATTACGCCGGCGCGC‑AGCAGGTCATGCATGTTCAGGGCGCCGGTCGGGCGGTCTTCACGGTCCACCACCACCAGGGCGTTGATCTTGTGGT‑CTTCCAtt                                                                                                                                          <  5:353997/149‑1 (MQ=255)
         tGCCTTTGCCGCGTTGCATGAGGTCCTGGTTCATTCCATCGCTCCTTTACATTACGCCGGCGCGC‑AGCAGGTCATGCATGTTCAGGGCGCCGGTCGGGCGGTCTTCACGGTCCACCACCACCAGGGCGTTGATCTTGTGGT‑CTTCCAtt                                                                                                                                          <  6:14571/149‑1 (MQ=255)
          gCCTTTGCCGCGTTGCATGAGGTCCTGGTTCATTCCATCGCTCCTTTACATTACGCCGGCGCGC‑AGCAGGTCATGCATGTTCAGGGCGCCGGTCGGGCGGTCTTCACGGTCCACCACCACCAGGGCGTTGATCTTGTGGT‑CTTCCATTa                                                                                                                                         >  4:164523/1‑149 (MQ=255)
          gCCTTTGCCGCGTTGCATGAGGTCCTGGTTCATTCCATCGCTCCTTTACATTACGCCGGCGCGC‑AGCAGGTCATGCATGTTCAGGGCGCCGGTCGGGCGGTCTTCACGGT‑‑‑CCACCACCAGGGCGTTGATCTtgtg                                                                                                                                                     >  1:319204/1‑135 (MQ=255)
          gCCTTTGCCGCGTTGCATGAGGTCCTGGTTCATTCCATCGCTCCTTTACATTACGCCGGCGCGC‑AGCAGGTCATGCATGTTCAGGGCGCCGGTCGGGCGGTCTTCACGGT‑‑‑CCACCACCAGGGCGTTGATCTtgtg                                                                                                                                                     <  2:319204/135‑1 (MQ=255)
             tttGCCGCGTTGCATGAGGTCCTGGTTCATTCCATCGCTCCTTTACATTACGCCGGCGCGC‑AGCAGGTCATGCATGTTCAGGGCGCCGGTCGGGCGGTCTTCACGGTCCACCACCACCAGGGCGTTGATCTTGTGGT‑CTTCCATTATCt                                                                                                                                      >  6:198896/1‑149 (MQ=255)
             tttGCCGCGTTGCATGAGGTCCTGGTTCATTCCATCGCTCCTTTACATTACGCCGGCGCGC‑AGCAGGTCATGCATGTTCAGGGCGCCGGTCGGGCGGTCTTCACGGTCCACCACCACCAGGGCGTTGATCTTGTGGT‑CTTCCATTATCt                                                                                                                                      >  1:115036/1‑149 (MQ=255)
                 ccGCGTTGCATGAGGTCCTGGTTCATTCCATCGCTCCTTTACATTACGCCGGCGCGC‑AGCAGGTCATGCATGTTCAGGGCGCCGGTCGGGCGGTCTTCACGGTCCACCACCACCAGGGCGTTGATCTTGTGGT‑CTTCCATTATCTTCAg                                                                                                                                  <  3:164523/149‑1 (MQ=255)
                        gCATGAGGTCCTGGTTCATTCCATCGCTCCTTTACATTACGCCGGCGCGC‑AGCAGGTCATGCATGTTCAGGGCGCCGGTCGGGCGGTCTTCACGGTCCACCACCACCAGGGCGTTGATCTTGTGGT‑CTTCCATTATCTTCAGGGCTTcg                                                                                                                           >  3:105138/1‑149 (MQ=255)
                           tGAGGTCCTGGTTCATTCCATCGCTCCTTTACATTACGCCGGCGCGC‑AGCAGGTCATGCATGTTCAGGGCGCCGGTCGGGCGGTCTTCACGGTCCACCACCACCAGGGCGTTGATCTTGTGGT‑CTTCCATTATCTTCAGGGCTTcggcg                                                                                                                        <  8:72836/149‑1 (MQ=255)
                           tGAGGTCCTGGTTCATTCCATCGCTCCTTTACATTACGCCGGCGCGC‑AGCAGGTCATGCATGTTCAGGGCGCCGGTCGGGCGGTCTTCACGGTCCACCACCACCAGGGCGTTGATCTTGTGGT‑CTTACATTATCTTCAGGGCTTcggcg                                                                                                                        >  8:113509/1‑149 (MQ=255)
                             aGGTCCTGGTTCATTCCATCGCTCCTTTACATTACGCCGGCGCGC‑AGCAGGTCATGCATGTTCAGGGCGCCGGTCGGGCGGTCTTCACGGTCCACCACCACCAGGGCGTTGATCTTGTGGT‑CTTCCATTATCTTCAGGGCTTCggcggc                                                                                                                      <  8:195850/149‑1 (MQ=255)
                                  cTGGTTCATTCCATCGCTCCTTTACATTACGCCGGCGCGC‑AGCAGGTCATGCATGTTCAGGGCGCCGGTCGGGCGGTCTTCACGGTCCACCACCACCAGGGCGTTGATCTTGTGGT‑CTTCCATTATCTTCAGGGCTTCGGCGGCGAGCa                                                                                                                 >  3:252970/1‑149 (MQ=255)
                                     gTTCATTCCATCGCTCCTTTACATTACGCCGGCGCGC‑AGCAGGTCATGCATGTTCAGGGCGCCGGTCGGGCGGTCTTCACGGTCCACCACCACCAGGGCGTTGATCTTGTGGT‑CTTCCATTATCTTCAGGGCTTCGGCGGCGAGCATGt                                                                                                              >  5:189884/1‑149 (MQ=255)
                                     gTTCATTCCATCGCTCCTTTACATTACGCCGGCGCGC‑AGCAGGTCATGCATGTTCAGGGCGCCGGTCGGGCGGTCTTCACGGTCCACCACCACCAGGGCGTTGATCTTGTGGT‑CTTCCATTATCTTCAGGGCTTCGGCGGCGAGCATGt                                                                                                              >  6:75667/1‑149 (MQ=255)
                                     gTTCATTCCATCGCTCCTTTACATTACGCCGGCGCGC‑AGCAGGTCATGCATGTTCAGGGCGCCGGTCGGGCGGTCTTCACGGTCCACCACCACCAGGGCGTTGATCTTATGGT‑CTTCCATTATCTTCAGGGCTTCGGCGGCGAGCATGt                                                                                                              >  2:175294/1‑149 (MQ=255)
                                     gTTCATTCCATCGCTCCTTTACATTACGCCGGCGCGC‑AGCAGGTCATGCATGTTCAGGGCGCCGGTCGGGCGGTCTTCACGGTCCACCACCACCAGGGCGTTGATCTTATGGT‑CTTCCATTATCTTCAGGGCTTCGGCGGCGAGCATGt                                                                                                              >  2:169601/1‑149 (MQ=255)
                                          ttCCATCGCTCCTTTACATTACGCCGGCGCGC‑AGCAGGTCATGCATGTTCAGGGCGCCGGTCGGGCGGTCTTCACGGTCCACCACCACCAGGGCGTTGATCTTGTGGT‑CTTCCATTATCTTCAgg                                                                                                                                 <  8:307894/125‑1 (MQ=255)
                                          ttCCATCGCTCCTTTACATTACGCCGGCGCGC‑AGCAGGTCATGCATGTTCAGGGCGCCGGTCGGGCGGTCTTCACGGTCCACCACCACCAGGGCGTTGATCTTGTGGT‑CTTCCATTATCTTCAgg                                                                                                                                 >  7:307894/1‑125 (MQ=255)
                                            ccATCGCTCCTTTACATTACGCCGGCGCGC‑AGCAGGTCATGCATGTTCAGGGCGCCGGTCGGGCGGTCTTCACGGTCCACCACCACCAGGGCGTTGATCTTGTGGT‑CTTCCATTATCTTCAGGGCTTCGGCGGCGAGCATGTCGGCAcg                                                                                                       >  5:48859/1‑149 (MQ=255)
                                             cATCGCTCCTTTACATTACGCCGGCGCGCGAGCAGGTCATGCATGTTCAGGGCGCCGGTCGGGCGGTCTTCACGGTCCACCACCACCAGGGCGTTGATCTTGTGGT‑CTTCCATTATCTTCAGGGCTTCGGCGGCGAGCATGTCGGCAcg                                                                                                       >  1:218537/1‑149 (MQ=255)
                                              aTCGCTCCTTTACATTACGCCGGCGCGC‑AGCAGGTCATGCATGTTCAGGGCGCCGGTCGGGCGGTCTTCACGGTCCACCACCACCAGGGCGTTGATCTTGTGGT‑CTTCCATTATCTTCAGGGCTTCGGCGGCGAGCATGTCGGCAcgcg                                                                                                     <  4:252970/149‑1 (MQ=255)
                                                  cTCCTTTACATTACGCCGGCGCGC‑AGCAGGTCATGCATGTTCAGGGCGCCGGTCGGGCGGTCTTCACGGTCCACCACCACCAGGGCGTTGATCTTGTGGT‑CTTCCATTATCTTCAGGGCTTCGGCGGCGAGCATGTCg                                                                                                            <  1:235015/138‑1 (MQ=255)
                                                  cTCCTTTACATTACGCCGGCGCGC‑AGCAGGTCATGCATGTTCAGGGCGCCGGTCGGGCGGTCTTCACGGTCCACCACCACCAGGGCGTTGATCTTGTGGT‑CTTCCATTATCTTCAGGGCTTCGGCGGCGAGCATGTCg                                                                                                            >  2:235015/1‑138 (MQ=255)
                                                  cTCCTTTACATTACGCCGGCGCGC‑AGCAGGTCATGCATGTTCAGGGCGCCGGTCGGGCGGTCTTCACGGTCCACCACCACCAGGGCGTTGATCTTGTGGT‑CTTCCATTATCTTCAGGGCTTCGGCGGCGAGCATGTCGGCACGCGCAGt                                                                                                 >  2:274360/1‑149 (MQ=255)
                                                    ccTTTACATTACGCCGGCGCGC‑AGCAGGTCATGCATGTTCAGGGCGCCGGTCGGGCGGTCTTCACGGTCCACCACCACCAGGGCGTTGATCTTGTGGT‑CTTCCATTATCTTCAGGGCTTCGGCGGCGAGCATGTCGGCACGCGCAGTCt                                                                                               >  8:318766/1‑149 (MQ=255)
                                                     cTTTACATTACGCCGGCGCGC‑AGCAGGTCATGCATGTTCAGGGCGCCGGTCGGGCGGTCTTCACGGTCCACCACCACCAGGGCGTTGATCTTGTGGT‑CTTCCATTATCTTCAGGGCTTCGGCGGCGAGCATGTCGGCACGCGCAGTCtt                                                                                              >  1:281084/1‑149 (MQ=255)
                                                     cTTTACATTACGCCGGCGCGC‑AGCAGGTCATGCATGTTCAGGGCGCCGGTCGGGCGGTCTTCACGGTCCACCACCACCAGGGCGTTGATCTTGTGGT‑CTTCCATTATCTTCAGGGCTTCGGCGGCGAGCATGTCGGCACGCGCAGTCtt                                                                                              <  7:24782/149‑1 (MQ=255)
                                                     cTTTACATTACGCCGGCGCGC‑AGCAGGTCATGCATGTTCAGGGCGCCGGTCGGGCGGTCTTCACGGTCCACCACCACCAGGGCGTTGATCTTGTGGT‑CTTCCATTATCTTCAGGGCTTCGGCGGCGAGCATGTCGGCACGCGCAGTCtt                                                                                              <  1:274360/149‑1 (MQ=255)
                                                        tACATTACGCCGGCGCGC‑AGCAGGTCATGCATGTTCAGGGCGCCGGTCGGGCGGTCTTCACGGTCCACCACCACCAGGGCGTTGATCTTGTGGT‑CTTCCATTATCTTCAGGGCTTCGGCGGCGAGCATGTCGGCACGCGCAGTCTTg                                                                                             <  2:218537/147‑1 (MQ=255)
                                                          cATTACGCCGGCGCGC‑AGCAGGTCATGCATGTTCAGGGCGCCGGTCGGGCGGTCTTCACGGTCCACCACCACCAGGGCGTTGATCTTGTGGT‑CTTCCATTATCTTCAGGGCTTCGGCGGCGAGCATGTCGGCACGCGCAGTCTTGCCAt                                                                                         <  4:121652/149‑1 (MQ=255)
                                                          cATTACGCCGGCGCGC‑AGCAGGTCATGCATGTTCAGGGCGCCGGTCGGGCGGTCTTCACGGTCCACCACCACCAGGGCGTTGATCTTGTGGT‑CTTCCATTATCTTCAGGGCTTCGGCGGCGAGCATGTCGGCACGCGCAGTCTTGCCAt                                                                                         <  4:105138/149‑1 (MQ=255)
                                                           aTTACGCCGGCGCGC‑AGCAGGTCATGCATGTTCAGGGCGCCGGTCGGGCGGTCTTCACGGTCCACCACCACCAGGGCGTTGATCTTGTGGT‑CTTCCATTATCTTCAGGGCTTCGGCGGCGAGCATGTCGGCACGCGCAGTCTTGCCAt                                                                                         >  2:72881/1‑148 (MQ=255)
                                                                  cGGCGCGC‑AGCAGGTCATGCATGTTCAGGGCGCCGGTCGGGCGGTCTTCACGGTCCACCACCACCAGGGCGTTGATCTTGTGGT‑CTTCCATTATCTTCAGGGCTTCGGCGGCGAGCATGTCGGCACGCGCAGTCTTGCCATGACCGGt                                                                                  >  7:169264/1‑148 (MQ=255)
                                                                  cGGCGCGC‑AGCAGGTCATGCATGTTCAGGGCGCCGGTCGGGCGGTCTTCACGGTCCACCACCACCAGGGCGTTGATCTTGTGGT‑CTTCCATTATCTTCAGGGCTTCGGCGGCGAGCATGTCGGCACGCGCAGTCTTGCCATGAACGGTc                                                                                 <  7:318766/149‑1 (MQ=255)
                                                                   ggCGCGC‑AGCAGGTCATGCATGTTCAGGGCGCCGGTCGGGCGGTCTTCACGGTCCACCACCACCAGGGCGTTGATCTTGTGGT‑CTTCCATTATCTTCAGGGCTTCGGCGGCGAGCATGTCGGCACGCGCAGTCTTGCCATGAACGGTCa                                                                                >  5:136124/1‑149 (MQ=255)
                                                                            gcagGTCATGCATGTTCAGGGCGCCGGTCGGGCGGTCTTCACGGTCCACCACCACCAGGGCGTTGATCTTGTGGT‑CTTCCATTATCTTCAGGGCTTCGGCGGCGAGCATGTCGGCACGCGCAGTCTTGCCATGAACGGTCATGACCTgg                                                                        <  6:48859/149‑1 (MQ=255)
                                                                             cagGTCATGCATGTTCAGGGCGCCGGTCGGGCGGTCTTCACGGTCCACCACCACCAGGGCGTTGATCTTGTGGT‑CTTCCATTATCTTCAGGGCTTCGGCGGCGAGCATGTCGGCACGCGCAGTCTTGCCATGAACGGTCATGACCTGGt                                                                       <  5:198896/149‑1 (MQ=255)
                                                                               ggTCATGCATGTTCAGGGCGCCGGTCGGGCGGTCTTCACGGTCCACCACCACCAGGGCGTTGATCTTGTGGT‑CTTCCATTATCTTCAGGGCTTCGGCGGCGAGCATGTCGGCACGCGCAGTCTTGCCATGAACGGTCATGACCTGGTc                                                                      >  8:6258/1‑148 (MQ=255)
                                                                                 tCATGCATGTTCAGGGCGCCGGTCGGGCGGTCTTCACGGTCCACCACCACCAGGGCGTTGATCTTGTGGT‑CTTCCATTATCTTCAGGGCTTCGGCGGCGAGCATGTCGGCCCGCGCAGTCTTGCCATGAACGGTCATGACCTGGTCGAt                                                                   >  8:206932/1‑149 (MQ=255)
                                                                                             aGGGCGCCGGTCGGGCGGTCTTCACGGTCCACCACCACCAGGGCGTTGATCTTGTGGT‑CTTCCATTATCTTCAGGGCTTCGGCGGCGAGCATGTCGGCACGCGCAGTCTTGCCATGAACGGTCATGACCTGGTCGATCAGGGTGgtgt                                                        <  7:206932/148‑1 (MQ=255)
                                                                                                gcgcCGGTCGGGCGGTCTTCACGGTCCACCACCACCAGGGCGTTGATCTTGTGGT‑CTTCCATTATCTTCAGGGCTTCGGCGGCGAGCATGTCGGCACGCGCAGTCTTGCCATGAACGGTCATGACATGGTCGATCAGGGTGGTGTGca                                                     >  8:4413/1‑148 (MQ=255)
                                                                                                gcgcCGGTCGGGCGGTCTTCACGGTCCACCACCACCAGGGCGTTGATCTGGTGGT‑CTTCCATTATCTTCAGGGCTTCGGCGGCGAGCATGTCGGCACGCGCAGTCTTGCCATGAACGGTCATGACCTGGTCGATCAGGGTGGTGTGcac                                                    <  4:135335/149‑1 (MQ=255)
                                                                                                         gggCGGTCTTCACGGTCCACCACCACCAGGGCGTTGATCTTGTGGT‑CTTCCATTATCTTCAGGGCTTCGGCGGCGAGCATGTCGGCACGCGCAGTCTTGCCATGAACGGTCATGACCTGGTCGATCAGGGTGGTGTGCACATCGATAtt                                           <  6:89268/149‑1 (MQ=255)
                                                                                                          ggCGGTCTTCACGGTCCACCACCACCAGGGCGTTGATCTTGTGGTGCTTCCATTATCTTCAGGGCTTCGGCGGCGAGCATGTCGGCACGCGCAGTCTTGCCATGAACGGTCATGACCTGGTCGATCAGGGTGGTGTGCACATCGATAtt                                           >  1:76982/1‑149 (MQ=255)
                                                                                                            cGGTCTTCACGGTCCACCACCACCAGGGCGTTGATCTTGTGGT‑CTTCCATTATCTTCAGGGCTTCGGCGGCGAGCATGTCGGCACGCGCAGTCTTGCCATGAACGGTCATGACCTGGTCGATCAGGGTGGTGTGCACATCGATATTGCg                                        <  7:113509/149‑1 (MQ=255)
                                                                                                            cGGTCTTCACGGTCCACCACCACCAGGGCGTTGATCTTGTGGT‑CTTCCATTATCTTCAGGGCTTCGGCGGCGAGCATGTCGGCACGCGCAGTCTTGCCATGAACGGTCATGACCTGGTCGATCAGGGTGGTGTGCACATCGATATTGCg                                        <  5:75667/149‑1 (MQ=255)
                                                                                                               tCTTCACGGTCCACCACCACCAGGGCGTTGATCTTGTGGT‑CTTCCATTATCTTCAGGGCTTCGGCGGCGAGCATGTCGGCACGCGCAGTCTTGCCATGAACGGTCATGACCTGGTCGATCAGGGTGGTGTGCACATCGATATTGCGGTc                                     >  6:274834/1‑149 (MQ=255)
                                                                                                                   cACGGTCCACCACCACCAGGGCGTTGATCTTGTGGT‑CTTCCATTATCTTCAGGGCTTCGGCGGCGAGCATGTCGGCACGCGCAGTCTTGCCATGAACGGTCATGACCTGGTCGATCAGGGTGGTGTGCACATCGATATTg                                          >  1:35911/1‑140 (MQ=255)
                                                                                                                   cACGGTCCACCACCACCAGGGCGTTGATCTTGTGGT‑CTTCCATTATCTTCAGGGCTTCGGCGGCGAGCATGTCGGCACGCGCAGTCTTGCCATGAACGGTCATGACCTGGTCGATCAGGGTGGTGTGCACATCGATATTg                                          <  2:35911/140‑1 (MQ=255)
                                                                                                                   cACGGTCCACCACCACCAGGGCGTTGATCTTATGGT‑CTTCCATTATCTTCAGGGCTTCGGCGGCGAGCATGTCGGCACGCGCAGTCTTGCCATGAACGGTCATGACCTGGTCGATCAGGGTGGTGTGCACATCGATATTGCGGTCCa                                   <  1:169601/147‑1 (MQ=255)
                                                                                                                   cACGGTCCACCACCACCAGGGCGTTGATCTTATGGT‑CTTCCATTATCTTCAGGGCTTCGGCGGCGAGCATGTCGGCACGCGCAGTCTTGCCATGAACGGTCATGACCTGGTCGATCAGGGTGGTGTGCACATCGATATTGCGGTCCa                                   <  1:175294/147‑1 (MQ=255)
                                                                                                                           accaccaccaGGGCGTTGATCTTGTGGT‑CTTCCATTATCTTCAGGGCTTCGGCGGCGAGCATGTCGGCACGCGCAGTCTTGCCATGAACGGTCATGACCTGGTCGATCAGGGTGGTGTGCACATCGATATTGCGGTCCAGGCTgcggcg                         >  2:155235/1‑149 (MQ=255)
                                                                                                                            ccaccaccaGGGCGTTGATCTTGTGGT‑CTTCCATTATCTTCAGGGCTTCGGCGGCGAGCATGTCGGCACGCGCAGTCTTGCCATGAACGGTCATGACCTGGTCGATCAGGGTGGTGTGCACATAGATATTGCGGTCCAGGCTGCGgcgc                        <  6:136124/149‑1 (MQ=255)
                                                                                                                             caccaccaGGGCGTTGATCTTGTGGT‑CTTCCATTATCTTCAGGGCTTCGGCGGCGAGCATGTCGGCACGCGCAGTCTTGCCATGAACGGTCATGACCTGGTCGATCAGGGTGGTGTGCACATCGATATTGCGGTCCAGGCTGCGGCGCa                       <  7:6258/149‑1 (MQ=255)
                                                                                                                              agcaccaGGGCGTTGAGCTTGGTGG‑CATCCATTATCTTCAGGGGTTCGGCGGCGAGCATGTCGGCACGCGCAGGCTTGCCATGAACGGTCATGACCTGGTCGATGAGTGTGGGGTGCACAGCGTTATTGCGGTCCAGGCTGCGGCGCAg                      <  4:93148/147‑1 (MQ=255)
                                                                                                                               ccaccaGGGCGTTGATCTTGTGGT‑CTTCCATTATCTTCAGGGCTTCGGCGGCGAGCATGTCGGCACGCGCAGTCTTGCCATGAACGGTCATGACCTGGTCGAACAGGGTGGTGGGCACATCGATATTGCGGTCCAGGCTGCGGCGCAgg                     >  4:279389/1‑149 (MQ=255)
                                                                                                                                    aGGGCGTTGATCTTGTGGT‑CTTCCATTATCTTCAGGGCTTCGGCGGCGAGCATGTCGGCACGCGCAGTCTTGCCATGAACGGTCATGACCTGGGCGATCAGGGTGGTGTGCACATCGATATTGCGGTCCAGGCTGCGGCGCAGGTCAcc                >  4:284668/1‑149 (MQ=255)
                                                                                                                                             aTCTTGTGGT‑CTTCCATTATCTTCAGGGCTTCGGCGGCGAGCATGTCGGCACGCGCAGTCTTGCCATGAACGGTCATGACCTGGTCGATCAGGGTGGTGTGCACATCGATATTGCGGTCCAGGCTGCGGCGCAGGTCACCGTCGGTGa        <  7:339135/148‑1 (MQ=255)
                                                                                                                                                  gtgGT‑CTTCCATTATCTTCAGGGCTTCGGCGGCGAGCATGTCGGCACGCGCAGTCTTGCCATGAACGGTCATGACCTGGTCGATCAGGGTGGTGTGCACATCGATATTGCGGTCCAGGCTGCGGCGCAGGTCACCGTCGGTGAAg      <  5:274834/145‑1 (MQ=255)
                                                                                                                                                  gtgGT‑CTTCCATTATCTTCAGGGCTTCGGCGGCGAGCATGTCGGCACGCGCAGTCTTGCCATGAACGGTCATGACCTGGTCGATCAGGGTGGTGTGCACATCGATATTGCGGTCCAGGCTGCGGCGCAGGTCACCGTCGGTGAAGATcc  >  6:252425/1‑149 (MQ=255)
                                                                                                                                                  gtgGT‑CTTCCATTATCTTCAGGGCTTCGGCGGCGAGCATGTCGGCACGCGCAGTCTTGCCATGAACGGTCATGACCTGGTCGATCAGGGTGGTGTGCACATCGATATTGCGGTCCAGGCTGCGGCGCAGGTCACCGTCGGTGAAGATcc  >  1:281655/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
CCAGCTTGATGCCTTTGCCGCGTTGCATGAGGTCCTGGTTCATTCCATCGCTCCTTTACATTACGCCGGCGCGC‑AGCAGGTCATGCATGTTCAGGGCGCCGGTCGGGCGGTCTTCACGGTCCACCACCACCAGGGCGTTGATCTTGTGGT‑CTTCCATTATCTTCAGGGCTTCGGCGGCGAGCATGTCGGCACGCGCAGTCTTGCCATGAACGGTCATGACCTGGTCGATCAGGGTGGTGTGCACATCGATATTGCGGTCCAGGCTGCGGCGCAGGTCACCGTCGGTGAAGATCC  >  NC_002947/1098536‑1098829

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