Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,094,517 C→A 6.3% G252C (GGC→TGC)  pycB ← pyruvate carboxylase subunit B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,094,5170CA 6.3% 196.3 / 2.7 63G252C (GGC→TGC) pycBpyruvate carboxylase subunit B
Reads supporting (aligned to +/- strand):  ref base C (33/26);  new base A (2/2);  total (35/28)
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.00e-01

TGCCGCCCGGTACCTGGTTGACCTGCACCCGGGTGTCTACCGCGGTGAACTCGCTTTCGAACTGGTGGTACTTCTTGCGCACAGCGTAGAAGTACAAGCCGATCTCCTGCAGAAGCTCCAGGTCCAGGCCGGTGTCGAACTCGCTGCCCTTGAGTGCAGCAACCATCGACTCGGTACCCGGATGGCTGGTGCCCCAGGCGAAGCTGGAGATGGCGGTGTCGATGTGGTCCGCGCCATTTTCCACGGCCTTCAGCTGGCACATGGCTGCCAGGCCGGCGGTATCGTGGGAGTGAATG  >  NC_002947/6094370‑6094665
                                                                                                                                                    |                                                                                                                                                     
tGCCGCCCGGTACCTGGTTGACCTGCACCGGGGTGTCTACCGCGGTGAACTCGCTTTCGAACTGGTGGTACTTCTTGCGCACAGCGTAGAAGTACAAGCCGATCTCCTGCAGAAGCTCCAGGTCCAGGCCGGTGTCGAACTCGCTGccc                                                                                                                                                      <  3:205343/149‑1 (MQ=255)
tGCCGCCCGGTACCTGGTTGACCTGCACCCGGGTGTCTACCGCGGTGAACTCGCTTTCGAACTGGTGGTACTTCTTGCGCACAGCGTAGAAGTACAAGCCGATCTCCTGCAGAAGCTCCAGGTCCAGGCCGGTGTCGAACTCGCTGccc                                                                                                                                                      <  8:13224/149‑1 (MQ=255)
 gccgccCGGTACCTGGTTGACCTGCACCCGGGTGTCTACCGCGGTGAACTCGCTTTCGAACTGGTGGTACTTCTTGCGCACAGCGTAGAAGTACAAGCCGATCTCCTGCAGAAGCTCCAGGTCCAGGCCGGTGTCGAACTCGCTGCCCt                                                                                                                                                     >  5:294911/1‑149 (MQ=255)
      ccGGTACCTGGTTGACCTGCACCCGGGTGTCTACCGCGGTGAACTCGCTTTCGAACTGGTGGTACTTCTTGCGCACAGCGTAGAAGTACAAGCCGATCTCCTGCAGAAGCTCCAGGTCCAGGCCGGTGTCGAACTCGCTGCCCTTGAGt                                                                                                                                                <  7:20705/149‑1 (MQ=255)
              tGGTTGACCTGCACCCGGGTGTCTACCGCGGTGAACTCGCTTTCGAACTGGTGGTACTTCTTGCGCACAGCGTAGAAGTACAAGCCGATCTCCTGCAGAAGCTCCAGGTCCAGGCCGGTGTCGAACTCGCTGCCCTTGAGTGCAGCAAc                                                                                                                                        >  2:294510/1‑149 (MQ=255)
                gTTGACCTGCACCCGGGTGTCTACCGCGGTGAACTCGCTTTCGAACTGGTGGTACTTCTTGCGCACAGCGTAGAAGTACAAGCCGATCTCCTGCAGAAGCTCCAGGTCCAGGCCGGTGTCGAACTCGCTGCACTCAAGGgcagcgagtt                                                                                                                                      <  2:328782/149‑6 (MQ=255)
                         cACCCGGGTGTCTACCGCGGTGAACTCGCTTTCGAACTGGTGGTACTTCTTGCGCACAGCGTAGAAGTACAAGCCGATCTCCTGCAGAAGCTCCAGGTCCAGGCCGGTGTCGAACTCGCTGCCCTTGAGTGCAGCAACCATCGACTCg                                                                                                                              <  8:75983/148‑1 (MQ=255)
                            ccGGGTGTCTACCGCGGTGAACTCGCTTTCGAACTGGTGGTACTTCTTGCGCACAGCGTAGAAGTACAAGCCGATCTCCTGCAGAAGCTCCAGGTCCAGGCCGGTGTCGAACTCGCTGCCCTTGAGTGCAGCAACCATCGACTCGGTAc                                                                                                                          >  7:222266/1‑149 (MQ=255)
                                   tCTACCGCGGTGAACTCGCTTTCGAACTGGGGGTACTTCTTGCGGACAGCGTAGAAGTACAAGCCGATCTCCTGCATAAGCTCCAGGTCCAGGCCGGTGTCGAACTCGCTGCCCTTGAGTGCAGCAACCATCGACTCGGTACCCGGATg                                                                                                                  <  7:241651/149‑1 (MQ=255)
                                     tACCGCGGTGAACTCGCTTTCGAACTGGTGGTACTTCTTGCGCACAGCGTAGAAGTACAAGCCGATCTCCTGCAGAAGCTCCAGGTCCAGGCCGGTGTCGAACTCGCTGCCCTTGAGTGCAGCAACCATCGACTCGGTACCCGGATGGc                                                                                                                >  7:197789/1‑149 (MQ=255)
                                      aCCGCGGTGAACTCGCTTTCGAACTGGTGGTACTTCTTGCGCACAGCGTAGAAGTACAAGCCGATCTCCTGCAGAAGCTCCAGGTCCAGGCCGGTGTCGAACTCGCTGCCCTTGAGTGCAGCAACCATCGACTCGGTACCCGGATGGCt                                                                                                               >  6:364279/1‑149 (MQ=255)
                                       ccGCGGTGAACTCGCTTTCGAACTGGTGGTACTTCTTGCGCACAGCGTAGAAGTACAAGCCGATCTCCTGCAGAAGCTCCAGGTCCAGGCCGGTGTCGAACTCGCTGCCCTTGAGTGCAGCAACCATCGACTCGGTACCCgg                                                                                                                      <  1:256467/142‑1 (MQ=255)
                                       ccGCGGTGAACTCGCTTTCGAACTGGTGGTACTTCTTGCGCACAGCGTAGAAGTACAAGCCGATCTCCTGCAGAAGCTCCAGGTCCAGGCCGGTGTCGAACTCGCTGCCCTTGAGTGCAGCAACCATCGACTCGGTACCCgg                                                                                                                      >  2:256467/1‑142 (MQ=255)
                                         gcgGTGAACTCGCTTTCGAACTGGTGGTACTTCTTGCGCACAGCGTAGAAGTACAAGCCGATCTCCTGCAGAAGCTCCAGGTCCAGGCCGGTGTCGAACTCGCTGCCCTTGAGTGCAGCAACCATCGACTCGGTACCCGGATGGCTGGt                                                                                                            <  6:230145/149‑1 (MQ=255)
                                          cgGTGAACTCGCTTTCGAACTGGTGGTACTTCTTGCGCACAGCGTAGAAGTACAAGCCGATCTCCTGAAGAAGCTCCAGGTCCAGGCCGGTGTCGAACTCGCTGCCCTTGAGTGCAGCAACCATCGACTCGGTACACGGATGGCTGGTg                                                                                                           >  4:270982/1‑149 (MQ=255)
                                            gTGAACTCGCTTTCGAACTGGTGGTACTTCTTGCGCACAGCGTAGAAGTACAAGCCGATCTCCTGCAGAAGCTCCAGGTCGAGGCCGGTGTCGAGCTCGCTGCAGTTGAGTGCAGCAACAATCGACTCGGTACTCGGATGGATGGTgaa                                                                                                         >  3:211224/1‑147 (MQ=255)
                                             tGAACTCGCTTTCGAACTGGTGGTACTTCTTGCGCACAGCGTAGAAGTACAAGCCGATCTCCTGCAGAAGCTCCAGGTCCAGGCCGGTGTCGAACTCGCTGCCCTTGAGTGCAGCAACCATCGACTCGGTACCCGGATGGCTGGTGccc                                                                                                        <  3:270982/149‑1 (MQ=255)
                                               aaCTCGCTTTCGAACTGGTGGTACTTCTTGCGCACAGCGTAGAAGTACAAGCCGATCTCCTGCAGAAGCTCCAGGTCCAGGCCGGTGTCGAACTCGCTGCCCTTGAGTGCAGCAACCATCGACTCGGTACCCGGATGGCTGGTGcccc                                                                                                       <  8:222266/148‑1 (MQ=255)
                                                  tCGCTTTCGAACTGGTGGTACTTCTTGCGCACAGCGTAGAAGTACAAGCCGATCTCCTGCAGAAGCTCCAGGTCCAGGCCGGTGTCGAACTCGCTGCCCTTGAGTGCAGCAACCATCGACTCGGTACCCGGATGGCTGGTGCCCCAgg                                                                                                    >  1:242948/1‑148 (MQ=255)
                                                    gCTTTCGAACTGGTGGTACTTCTTGCGCACAGCGTAGAAGTACAAGCCGATCTCCTGCAGAAGCTCCAGGTCCAGGCCGGTGTCGAACTCGCTGCCCTTGAGTGCAGCAACCATCGACTCGGTACCCGGATGGCTGGTGCCCCAgg                                                                                                    <  2:242948/146‑1 (MQ=255)
                                                    gCTTTCGAACTGGTGGTACTTCTTGCGCACAGCGTAGAAGTACAAGCCGATCTCCTGCAGAAGCTCCAGGTCCAGGCCGGTGTCGAACTCGCTGCCCTTGAGTGCAGCAACCATCGACTCGGTACCCGGATGGCTGGTGCCCCAGGCg                                                                                                  >  8:118609/1‑148 (MQ=255)
                                                            aCTGGTGGTACTTCTTGCGCACAGCGTAGAAGTACAAGCCGATCTCCTGCAGAAGCTCCAGGTCCAGGCCGGTGTCGAACTCGCTGCCCTTGAGTGCAGCAACCATCGACTCGGTACCCGGATGGCTGGTGCCCCAg                                                                                                     >  3:71746/1‑137 (MQ=255)
                                                            aCTGGTGGTACTTCTTGCGCACAGCGTAGAAGTACAAGCCGATCTCCTGCAGAAGCTCCAGGTCCAGGCCGGTGTCGAACTCGCTGCCCTTGAGTGCAGCAACCATCGACTCGGTACCCGGATGGCTGGTGCCCCAg                                                                                                     <  4:71746/137‑1 (MQ=255)
                                                              tggtggTACTTCTTGCGCACAGCGTAGAGAGTACAAGCCGATCTCCTGCAGAAGCTCCAGGTCCAGGCCGGTGTCGAACTCGCTGCCCTTGAGTGCAGCAACCATCGACTCGGTACCCGGATGGCTGGTGCCCCAGGCGAAGCTg                                                                                            >  1:20952/1‑145 (MQ=255)
                                                                gtggtACTTCTTGCGCACAGCGTAGAAGTACAAGCCGATCTCCTGCAGAAGCTCCAGGTCCAGGCCGGTGTCGAACTCGCTGCCCTTGAGTGCAGCAACCATc                                                                                                                                    >  6:309829/1‑103 (MQ=255)
                                                                gtggtACTTCTTGCGCACAGCGTAGAAGTACAAGCCGATCTCCTGCAGAAGCTCCAGGTCCAGGCCGGTGTCGAACTCGCTGCCCTTGAGTGCAGCAACCATc                                                                                                                                    <  5:309829/103‑1 (MQ=255)
                                                                gtggtACTTCTTGCGCACAGCGTAGAAGTACAAGCCGATCTCCTGCAGAAGCTCCAGGTCCAGGCCGGTGTCGAACTCGCTGCCCTTGAGTGCAGCAACCATCGACTCGGTACCCGGATGGCTGGTGCCCCAGGCGAAGCTGGAGATgg                                                                                     <  7:118609/149‑1 (MQ=255)
                                                                    tACTTCTTGCGCACAGCGTAGAAGTACAAGCCGATCTCCTGCAGAAGCTCCAGGTCCAGGCCGGTGTCGAACTCGCTGCCCTTGAGTGCAGCAACCATCGACTCGGTACCCGGATGGCTGGTGCCCCa                                                                                                      <  7:191812/128‑1 (MQ=255)
                                                                    tACTTCTTGCGCACAGCGTAGAAGTACAAGCCGATCTCCTGCAGAAGCTCCAGGTCCAGGCCGGTGTCGAACTCGCTGCCCTTGAGTGCAGCAACCATCGACTCGGTACCCGGATGGCTGGTGCCCCa                                                                                                      >  8:191812/1‑128 (MQ=255)
                                                                          ttGCGCACAGCGTAGAAGTACAAGCCGATCTCCTGCAGAAGCTCCAGGTCCAGGCCGGTGTCGAACTCGCTGCCCTTGAGTGCAGCAACCATCGACTCGGTACCCGGATGGCTGGTGCCCCAGGCGAAGCTGGAGATGGCCGTGTCGAt                                                                           >  3:231962/1‑149 (MQ=255)
                                                                            gcgcACAGCGTAGAAGTACAAGCCGATCTCCTGCAGAAGCTCCAGGTCCAGGCCGGTTTCGAACTCGCTGCCCTTGAGTGCAGCAACCATCGACTCGGTACCCGGATGGCTGGTGCCCCAGGCGAAGCTGGAGATGGCGGTGTCGa                                                                            <  6:294911/146‑1 (MQ=255)
                                                                            gcgcACAGCGTAGAAGTACAAGCCGATCTCCTGCAGAAGCTCCAGGTCCAGGCCGGTGTCGAACTCGCTGCCCTTGAGTGCAGCAACCATCGACTCGGTACCCGGATGGCTGGTGCCCCAGGCGAAGCTGGAGATGGCGGTGTCGAtg                                                                          <  2:20952/148‑1 (MQ=255)
                                                                              gcACAGCGTAGAAGTACAAGCCGATCTCCTGCAGAAGCTCCAGGTCCAGGCCGGTGTCGAACTCGCTGCCCTTGAGTGCAGCAACCATCGACTCGGTACCCGGATGGCTGGTGCCCCAGGCGAAGCTGGAGATGGCGGTGTCGATGTgg                                                                       >  4:150363/1‑149 (MQ=255)
                                                                                      tAGAAGTACAAGCCGATCTCCTGCAGAAGCTCCAGGTCCAGGCCGGTGTCGAACTCGCTGCCCTTGAGTGCAGCAACCATCGACTCGGTACCCGGATGGCTGGTGCCCCAGGCGAAGCTGGAGATGGCGGTGTCGATGTGGTCCgcgc                                                                >  5:221731/1‑148 (MQ=255)
                                                                                             tACAAGCCGATCTCCTGCAGAAGCTCCAGGTCCAGGCCGGTGTCGAACTCGCTGCCCTTGAGTGCAGCAACCATCGACTCGGTACCCGGATGGCTGGTGCCCCAGGCGAAGCTGGAGATGGCGGTGTCGATGTGGTCCGCGCCATTTTc                                                         >  1:137622/1‑149 (MQ=255)
                                                                                              aCAAGCCGATCTCCTGCAGAAGCTCCAGGTCCAGGCCGGTGTCGAACTCGCTGCCCTTGAGTGCAGCAACCATCGACTCGGTACCCGGATGGCTGGTGCCCCAGGCGAAGCTGGAGATg                                                                                      <  4:294110/119‑1 (MQ=255)
                                                                                              aCAAGCCGATCTCCTGCAGAAGCTCCAGGTCCAGGCCGGTGTCGAACTCGCTGCCCTTGAGTGCAGCAACCATCGACTCGGTACCCGGATGGCTGGTGCCCCAGGCGAAGCTGGAGATg                                                                                      >  3:294110/1‑119 (MQ=255)
                                                                                               cAAGCCGATCTCCTGCAGAAGCTCCAGGTCCAGGCCGGTGTCGAACTCGCTGCCCTTGAGTGCAGCAACCATCGACTCGGTACCCGGATGGCTGGTGCCCCAGGCGAAGCTgg                                                                                           >  1:135629/1‑113 (MQ=255)
                                                                                               cAAGCCGATCTCCTGCAGAAGCTCCAGGTCCAGGCCGGTGTCGAACTCGCTGCCCTTGAGTGCAGCAACCATCGACTCGGTACCCGGATGGCTGGTGCCCCAGGCGAAGCTgg                                                                                           <  2:135629/113‑1 (MQ=255)
                                                                                                aaGCGGATGTCGTGCATAATCTCCAGGTCGAGGCCGGTGTCGAAGTCGCTGCCCTTGAGTGCAGGAACCATCGACTCGGTTCCCGGATGGCTGGTGCCGGAGGCGAAGCTGGAGATGGCGGTGTCGATGTGGTCCGCGCCATTTTCCAc                                                      <  4:231962/149‑1 (MQ=255)
                                                                                                aaGCCGATCTCCTGCAGAAGCTCCAGGTCCAGGCCGGTGTCGAACTCGCTGCCCTTCAGTGCAGCAACCATCGACTCGGTACCCGGATGGCTGGTGCCCCAGGCGAAGCTGGAGATGGCGGTGGCGATGTGGTCCGCGCCATTTTCCAc                                                      <  7:232303/149‑1 (MQ=255)
                                                                                                  gCCGATCTCCTGCAGAAGCTCCAGGTCCAGGCCGGTGTCGAACTCGCTGCCCTTGAGCGCAGCAACCATCGACTCGGTACCCGGATGGCTGGTGCCCCAGGCGAAGCTGGAGATGGCGGAGTCGATGTGGTCCGCGCCATTTTCCACgg                                                    >  8:11321/1‑149 (MQ=255)
                                                                                                   ccGATCTCCTGCAGAAGCTCCAGGTCCAGGCCGGTGTCGAACTCGCTGCCCTTGAGTGCAGCAACCATCGACTCGGTACCCGGATGGCTGGTGCCCCAGGCGAAGCTGGAGATGGCGGTGTCGATGTGGTCCGCGCCATTTTCCACGGc                                                   <  5:249026/149‑1 (MQ=255)
                                                                                                      aTCTCCTGCAGAAGCTCCAGGTCCAGGCCGGTGTCGACCTCGCTGCCCTTGAGTGCAGCAACCATCGACTCGGTACCCGGATGGCTGGTGCCCCAGGCGAAGCTGGAGATGGCGGTGTCGATGTGGTCCGCGCCATTTTCCACGGCCtt                                                >  5:219351/1‑149 (MQ=255)
                                                                                                        ctcCTGCAGAAGCTCCAGGTCCAGGCCGGTGTCGAACTCGCTGCCCTTGAGTGCAGCAACCATCGACTCGGTACCCGGATGGCTGGTGCCCCAGGCGAAGCTGGAGATGGCGGTGTCGATGTGGTCCGCGCCATTTTCCACGGCCtt                                                >  7:36599/1‑147 (MQ=255)
                                                                                                           cTGCAGAAGCTCCAGGTCCAGGCCGGTGTCGAACTCGCTGCCCTTGAGTGCAGCAACCATCGACTCGGTACCCGGATGGCTGGTGCCCCAGGCGAAGCTgg                                                                                           <  2:1902/101‑1 (MQ=255)
                                                                                                           cTGCAGAAGCTCCAGGTCCAGGCCGGTGTCGAACTCGCTGCCCTTGAGTGCAGCAACCATCGACTCGGTACCCGGATGGCTGGTGCCCCAGGCGAAGCTgg                                                                                           >  1:1902/1‑101 (MQ=255)
                                                                                                            tGCAGAAGCTCCAGGTCCAGGCCGGTGTCGAACTCGCTGCCCTTGAGTGCAGCAACCATCGACTCGGTACCCGGATGGCTGGTGCCCCAGGCGAAGCTGGAGATGGCGGTGTCGATGTGGTCCGCGCCATTTTCCACGGCCTTCAGCTg                                          >  4:138323/1‑149 (MQ=255)
                                                                                                              ctgAAGCTCCAGGTCCAGGCCGGTGTCGAACTCGCTGCCCTTGAGTGCAGCAACCATCGACTCGGTACCCGGATGGCTGGTGCCCCAGGCGAAGCTGGAGATGGCGGTGTCGATGTGGTCCGCGCCATTTTCCACGGCCTTCAGCTgg                                         >  4:26925/3‑148 (MQ=255)
                                                                                                                gAAGCTCCAGGTCCAGGCCGGTGTCGAACTCGCTGCCCTTGAGTGCAGCAACCATCGACTCGGTACCCGGATGGCTGGTGCCCCAGGCGAAGCTGGAGATGGCGGTGTCGATGTGGTCCGCGCCATTTTCCACGGCCTTCAGCTGGca                                       >  2:86765/1‑148 (MQ=255)
                                                                                                                       caggtccaggCCGGTGTCGAACTCGCTGCCCTTGAGTGCAGCAACCATCGACTCGGTACCCGGATGGCTGGTGCCCCAGGCGAAGCTGGAGATGGCGGTGTCGATGTGGTCCGCGCCATTTTCCACGGCCTTCAGCTGGCACATGGCTg                               >  6:302480/1‑149 (MQ=255)
                                                                                                                            ccaggCCGGTGTCGAACTCGCTGCCCTTGAGTGCAGCAACCATCGACTCGGTACCCGGATGGCTGGTGCCCCAGGCGAAGCTGGAGATGGCGGTGTCGATGTGGTCCGCGCCATTTTCCACGGCCTTCAGCTGGCACATGGCTGCCAgg                          >  5:13421/1‑149 (MQ=255)
                                                                                                                             caggCCGGTGTCGAACTCGCTGCCCTTGAGTGCAGCAACCATCGACTCGGTACCCGGCTTGGCTGGTGCCCCAGTTGAAGCTGGAGATGGCGGTGTCGATGTGGTCCGCGCCATTTTCCACGGCCTTCAGCTGGCACATGGCTGCCAg                           <  1:294510/148‑1 (MQ=255)
                                                                                                                                         gAACTCGCTGCCCTTGAGTGCAGCAACCATCGACTCGGTACCCGGATGGCTGGTGCCCCAGGCGAAGCTGGAGATGGCGGTGTCGATGTGGTCCGCGCCATTTTCCACGGCCTTGAGCTTGCACATGGCTGCCAGGCCGGCGGTATCGt             <  2:137622/149‑1 (MQ=255)
                                                                                                                                         gAACTCGCTGCCCTTGAGTGCAGCAACCATCGACTCGGTACCCGGATGGCTGGTGCCCCAGGCGAAGCTGGAGATGGCGGTGTCGATGTGGTCCGCGCCATTTTCCACGGCCTTCAGCTGGCACATGGCTGCCAGGCCGGCGGTATCGt             <  8:36599/149‑1 (MQ=255)
                                                                                                                                         gAACTCGCTGCCCTTGAGTGCAGCAACCATCGACTCGGTACCCGGATGACAGGTGCCCCAGGCGAAGCTGGAGAAGGCGGTGTCCATGTGGTACGCGCCATTTTCCACGGCCTTCAGCTGGCACATGGCATCCAGGCCTGCGGTTTCGt             >  1:84971/1‑149 (MQ=255)
                                                                                                                                          aaCTCGCTGCCCTTGAGTGCAGCAACCATCGACTCGGTACCCGGATGGCTGGTGCCCCAGGCGAAGCTGGAGATGGCGGTGTCGATGTGGTCCGCGCCATTTTCCACGGCCTTCAGCTGGCACATGGCTGCCAGGCCGGCGGTATCGTg            >  4:357558/1‑149 (MQ=255)
                                                                                                                                          aaCTCGCTGCCCTTGAGTGCAGCAACCATCGACTCGGTACCCGGATGGCTGGTGCCCCAGGCGAAGCTGGAGATGGCGGTGTCGATGTGGTCCGCGCCATTTTCCACGGCCTTCAGCTGGCACATGGCTGCCAGGCCGGCGGTATCGTg            <  8:322143/149‑1 (MQ=255)
                                                                                                                                           aCTCGCTGCCCTTGAGTGCAGCAACCATCGACTCGGTACCCGGATGGCTGGTGCCCCAGGCGAAGCTGGAGATGGCGGTGTCGATGTGGTCCGCGCCATTTTCCACGGCCTTCAGCTGGCACATGGCTGCCAGGCCGGCGGTATCGTgg           <  5:150885/149‑1 (MQ=255)
                                                                                                                                            cTCGCTGCCCTTGAGTGCAGCAACCATCGACTCGGTACCCGGATGGCTGGTGCCCCAGGCGAAGCTGGAGATGGCGGTGTCGATGTGGTCCGCGCCATTTTCCACGGCCTTCAGCTGGCACATGGCTGCCAGGCAGGCGGTATCGTgg           >  1:52965/1‑148 (MQ=255)
                                                                                                                                             ttgCTGCACTCAAGGGCAGCAACCATCGACTCGGTACCCGGATGGCTGGTGCCCCAGGCGAAGCTGGAGATGGCGGTGTCGAt                                                                           <  8:14806/81‑1 (MQ=255)
                                                                                                                                             ttgCTGCACTCAAGGGCAGCAACCATCGACTCGGTACCCGGATGGCTGGTGCCCCAGGCGAAGCTGGAGATGGCGGTGTCGAt                                                                           >  7:14806/3‑83 (MQ=255)
                                                                                                                                               gCTGCCCTTGAGTGCAGCAACCATCGACTCGGTACCCGGATGGCTGGTGCCCCAGGCGAAGCTGGAGATGGCGGTGTCGATGTGGTCCGCGCCATTTTCCACGGCCTTCAGCTGGCACATg                                   <  1:330768/121‑1 (MQ=255)
                                                                                                                                               gCTGCCCTTGAGTGCAGCAACCATCGACTCGGTACCCGGATGGCTGGTGCCCCAGGCGAAGCTGGAGATGGCGGTGTCGATGTGGTCCGCGCCATTTTCCACGGCCTTCAGCTGGCACATg                                   >  2:330768/1‑121 (MQ=255)
                                                                                                                                                 tGCCCTTGAGTGCAGCAACCATCGACTCGGTACCCGGATGGCTGGTGCCCCAGGCGAAGCTGGAGATGGCGGTGTCGATGTGGTCCGCGCCATTTTCCACGGCCTTCAGCTGGCACATGGCTGCCAGGCCGGCGGTATCGTGGGAGtg      >  5:130150/1‑148 (MQ=255)
                                                                                                                                                   cccTTGAGTGCAGCAACCATCGACTCGGTACCCGGATGGCTGGTGCCCCAGGCGAAGCTGGAGATGGCGGTGTCGATGTGGTCCGCGCCATTTTCCACGGCCTTCAGCTGGCACATGGCTGCCAGGCCGGCGGTATCGTGGGAGtgaat   <  3:357558/149‑1 (MQ=255)
                                                                                                                                                    ccTTGAGTGCAGCAACCATCGACTCGGTACCCGGATGGCTGGTGCCCCAGGCGAAGCTGGAGATGGCGGTGTCGATGTGGTCCGCGCCATTTTCCACGGCCTTCAGCTGGCACATGGCTGCCAGGCCGGCGGTATCGTGGGAGtgaatg  <  1:86765/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TGCCGCCCGGTACCTGGTTGACCTGCACCCGGGTGTCTACCGCGGTGAACTCGCTTTCGAACTGGTGGTACTTCTTGCGCACAGCGTAGAAGTACAAGCCGATCTCCTGCAGAAGCTCCAGGTCCAGGCCGGTGTCGAACTCGCTGCCCTTGAGTGCAGCAACCATCGACTCGGTACCCGGATGGCTGGTGCCCCAGGCGAAGCTGGAGATGGCGGTGTCGATGTGGTCCGCGCCATTTTCCACGGCCTTCAGCTGGCACATGGCTGCCAGGCCGGCGGTATCGTGGGAGTGAATG  >  NC_002947/6094370‑6094665

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: