Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	71,801	A→T	10.1%	L35Q (CTG→CAG)	glyS ←	glycine‑tRNA ligase subunit beta

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	71,801	0	A	T	10.1%	189.4 / 10.3	69	L35Q (CTG→CAG)	glyS	glycine‑tRNA ligase subunit beta
Reads supporting (aligned to +/- strand):  ref base A (34/28);  new base T (5/2);  total (39/30)
Fisher's exact test for biased strand distribution p-value = 6.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-01

AAGGCAGCCTGCAGGGGCGGGCCGTCGATGTTGATGCTGCGGTCAGGTTGCTGCACGTCCAGCTGGCGAATCAGCACGGCCAGACGACGCGGCGCGGCGTACACCTGCTTGCCGGTGTAGTTCAGGCCAGCAGCCTGCAGGCCTTTTTCGATACCGGCCAGGAAGGCGTCGCCAAGGCTGGCGAGGGCCTTGGGTGGCAGCTCTTCGGTGCCCAGTTCTACCAGGAAATCTTGAGCACTCATTGTGCAGCCTCCAGCTTAGCCAACACTTCGTCACGCAGTTCAGGG  >  NC_002947/71663‑71949                                                                                                                                           |                                                                                                                                                     aaGGCAGCCTGCAGGGGCGGGCCGTCGATGTTGATGCTGCGGTCAGGTTGCTGCACGTCCAGCTGGCGAATCAGCACGGCCAGACGACGCGGCGCGGCGTACACCTGCTTGCCGGTGTAGTTCAGGCCAGCAGCCTGCAGGCCttttt                                                                                                                                             <  7:304937/148‑1 (MQ=255)    gCAGCCTGCAGGGGCGGGCCGTCGATGTTGATGCTGCGGTCAGGTTGCTGCACGTCCAGCTGGCGAATCAGCACGGCCAGACGACGCGGCGCGGCGTACACCTGCTTGCCGGTGTAGTTCAGGCCTGCAGGCTGCTGGCCt                                                                                                                                                 <  4:203101/141‑1 (MQ=255)    gCAGCCTGCAGGGGCGGGCCGTCGATGTTGATGCTGCGGTCAGGTTGCTGCACGTCCAGCTGGCGAATCAGCACGGCCAGACGACGCGGCGCGGCGTACACCTGCTTGCCGGTGTAGTTCAGGCCTGCAGGCTGCTGGCCt                                                                                                                                                 >  3:203101/1‑141 (MQ=255)           gCAGGGGCGGGCCGTCGATGTTGATGCTGCGGTCAGGTTGCTGCACGTCCAGCTGGCGAATCAGCACGGCCAGACGACGCGGCGCGGCGTACACCTGCTTGCCGGTGTAGTTCAGGCCAGCAGCCTGCAGGCCTTTTTCGATACCGGc                                                                                                                                   >  4:280568/1‑148 (MQ=255)              ggggCGGGCCGTCGATGTTGATGCTGCGGTCAGGTTGCTGCACGTCCAGCTGGCGAATCAGCACGGCCAGACGACGCGGCGCGGCGTACACCTGCTTGCCGGTGTAGTTCAGGCCAGCAGCCTGCAGGCCTTTTTCGATACCGGCCAg                                                                                                                                >  8:28812/1‑148 (MQ=255)               gggcgggcCGTCGATGTTGATGCTGCGGTCAGGTTGCTGCACGTCCAGCTGGCGAATCAGCACGGCCAGACGACGCGGCGCGGCGTACACCTGCTTGCCGGTGTAGTTCAGGCCAGCAGCCTGCAGGCCTTTTTCGATACCGGCCAgg                                                                                                                               >  5:283228/1‑148 (MQ=255)                ggcgggcCGTCGATGTTGATGCTGCGGTCAGGTTGCTGCACGTCCAGCTGGCGAATCAGCACGGCCAGACGACGCGGCGCGGCGTACACCTGCTTGCCGGTGTAGTTCAGGCCAGCAGCCTGCAGGCCTTTTTCGATACCGGCCAGGaa                                                                                                                             <  7:216321/149‑1 (MQ=255)                    ggcCGTCGATGTAGATGCTGCGGTCAGGTTGCTGCACGTCCAGCTGGCGAATCAGCACGGCCAGACGACGCGGCGCGCCGTACACCTGCTTGCCGGTGTAGTTCAGGCCAGCAGCCTGCAGGCCTTTTTCGATACCGGCCAGGAAGGCg                                                                                                                         >  1:112195/1‑149 (MQ=255)                        gTCGATGTTGATGCTGCGGTCAGGTTGCTGCACGTCCAGCTGGCGAATCAGCACGGCCAGACGACGCGGCGCGGCGTACACCTGCTTGCCGGTGTAGTTCAGGCCAGCAGCCTGCAGGCCTTTTTCGATACCGGCCAGGAAGGCGTCGc                                                                                                                     >  8:229166/1‑149 (MQ=255)                               ttGATGCTGCGGTCAGGTTGCTGCACGTCCAGCTGGCGAATCAGCACGGCCAGACGACGCGGCGCGGCGTACACCTGCTTGCCGGTGTAGTTCAGGCCAGCAGCCTGCAGGCCTTTTTCGATACCGGCCAGGAAGGCGTCGCCagcctt                                                                                                              >  3:30986/1‑144 (MQ=255)                                tGATGCTGCGGTCAGGTTGCTGCACGTCCAGCTGGCGAATCAGCACGGCCAGACGACGCGGCGCGGCGTACACCTGCTTGCCGGTGTAGTTCAGGCCAGCAGCCTGCAGGCCTTTTTCGATACCGGCCAGGAAGGCGTCGCCAAGGCTg                                                                                                             >  1:1911/1‑149 (MQ=255)                                      tgcGGTCAGGTTGCTGCACGTCCAGCTGGCGAATCAGCACGTCCAGACGACGCGGCGCGGCGTACACCTGCTTGCCGGTGTAGTTCAGGCCAGCATCCTGCAGGCCTTTTTCGATAACGTCCAGGAAGGCGTCTCCAAGGCTGGCGAgg                                                                                                       <  1:320641/149‑1 (MQ=255)                                       gcGGTCAGGTTGCTGCACGTCCAGCTGGCGAATCAGCACGGCCAGACGACGCGGCGCGGCGTACACCTGCTTGCCGGTGTAGTTCAGGCCAGCAGCCTGCAGGCCTGTTTCGATat                                                                                                                                       >  4:167952/1‑115 (MQ=255)                                        cGGTCAGGTTGCTGCACGTCCAGCTGGCGAATCAGCACGGCCAGACGACGCGGCGCGGCGTACACCTGCTTGCCGGTGTAGTTCAGGCCAGCAGCCTGCAGGCCTTTTTCGATACCGGCCAGGAAGGCGTCGCCAAGGCTGGCGAggg                                                                                                      <  6:283228/148‑1 (MQ=255)                                            cAGGTTGCTGCACGTCCAGCTGGCGAATCAGCACGGCCAGACGACGCGGCGCGGCGTACACCTGCTTGCCGGTGTAGTTCAGGCCAGCAGCCTGCAGGCCTTTTTCGATACCGGCCAGGAAGGCGTCGCCAAGGCTGGCGAGGGCCTTg                                                                                                 >  7:111054/1‑149 (MQ=255)                                               gTTGCTGCACGTCCAGCTGGCGAATCAGCACGGCCAGACGACGCGGCGCGGCGTACACCTGCTTGCCGGTGTAGTTCAGGCCAGCAGCCTGCAGGCCTTTTTCGATACCGGCCAGGAAGGCGTCGCCAAGGCTGGCGAGGGCCTTgg                                                                                                >  7:312763/1‑147 (MQ=255)                                               gTTGCTGCACGTCCAGCTGGCGAATCAGCACGGCCAGACGACGCGGCGCGGCGTACACCTGCTTGCCGGTGTAGTTCAGGCCAGCAGCCTGCAGGCCTTTTTCGATACCGGCCAGGAAGGCGTCGCCAAGGCTGGCGAGGGCCTTGGGt                                                                                              <  2:112195/149‑1 (MQ=255)                                               gTTGCTGCACGTCCAGCTGGCGAATCAGCACGGCCAGACGACGCGGCGCGGCGTACACCTGCTTGCCGGTGTAGTTCAGGCCAGCAGCCTGCAGGCCTTTTTCGATACCGGCCAGGAAGGCGTCGCCAAGGCTGGCGAGGGCCTTGGGt                                                                                              <  1:289490/149‑1 (MQ=255)                                                  gctgcACGTCCAGCTGGCGAATCAGCACGGCCAGACGACGCGGCGCGGCGTACACCTGCTTGCCGGTCTAGTTCAGGCCAGCAGCCTGCAGGCCTTTTTCGATACCGGCCAGGAAGGCGTCGCCACGGCTGGCGAGGGCCTTGGGTgg                                                                                            >  7:301710/1‑148 (MQ=255)                                                     gcACGTCCAGCTGGCGAATCAGCACGGCCAGACGACGCGGCGCGGCGTACACCTGCTTGCCGGTGTAGTTCAGGCCTGCAGCCTGCAGGCCTTTTTCGATACCGGCCAGGAAGGCGTCGc                                                                                                                     >  7:257186/1‑120 (MQ=255)                                                     gcACGTCCAGCTGGCGAATCAGCACGGCCAGACGACGCGGCGCGGCGTACACCTGCTTGCCGGTGTAGTTCAGGCCTGCAGCCTGCAGGCCTTTTTCGATACCGGCCAGGAAGGCGTCGc                                                                                                                     <  8:257186/120‑1 (MQ=255)                                                     gcACGTCCAGCTGGCGAATCAGCACGGCCAGACGACGCGGCGCGGCGTACACCTGCTTGCCGGTGTAGTTCAGGCCAGCAGCCTGCAGGCCTTTTTCGATACCGGCCAGGAAGGCGTCGCCAAGGCTGGCGAGGGCCTTGGGTGGCAGc                                                                                        >  7:65731/1‑149 (MQ=255)                                                      cACGTCCAGCTGGCGAATCAGCACGGCCAGACGACGCGGCGCGGCGTACACCTGCTTGCCGGTGTAGTTCAGGCCAGCAGCCTGCAGGCCTTTTTCGATACCGGCCAGGAAGGCGTCGCCAAGGCTGGCGAGGGCCTTGGGTGGCAGc                                                                                        <  7:23976/148‑1 (MQ=255)                                                          tCCAGCTGGCGACTCAGCACGGCCAGACGACGCGGCGCGGCGTACACCTGCTTGCCGGTGTAGTTCAGGCCAGCAGCCTGCAGGCCTTTTTCGATACCGGCCAGGAAGGCGTCGCCAAGGCTGGCGAGGGCCTTGGGTGGCAGCTCTTc                                                                                   <  2:1911/149‑1 (MQ=255)                                                            cAGCTGGCGAATCAGCACGGCCAGACGACGCGGCGCGGCGTACACCTGCTTGCCGGTGTAGTTCAGGCCAGCAGCCTGCAGGCCTTTTTCGATACCGGCCAGGAAGGCGTCGCCAAGGCTGGCGAGGGCCTTGGGTGGCAGCTCTTCg                                                                                  <  8:111054/148‑1 (MQ=255)                                                             aGCTGGCGAATCAGCACGGCCAGACGACGCGGCGCGGCGTACACCTGCTTGCCGGTGTAGTTCAGGCCAGCAGCCTGCAGGCCTTTTTCGATACCGGCCAGGAAGGCGTCGCCAAGGCTGGCGAGGGCCTTGGGTGGCAGCTCTTCGGt                                                                                <  2:244808/149‑1 (MQ=255)                                                             aGCTGGCGAATCAGCACGGCCAGACGACGCGGCGCGGCGTACACCTGCTTGCCGGTGTAGTTCAGGCCAGCAGCCTGCAGGCCTTTTTCGATACCGGCCAGGAAGGCGTCGCCAAGGCTGGCGAGGGCCTTGGGTGGCAGCTCTTCGGt                                                                                >  5:27377/1‑149 (MQ=255)                                                             aGCTGGCGAATCAGCACGGCCAGACGACGCGGCGCGGCGTACACCTGCTTGCCGGTGTAGTTCAGGCCAGCAGCCTGCAGGCCTTTTTCGATACCGGCCAGGAAGGCGTCGCCAAGGCTGGCGAGGGCCTTGGGTGGCAGCTCTTCGGt                                                                                >  5:51400/1‑149 (MQ=255)                                                                  gCGAATCAGCACGGCCAGACGACGCGGCGCGGCGTACACCTGCTTGCCGGTGTAGTTCAGGCCAGCAGCCTGCAGGCCTTTTTCGATACCGGCCAGGAAGGCGTCGCCAAGGCTGGCGAGGGCCTTGGGAGGCAGCTCTTCGGTGccc                                                                            <  8:312763/148‑1 (MQ=255)                                                                   cGAATCAGCACGGCCAGACGACGCGGCGCGGCGTACACCTGCTTGCCGGTGTAGTTCAGGCCAGCAGCCTGCAGGCCTTTTTCGATACCGGCCAGGAAGGCGTCGCCAAGGCTGGCGAGGGCCTTGGGTGGCAGCTCTTCGGTGCCCAg                                                                          >  7:238893/1‑149 (MQ=255)                                                                    gAATCAGCACGGCCAGACGACGCGGCGCGGCGTACACCTGCTTGCCGGTGTAGTTCAGGCCAGCAGCCTGCAGGCCTTTTTCGATACCGGCCAGGAAGGCGTCGCCAAGGCTGGCGAGGGCCTTGGGTGGCAGCTCTTCGGTGCCCAGt                                                                         >  4:211033/1‑149 (MQ=255)                                                                     aaTCAGCACGGCCAGACGACGCGGCGCGGCGTACACCTGCTTGCCGGTGTAGTTCAGGCCAGCAGCCTGCAGGCCTTTTTCGATACCGGCCAGGAAGGCGTCGCCAAGGCTGGCGAGGGCCTTGGGTGGCAGCTCTTCGGTGcc                                                                             <  5:81497/144‑1 (MQ=255)                                                                     aaTCAGCACGGCCAGACGACGCGGCGCGGCGTACACCTGCTTGCCGGTGTAGTTCAGGCCAGCAGCCTGCAGGCCTTTTTCGATACCGGCCAGGAAGGCGTCGCCAAGGCTGGCGAGGGCCTTGGGTGGCAGCTCTTCGGTGcc                                                                             >  6:81497/1‑144 (MQ=255)                                                                       tCAGCACGGCCAGACGACGCGGCGCGGCGTACACCTGCTTGCCGGTGTAGTTCAGGCCAGCAGCCTGCAGGCCTTTTTCGATACCGGCCAGGAAGGCGTCGCCAAGGCTGGCGAGGGCCTTGGGTGGCAGCTCTTCGGTGCCCAGTTc                                                                       >  4:294637/1‑148 (MQ=255)                                                                         aGCACGGCCAGACGACGCGGCGCGGCGTACACCTGCTTGCCGGTGTAGTTCAGGCCAGCAGCCTGCAGGCCTTTTTCGATACCGGCCAGGAAGGCGTCGCCAAGGCTGGCGAGGGCCTTGGGTGGCAGCTCTTCGGTGCCCAGTTCTAc                                                                    <  7:229166/149‑1 (MQ=255)                                                                         aGCACGGCCAGACGACGCGGCGCGGCGTACACCTGCTTGCCGGTGTAGTTCAGGCCAGCAGCCTGCAGGCCTTTTTCGATACCGGCCAGGAAGGCGTCGCCAAGGCTGGCGAGGGCCTTGGGTGGCAGCTCTTCGGTGCCCAGTTCTAc                                                                    <  6:27377/149‑1 (MQ=255)                                                                                   gacgacGCGGCGCGGCGTACACCTGCTTGCCGGTGTAGTTCAGGCCAGCAGCCTGCAGGCCTTTTTCGATACCGGCCAGGAAGGCGTCGCCAAGGCTGGCGAGGGCCTTGGGTGGCAGCTCTTCGGTGCCCAGTTCTACCAGGAAATCt                                                          <  3:280568/149‑1 (MQ=255)                                                                                    acgacgCGGCGCGGCGTACACCTGCTTGCCGGTGTAGTTCAGGCCAGCAGCCTGCAGGCCTTTTTCGATACCGGCCAGGAAGGCGTCGCCAAGGCTGGCGAGGGCCTTGGGTGGCAGCTCTTCGGTGCCCAGTTCTACCAGGAAATCtt                                                         >  8:150343/1‑149 (MQ=255)                                                                                         gcggcgcggcGTACACCTGCTTGCGGGTGTAGTTCAGGCCAGCAGCCTGCAGGCCTTTTTCGATACCGGCCAGGAAGGCGTCGCCAAGGCTGGCGAGGGCCTTGGGTGGCAGCTCTTCGGTGCCCAGTTCTACCAGGAAATCTTGAGCa                                                    >  5:196917/1‑149 (MQ=255)                                                                                             cgcggcgTACACCTGCTTGCCGGTGTAGTTCAGGCCAGCAGCCTGCAGGCCTTTTTCGATACCGGCCAGGAAGGCGTCGCCAAGGCTGGCGAGGGCCTTGGGTGGCAGCTCTTCGGTGCCCAGTTCTACCAGGAAATCTTGAGCACTCa                                                <  7:28812/149‑1 (MQ=255)                                                                                              gcggcgTACACCTGCTTGCCGGTGTAGTTCAGGCCAGCAGCCTGCAGGCCTTTTTCGATACCGGCCAGGAAGGCGTCGCCAAGGCTGGCGAGGGCCTTGGGTGGCAGCTCTTCGGTGCCCAGTTCTACCAGGAAATCTTGAGCACTCAt                                               >  6:236967/1‑149 (MQ=255)                                                                                                     acacCTGCTTGCCGGTGTAGTTCAGGCCAGCAGCCTGCAGGCCTTTTTCGATACCGGCCAGGAAGGCGTCGCCAAGGCTGGCGAGGGCCTTGGGTGGCAGCTCTTCGGTGCCCAGTTCTACCAGGAAATCTTGAGCACTCATTGTGCAg                                        >  7:79532/1‑149 (MQ=255)                                                                                                      cacCTGCTTGCCGGTGTAGTTCAGGCCAGCAGCCTGCAGGCCTTTTTCGATACCGGCCAGGAAGGCGTCGCCAAGGCTGGCGAGGGCCTTGGGTGGCAGCTCTTCGGTGCCCAGTTCTACCAGGAAATCTTGAGCACTCATTGTGCAGc                                       <  3:294637/149‑1 (MQ=255)                                                                                                      cacCTGCTTGCCGGTGTAGTTCAGGCCAGCAGCCTGCAGGCCTTTTTCGATACCGGCCAGGAAGGCGTCGCCAAGGCTGGCGAGGGCCTTGGGTGGCAGCTCTTCGGTGCCCAGTTCTACCAGGAAATCTTGAGCACTCATTGTGCAGc                                       <  8:238893/149‑1 (MQ=255)                                                                                                      cacCTGCTTGCCGGTGTAGTTCAGGCCAGCAGCCTGCAGGCCTTTTTCGATACCGGCCAGGAAGGCGTCGCCAAGGCTGGCGAGGGCCTTGGGTGGCAGCTCTTCGGTGCCCAGTTCTACCAGGAAATCTTGAGCACTCATTGTGCAGc                                       <  5:236967/149‑1 (MQ=255)                                                                                                       acCTGCTTGCCGGTGTAGTTCAGGCCAGCAGCCTGCAGGCCTTTTTCGATACCGGCCAGGAAGGCGTCGCCAAGGCTGGCGAGGGCCTTgg                                                                                                >  6:100919/1‑91 (MQ=255)                                                                                                       acCTGCTTGCCGGTGTAGTTCAGGCCAGCAGCCTGCAGGCCTTTTTCGATACCGGCCAGGAAGGCGTCGCCAAGGCTGGCGAGGGCCTTgg                                                                                                <  5:100919/91‑1 (MQ=255)                                                                                                       acCTGCTTGCCGGTGTAGTTCAGGCCAGCAGCCTGCAGGCCTTTTTCGATACCGGCCAGGAAGGCGTCGCCAAGGCTGGCGAGGGCCTTGGGTGGCAGCTCTTCGGTGCCCAGTTCTACCAGGAAATCTTGAGCACTCATTGTGCAGcc                                      >  5:191425/1‑149 (MQ=255)                                                                                                         cTGCTTGCCGGTGTAGTTCAGGCCAGCAGCCTGCAGGCCTTTTTCGATACCGGCCAGGAAGGCGTCGCCAAGGCTGGCGAGGGCCTTGGGTGGCAGCTCTTCGGTGCCCAGTTCTACCAGGAAATCTTGAGCACTCATTGTGCAGCCTc                                    >  6:276921/1‑149 (MQ=255)                                                                                                           gCTTGCCGGTGTAGTTCAGGCCAGCAGCCTGCAGGCCTTTTTCGATACCGGCCAGGAAGGCGTCGCCAAGGCTGGCGAGGGCCTTGGGTGGCAGCTCTTCGGTGCCCAGTTCTACCAGGAAATCTTGa                                                       >  3:212697/1‑128 (MQ=255)                                                                                                           gCTTGCCGGTGTAGTTCAGGCCAGCAGCCTGCAGGCCTTTTTCGATACCGGCCAGGAAGGCGTCGCCAAGGCTGGCGAGGGCCTTGGGTGGCAGCTCTTCGGTGCCCAGTTCTACCAGGAAATCTTGa                                                       <  4:212697/128‑1 (MQ=255)                                                                                                             ttGCCGGTGTAGTTCAGGCCAGCAGCCTGCAGGCCTTTTTCGATACCGGCCAGGAAGGCGTCGCCAAGGCTGGCGAGGGCCTTGGGTGGCAGCTCTTCGGTGCCCAGTTCTACCAGGAAATCTTGAGCACTCATTGTGCAGCCTCCAGc                                >  4:258158/1‑149 (MQ=255)                                                                                                               gCCGGTGTAGTTCAGGCCAGCAGCCTGCAGGCCTTTTTCGATACCGGCCAGGAAGGCGTCGCCAAGGCTGGCGAGGGCCTTGGGTGGCAGCTCTTCGGTGCCCAGTTCTACCAGGAAATCTTGAGCACTCATTGTGCAGCCTCCAGc                                <  3:211033/147‑1 (MQ=255)                                                                                                                ccGGTGTAGTTCAGGCCAGCAGCCTGCAGGCCTTTTTCGATACCGGCCAGGAAGGCGTCGCCAAGGCTGGCGAGGGCCTTGGGTGGCAGCTCTTCGGTGCCCAGTTCTACCAGGAAATCTTGAGCACTCATTGTGCAGCCTCCAGCTTa                             <  6:119969/149‑1 (MQ=255)                                                                                                                      tAGTTCAGGCCAGCAGCCTGCAGGCCTTTTTCGATACCGGCCAGGAAGGCGTCGCCAAGGCTGGCGAGGGCCTTGGGTGGCAGCTCTTCGGTGCCCAGTTCTACCAGGAAATCTTGAGCACTCATTGTGCAGCCTCCAGCTTAGCCaa                        >  5:163964/1‑148 (MQ=255)                                                                                                                      tAGTTCAGGCCAGCAGCCTGCAGGCCTTTTTCGATACCGGCCAGGAAGGCGTCGCCAAGGCTGGCGAGGGCCTTGGGTGGCAGCTCTTCGGTGCCCAGTTCTACCAGGAAATCTTGAGCACTCATTGTGCAGCCTCCAGCTTAGCCAac                       <  8:79532/149‑1 (MQ=255)                                                                                                                       gaaaaaGGCCTGCAGGCTGCTGGCCTTTTTCGATACCGGCCAGGAAGGCGTCGCCAAGGCTGGCGAGGGCCTTGGGTGGCAGCTCTTCGGTGCCCAGTTCTACCAGGAAATCTTGAGCACTCAtt                                              <  4:48989/120‑1 (MQ=255)                                                                                                                       gaaaaaGGCCTGCAGGCTGCTGGCCTTTTTCGATACCGGCCAGGAAGGCGTCGCCAAGGCTGGCGAGGGCCTTGGGTGGCAGCTCTTCGGTGCCCAGTTCTACCAGGAAATCTTGAGCACTCATTGTGCAGCCTCCAGCTTAGCCAaca                      >  1:44653/6‑149 (MQ=255)                                                                                                                       gaaaaaGGCCTGCAGGCTGCTGGCCTTTTTCGATACCGGCCAGGAAGGCGTCGCCAAGGCTGGCGAGGGCCGTGGGTGGCAGCTCTTCGGTGCCCAGTTCTACCAGGAAATCTTGAGCACTCAtt                                              >  3:48989/6‑125 (MQ=255)                                                                                                                        aaaaaGGCCTGCAGGCTGCTGGCCTTTTTCGATACCGGCCAGGAAGGCGTCGCCAAGGCTGGCGAGGGCCTTGGGTGGCAGCTCTTCGGTGCCCAGTTCTACCAGGAAATCTTGAGCACTCATTGTGCAGCCTCCAGCTTAGCCAacac                     >  4:16572/5‑149 (MQ=255)                                                                                                                         ttCAGGCCAGCAGCCTGCAGGCCTTTTTCGATACCGGCCAGGAAGGCGTCGCCAAGGCTGGCGAGGGCCTTGGGTGGCAGCTCTTCGGTGCCCAGTTCTACCAGGAAATCTTGAGCACTCAtt                                              >  2:99630/1‑123 (MQ=255)                                                                                                                         ttCAGGCCAGCAGCCTGCAGGCCTTTTTCGATACCGGCCAGGAAGGCGTCGCCAAGGCTGGCGAGGGCCTTGGGTGGCAGCTCTTCGGTGCCCAGTTCTACCAGGAAATCTTGAGCACTCAtt                                              <  1:99630/123‑1 (MQ=255)                                                                                                                         ttCAGGACAGCAGCCTGCAGGCCTTTTTCGATCCCGGCCAGGAAGGCGTCGCCAACGCTGGCGAGGGCCTTGGGTGGCAGCTCTTCGGTGCCCAGTTCTAACAGGAAATCTTGAGCGCTCATTGTGCAGCCTCAAGCTTAGACAaccc                     >  2:251609/1‑146 (MQ=255)                                                                                                                             ggCCTGCAGGCTGCTGGCCTTTTTCGATACCGGCCAGGAAGGCGTCGCCAAGGCTGGCGAGGGCCTTGGGTGGCAGCTCTTCGGTGCCCAGTTCTACCAGGAAATCTTGAGCACTCATTGTGCAGCCTCCAGCTTAGCCAACACTTCGt                >  4:259213/1‑149 (MQ=255)                                                                                                                                  gcagcCTGCAGGCCTTTTTCGATACCGGCCAGGAAGGCGTCGCCAAGGCTGGCGAGGGCCTTGGGTGGCAGCTCTTCGGTGCCCAGTTCTACCAGGAAATCtt                                                         >  6:291614/1‑103 (MQ=255)                                                                                                                                  gcagcCTGCAGGCCTTTTTCGATACCGGCCAGGAAGGCGTCGCCAAGGCTGGCGAGGGCCTTGGGTGGCAGCTCTTCGGTGCCCAGTTCTACCAGGAAATCtt                                                         <  5:291614/103‑1 (MQ=255)                                                                                                                                   cagcCTGCAGGCCTTTTTCGATACCGGCCAGGAAGGCGTCGCCAAGGCTGGCGAGGGCCTTGGGTGGCAGCTCTTCGGTGCCCAGTTCTACCAGGAAATCTTGAGCACTCATTGTGCAGCCTCCAGCTTAGCCAACACTTCGTCACg            >  7:235/1‑147 (MQ=255)                                                                                                                                      ccTGCAGGCCTTTTTCGATACCGGCCAGGAAGGCGTCGCCAAGGCTGGCGAGGGCCTTGGGTGGCAGCTCTTCGGTGCCCAGTTCTACCAGGAAATCTTGAGCACTCATTGTGCAGCCTCCAGCTTAGCCAACACTTCGTCACGCAGtt       <  5:276921/149‑1 (MQ=255)                                                                                                                                        tGCAGGCCTTTTTCGATACCGGCCAGGAAGGCGTCGCCAAGGCTGGCGAGGGCCTTGGGTGGCAGCTCTTCGGTGCCCAGTTCTACCAGGAAATCTTGAGCACTCATTGTGCAGCCTCCAGCTTAGCCAACACTTCGTCACGCAGtt       >  3:120860/1‑147 (MQ=255)                                                                                                                                           aGGCCTTTTTCGATACCGGCCAGGAAGGCGTCGCCAAGGCTGGCGAGGGCCTTGGGTGGCAGCTCTTCGGTGCCCAGTTCTACCAGGAAATCTTGAGCACTCATTGTGCAGCCTCCAGCTTAGCCAACACTTCGTCACGCAGTTCAggg  <  6:196917/149‑1 (MQ=255)                                                                                                                                           |                                                                                                                                                     AAGGCAGCCTGCAGGGGCGGGCCGTCGATGTTGATGCTGCGGTCAGGTTGCTGCACGTCCAGCTGGCGAATCAGCACGGCCAGACGACGCGGCGCGGCGTACACCTGCTTGCCGGTGTAGTTCAGGCCAGCAGCCTGCAGGCCTTTTTCGATACCGGCCAGGAAGGCGTCGCCAAGGCTGGCGAGGGCCTTGGGTGGCAGCTCTTCGGTGCCCAGTTCTACCAGGAAATCTTGAGCACTCATTGTGCAGCCTCCAGCTTAGCCAACACTTCGTCACGCAGTTCAGGG  >  NC_002947/71663‑71949

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