Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 70,384 G→A 6.0% G507G (GGC→GGT)  glyS ← glycine‑tRNA ligase subunit beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00294770,3840GA 6.0% 197.1 / 2.2 66G507G (GGC→GGT) glySglycine‑tRNA ligase subunit beta
Reads supporting (aligned to +/- strand):  ref base G (28/34);  new base A (2/2);  total (30/36)
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.00e-01

AAGTCCAGGGCGGAGCCTGGCTGCAGTGCGCGTACCGACAGGTACGTGGCGACGTCGATGCCTTCGTCTTCGTAACGCGCGCGCAGGCGGTCGAAAATGAACTCCAGTACTTGCTCGGCCAGGCC‑GGGGGCCTTGATCTTGGCGCCGAACTG‑CTTGACC‑GCAAACTCAACCGCGCCGGTCAGGTCCAGGTCCAGTTGCTTCTCGATCAGGATACGCAGCACGCCCAGAGCAGCACGGCGCAGGGCGTACGGGTCCTTGCTGCCGGTAGGCAGCATGCCGAT  >  NC_002947/70241‑70521
                                                                                                                                                |                                                                                                                                           
aaGTCCAGGGCGGAGCCTGGCTGCAGTGCGCGTACCGACAGGTACGTGGCGACGTCGATGCCTTCGTCTTCGTAACGCGCGCGCAGGCGGTCGAAAATGAACTCCAGTACTTGCTCGGCCAGGCC‑GGGGGCCTTGATCTTGGCGCCGaa                                                                                                                                        >  3:327586/1‑149 (MQ=255)
aaGTCCAGGGCGGAGCCTGGCTGCAGTGCGCGTACCGACAGGTACGTGGCGACGTCGATGCCTTCGTCTTCGTAACGCGCGCGCAGGCGGTCGAAAATGAACTCCAGTACTTGCTCGGCCAGGCC‑GGGGGCCTTGATCTTGGCGCCGaa                                                                                                                                        >  3:161135/1‑149 (MQ=255)
      aGGGCGGAGCCTGGCTGCAGTGCGCGTACCGACAGGTACGTGGCGACGTCGATGCCTTCGTCTTCGTAACGCGCGCGCAGGCGGTCGAAAATGAACTCCAGTACTTGCTCGGCCAGGCC‑GGGGGCCTTGATCTTGGCGCCGAACTG‑Ctt                                                                                                                                 <  6:2233/149‑1 (MQ=255)
          cGGAGCCTGGCTGCAGTGCGCGTACCGACAGGTACGTGGCGACGTCGATGCCTTCGTCTTCGTAACGCGCGCGCAGGCGGTCGAAAATGAACTCCAGTACTTGCTCGGCCAGGCC‑GGGGGCCTTGATCTTGGCGCCGAACTG‑CTTGAcc                                                                                                                             <  4:327586/149‑1 (MQ=255)
              gCCTGGCTGCAGTGCGCGTACCGACAGGTACGTGGCGACGTCGATGCCTTCGTCTTCGTAACGCGCGCGCAGGCGGTCGAAAATGAACTCCAGTACTTGCTCGGCCAGGCC‑GGGGGCCTTGATCTTGGCGCCGAACTG‑CTTGACC‑GCaa                                                                                                                        <  3:114140/149‑1 (MQ=255)
                cTGGCTGCAGTGCGCGTACCGACAGGTACGTTGCGACGTCGATGCCTTCGTCTTCGTAACGCGCGCGCAGGCGGTCGAAAATGAACTCCAGTACTTGCTCGGCCAGGCC‑GGGGGCCTTGATCTTGGCGCCGAACTG‑CTTGACC‑GCAAAc                                                                                                                      <  8:2911/149‑1 (MQ=255)
                cTGGCTGCAGTGCGCGTACCGACAGGTACGTGGCGACGTCGATGCCTTCGTCTTCGTAACGCGCGCGCAGGCGGTCGAAAATGAACTCCAGTACTTGCTCGGCCAGGCC‑GGGGGCCTTGATCTTGGCGCCGAACTG‑CTTGACC‑GCAAAc                                                                                                                      >  1:131120/1‑149 (MQ=255)
                cTGGCTGCAGTGCGCGTACCGACAGGTACGTGGCGACGTCGATGCCTTCGTCTTCGTAACGCGCGCGCAGGCGGTCGAAAATGAACTCCAGTACTTGCTCGGCCAGGCC‑GGGGGCCTTGATCTTGGCGCCGAACTG‑CTTGACC‑GCAAAc                                                                                                                      <  4:165565/149‑1 (MQ=255)
                               gTACCGACAGGTACGTGGCGACGTCGATGCCTTCGTCTTCGTAACGCGCGCGCAGGCGGTCGAAAATGAACTCCAGTACTTGCTCGGCCAGGCC‑GGGGGCCTTGATCTTGGCGCCGAACTG‑CTTGACC‑GCAAACTCAACCGCGCCGGt                                                                                                        >  5:281725/1‑148 (MQ=255)
                               gTACCGACAGGTACGTGGAGACGTCGATGCCTTCGTCTTAGTAACTCGCGCGCAGGCGGTCGAAAATGAACTCCAGTACTTGCTCGGCCAGGCC‑GGGGTCCTTGATCTTGGCGCCGAACTG‑CTTGACC‑GCAAACTCAACCGCGCCGGTc                                                                                                       <  3:79948/149‑1 (MQ=255)
                                   cGACAGGTACGTGGCGACGTCGATGCCTTCGTCTTCGTAACGCGCGCGCAGGCGGTCGAAAATGAACTCCAGTACTTGCTCGGCCAGGCC‑GGGGGCCTTGATCTTGGCGCCGAACTG‑CTTGACC‑GCAAACTCAACCGCGCCGGTcaggt                                                                                                   >  2:173879/1‑149 (MQ=255)
                                     aCAGGTACGTGGCGACGTCGATGCCTTCGTCTTCGTAACGCGCGCGCAGGCGGTCGAAAATGAACTCCAGTACTTGCTCAGCCAGGCC‑GGGGGCCTTGATCTTGGCGCCGAACTG‑CTTGACC‑GCAAACTCAACCGCGCCGGTcaggtcc                                                                                                 <  1:325729/149‑1 (MQ=255)
                                       aGGTACGTGGCGACGTCGATGCCTTCGTCTTCGTAACGCGCGCGCAGGCGGTCGAAAATGAACTCCAGTACTTGCTCGGCCAGGCC‑GGGGGCCTTGATCTTGGCGCCGAACTG‑CTTGACC‑GCAAACTCAACCGCGCCGGTcaggtgcag                                                                                               <  2:131120/149‑1 (MQ=255)
                                       aGGTACGTGGCGACGTCGATGCCTTCGTCTTCGTAACGCGCGCGCAGGCGGTCGAAAATGAACTCCAGTACTTGCTCGGCCAGGCC‑GGGGGCCTTGATCTTGGCGCCGAACTG‑CTTGACC‑GCAAACTCAACCGCGCCGGTcaggtccag                                                                                               <  5:285867/149‑1 (MQ=255)
                                         gTACGTGGCGACGTCGATGCCTTCGTCTTCGTAACGCGCGCGCAGGCGGTCGAAAATGAACTCCAGTACTTGCTCGGCCAGGCC‑GGGGGCCTTGATCTTGGCGCCGAACTG‑CTTGACC‑GCaaa                                                                                                                       <  2:249136/123‑1 (MQ=255)
                                         gTACGTGGCGACGTCGATGCCTTCGTCTTCGTAACGCGCGCGCAGGCGGTCGAAAATGAACTCCAGTACTTGCTCGGCCAGGCC‑GGGGGCCTTGATCTTGGCGCCGAACTG‑CTTGACC‑GCaaa                                                                                                                       >  5:38745/1‑123 (MQ=255)
                                         gTACGTGGCGACGTCGATGCCTTCGTCTTCGTAACGCGCGCGCAGGCGGTCGAAAATGAACTCCAGTACTTGCTCGGCCAGGCC‑GGGGGCCTTGATCTTGGCGCCGAACTG‑CTTGACC‑GCaaa                                                                                                                       <  6:38745/123‑1 (MQ=255)
                                         gTACGTGGCGACGTCGATGCCTTCGTCTTCGTAACGCGCGCGCAGGCGGTCGAAAATGAACTCCAGTACTTGCTCGGCCAGGCC‑GGGGGCCTTGATCTTGGCGCCGAACTG‑CTTGACC‑GCAAACTCAACCGCGCCGGTcaggtccaggt                                                                                             >  5:197741/1‑149 (MQ=255)
                                         gTACGTGGCGACGTCGATGCCTTCGTCTTCGTAACGCGCGCGCAGGCGGACGAAAATGAACTCCAGTACTTGCTCTGCCAGGCC‑GGGGGCCTTGATCTTGGCGCCGAACTG‑CTTGACC‑GCaaa                                                                                                                       >  1:249136/1‑123 (MQ=255)
                                          tACGTGGCGACGTCGATGCCTTCGTCTTCGTAACGCGCGCGCAGGCGGTCGAAAATGAACTCCAGTACTTGCTCGGCCAGGCC‑GGGGGCCTTGATCTTGGCGCCGAACTG‑CTTGACC‑GCAAACTCAACCGCGCCGGTcaggtccaggtc                                                                                            <  1:173879/149‑1 (MQ=255)
                                           aCGTGGCGACGTCGATGCCTTCGTCTTCGTAACGCGCGCGCAGGCGGTCGAAAATGAACTCCAGTACTTGCTCGGCCAGGCC‑GGGGGCCTTGATCTTGGCGCCGAACTG‑CTTGACC‑GCAAACTCAACCGCGCCGGTcaggtccaggtc                                                                                            >  3:210761/1‑148 (MQ=255)
                                             gTGGCGACGTCGATGCCTTCGTCTTCGTAACGCGCGCGCAGGCGGTCGAAAATGAACTCCAGTACTTGCTCGGCCATGCC‑GGGGGCCTTGATCTTGGCGCCGAACTG‑CTTGACC‑GCAAACTCAACCGCGCCGTTCAGgtccatgtccag                                                                                         <  3:289223/149‑1 (MQ=255)
                                             gTGGCGACGTCGATGCCTTCGTCTTCGTAACGCGCGCGCAGGCGGTCGAAAATGAACTCCAGTACTTGCTCGGCCAGGCC‑GGGGGCCTTGATCTTGGCGCCGAACTG‑CTTg                                                                                                                                <  7:302112/111‑1 (MQ=255)
                                             gTGGCGACGTCGATGCCTTCGTCTTCGTAACGCGCGCGCAGGCGGTCGAAAATGAACTCCAGTACTTGCTCGGCCAGGCC‑GGGGGCCTTGATCTTGGCGCCGAACTG‑CTTg                                                                                                                                >  8:302112/1‑111 (MQ=255)
                                                 cGACGTCGATGCCTTCGTCTTCGTAACGCGCGCGCAGGCGGTCGAAAATGAACTCCAGTACTTGCTCGGCCAGGCC‑GGGGGCCTTGATCTTGGCGCCGAACTG‑CTTGACC‑GCAAACTCAACCGCGCCGGTCAGGTCCAGGTCCAGTTGc                                                                                     <  2:227550/149‑1 (MQ=255)
                                                       cGATGCCTTCGTCTTCGTAACGCGCGCGCAGGCGGTCGAAAATGAACTCCAGTACTTGCTCGGCCAGGCA‑GGGGGCCTTGATCTTGGCGCCGAACTG‑CTTGACC‑GCAAACTCAACCGCGCCGGTCAGGTCCAGGTCCAGTTGCTtctc                                                                                >  2:307322/1‑148 (MQ=255)
                                                        gATGCCTTCGTCTTCGTAACGCGCGCGCAGGCGGTCGAAAATGAACTCCAGTACTTGCTCGGCCAGGCC‑GGGGGCCTTGATCTTGGCGCCGAACTG‑CTTGACC‑GCAAACTCAACCGCGCCGGTCAGGTCCAGGTCCAGTTGCTTCTCGa                                                                              >  7:136712/1‑149 (MQ=255)
                                                         aTGCCTTCGTCTTCGTAACGCGCGCGCAGGCGGTCGAAAATGAACTCCAGTACTTGCTCGGCCAGGCC‑GGGGGCCTTGATCTTGGCGCCGAACTG‑CTTGACC‑GCAAACTCAACCGCGCCGGTCAGGTCCAGGTCCAGTTGCTTCTCGa                                                                              >  1:168005/1‑148 (MQ=255)
                                                          tGCCTTCGTCTTCGTAACGCGCGCGCAGGCGGTCGAAAATGAACTCCAGTACTTGCTCGGCCAGGCC‑GGGGGCCTTGATCTTGGCGCCGAACTG‑CTTGACC‑GCAAACTCAACCGCGCCGGTCAGGTCCAGGTCCAGTTGCTTCTCGATc                                                                            <  2:168005/149‑1 (MQ=255)
                                                            ccTTCGTCTTCGTAACGCTCGCGCAGGCGGTCGAAAATGAACTCCAGTACTTGCTCGGCGAGGCC‑GGGGGCATTGATCTTGGCGACGAACTG‑CTTGACC‑GCAAACTCAAGCGCGCCGGTCAGGTCCAGGTCCTGTTGCTTCTCGATCat                                                                          <  4:200065/149‑2 (MQ=255)
                                                            ccTTCGTCTTCGTAACGCGCGCGCAGGCGGTCGAAAATGAACTCCAGTACTTGCTCGGCCAGGCC‑GGTGGCCTTGATCTTGGCGCCGTACTG‑CTTGACC‑GCAAACTCAACCGCGCCGGTCAGGTCCAGGTCCAGTTGCTTCTCGATCAg                                                                          <  7:69312/149‑1 (MQ=255)
                                                            ccTTCGTCTTCGTAACGCGCGCGCAGGCGGTCGAAAATGAACTCCAGTACTTGCTCGGCCAGGCC‑GGGGGCCTTGATCTTGGCGCCGAACTG‑CTTGACC‑GCAAACTCAACCGCGCCGGTCAGGTCCAGGTCCAGTTGCTTCTCGATCa                                                                           <  4:161135/148‑1 (MQ=255)
                                                            ccTTCGTCTTCGTAACGCGCGCGCAGGCGGTCGAAAATGAACTCCAGTACTTGCTCGGCCAGGCC‑GGGGGCCTTGATCTTGGCGCCGAACTG‑CTTGACC‑GCAAACTCAACCGCGCCGGTCAGGTCCAGGTCCAGTTGCTTCTCGATCAg                                                                          <  7:106304/149‑1 (MQ=255)
                                                              ttcgtcttcgtAACGCGCGCGCAGGCGGTCGAAAATGAACTCCAGTACTTGCTCGGCCAGGCC‑GGGGGCCTTGATCTTGGCGCCGAACTG‑CTTGACC‑GCAAACTCAACCGCGCCGGTCAGGTCCAGGTCCAGTTGATTCTCTATCAGGa                                                                        >  4:235331/1‑149 (MQ=255)
                                                                  tcttcgtAACGCGCGCGCAGGCGGTCGAAAATGAACTCCAGTACTTGCTCGGCCAGGCC‑GGGGGCCTTGATCTTGGCGCCGAACTG‑CTTGACC‑GCAAACTCACCCGCGCCGGTCAGGTCCAGGTCCAGTTGCTTCTCGATCAGGATACg                                                                    >  4:40052/1‑149 (MQ=255)
                                                                  tcttcgtAACGCGCGCGCAGGCGGTCGAAAATGAACTCCAGTACTTGCTCGGCCAGGCC‑GGGGGCCTTGATCTTGGCGCCGAACTG‑CTTGACC‑GCAAACTCAACCGCGCCGGTCAGGTCCAGGTCCAGTTGCTTCTCGATCAGGATACg                                                                    >  6:140709/1‑149 (MQ=255)
                                                                       gtAACGCGCGCGCAGGCGGTCGAAAATGAACTCCAGTACTTGCTCGGCCAGGCC‑GGGGGCCTTGATCTTGGCGCCGAACTG‑CTTGACC‑GCAAACTCAACCGCGCCGGTCAGGTCCAGGTCCAGTTGCTTCTCGATCAGGATACgcagc                                                                <  8:136712/148‑1 (MQ=255)
                                                                             cgcgcgcAGGCGGTCGAAAATGAACTCCAGTACTTGCTCGGCCAGGCC‑GGGGGCCTTGATCTTGGCGCCGAACTG‑CTTGACC‑GCAAACTCAACCGCGCCGGTCAGGTCCAGGTCCAGTTGCTTCTCGATCAGGATACGCAGCACGCCCa                                                         <  6:281725/149‑1 (MQ=255)
                                                                             cgcgcgcAGGCGGTCGAAAATGAACTCCAGTACTTGCTCGGCCAGGCC‑GGGGGCCTTGATCTTGGCGCCGAACTG‑CTTGACC‑GCAAACTCAACCGCGCCGGTCAGGTCCAGGTCCAGTTGCTTCTCGATCAGGATACGCAGCACGCCCa                                                         <  5:140709/149‑1 (MQ=255)
                                                                                  gcAGGCGGTCGAAAATGAACTCCAGTACTTGCTCGGCCAGGCC‑GGGGGCCTTGATCTTGGCGCCGAACTG‑CTTGACC‑GCAAACTCAACCGCGCCGGTCAGGTCCAGGTCCAGTTGCTTCTCGATCAGGATACGCAGCACGCCCAGagca                                                    <  5:90965/149‑1 (MQ=255)
                                                                                             aaaaTGAACTCCAGTACTTGCTCGGCCAGGCCGGGGGGCCTTGATCTTGGCGCCGAACTG‑CTTGACC‑GCAAACTCAACCGCGCCGGTCAGGTCCAGGTCCAGTTGCTTCTCGATCAGGATACGCAGCACGCCCAGAGCAGCACGgcg                                            >  1:41338/1‑147 (MQ=255)
                                                                                             aaaaTGAACTCCAGTACTTGCTCGGCCAGGCC‑GGGGGCCTTGATCTTGGCGCCGAACTG‑CTTGACC‑GCAAACTCAACCGCGCCGGTCAGGTCCAGGTCCAGTTGCTTCTCGATCAGGATACGCAGCACGCCCAGAGCAGCACGGCGCAg                                         >  7:261352/1‑149 (MQ=255)
                                                                                              aaaTGAACTCCAGTACTTGCTCGGCCAGGCC‑GGGGGCCTTGATCTTGGCGCCGAACTG‑CTTGACC‑GCAAACTCAACCGCGCCGGTCAGGTCCAGGTCCAGTTGCTTCTCGATCAGGATACGCAGCACGCCCAGAGCAGCACGGCGCAgg                                        <  1:165347/149‑1 (MQ=255)
                                                                                                    aCTCCAGTACTTGCTCGGCCAGGCC‑GGGGGCCTTGATCTTGGCGCCGAACTG‑CTTGACC‑GCAAACTCAACCGCGCCGGTCAGGTCCAGGTCCAGTTGCTTCTCGATCAGGATACGCAGCACGCCCAGAGCAGCACGGCGCAGGGCGTAc                                  <  3:192410/149‑1 (MQ=255)
                                                                                                      tCCAGTACTTGCTCGGCCAGGCC‑GGGGGCCTTGATCTTGGCGCCGAACTGAATTGCCC‑GCAAACTCAATCGCGCCGGTCAGGTCATGGTCCAGTTGCTTCTCGATCAGGATACGCAGCACGCCAAGCGCAGCACGGCGCTGGGCGTACg                                 >  1:254503/1‑149 (MQ=255)
                                                                                                      tCCAGTACTTGCTCGGCCAGGCC‑GGGGGCCTTGATCTTGGCGCCGAACTG‑CTTGACC‑GCAAACTCAACCGCGCCGGTCAGGTCCAGGTCCAGTTGCTTCTCGATCAGGATACGCAGCACGCCCAGAGCAGCACGGCGCAGGGCGTACgg                                >  3:55793/1‑149 (MQ=255)
                                                                                                      tCCAGGACTTGCTCGGCCAGGCC‑GGGGGCCTTGATCTTGGCGCCGAACTG‑CTTGACC‑GCAAACTCAACCGCGCCGGTCAGGTCCAGGTCCAGTTGCTTCTCGATCAGGATACGCAGCACGCCCAGAGCAGCACGGCGCAGGGCGTACgg                                >  7:303105/1‑149 (MQ=255)
                                                                                                        cAGTACTTGCTCGGCCAGGCC‑GGGGGCCTTGATCTTGGCGCCGAACTG‑CTTGACC‑GCAAACTCAACCGCGCCGGTCAGGTCCAGGTCCAGTTGCTTCTCGATCAGGATACGCAGCACGCCCAGAGCAGCACGGCGCAGGGCGTACGGGt                              >  7:312023/1‑149 (MQ=255)
                                                                                                        cAGTACTTGCTCGGCCAGGCC‑GGGGGCCTTGATCTTGGCGCCGAACTG‑CTTGACC‑GCAAACTCAACCGCGCCGGTCAGGTCCAGGTCCAGTTGCTTCTCGATCAGGATACGCAGCACGCCCAGAGCAGCACGGCACAGGGCGTACgg                                <  4:210761/147‑1 (MQ=255)
                                                                                                             cTTGCTCGGCCAGGCC‑GGGGGCCTTGATCTTGGCGCCGAACTG‑CTTGACC‑GCAAACTCAACCGCGCCGGTCAGGTCCAGGTCCAGTTGCTTCTCGATCAGGATACGCAGCACGCCCAGAGCAGCACGGCGCAGGGCGTACGGGTCCtt                          <  1:307322/148‑1 (MQ=255)
                                                                                                              ttGCTCGGCCAGGCC‑GGGGGCCTTGATCTTGGCGCCGAACTG‑CTTGACC‑GCAAACTCAACCGCGCCGGTCAGGTCCAGGTCCAGTTGCTTCTCTATCAGGATACGCAGCACGCCCAGAGCAGCACGGCGCAGGGCGTACGGGTCCTtgc                        >  2:167685/1‑149 (MQ=255)
                                                                                                                  tCGGCCAGGCC‑GGGGGCCTTGATCTTGGCGCCGAACTG‑CTTGACC‑GCAAACTCAACCGCGCCGGTCAGGTCCAGGTCCAGTTGCTTCTCGATCAGGATACGCAGCACGCCCAGAGCAGCACGGCGCAGGGCGTACGGGTCCTtgctgc                     >  3:191152/1‑148 (MQ=255)
                                                                                                                  tCGGCCAGGCC‑GGGGGCCTTGATCTTGGCGCCGAACTG‑CTTGACC‑GCAAACTCAACCGCGCCGGTCAGGTCCAGGTCCAGTTGCTTCTCGATCAGGATACGCAGCACGCCCAGAGCAGCACGGCGCAGGGCGTACGGGTCCTTGCTGcc                    >  8:255566/1‑149 (MQ=255)
                                                                                                                    ggCCAGGCC‑GGGGGCCTTGATCTTGGCGCCGAACTG‑CTTGACC‑GCAAACTCAACCGCGCCGGTCAGGTCCAGGTCCAGTTGCTTCTCGATCAGGATACGCAGCACGCCCAGAGCAGCACGGCGCAGGGCGTACGGGTCCTTGCTGCCgg                  <  8:303105/149‑1 (MQ=255)
                                                                                                                    ggCCAGGCC‑GGGGGCCTTGATCTTGGCGCCGAACTG‑CTTGACC‑GCAAACTCAACCGCGCCGGTCAGGTCCAGGTCCAGTTGCTTCTCGATCAGGATACGCAGCACGCCCAGAGCAGCACGGCGCAGGGCGTACGGGTCCTTGCTGCCg                   >  6:198868/1‑148 (MQ=255)
                                                                                                                    ggCCAGGCC‑GGGGGCCTTGATCTTGGCGCCGAACTG‑CTTGACC‑GCAAACTCAACCGCGCCGGTCAGGTCCAGGTCCAGTTGCTTCTCGATCAGGATACGCAGCACGCCCAGAGCAGCACGGCGCAGGGCGTACGGGTCCTTGCTGCCg                   <  5:198868/148‑1 (MQ=255)
                                                                                                                         ggCC‑GGGGGCCTTGATCTTGGCGCCGAACTG‑CTTGACC‑GCAAACTCAACCGCGCCGGTCAGGTCCAGGTCCAGTTGCTTCTCGATCAGGATACGCAGCACGCCCAGAGCAGCACGGCGCAGGGCGTACGGGTCCTTGCTGCCGGTAGgc             <  4:55793/149‑1 (MQ=255)
                                                                                                                           cc‑GGGGGCCTTGATCTTGGCGCCGAACTG‑CTTGACC‑GCAAACTCAACCGCGCCGGTCAGGTCCAGGTCCAGTTGCTTCTCGATCAGGATACGCAGCACGCCCAGAGCAGCACGGCGCAGGGCGTACGGGTCCTTGCTGCCAGTAGgcag           <  4:191152/149‑1 (MQ=255)
                                                                                                                            c‑GGGGGCCTTGATCTTGGCGCCGAACTG‑CTTGACC‑GCAAACTCAACCGCGCCGGTCAGGTCCAGGTCCAGTTGCTTCTCGATCAGGATACGCAGCACGCCCAGAGCAGCACGGCGCAGGGCGTACGGGTCCTTGCTGCCGGTAGgcagc          <  6:197741/149‑1 (MQ=255)
                                                                                                                               ggggCCTTGATCTTGGCGCCGAACTG‑CTTGACC‑GCAAACTCAACCGCGCCGGTCAGGTCCAGGTCCAGTTGCTTCTCGATCAGGATACGCAGCACGCCCAGAGCAGCACGGCGCAGGGCGTACGGGTCCTTGCTGCCGGTAGGCAGCAt        <  2:41338/149‑1 (MQ=255)
                                                                                                                                  gCCTTGATCTTGGCACCGAACTG‑CTTGACC‑GCAAACTCAACCGCGCCGGTCAGGTCCAGGTCCAGTTGCTTCTCGATCAGGATACGCAGCACGCCCAGAGCAGCACGGCGCAGGGCGTACGGGTCCTtgctgc                     >  6:250630/1‑133 (MQ=255)
                                                                                                                                  gCCTTGATCTTGGCACCGAACTG‑CTTGACC‑GCAAACTCAACCGCGCCGGTCAGGTCCAGGTCCAGTTGCTTCTCGATCAGGATACGCAGCACGCCCAGAGCAGCACGGCGCAGGGCGTACGGGTCCTtgctgc                     >  6:245967/1‑133 (MQ=255)
                                                                                                                                  gCCTTGATCTTGGCACCGAACTG‑CTTGACC‑GCAAACTCAACCGCGCCGGTCAGGTCCAGGTCCAGTTGCTTCTCGATCAGGATACGCAGCACGCCCAGAGCAGCACGGCGCAGGGCGTACGGGTCCTtgctgc                     <  5:245967/133‑1 (MQ=255)
                                                                                                                                  gCCTTGATCTTGGCACCGAACTG‑CTTGACC‑GCAAACTCAACCGCGCCGGTCAGGTCCAGGTCCAGTTGCTTCTCGATCAGGATACGCAGCACGCCCAGAGCAGCACGGCGCAGGGCGTACGGGTCCTtgctgc                     <  5:250630/133‑1 (MQ=255)
                                                                                                                                     ttGATCTTGGCGCCGAACTG‑CTTGACCGGCAAACTCAACCGCGCCGGTCAGGTCCAGGTCCAGTTGCTTCTCGATCAGGATACGCAGCACGCCCAGAGCAGCACGGCGCAGGGCGTACGGGTCCTTGCTGCCGGTAGGCAGCATGCCGa   >  1:19683/1‑149 (MQ=255)
                                                                                                                                     ttGATCTTGGCGCCGAACTG‑CTTGACC‑GCAAACTCAACCGCGCCGGTCAGGTCCAGGTCCAGTTGCTTCTCGATCAGGATACGCAGCACGCCCAGAGCAGCACGGCGCAGGGCGTACGGGTCCTTGCTGCCGGTAGGCAGCATGCCGAt  >  5:29157/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                           
AAGTCCAGGGCGGAGCCTGGCTGCAGTGCGCGTACCGACAGGTACGTGGCGACGTCGATGCCTTCGTCTTCGTAACGCGCGCGCAGGCGGTCGAAAATGAACTCCAGTACTTGCTCGGCCAGGCC‑GGGGGCCTTGATCTTGGCGCCGAACTG‑CTTGACC‑GCAAACTCAACCGCGCCGGTCAGGTCCAGGTCCAGTTGCTTCTCGATCAGGATACGCAGCACGCCCAGAGCAGCACGGCGCAGGGCGTACGGGTCCTTGCTGCCGGTAGGCAGCATGCCGAT  >  NC_002947/70241‑70521

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