Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 69,935 G→C 5.2% A657G (GCC→GGC)  glyS ← glycine‑tRNA ligase subunit beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00294769,9350GC 5.2% 251.8 / 3.0 77A657G (GCC→GGC) glySglycine‑tRNA ligase subunit beta
Reads supporting (aligned to +/- strand):  ref base G (32/41);  new base C (2/2);  total (34/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.26e-01

ATGTACGCATCGGAGTCGTAGTTGATCACCCCGTCACGGTCGAGAATCAGCAGTTTCAAGGCTTACCCCAGCAGCGAAATGTCAGCCACGCCCAGGAACAGGCCGCGCAGACGGCTGAGCAGGGCATAACGGTTGGCACGTACCTTGGCGTCTTCGGCGTTGACCATGACCGCCTCGAAGAAGGCATCGACCGGGTCACGCAGGGCCGCCAGGCGGGCCAGCGATTCGCTGTACTGACGTGCTGCAGCCATCGGTTGCACGGCCTGGTCGGCCTGCTGGATGGCCGAGTA  >  NC_002947/69788‑70077
                                                                                                                                                     |                                                                                                                                               
aTGTACGCATCGGAGTCGTAGTTGATCACCCCGTCACGGTCGAGAATCAGCAGTTTCAAGGCTTACCCCAGCAGCGAAATGTCAGCCACGCCCAGGAACAGGCCGCGCAGACGGCTGAGCAGGGCATAACGGTTGGCACGTACCTTgg                                                                                                                                                 <  4:52344/148‑1 (MQ=255)
 tGTACGCATCGGAGTCGTAGTTGATCACCCCGTCACGGTCGAGAATCAGCAGTTTCAAGGCTTACCCCAGCAGCGAAATGTCAGCCACGCCCAGGAACAGGCCGCGCAGACGGCTGAGCAGGGCATAACGGTTGGCACGTACCTTGGc                                                                                                                                                <  1:137912/148‑1 (MQ=255)
 tGTACGCATCGGAGTCGTAGTTGATCACCCCGTCACGGTCGAGAATCAGCAGTTTCAAGGCTTACCCCAGCAGCGAAATGTCAGCCACGCCCAGGAACAGGCCGCGCAGACGGCTGAGCAGGGCATAACGGTTGGCACGTACCTTGGc                                                                                                                                                <  8:175322/148‑1 (MQ=255)
  gTACGCATCGGAGTCGTAGTTGATCACCCCGTCACGGTCGAGAATCAGCAGTTTCAAGGCATACCCCAGCAGCGAAATGTCAGCCACGCCCACGAACAGGCCGCGCAGACGGCTGAGCAGGGCATAACGGTTGGCACGTACCTTGGCGt                                                                                                                                              >  2:313253/1‑149 (MQ=255)
      gCATCGGAGTCGTAGTTGATCACCCCGTCACGGTCGAGAATCAGCAGTTTCAAGGCTTACCCCAGCAGCGAAATGTCAGCCACGCCCAGGAACAGGCCGCGCAGACGGCTGAGCAGGGCATAACGGTTGGCACGTACCTTGGCGTCTTc                                                                                                                                          <  2:236412/149‑1 (MQ=255)
        aTCGGAGTCGTAGTTGATCACCCCGTCACGGTCGAGAATCAGCAGTTTCAAGGCTTACCCCAGCAGCGAAATGTCAGCCACGCCCAGGAACAGGCCGCGCAGACGGCTGAGCAGGGCATAACGGTTGGCACGTACCTTGGCGTCTTCgg                                                                                                                                        <  1:303637/149‑1 (MQ=255)
          cGGAGTCGTAGTTGATCACCCCGTCACGGTCGAGAATCAGCAGTTTCAAGGCTTACCCCAGCAGCGAAATGTCAGCCACGCCCAGGAACAGGCCGCGCAGACGGCTGAGCAGGGCATAACGGTTGGCACGTACCTTGGCGTCTTCGGCg                                                                                                                                      >  8:253973/1‑149 (MQ=255)
            gAGTCGTAGTTGATCACCCCGTCACGGTCGAGAATCAGCAGTTTCAAGGCTTACCCCAGCAGCGAAATGTCAGCCACGCCCAGGAACAGGCCGCGCAGACGGCTGAGCAGGGCATAACGGTTGGCACGTACCTTGGCGTCTTCGGCGtt                                                                                                                                    <  4:391371/149‑1 (MQ=255)
                 gTAGTTGATCACCCCGTCACGGTCGAGAATCAGCAGTTTCAAGGCTTACCCCAGCAGCGAAATGTCAGCCACGCCCAGGAACAGGCCGCGCAGACGGCTGAGCAGGGCATGACGGTTGGCACGTACCTTGGCGTCTTCGGCGTTGACCa                                                                                                                               >  4:215763/1‑149 (MQ=255)
                  tAGTTGATCACCCCGTCACGGTCGAGAATCAGCAGTTTCAAGGCTTACCCCAGCAGCGAAATGTCAGCCACGCCCAGGAACAGGCCGCGCAGACGGCTGAGCAGGGCATAACGGTTGGCACGTACCTTGGCGTCTTCGGCGTTGAc                                                                                                                                 >  1:207100/1‑146 (MQ=255)
                  tAGTTGATCACCCCGTCACGGTCGAGAATCAGCAGTTTCAAGGCTTACCCCAGCAGCGAAATGTCAGCCACGCCCAGGAACAGGCCGCGCAGACGGCTGAGCAGGGCATAACGGTTGGCACGTACCTTGGCGTCTTCGGCGTTGAc                                                                                                                                 <  2:207100/146‑1 (MQ=255)
                    gTTGATCACCCCGTCACGGTCGAGAATCAGCAGTTTCAAGGCTTACCCCAGCAGCGAAATGTCAGCCACGCCCAGGAACAGGCCGCGCAGACGGCCTGAGCAGGGCATAACGGTTGGCACGTACCTTGGCGTCTTCGGCGTTGACCAt                                                                                                                              <  1:313253/148‑1 (MQ=255)
                    gTTGATCACCCCGTCACGGTCGAGAATCAGCAGTTTCAAGGCTTACCCCAGCAGCGAAATGTCAGCCACGCCCAGGAACAGGCCGCGCAGACGGCTGAGCAGGGCATAACGGTTGGCACGTACCTTGGCGTCTTCGGCGTTGACCATGa                                                                                                                            >  8:42836/1‑149 (MQ=255)
                         tCACCCCGTCACGGTCGAGAATCAGCAGTTTCAAGGCTTACCCCAGCAGCGAAATGTCAGCCACGCCCAGGAACAGGCCGCGCAGACGGCTGAGCAGGGCATAACGGTTGGCACGTACCTTGGCGTCTTCGGCGTTGACCATGACCGcc                                                                                                                       <  6:91843/149‑1 (MQ=255)
                           aCCCCGTCACGGTCGAGAATCAGCAGTTTCAAGGGCTTACCCCAGCAGCGAAATGTCAGCCACGCCCAGGAACAGGCCGCGCAGACGGCTGAGCAGGGCATAACGGTTGGCACGTACCTTGGCGTCTTCGGCGTTGACCATGACCGCCt                                                                                                                      >  1:256662/1‑149 (MQ=255)
                            ccccGTCACGGTCGAGAATCAGCAGTTTCAAGGCTTACCCCAGCAGCGAAATGTCAGCCACGCCCAGGAACAGGCCGCGCAGACGGCTGAGCAGGGCATAACGGTTGGCACGTACCTTGGCGTCTTCGGCGtt                                                                                                                                    >  5:36991/1‑133 (MQ=255)
                            ccccGTCACGGTCGAGAATCAGCAGTTTCAAGGCTTACCCCAGCAGCGAAATGTCAGCCACGCCCAGGAACAGGCCGCGCAGACGGCTGAGCAGGGCATAACGGTTGGCACGTACCTTGGCGTCTTCGGCGtt                                                                                                                                    <  6:36991/133‑1 (MQ=255)
                                gTCACGGTCGAGAATCAGCAGTTTCAAGGCTTACCCCAGCAGCGAAATGTCAGCCACGCCCAGGAACAGGCCGCGCAGACGGCTGAGCAGGGCATAACGGTTGGCACGTACCTTGCCGTCTTCGGCGTTGACCATGACCGCCTc                                                                                                                     >  2:228331/1‑144 (MQ=255)
                                gTCACGGTCGAGAATCAGCAGTTTCAAGGCTTACCCCAGCAGCGAAATGTCAGCCACGCCCAGGAACAGGCCGCGCAGACGGCTGAGCAGGGCATAACGGTTGGCACGTACCTTGCCGTCTTCGGCGTTGACCATGACCGCCTc                                                                                                                     <  1:228331/144‑1 (MQ=255)
                                 tCACGGTCGAGAATCAGCAGTTTCAAGGCTTACCCCAGCAGCGAAATGTCAGCCACGCCCAGGAACAGGCCGCGCAGACGGCTGAGCAGGGCATAACGGTTGGCACGTACCTTGGCGTCTTCGGCGTTGACCATGACCGCCTCgaagaa                                                                                                               <  8:362611/149‑1 (MQ=255)
                                    cGGTCGACAAGCAGCAGCCTCCACGCTTACCCCAGCAGCGAAATGCCAGCCACGCCCAGGAACAGGCCGCGCAGACGGCTGAGCAGGGCATAACGGGTGGCACGTACCTTGGCGTCTTCTGCGCTGACCATGACCGCCTCGAAGAAGGc                                                                                                            <  6:291658/149‑1 (MQ=255)
                                       tCGAGAATCAGCAGATTCAAGGCTTACCCCAGCAGCGAAATGTCAGACACGCCCAGGAACAGGCCGCGCAGACGGCTGAGCAGTGCATAACGGTTGGCACGTACCTTGGTGTCTTCGGCGTTGACCATGACCGCCTAGAAGAAGGCATc                                                                                                         >  8:18430/1‑149 (MQ=255)
                                         gagaATCAGCAGTTTCAAGGCTTACCCCAGCAGCGAAATGTCAGCCACGCCCAGGAACAGGCCGCGCAGACGGCTGAGCAGGGCATAACGGTTGGCACGTACCTTGGCGTCTTCGGCGTTGACCATGACCGCCTCGAAGAAGGCATCGa                                                                                                       >  8:224809/1‑149 (MQ=255)
                                          agaATCAGCAGTTTCAAGGCTTACCCCAGCAGCGAAATGTCAGCCACGCCCAGGAACAGGCCGCGCAGACGGCTGAGCAGGGCATAACGGTTGGCACGTACCTTGGCGTCTTCGGCGTTGACCATGACCGCCTCGAAGAAGGCATCGAc                                                                                                      >  2:235269/1‑149 (MQ=255)
                                           gaATCAGCAGTTTCAAGGCTTACCCCAGCAGCGAAATGTCAGCCACGCCCAGGAACAGGCCGCGCAGACGGCTGAGCAGGGCATAACGGTTGGCACGTACCTTGGCGTCTTCGGCGTTGACCATGACCGCCTCGAAGAAGGCATCGAcc                                                                                                     <  3:91436/149‑1 (MQ=255)
                                            aaTCAGCAGTTTCAAGGCTTACCCCAGCAGCGAAATGTCAGCCACGCCCAGGAACAGGCCGCGCAGACGGCTGAGCAGGGCATAACGGTTGGCACGTACCTTGGCGTCTTCGGCGTTGACCATGACCGCCTCGAAGAAGGCATCGACCg                                                                                                    >  7:198928/1‑149 (MQ=255)
                                            aaTCAGCAGTTTCAAGGCTTACCCCAGCAGCGAAATGTCAGCCACGCCCAGGAACAGGCCGCGCAGACGGCTGAGCAGGGCATAACGGTTGGCACGTACCTTGGCGTCTTCGGCGTTGACCATGACCGCCTCGAAGAAGGCATCGACCg                                                                                                    >  3:323979/1‑149 (MQ=255)
                                            aaTCAGCAGTTTCAAGGCTTACCCCAGCAGCGAAATGTCAGCCACGCCCAGGAACAGGCCGCGCAGACGGCTGAGCAGGGCATAACGGTTGGCACGTACCTTGGCGTCTTCGGCGTTGACCATGACCGCCTCGAAGAAGGCATCGACCg                                                                                                    <  7:42836/149‑1 (MQ=255)
                                             aTCAGCAGTTTCAAGGCTTACCCCAGCAGCGAAATGTCAGCCACGCCCAGGAACAGGCCGCGCAGACGGCTGAGCAGGGCATAACGGTTGGCACGTACCTTGGCGTCTTCGGCGTTGACCATGACCGCCTCGAAGAAGGCATCGACCgg                                                                                                   <  4:323979/149‑1 (MQ=255)
                                              tCAGCAGTTTCAAGGCTTACCCCAGCAGCGAAATGTCAGCCACGCCCAGGAACAGGCCGCGCAGACGGCTGAGCAGGGCATAACGGTTGGCACGTACCTTGGCGTCTTCGGCGTTGACCATGACCGCCt                                                                                                                      >  6:238675/1‑129 (MQ=255)
                                              tCAGCAGTTTCAAGGCTTACCCCAGCAGCGAAATGTCAGCCACGCCCAGGAACAGGCCGCGCAGACGGCTGAGCAGGGCATAACGGTTGGCACGTACCTTGGCGTCTTCGGCGTTGACCATGACCGCCt                                                                                                                      <  5:238675/129‑1 (MQ=255)
                                                 gcagTTTCAAGGCTTACCCCAGCAGCGAAATGTCAGCCACGCCCAGGAACAGGCCGCGCAGACGGCTGAGCAGGGCATAACGGTTGGCACGTACCTTGGCGTCTTCGGCGTTGACCATGACCGCCTCGAAGAAGGCATCGACCGGGTCa                                                                                               >  5:46371/1‑149 (MQ=255)
                                                 gcagTTTCAAGGCTTACCCCAGCAGCGAAATGTCAGCCACGCCCAGGAACAGGCCGCGCAGACGGCTGAGCAGGGCATAACGGTTGGCACGTACCTTGGCGTCTTCGGCGTTGACCATGACCGCCTCGAAGAAGGCATCGACCGGGTCa                                                                                               >  5:412923/1‑149 (MQ=255)
                                                 gcagTTTCAAGGCTTACCCCAGCAGCGAAATGTCAGCCACGCCCAGGAACAGGCCGCGCAGACGGCTGAGCAGGGCATAACGGTTGGCACGTACCTTGGCGTCTTCGGCGTTGACCATGACCGCCTCGAAGAAGGCATCGACCGGGTCa                                                                                               <  4:223140/149‑1 (MQ=255)
                                                 gcagTTTCAAGGCTTACCCCAGCAGCGAAATGTCAGCCACGCCCAGGAACAGGCAGCGCAGACGGCTGAGCAGGGCATAACCGTTGGCACGTACCTTGGCGTCTTCGGCGTTGACCATGACCGCCTCGAAGAAGGCATCGACCGGGTCa                                                                                               >  6:396591/1‑149 (MQ=255)
                                                    gTTTCAAGGCTTACCCCAGCAGCGAAATGTCAGCCACGCCCAGGAACAGGCCGCGCAGACGGCTGAGCAGGGCATCACGGTTGGCACGTACCTTGGCGTCTTCGGCGTTGACCATGACCGCCTCGAAGAAGGCATCGACCGGGTCACGc                                                                                            <  3:215763/149‑1 (MQ=255)
                                                       tCAAGGCTTACCCCAGCAGCGAAATGTCAGCCACGCCCAGGAACAGGCCGCGCAGACGGCTGAGCAGGGCATAACGGTTGGCACGTACCTTGGCGTCTTCGGCGTTGACCATGACCGCCTCGAAGAAGGCATCGACCGGGTCACGCAg                                                                                          <  7:253973/148‑1 (MQ=255)
                                                             gCTTACCCCAGCAGCGAAATGTCAGCCACGCCCAGGAACAGGCCGCGCAGACGGCTGAGCAGGGCATAACGGTTGGCACGTACCTTGGCGTCTTCGGCGTTGACCATGACCGCCTCGAAGAAGGCATCGACCGGGTCACGCAGGgccgc                                                                                   >  5:426015/1‑149 (MQ=255)
                                                                 aCCCCAGCAGCGAAATGTCAGCCACGCCCAGGAACAGGCCGCGCAGACGGCTGAGCAGGGCATAACGGGTGGCACGTACCTTGGCGTCTTCGGCGTTGACCATGACCGCCTCGAAGAAGGCATCGACCGGGTCACGCAGGGCCGCCAg                                                                                <  8:30448/148‑1 (MQ=255)
                                                                  ccccAGCAGCGAAATGTCAGCCACGCCCAGGAACAGGCCGCGCAGACGGCTGAGCAGGGCATAACGGTTGGCACGTACCTTGGCGTCTTCGGCGTTGACCATGACCGCCTCGAAGAAGGCATCGACCGGGTCACGCAGGGCCGCCAGGc                                                                              <  6:345832/149‑1 (MQ=255)
                                                                     cagcagCGAAATGTCAGCCACGCCCAGGAACAGGCCGCGCAGACGGCTGAGCAGGGCATAACGGTTGGCACGTACCTTGGCGTCTTCGGCGTTGACCATGACCGCCTCgaagaa                                                                                                               <  4:87001/114‑1 (MQ=255)
                                                                     cagcagCGAAATGTCAGCCACGCCCAGGAACAGGCCGCGCAGACGGCTGAGCAGGGCATAACGGTTGGCACGTACCTTGGCGTCTTCGGCGTTGACCATGACCGCCTCgaagaa                                                                                                               >  3:87001/1‑114 (MQ=255)
                                                                           cGAAATGTCAGCCACGCCCAGGAACAGGCCGCGCAGACGGCTGAGCAGGGCATAACGGTTGGCACGTACCTTGGCGTCTTCGGCGTTGACCATGACCGCCTCGAAGAAGGCATCGACCGGGTCACGCAGGGCCGCCAGGCGGGCCAGc                                                                      >  2:232637/1‑148 (MQ=255)
                                                                             aaaTGTCAGCCACGCCCAGGAACAGGCCGCGCAGACGGCTGAGCAGGGCATAACGGTTGGCACGTACCTTGGCGTCTTCGGCGTTGACCATGACCGCCTCGAAGAAGGCATCGACCGGGTCACGCAGGGCCGCCAGGCGGGCCAGCGAt                                                                   >  6:90111/1‑149 (MQ=255)
                                                                              aaTGTCAGCCACGCCCAGGAACAGGCCGCGCAGACGGCTGAGCAGGGCATAACGGTTGGCACGTACCTTGGCGTCTTCGGCGTTGACCATGACCGCCTCGAAGAAGGCATCGACCGGGTCACGCAGGGCCGCCAGGCGGGCCAGCGAtt                                                                  <  6:412923/149‑1 (MQ=255)
                                                                                tGTCAGCCACGCCCAGGAACAGGCCGCGCAGACGGCTGAGCAGGGCATAACGGTTGGCACGTACCTTGGCGTCTTCGGCGTTGACCATGACCGCCTCGAAGAAGGCATCGACCGGGTCACGCAGGGCCGCCAGGCGGGCCAGCGATTCg                                                                >  6:86617/1‑149 (MQ=255)
                                                                                  tCAGCCACGCCCAGGAACAGGCCGCGCAGACGGCTGAGCAGGGCATAACGGTTGGCACGTACCTTGGCGTCTTCGGCGTTGACCATGACCGCCTCGAAGAAGGCATCGACCGGGTCACGCAGGGCCGCCAGGCGGGCCAGCGATTCGCt                                                              <  2:256662/149‑1 (MQ=255)
                                                                                    aGCCACGCCCAGGAACAGGCCGCGCAGACGGCTGAGCAGGGCATAACGGTTGGCACGTACCTTGGCGTCTTCGGCGTTGACCATGACCGCCTCGAAGAAGGCATCGACCGGGTCACGCAGCGGCCGCCAGGCGGGCCAGCGATTCGCTg                                                             <  1:232637/149‑1 (MQ=255)
                                                                                       cACGCCCAGGAACAGGCCGCGCAGACGGCTGAGCAGGGCATAACGGTTGGCACGTACCTTGGCGTCTTCGGCGTTGACCATGACCGCCTCGAAGAAGGCATCGACCGGGTCACGCAGGgc                                                                                      >  5:232035/1‑120 (MQ=255)
                                                                                       cACGCCCAGGAACAGGCCGCGCAGACGGCTGAGCAGGGCATAACGGTTGGCACGTACCTTGGCGTCTTCGGCGTTGACCATGACCGCCTCGAAGAAGGCATCGACCGGGTCACGCAGGgc                                                                                      <  6:232035/120‑1 (MQ=255)
                                                                                            ccAGGAACAGGCCGCGCAGACGGCTGAGCAGGGCATAACGGTTGGCACGTACCTTGGCGTCTTCGGCGTTGACCATGACCGCCTCGAAGAAGGCATCGACCGGGTCACGCAGGGCCGCCAGGCGGGCCAGCGATTCGCTGTACTGACGt                                                    >  5:120099/1‑149 (MQ=255)
                                                                                            ccAGGAACAGGCCGCGCAGACGGCTGAGCAGGGCATAACGGTTGGCACGTACCTTGGCGTCTTCGGCGTTGACCATGACCGCCTCGAAGAAGGCATCGACCGGGTCACGCAGGGCCGCCAGGCGGGCCAGCGATTCGCTGTACTGACGt                                                    >  3:4027/1‑149 (MQ=255)
                                                                                               ggAACAGGCCGCGCAGACGGCTGAGCAGGGCATAACGGTTGGCACGTACCTTGGCGTCTTCGGCGTTGACCATGACCGCCTCGAAGAAGGCATCGACCGGGTCACGCAGGGCCGCCAGGCGGGCCAGCGATTCg                                                                <  3:189608/134‑1 (MQ=255)
                                                                                               ggAACAGGCCGCGCAGACGGCTGAGCAGGGCATAACGGTTGGCACGTACCTTGGCGTCTTCGGCGTTGACCATGACCGCCTCGAAGAAGGCATCGACCGGGTCACGCAGGGCCGCCAGGCGGGCCAGCGATTCg                                                                >  4:189608/1‑134 (MQ=255)
                                                                                                   cAGGCCGCGCAGACGGCTGAGCAGGGCATAACGGTTGGCACGTACCTTGGCGTCTTCGGCGTTGACCATGACCGCCTCGAAGAAGGCATCGACCGGGTCACGCAGGGCCGCCAGGCGGGCCAGCGATTCGCTGTACTGACGTGCTGCAg                                             <  7:224809/149‑1 (MQ=255)
                                                                                                           gcAGACGGCTGAGCAGGGCATAACGGTTGGCACGTACCTTGGCGTCTTCGGCGTTGACCATGACCGCCTCGAAGAAGGCATCGACCGGGTCACGCAGGGCCGCCAGGCGGGCGAGCGATTCGCTGTACTGACGTGCTGCAGCCATCGGt                                     <  5:396591/149‑1 (MQ=255)
                                                                                                           gcAGACGGCTGAGCAGGGCATAACGGTTGGCACGTACCTTGGCGTCTTCGGCGTTGACCATGACCGCCTCGAAGAAGGCATCGACCGGGTCACGCAGGGCCGCCAGGCGGGCCAGCGATTCGCTGTACTGACGTGCTGCAGCCATCGGt                                     <  6:46371/149‑1 (MQ=255)
                                                                                                                cggcggAGCAGGGCATACCGGTTGGCCCGTACCTTGGCGTCTGCGGCGTTGACCATGCCCGCCTCGAAGAAGGCATCGACCGGGTCACGCAGGGCCGCCAGGCGGGCCAGCGATTCGCTGTACTGACGTGCTGCAGCCATCGGTTGCAc                                <  6:410188/149‑1 (MQ=255)
                                                                                                                     tGAGCAGGGCATAACGGTTGGCACGTACCTTGGCGTCTTCGGCGTTGACCATGACCGCCTCGAACAAGGCATCGACCGGGTCACGCAGGGCCGCCAGGCGGGCCAGCGATTCGCTGTCCTGACGTCCTGCAGCCATCGGTTGCACGGcc                            >  3:12082/1‑149 (MQ=255)
                                                                                                                           gggCATAACGGTTGGCACGTACGTTGGTTTCTTCGGCGTTGACCATGACCGCCTCGAATAAGGCATCGACCGGGTCAGGCAGGGCCGCCAGGCGGGCCAGCGATTCGCTGTACTGACGTGCTTCAGCCATCGGTTTCACGGCCTGGTCg                      <  7:18430/149‑1 (MQ=255)
                                                                                                                            ggCATAACGGTTGGCACGTACCTTGGCGTCTTCGGCGTTGACCATGACCGCCTCGAAGAAGGCATCGACCGGGTCACGCAGGGCCGCCAGGCGGGCCAGCGATTCGCTGTACTGACGTGCTGCAGCCATCGGTTg                                   <  4:305102/135‑1 (MQ=255)
                                                                                                                            ggCATAACGGTTGGCACGTACCTTGGCGTCTTCGGCGTTGACCATGACCGCCTCGAAGAAGGCATCGACCGGGTCACGCAGGGCCGCCAGGCGGGCCAGCGATTCGCTGTACTGACGTGCTGCAGCCATCGGTTg                                   >  3:305102/1‑135 (MQ=255)
                                                                                                                            ggCATAACGGTTGGCACGTACCTTGGCGTCTTCGGCGTTGACCATGACCGCCTCGAAGAAGGCATCGACCGGGTCACGCAGGGCCGCCAGGCGGGCCAGCGATTCGCTGTACTGACGTGCTGCAGCCATCGGTTGCACGGCCTGGTCg                      <  6:120099/148‑1 (MQ=255)
                                                                                                                                tAACGGTTGGCACGTACCTTGGCGTCTTCGGCGTTGACCATGACCGCCTCGAAGAAGGCATCGACCGGGTCACGCAGGGCCGCCAGGCGGGCCAGCGATTCGCTGTACTGACGTGCTGCAGCCATCGGtt                                    <  6:98824/130‑1 (MQ=255)
                                                                                                                                tAACGGTTGGCACGTACCTTGGCGTCTTCGGCGTTGACCATGACCGCCTCGAAGAAGGCATCGACCGGGTCACGCAGGGCCGCCAGGCGGGCCAGCGATTCGCTGTACTGACGTGCTGCAGCCATCGGtt                                    >  5:98824/1‑130 (MQ=255)
                                                                                                                                tAACGGTTGGCACGTACCTTGGCGTCTTCGGCGTTGACCATGACCGCCTCGAAGAAGGCATCGACCGGGTCACGCAGGGCCGCCAGGCGGGCCAGCGATTCGCTGTACTGACGTGCTGCAGCCATCGGTTGCACGGCCTGGTCGGc                    >  4:113594/1‑146 (MQ=255)
                                                                                                                                tAACGGTTGGCACGTACCTTGGCGTCTTCGGCGTTGACCATGACCGCCTCGAAGAAGGCATCGACCGGGTCACGCAGGGCCGCCAGGCGGGCCAGCGATTCGCTGTACTGACGTGCTGCAGCCATCGGTTGCACGGCCTGGTCGGc                    <  3:113594/146‑1 (MQ=255)
                                                                                                                                  aCGGTTGGCACGTACCTTGGCGTCTTCGGCGTTGACCATGACCGCCTCGAAGAAGGCATCGACCGGGTCACGCAGGGCCGCCAGGCGGGCCAGCGATTCGCTGTACTGACGTGCTGCAGCCATCGGTTGCACGGCCTGGTCGGCctgct               <  5:90111/149‑1 (MQ=255)
                                                                                                                                         gCACGTACCTTGGCGTCTTCGGCGTTGACCATGACCGCCTCGAAGAAGGCATCGACCGGGTCACGCAGGGCCGCCAGGCGGGCCAGCGATTCGCTGTACTGACGTGCTGCAGCCATCGGTTGCACGGCCTGGTCGGCCTGCTGGATgg         <  3:16104/148‑1 (MQ=255)
                                                                                                                                           aCGTACCTTGGCGTCTTCGGCGTTGACCATGACCGCCTCGAAGAAGGCATCGACCGGGTCACGCAGGGCCGCCAGGCGGGCCAGCGATTCGCTGTACTGAc                                                      >  2:144439/1‑101 (MQ=255)
                                                                                                                                           aCGTACCTTGGCGTCTTCGGCGTTGACCATGACCGCCTCGAAGAAGGCATCGACCGGGTCACGCAGGGCCGCCAGGCGGGCCAGCGATTCGCTGTACTGAc                                                      <  1:144439/101‑1 (MQ=255)
                                                                                                                                            cGTACCTTGGCGTCTTCGGCGTTGACCATGACCGCCTCgaagaa                                                                                                               <  2:146805/44‑1 (MQ=255)
                                                                                                                                            cGTACCTTGGCGTCTTCGGCGTTGACCATGACCGCCTCgaagaa                                                                                                               >  1:146805/1‑44 (MQ=255)
                                                                                                                                            cGTACCTTGGCGTCTTCGGCGTTGACCATGACCGCCTCGAAGAAGGCATCGACCGGGTCACGCAGGGCCGCCAGGCGGGCCAGCGATTCGCTGTACTGACGTGCTGCAGCCATCGGTTGCACGGCCTGGTCGGCCTGCTGGATGGCCGa     <  8:198928/149‑1 (MQ=255)
                                                                                                                                                ccTTGGCGTCTTCGGCGTTGACCATGACCGCCTCGAAGAAGGCATCGACCGGGTCACGCAGGGCCGCCAGGCGGGCCAGCGATTCGCTGTACTGACGTGCTGCAGCCATCGGTTGCACGGCCTGGTCGGCCTGCTGGATGGCCGAGTa  >  3:276299/1‑148 (MQ=255)
                                                                                                                                                  ttGGCGTGTTCGGCGTTGACCATGACCGCCTCGAAGAAGGCATCGACCGGGTCACGCAGGGCCGCCAGGCGGGCCAGCGATTCGCTGTACTGACGTGCTGCAGCCATCGGTTGCACGGCCTGGTCGGCctgct               <  6:16052/133‑1 (MQ=255)
                                                                                                                                                  ttGCCGTCTTCGGCGTTGACCATGACCGCCTCGAAGAAGGCATCGACCGGGTCACGCAGGgccgc                                                                                   <  3:205413/65‑1 (MQ=255)
                                                                                                                                                  ttGCCGTCTTCGGCGTTGACCATGACCGCCTCGAAGAAGGCATCGACCGGGTCACGCAGGgccgc                                                                                   >  4:205413/1‑65 (MQ=255)
                                                                                                                                                     |                                                                                                                                               
ATGTACGCATCGGAGTCGTAGTTGATCACCCCGTCACGGTCGAGAATCAGCAGTTTCAAGGCTTACCCCAGCAGCGAAATGTCAGCCACGCCCAGGAACAGGCCGCGCAGACGGCTGAGCAGGGCATAACGGTTGGCACGTACCTTGGCGTCTTCGGCGTTGACCATGACCGCCTCGAAGAAGGCATCGACCGGGTCACGCAGGGCCGCCAGGCGGGCCAGCGATTCGCTGTACTGACGTGCTGCAGCCATCGGTTGCACGGCCTGGTCGGCCTGCTGGATGGCCGAGTA  >  NC_002947/69788‑70077

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: