Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,026,750	A→T	12.5%	T20T (ACA→ACT)	ocd →	ornithine cyclodeaminase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,026,750	0	A	T	12.5%	146.7 / 11.7	56	T20T (ACA→ACT)	ocd	ornithine cyclodeaminase
Reads supporting (aligned to +/- strand):  ref base A (28/21);  new base T (2/5);  total (30/26)
Fisher's exact test for biased strand distribution p-value = 2.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.39e-01

GGTGGGAATGTCGACATGGCTCAACACCAACGCATCATCCTGGAAGCTGTGCGCGACCGCCGCCCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGTCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTCCGAGCACAATGGCGAGGTGGACGTCAAAACCGCTTACCTGCCAGACCTTGAGAATTTCG  >  NC_002947/5026602‑5026895                                                                                                                                                     |                                                                                                                                                  ggTGGGAATGTCGACATGGCTCAACACCAACGCATCATCCTGGAAGCTGTGCGCGACCGCCGCCCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGTCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATaca                                                                                                                                                   <  2:267839/149‑1 (MQ=255) ggTGGGAATGTCGACATGGCTCAACACCAACGCATCATCCTGGAAGCTGTGCGCGACCGCCGCCCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGTCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATaca                                                                                                                                                   <  4:150302/149‑1 (MQ=255)         tGTCGACATGGCTCAACACCAACGCATCATCCTGGAAGCTGTGCGCGACCGCCGCCCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGTCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTAt                                                                                                                                           <  8:231848/149‑1 (MQ=255)             gACATGGCTCATCACCAACGCATCATCCTGGAAGCTGTGCGCGACCGCCGCCCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGTCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGAt                                                                                                                                       <  4:172156/149‑1 (MQ=255)             gACATGGCTCAACACCAACGCATCATCCTGGAAGCTGTGCGCGACCGCCGCCCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGTCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACTGTCTGTATCGa                                                                                                                                        <  5:100979/148‑1 (MQ=255)             gACATGGCTCAACACCAACGCATCATCCTGGAAGCTGTGCGCGACCGCCGCCCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGTCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACTGTCTGTATCGa                                                                                                                                        <  2:161553/148‑1 (MQ=255)              aCATGGCTCAACACCAACGCATCATCCTGGAAGCTGTGCGCGACCGCCGCCCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGTCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACTGTCTGTATCGagc                                                                                                                                      <  8:46853/149‑3 (MQ=255)                 tGGCTCAACACCAACGCATCATCCTGGAAGCTGTGCGCGACCGCCGCCCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGTCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAg                                                                                                                                   >  6:62027/1‑149 (MQ=255)                     tCAACACCAACGCATCATCCTGGAAGCTGTGCGCGACCGCCGCCCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGTCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACTGTCTGTATCGa                                                                                                                                        >  4:166111/1‑140 (MQ=255)                     tCAACACCAACGCATCATCCTGGAAGCTGTGCGCGACCGCCGCCCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGTCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACTGTCTGTATCGa                                                                                                                                        <  3:166111/140‑1 (MQ=255)                          acCAACGCATCATCCTGGAAGCTGTGCGCGCCCGCCGCCCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGTCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGCTGTAGTGGAAAACg                                                                                                                          >  5:131660/1‑149 (MQ=255)                          acCAACGCATCATCCTGGAAGCTGTGCGCGACCGCCGCCCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGTCCTTCTGCTTAACCACCGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACg                                                                                                                          >  7:84196/1‑149 (MQ=255)                            cAACGCATCATCCTGGAAGCTGTGCGCGACCGCCGCCCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGTCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATacagac                                                                                                                                                >  5:55876/1‑125 (MQ=255)                            cAACGCATCATCCTGGAAGCTGTGCGCGACCGCCGCCCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGTCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATacagac                                                                                                                                                <  6:55876/125‑1 (MQ=255)                              aCGCATCATCCTGGAAGCTGTGCGCGACCGCCGCCCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGTCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATacagaca                                                                                                                                               >  5:152829/1‑124 (MQ=255)                              aCGCATCATCCTGGAAGCTGTGCGCGACCGCCGCCCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGTCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATacagaca                                                                                                                                               <  6:152829/124‑1 (MQ=255)                                   tcatcCTGGAAGCTGTGCGCGACCGCCGCCCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGTCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCAt                                                                                                                 <  1:327817/149‑1 (MQ=255)                                     atcCTGGAAGCTGTGCGCGACCGCCGCCCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGTCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCAttt                                                                                                                     >  4:138509/1‑143 (MQ=255)                                     atcCTGGAAGCTGTGCGCGACCGCCGCCCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGTCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCAttt                                                                                                                     <  3:138509/143‑1 (MQ=255)                                         tGGAAGCTGTGCGCGACCGCCGCCCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGTCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTgg                                                                                                           >  8:326784/1‑149 (MQ=255)                                            aaGCTGTGCGCGACCGCCGCCCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGTCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCa                                                                                                        <  7:326784/149‑1 (MQ=255)                                                  agcgcGACCGCCGCCCAGCACGTGCTTAAAGCGGTGCCATGTGAGCGAAGTCCTTCTGGTTAAGCAACGCGAAATCCGTGAGTGCGTAACGCGCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCa                                                                                                  <  6:131660/148‑1 (MQ=255)                                                        aCCGCCGCCCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGTCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTgg                                                                                            <  6:221081/149‑1 (MQ=255)                                                        aCCGCCGCCCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGTCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTgg                                                                                            >  5:221081/1‑149 (MQ=255)                                                            ccgccCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGTCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCAt                                                                                        >  3:70262/1‑149 (MQ=255)                                                            ccgccCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGTCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCAt                                                                                        >  3:157381/1‑149 (MQ=255)                                                                ccAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGTCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGATCATGCCg                                                                                    >  8:41197/1‑149 (MQ=255)                                                                 cAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGTCCTTCTGCTTAACCACCGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGt                                                                                   <  8:84196/149‑1 (MQ=255)                                                                              aaGAGGTGCCATGTGAGCGAAGTCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGt                                                                       >  5:368044/1‑148 (MQ=255)                                                                                  ggTGCCATGTGAGCGAAGTCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGATCATGCCGTCAATCCTCCGGTTAGAc                                                                  <  7:41197/149‑1 (MQ=255)                                                                                      ccATGTGAGCGAAGTCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACAtct                                                              >  7:61630/1‑149 (MQ=255)                                                                                         tgtgAGCGAAGTCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTcc                                                            >  4:323919/1‑148 (MQ=255)                                                                                         tgtgAGCGAAGTCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTcc                                                            >  8:120267/1‑148 (MQ=255)                                                                                           tgAGCGAAGTCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTCCGAg                                                         <  8:61630/149‑1 (MQ=255)                                                                                                  aGTCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTCCGAGCACAATg                                                  <  4:157381/149‑1 (MQ=255)                                                                                                    tCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTCCGAGCACAATGGc                                                >  1:364865/1‑149 (MQ=255)                                                                                                       ttCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACAtct                                                              >  1:195756/1‑132 (MQ=255)                                                                                                       ttCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACAtct                                                              <  2:195756/132‑1 (MQ=255)                                                                                                           gCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTCCGAGCACAATGGCGAGGTg                                          >  5:332034/1‑148 (MQ=255)                                                                                                             ttAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTCCGAGCACAATGGCGAGGTGGAc                                       >  7:120737/1‑149 (MQ=255)                                                                                                               aaCCAACGCGAAATCCGTGAGTGCTTAAAGCTCGATACAGACTGTATCGATGTAGTGGAAAACGCATTTGCATCGCTTGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTCCGAGCACAATGGCGAGGTGGAc                                       <  7:253134/147‑1 (MQ=255)                                                                                                                 ccAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTcc                                                            <  4:135734/124‑1 (MQ=255)                                                                                                                 ccAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTcc                                                            >  3:135734/1‑124 (MQ=255)                                                                                                                 ccAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTCCGAGCACAATGGCGAGGTGGACGTc                                    >  8:81162/1‑148 (MQ=255)                                                                                                                      gcgAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTCCGAGCACAATGGCGAGGTGGACGTCAAAAcc                              >  7:67564/1‑149 (MQ=255)                                                                                                                        gAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTCCGAGCACAATGGCGAGGTGGACGTCAAAACCGc                            <  5:62027/149‑1 (MQ=255)                                                                                                                        gAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGACGTCATGCCGTCAATCCTCCGGTTAGACATCTCCGAGCACAATGGCGAGGTGGACGTCAAAACCGc                            >  4:305830/1‑149 (MQ=255)                                                                                                                          aaTCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTCCGAGCACAATGGCGAGGTGGACGTCAAAACCGCt                           >  3:19299/1‑148 (MQ=255)                                                                                                                              cGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTCCGAGCACAATGGCGAGGTGGACGTCAAAACCGCTTACCt                      >  8:208887/1‑149 (MQ=255)                                                                                                                                tGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTCCGAGCACAATGGCGAGGTGGACGTCAAAACCGCTTACCTGc                    <  7:208887/149‑1 (MQ=255)                                                                                                                                  aGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTCCGAGCACAATGGCGAGGTGGACGTCAAAACCGCTTACCTGCCa                  <  8:120737/149‑1 (MQ=255)                                                                                                                                       cactacatCGATACTGTCTGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTCCGAGCACAATGGCGAGGTGGACGTCAAAACCGCTTACCTGCCAGACCt             >  3:29855/8‑149 (MQ=255)                                                                                                                                         ctacatCGATACTGTCTGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTCCGAGCACAATGGCGAGGTGGACGTCAAAACCGCTTACCTGCCAGACCTTg           <  4:29855/144‑1 (MQ=255)                                                                                                                                           cGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTCCGAGCACAATGGCGAGGTGGACGTCAAAACCGCTTACCTGCCAGACCtt            >  2:298278/1‑146 (MQ=255)                                                                                                                                           cGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTCCGAGCACAATGGCGAGGTGGACGTCAAAACCGCTTACCTGCCAGACCTTgag         <  4:19299/149‑1 (MQ=255)                                                                                                                                            gCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTCCGAGCACAATGGCGAGGTGGACGTCAAAACCGCTTACCTGCCAGACCTTgag         >  7:95747/1‑148 (MQ=255)                                                                                                                                              tCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTCCGAGCACAATGGCGAGGTGGACGTCAAAACCGCTTACCTGCCAGACCTTGAGAAt      >  3:52316/1‑149 (MQ=255)                                                                                                                                                gATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTCCGAGCACTATGGCGAGGTGGACGTCAAAACCGCTTACCTGCCAGACCTTGAGAAttt    >  8:105516/1‑149 (MQ=255)                                                                                                                                                gATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTCCGAGCACAATGGCGAGGTGGACGTCAAAACCGCt                           <  8:25830/126‑1 (MQ=255)                                                                                                                                                gATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTCCGAGCACAATGGCGAGGTGGACGTCAAAACCGCt                           >  7:25830/1‑126 (MQ=255)                                                                                                                                                 aTACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTCCGAGCACAATGGCGAGGTGGACGTCAAAACCGCTTACCTGCCAGACCTTGAGAATTTc   <  8:95747/149‑1 (MQ=255)                                                                                                                                                  tACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTCCGAGCACAATGGCGAGGTGGACGTCAAAACCGCTTACCTGCCAGACCTTGAGAATTTCg  <  3:323919/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                  GGTGGGAATGTCGACATGGCTCAACACCAACGCATCATCCTGGAAGCTGTGCGCGACCGCCGCCCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGTCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTCCGAGCACAATGGCGAGGTGGACGTCAAAACCGCTTACCTGCCAGACCTTGAGAATTTCG  >  NC_002947/5026602‑5026895

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