Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,026,752	A→T	5.9%	D21V (GAC→GTC)	ocd →	ornithine cyclodeaminase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,026,752	0	A	T	5.9%	204.5 / 2.5	68	D21V (GAC→GTC)	ocd	ornithine cyclodeaminase
Reads supporting (aligned to +/- strand):  ref base A (32/32);  new base T (2/2);  total (34/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.79e-01

TGGGAATGTCGACATGGCTCAACACCAACGCATCATCCTGGAAGCTGTGCGCGACCGCCGCCCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGTCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTCCGAGCACAATGGCGAGGTGGACGTCAAAACCGCTTACCTGCCAGACCTTGAGAATTTCGCTA  >  NC_002947/5026604‑5026898                                                                                                                                                     |                                                                                                                                                   tGGGAATGTCGACATGGCTCAACACCAACGCATCATCCTGGAAGCTGTGCGCGACCGCCGCCCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGTCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATacaga                                                                                                                                                    <  3:101278/149‑1 (MQ=255)           gACATGGCTCAACACCAACGCATCATCCTGGAAGCTGTGCGCGACCGCCGCCCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGTCCTTCTGCTTAACCAGCGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGAt                                                                                                                                          <  8:153421/149‑1 (MQ=255)             cATGGCTCAACACAAACGCATCATCCTGTAAGCTGTGCGCGACCGCCGCCCTGCAGGTGCTTAATGAGGTGCATTGTGAGCTATGTCTTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACTGTCTGTATCGagcg                                                                                                                                        <  1:115539/149‑4 (MQ=255)              aTGGCTCAACACCAACGCATCATCCTGGAAGCTGTGCGCGACCGCCGCCCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGTCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATg                                                                                                                                         >  7:215849/1‑147 (MQ=255)                 gCTCAACACCAACGCATCATCCTGGAAGCTGTGCGCGACCGCCGCCCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGTCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGt                                                                                                                                     <  6:312941/148‑1 (MQ=255)                 gCTCAACACCAACGCATCATCCTGGAAGCTGTGCGCGACCGCCGCCCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGTCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTg                                                                                                                                    <  3:187908/149‑1 (MQ=255)                   tCAACACCAACGCATCATCCTGGAAGCTGTGCGCGACCGCCGCCCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGTCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTgg                                                                                                                                   >  4:224089/1‑148 (MQ=255)                   tCAACACCAACGCATCATCCTGGAAGCTGTGCGCGACCGCCGCCCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGTCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTgg                                                                                                                                   >  5:145192/1‑148 (MQ=255)                     aaCACCAACGCATCATCCTGGAAGCTGTGCGCGACCGCCGCCCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGTCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGaa                                                                                                                                 >  5:159120/1‑148 (MQ=255)                     aaCACCAACGCATCATCCTGGAAGCTGTGCGCGACCGCCGCCCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGTCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGaa                                                                                                                                 >  1:195376/1‑148 (MQ=255)                     aaCACCAACGCATCATCCTGGAAGCTGTGCGCGACCGCCGCCCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGTCCTTCTGCTTAACCAACGCGAAATCAGTGAGTGCGTAACGCTCGATACAGAAAGTATCGATGTAGTGGaa                                                                                                                                 >  7:266154/1‑148 (MQ=255)                      acacCAACGCATCATCCTGGAAGCTGTGCGCGACCGCCGCCCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGTCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGaaaa                                                                                                                               >  5:372179/1‑149 (MQ=255)                        acCAACGCATCATCCTGGAAGCTGTGCGCGACCGCCGCCCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGTCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACg                                                                                                                             >  4:35821/1‑149 (MQ=255)                        acCAACGCATCATCCTGGAAGCTGTGCGCGACCGCCGCCCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGTCCTTCTGCTTAAAAAACGCGAAATCCGTGAGTGCGTAAAGCTCGATACAGACAGTATCGATGAAGTGGAAAACg                                                                                                                             >  3:327377/1‑149 (MQ=255)                         cgaaCGAATCATCCTGGAATCTGTTCGCGACCGCGGGCCATCTCGTGCTTATAGAGGTGCCATGTGATAGAAGTCGTTCTGGTTAACCAACGCGAAATCCGTGAGTGCGTAACGGTCGATACAGACAGTATCGATGTAGTGGAAAACga                                                                                                                            <  2:18720/147‑2 (MQ=255)                          cAACGCATCATCCTGGAAGCTGTGCGCGACCGCCGCCCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGTCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCa                                                                                                                           >  2:382277/1‑149 (MQ=255)                             cGCATCATCCTGGAAGCTGTGCGCGACCGCCGCCCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGTCCTTCTGCTTAACCACCGAGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCAttt                                                                                                                        <  4:69674/149‑1 (MQ=255)                             cGCATCATCCTGGAAGCTGTGCGCGACCGCCGCCCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGTCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCAttt                                                                                                                        <  6:372179/149‑1 (MQ=255)                             cGCATCATCCTGGAAGCTGTGCGCGACCGCCGCCCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGTCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCAttt                                                                                                                        <  5:416503/149‑1 (MQ=255)                             cGCATCATCCTGGAAGCTGGGCGCGACCGCCGCCCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGTCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCAttt                                                                                                                        <  4:75158/149‑1 (MQ=255)                                   atcCTGGAAGCTGTGCGCGACCGCCGCCCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGTCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCg                                                                                                                  >  8:230582/1‑149 (MQ=255)                                    tcCTGGAAGCTGTGCGCGACCGCCGCCCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGTCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACTGTCTGTATCGa                                                                                                                                           <  3:215203/123‑1 (MQ=255)                                    tcCTGGAAGCTGTGCGCGACCGCCGCCCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGTCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACTGTCTGTATCGa                                                                                                                                           >  4:215203/1‑123 (MQ=255)                                      cTGGAAGCTGTGCGCGACCGCCGCCCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGTCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTg                                                                                                               <  6:148499/149‑1 (MQ=255)                                      cTGGAAGCTGTGCGCGACCGCCGCCCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGTCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTg                                                                                                               <  8:215849/149‑1 (MQ=255)                                       tGGAAGCTGTGCGCGACCGCCGCCCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGTCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTgg                                                                                                              >  5:395740/1‑149 (MQ=255)                                         gAAGCTGTGCGCGACCGCCGCCCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGTCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGc                                                                                                             >  5:141307/1‑148 (MQ=255)                                             cTGTGCGCGACCGCCGCCCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGTCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCagag                                                                                                        <  1:382277/149‑1 (MQ=255)                                              tgtgCGCGACCGCCGCCCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGTCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAgg                                                                                                       >  2:151186/1‑149 (MQ=255)                                                   gcgACCGCCGCCCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGTCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAg                                                                                                  <  3:399988/149‑1 (MQ=255)                                                        cgccgcCCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGTCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTc                                                                                             <  6:395740/149‑1 (MQ=255)                                                          ccgccCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGTCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCAt                                                                                           >  7:82632/1‑149 (MQ=255)                                                          ccgccCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGTCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCAt                                                                                           >  2:268026/1‑149 (MQ=255)                                                          ccgccCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGTCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCAt                                                                                           <  7:230582/149‑1 (MQ=255)                                                           cgccCAGCAGGTGCTTAAAGTGGTGCCATGTGAGCGATGTCCTTCTTCTTAACAAACGCGAAATCCGTGAGTGCGTTACGCTCGATACAGATAGTATCGATGTAGTGGAAAACGCATTTTGATCGCTGTCGAGAGGTAAAGTGGTCATg                                                                                          <  8:60833/149‑1 (MQ=255)                                                            gccCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGTCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGc                                                                                         <  2:195376/149‑1 (MQ=255)                                                                     gTGCTTAAAGAGGTGCCATGTGAGCGAAGTCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAAtcc                                                                                >  8:43372/1‑149 (MQ=255)                                                                        cTTAAAGAGGTGCCATGTGAGCGAAGTCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAAtcctc                                                                              <  6:141307/148‑1 (MQ=255)                                                                        cTTAAAGAGGTGCCATGTGAGCGAAGTCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAAtcctc                                                                              <  8:82632/148‑1 (MQ=255)                                                                           aaaGAGGTGCCATGTGAGCGAAGTCCTTGTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCg                                                                            >  1:228943/1‑147 (MQ=255)                                                                           aaaGAGGTGCCATGTGAGCGAAGTCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCg                                                                            <  2:228943/147‑1 (MQ=255)                                                                           aaaGAGGTGCCATGTGAGCGAAGTCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCg                                                                            <  6:159120/147‑1 (MQ=255)                                                                           aaaGAGGTGCCATGTGAGCGAAGTCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGt                                                                          >  5:33434/1‑149 (MQ=255)                                                                             agagGTGCCATGTGAGCGAAGTCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTa                                                                        <  1:168174/149‑1 (MQ=255)                                                                               agGTGCCATGTGAGCGAAGTCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGATAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCAGTCAATCCTCCGGTTAGa                                                                      >  2:306347/1‑149 (MQ=255)                                                                                      aTGTGAGCGAAGTGCTTATGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCGCGATACAGACAGTATAGATGTAGTGGAAAACGCATTTGCATCGCTGGGCAGAGGCAAAGGGGTCATGCCGTCCATCCTCCGGGTAGACATGTcc                                                               >  6:327697/1‑149 (MQ=255)                                                                                       tgtgAGCGAAGTCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTcc                                                               >  4:209339/1‑148 (MQ=255)                                                                                        gtgAGCGAAGTCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCTATCCTCCGGTTACACATCTcc                                                               >  8:217836/1‑147 (MQ=255)                                                                                             cGAAGTCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTCCGAGcac                                                         >  5:261801/1‑148 (MQ=255)                                                                                              gAAGTCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTCCGAGCACaa                                                       >  4:47160/1‑149 (MQ=255)                                                                                                    cTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTCCGAGCACAATGGc                                                   >  1:294623/1‑147 (MQ=255)                                                                                                       cTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTCCGAGCACAATGGCGAg                                                <  4:62997/147‑1 (MQ=255)                                                                                                        tGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTCCGAGCACAATGGCGAGGTg                                             >  3:58113/1‑149 (MQ=255)                                                                                                                cAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTCCGAGCACAATGGCGAGGTGGACGTCa                                      >  1:311625/1‑148 (MQ=255)                                                                                                                              tGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTa                                                                        <  6:29148/100‑1 (MQ=255)                                                                                                                              tGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTa                                                                        >  5:29148/1‑100 (MQ=255)                                                                                                                              tGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTCCGAGCACAATGGCGAGGTGGACGTCAAAACCGCTTACCTGc                       <  1:21672/149‑1 (MQ=255)                                                                                                                              tGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTCCGAGCACAATGGCGAGGTGGACGTCAAAACCGCTTACCTGc                       <  2:294623/149‑1 (MQ=255)                                                                                                                              tGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTCCGAGCACAATGGCGAGGTGGACGTCAAAACCGCTTACCTGc                       <  3:35821/149‑1 (MQ=255)                                                                                                                                aGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTCCGAGCACAATGGCGAGGTGGACGTCAAAACCGCTTACCTGCCa                     <  1:151186/149‑1 (MQ=255)                                                                                                                                aGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTCCGAGCACAATGGCGAGGTGGACGTCAAAACCGCTTACCTGCCa                     <  7:43372/149‑1 (MQ=255)                                                                                                                                      actacatCGATACTGTCTGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTCCGAGCACAATGGCGAGGTGGACGTCAAAACCGCTTACCTGCCAGACCtt               >  2:279814/7‑149 (MQ=255)                                                                                                                                        aCGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTCCGAGCACAATGGCGAGGTGGACGTCAAAACCGCTTACCTGCCAGACCTTg              <  1:268026/148‑1 (MQ=255)                                                                                                                                         cGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTCCGAGCACAATGGCGAGGTGGACGTCAAAACCGCTTACCTGCCAGACCtt               >  2:353686/1‑146 (MQ=255)                                                                                                                                         cGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTCCGAGCACAATGGCGAGGTGGACGTCAAAACCGCTTACCTGCCAGACCTTga             <  1:353686/148‑1 (MQ=255)                                                                                                                                           cTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTCCGAGCACAATGGCGAGGTGGACGTCAAAACCGCTTACCTGCCAGACCTTGAGaa          <  3:224089/149‑1 (MQ=255)                                                                                                                                               aTACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTCCGAGCACAATGGCGAGGTGGACGTCAAAACCGCTTACCTGCCAGACCTTGAGAATTTc      <  3:209339/149‑1 (MQ=255)                                                                                                                                               aTACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTCCGAGCACAATGGCGAGGTGGACGTCAAAACCGCTTACCTGCCAGACCTTGAGAATTTc      >  3:336119/1‑149 (MQ=255)                                                                                                                                               aTACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTCCGAGCACAATGGCGAGGTGGACGTCAAAACCGCTTACCTGCCAGACCTTGAGAATTTc      >  3:259641/1‑149 (MQ=255)                                                                                                                                                   agacagTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTCCGAGCACAATGGCGAGGTGGACGTCAAAACCGCTTACCTGCCAGACCTTGAGAATTTCGCTa  >  6:411629/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                   TGGGAATGTCGACATGGCTCAACACCAACGCATCATCCTGGAAGCTGTGCGCGACCGCCGCCCAGCACGTGCTTAAAGAGGTGCCATGTGAGCGAAGTCCTTCTGCTTAACCAACGCGAAATCCGTGAGTGCGTAACGCTCGATACAGACAGTATCGATGTAGTGGAAAACGCATTTGCATCGCTGGCCAGAGGCAAAGTGGTCATGCCGTCAATCCTCCGGTTAGACATCTCCGAGCACAATGGCGAGGTGGACGTCAAAACCGCTTACCTGCCAGACCTTGAGAATTTCGCTA  >  NC_002947/5026604‑5026898

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