Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,747,155	G→T	14.8%	A9S (GCC→TCC)	ydhS →	putative oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,747,155	0	G	T	14.8%	124.5 / 2.9	62	A9S (GCC→TCC)	ydhS	putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (30/21);  new base T (2/9);  total (32/30)
Fisher's exact test for biased strand distribution p-value = 2.00e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGCCGGGCCTACGGGGATCTACACCTTATTCTCGCTTCTACAGCAACAAACTCCACTTTCTATTTCTATCTTCGAGCAGGCTGACGAGGCCGGTGTCGGGATGCCATACAGTGATGAGGAAAACTCAAAAATGATGCTGGCAAATATTGC  >  NC_000913/1747007‑1747303                                                                                                                                                     |                                                                                                                                                     aGAATTTAATCAAGGGCGGTTAGCGCCCTATTCATCACTGTCTGAAATTTCTCAAATACTAAAAATCTCAAACAAACTTATCTGATAACACTAAATTCGAAAGAATGCTTACAGGTAAGTAACAATGAAAAAAATTTCTATAGTGGGTg                                                                                                                                                      >  2:767936/1‑149 (MQ=255)          taaatgtcGGATATCGACCTTTTCATCCCTGTCGGAAATTTCTCAAATTCTAAAAATCTCAACCAGAGTTATATGATAAAAATAAATACCAAATAATGCGAAAATGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTTCAGGGCATa                                                                                                                                             >  2:548813/8‑149 (MQ=255)          tCAAGGGCGGTTAGCGCCCTTTTCATCACTGTCTGAAATTTCTCAAATTCTAAAAATATCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGCCGGGCCTa                                                                                                                                             >  1:1405698/1‑149 (MQ=255)            aaGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGCCGGGCCTACg                                                                                                                                           >  1:763407/1‑149 (MQ=255)             aGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGCCGGGCCTACgg                                                                                                                                          >  2:768016/1‑149 (MQ=255)              gggCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAGAATTGCTATTGTGGGTGCCGGGCCTACggg                                                                                                                                         >  1:791878/1‑149 (MQ=255)                 cGGTTCGCGCCCTTTTCATCACTCTCTGAAATTTCACAAATTCTAAAAATCTCAACCAAACTTATCAGATAACACTAAAATCGTAAGAATGCGTACAGGAAAGTAACAAAGAAAAAAAATCCTATTGGGGGTGCCGGGCCTACGGGGAt                                                                                                                                      <  2:839365/149‑1 (MQ=255)                  ggTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGCCGGGCCTACGGGGATc                                                                                                                                     <  1:506750/149‑1 (MQ=255)                  ggTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGCCGGGCCTACGGGGATc                                                                                                                                     <  1:1398904/149‑1 (MQ=255)                        cgcCCTTTTCATCCCTTTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGCCGGGCCTACGGGGATCTACAcc                                                                                                                               <  1:1194711/149‑1 (MQ=255)                         gcCCTTTTCATCCCTGTCTGAAATTTATCAAATTCTAAAAATCTCACCCAAACTTATCTGATAACACTAAATACGAAAGAAGGCGTACAGGGAAGTAACAATGAAAAAAATTGCTATTGTGGGTGCCGGGCCTACGGGGATCTACACCt                                                                                                                              <  1:1040952/149‑1 (MQ=255)                          cccTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGCCGGGCCTACGGGGATCTACACCtt                                                                                                                             >  1:655020/1‑149 (MQ=255)                           ccTTTTCATCCCAGGCGGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTTTTTTTTTTGTTTCGGGGACTACGGTGATCTACACCTTa                                                                                                                            <  2:1320840/149‑1 (MQ=255)                                    cccTGTCTGAAATTTCTCAAAATCTAAAAATCTCAACCAAACTTATCTGATAACACAAAATACGAAAGAATGCGTACACGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGCCGGGCCTACGGGGATCTACACCTTATTCTCGAtt                                                                                                                   >  2:1330305/1‑149 (MQ=255)                                    cccTGTATGAAATTTCTCAAATACTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATTCGTACATGTAAGTAACAATGAAAAATTTTTTTTTTTTTGTTTCGGGGACTACGGTGATCTACACCTTATTCTCGCtt                                                                                                                   <  2:1559609/149‑1 (MQ=255)                                      cTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGATACAATGAAAAAAATTCCTATTGTGGGTGCCGTGCCTACGCTGATCTACACCTTATTCTCGCTTCt                                                                                                                 >  2:1458745/1‑149 (MQ=255)                                       tGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTTTTTTTTTTGTTTCCGGGACTACGGTGATCTACACCTTATTCTCGCTTCTa                                                                                                                <  2:1539536/149‑1 (MQ=255)                                         tCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGCCGGGCCTACGGGGATCTACACCTTATTCTCGCTTCTACa                                                                                                              >  2:1349434/1‑149 (MQ=255)                                          cTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAATTTTTTTTTTTTGTTTCGGGGACTACGGTGATCTACACCTTATTCTCGCTTCTACAg                                                                                                             <  2:825139/149‑1 (MQ=255)                                           ttaaaTTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAAATGCTATTGTGGGTGCCGGGCCTACGGGGATCTACAAATTATTCTCGCTTCTACAGc                                                                                                            >  2:933005/3‑149 (MQ=255)                                           tGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGCCGGGCCTACGGGGATCTACACCTTATTCTCGCTTCTACAGc                                                                                                            >  1:704157/1‑149 (MQ=255)                                                tttCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGCCGGGCCTACGGGGATCTACACCTTATTCTCGCTTCTACAGcaacaa                                                                                                       >  2:404555/1‑149 (MQ=255)                                                    tcAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAATGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTTTTGTGGGTGCCGGGCCTACGGGGATCTACACCTTATTCTCGCTTCTACAGCAACAAACTc                                                                                                   <  1:1530086/149‑1 (MQ=255)                                                      aaaTTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGCCGGGCCTACggg                                                                                                                                         >  1:1190597/1‑109 (MQ=255)                                                      aaaTTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGCCGGGCCTACggg                                                                                                                                         <  2:1190597/109‑1 (MQ=255)                                                       aaTTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGCCGGGCCTACGGGGATCTACACCTTATTCTCGCTTCTACAGCAACAAACTCCAc                                                                                                <  1:1426896/149‑1 (MQ=255)                                                        aTTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGCCGGGCCTACGGGGATCTACACCTTATTCTCGCTTCTACAGCAACAAACTCCACt                                                                                               <  1:1353139/149‑1 (MQ=255)                                                             aaaaaTCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGCCGGGCCTACGGGGATCTACACCTTATTCTCGCTTCTACAGCAACAAACTCCACtttcta                                                                                          >  1:1304264/1‑149 (MQ=255)                                                             aaaaaTCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATCAAAAAAATTGCTATTGTGGGTTCCGGGCCTACGGGGATCTACACGTTATTCACGCTTCTAGAGCAACAAACTCCACtttcta                                                                                          >  2:530465/1‑149 (MQ=255)                                                             aaaaaTATCAACCAAACTTATGAGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACATTGAAAAAAATTGCTATTGTGGGTGCCGTGCCTACGGGGATCTACACCTTATTCTCGCTTCTACAGCAACAAACTCCACtttcta                                                                                          <  1:1266431/149‑1 (MQ=255)                                                              aaaaTCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGCCGGGCCTACGGGGATCTACACCTTATTCTCGCTTCTACAGCAACAAACTCCACtttctat                                                                                         >  2:39179/1‑149 (MQ=255)                                                                aaTCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAATTTTTTTTTGTTGTTTCCGGGCCTACGGTGATCTACACCTTATTCTCGCTTCTACAGCAACAAACTCCACtttctattt                                                                                       <  2:481138/149‑1 (MQ=255)                                                                aaTCTCAACCAAAATTATCTGATAACACTAAATTCTAAAGAATGCGTACATGTAAGTAACAATGAAAAAAATTGCTATTGTGTGTGCCGGGTCTACGGGGATCTACACCTTATTCTCGCTTCTACAGCAACAAACTCCACtttctattt                                                                                       <  1:114439/149‑1 (MQ=255)                                                                       aCCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGCCGGGCCTACGGGGATCTACACCTTATTCTCGCTTCTACAGCAACAAACTCCACTTTCTATTTCTATCtt                                                                                >  1:1481433/1‑149 (MQ=255)                                                                        ccAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGCCGGGCCTACGGGGATCTACACCTTATTCTCGCTTCTACAGCAACAAACTCCACTTTCTATTTCTATCTTc                                                                               <  1:778752/149‑1 (MQ=255)                                                                         cAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGCCGGGCCTACGGGGATCTACACCTTATTCTCGCTTCTACAGCAACAAACTCCACTTTCTATTTCTATCTTCg                                                                              >  1:816970/1‑149 (MQ=255)                                                                                   tGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGCCGGGCCTACGGGGATCTACACCTTATTCTCGATTCTACAGCAACAAACTCCACTTTCTATTTCTATCTTCGATCAGGCTGa                                                                    >  1:1041444/1‑149 (MQ=255)                                                                                     aTAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGCCGGGCCTACGGGGATCTACACCTTATTCTCGCTTCTACACCAACAAACTCCACTTTCTATTTCTATCTTCGAGCAGGCTGACg                                                                  >  2:821771/1‑149 (MQ=255)                                                                                       aaCACTAAATTCGAAAGAATGCGTACAGGTCATTAACAATTAAAATAATTTTTTTTTTTGTTTCGGGGACTACGGTGATCTACACCTTATTCTCGCTTCTACAGCAACAAACTCCACTTTATATTTCTATCTTCGAGCAGGCTGACGAg                                                                <  2:852880/149‑1 (MQ=255)                                                                                            tAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGCCGGGCCTACGGGGATCTACACCTTATTCTCGCTTCTACAGCAACAAACTCCACTTTCTATTTCTATCTTCGAGCAGGCTGACGAGGCCgg                                                           >  2:758730/1‑149 (MQ=255)                                                                                               aTTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTTATATTGTTGTTGCCGGGACTACGGGGATCTACACCTTATTCTCGCTTCTACAGCAACAAACTCCACTTTCTATTTCTATCTTCGAGCAGGCTGACGAGGCCGgtgt                                                        <  2:517379/149‑1 (MQ=255)                                                                                                ttCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGCCGGGCCTACGGGGATCTACACCTTATTCTCGCTTCTACAGCAAGAAACTCCACTTTCTATTTCTATCTTCGAGCAGGCTGACGAGGCCGGTGTc                                                       >  2:1166568/1‑149 (MQ=255)                                                                                                  cGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGCCGGGCCTACGGGGATCTACACCTTATTCTCGCTTCTACAGCAACAAACTCCACTTTCTATTTCTATATTCGAGCAGGCTGACGAGGCCGGTGTCgg                                                     >  2:1216982/1‑149 (MQ=255)                                                                                                  cGAAACAATGCGTACAGGTAAGTAACAATGTAAAAATTTTTTTTTTTGGTTTCGGGGACTACGGTGATCTACACCTTATTCTCGCTTCTACAGCAACAAACTCCACTTTCTATTTCTATCTTCGAGCAGGCTGACGAGGCCGGTGTCgg                                                     <  2:1502170/149‑1 (MQ=255)                                                                                                    aaaGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGCCGGGCCTACGGGGATCTACACCTTATTCTCGCTTCTACAGCAACAAACTCCACTTTCTATTTCTATCTTCGAGCAGGCTGACGAGGCCTGTGTCGGGa                                                   >  2:901226/1‑149 (MQ=255)                                                                                                      aGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTGGTTGCCGGGCCTACGGTGATCTACACCTTATTCTCGCTTCTACAGCAACAAACTCCACTTTCTATTTCTATCTTCGAGCAGGCTGACGAGGCCGGTGTCGGGATg                                                 <  2:762259/149‑1 (MQ=255)                                                                                                        aaTGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGCCGGGCCTACGGGGATCTACACCTTATTCTCGCTTCTACAGCAACAAACTCCACTTTCTATTTCTATCTTCGAGCAGGCTGACGAGGCCGGTGTCGGGATGcc                                               >  1:1396324/1‑149 (MQ=255)                                                                                                         aTGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGCCGGGCCTACGGGGATCTACACCTTATTCTCGCTTCTACAGCAACAAACTCCACTTTCTATTTCTATCTTCGAGCAGGCTGACGAGGCCGGTGTCGGGATGCCa                                              <  1:1146637/149‑1 (MQ=255)                                                                                                               aCAGGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGCCGGGCCTACGGGGATCTACACCTTATTCTCGCTTCTACAGCATCAAACTCCACTTTCTATTTCTATCTTCGAGCAGGCTGACGAGGCCGGTGTCGGGATGCCATACAGt                                        >  1:367562/1‑149 (MQ=255)                                                                                                                  ggTAAGTAACAATGAAAAAAATTGCTTTTTTGGTTGCCGGGCCTACGGTGATCTACACCTTATTCTCGCTTCTACAGCAACAAACTCCACTTTCTATTTCTATCTTCGAGCAGGCTGACGAGGCCGGTGTCGGGATGCCATACAGtgat                                     <  2:763407/149‑1 (MQ=255)                                                                                                                  ggTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGCCGGGCCTACGGGGATCTACACCTTATTCTCGCTTCTACAGCAACAAACTCCACTTTCTATTTCTATCTTCGAGCAGGCTGACGAGGCCGGTGTCGGGATGCCATACAGtgat                                     <  1:1349434/149‑1 (MQ=255)                                                                                                                   gtaagtaaCAATGAAAAAAATTGCTATTGTGGGTGCCGGGCCTACGGGGATCTACACCTTATTCTCGCTTCTACAGCAACAAACTCCACTTTCTATTTCTATCTTCGAGCAGGCTGACGAGGCCGGTGTCGGGATGCCATACAGtgatg                                    >  2:1523165/1‑149 (MQ=255)                                                                                                                   gtaagtaaCAATGAAAAAAATTGCTATTGTGGGTGCCGGGCCTACGGGGATCTACACCTTATTCTCGCTTCTACAGCAACAAACTCCACTTTCTATTTCTATCTTCGAGCAGGCTGACGAGGCCGGTGTCGGGATGCCATACAGtgatg                                    >  2:393848/1‑149 (MQ=255)                                                                                                                    taagtaaCAATGAAAAAAATTGCTATTGAGGGTGCCGGGACTACGGGGATCTACACCTTATTCTCGCTTCTACAGCAACAAACTCCACTTTCTATTTCTATCTTCGAGCAGGCTGACGAGGCCGGTGTCGGGATGCCATACAGtgatga                                   <  2:1126565/149‑1 (MQ=255)                                                                                                                     aagaaaCAATTAAAAAAATTGCTATTGTGGGTGCGGGGCGTACGGGGATCTACACCTTATTCTCGCTTCTACAGCAACATACTCCACTTTCTATTTCTATCTTCGAGCAGGCTGACGAGGCCGGTGTCGGGATGCCATACAGTGATGAg                                  <  1:185352/149‑1 (MQ=255)                                                                                                                      agtaaCAATGAAAAAAATTGCTATTGTGGGTGCCGGGCCTACGGGGATCTACACCTTATTCTCGCTTCTACAGCAACAAACTCCACTTTCTATTTATATCTTCGAGCAGGCTGACGAGGCCGGTGTCGGGATGCCATACAGTGATGAgg                                 >  2:1342162/1‑149 (MQ=255)                                                                                                                         aaCAATGAAAAAAATTGCTATTGTGGGTGCCGGGCCTACGGGGATCTACACCTTATTCTCGCTTCTACAGCAGCAAACTCCACTTTCTATTTCTATCTTCGAGCAGGCTGACGAGGCCGGTGTCGGGATGCCATACAGTGATGAGGaaa                              >  1:453586/1‑149 (MQ=255)                                                                                                                                 aaaaaaTTGCTATTGTGGGTGCCGGGCCTACGGGGATCTACACCTTATTCTCGCTTCTACAGCAACAAACTCCACTTTCTATTTCTATCTTCGAGCAGGCTGACGAGGCCGGTGTCGGGATGCCATACAGTGATGAGGAAAACTCaaaa                      >  2:1557998/1‑149 (MQ=255)                                                                                                                                 aaaaaaTTGCTATTGTGGGGTCCGCGCCTACGGGGATCTACACCTTATGCTCGCTTCTCCCGCAACAAACTCCCCTTTCTATTGCTATCTTCGAGGAGGCTGACGAGGCCGGTGTCGGGATGCCATACAGTGATGCGGAAAACTCaaaa                      <  1:748967/149‑1 (MQ=255)                                                                                                                                      aTTGCTATTATTGCTGCCGGGCCTACGGGGATGTACAGCTTATTCTCGCTTCTACAGCAACAAACTCCACTTTCTATTTCTATCTTGGAGCTGGCTGACGAGGACGGTGGCGGGATGCCATACAGTGATGATGAAAACTCAAAAatgat                 <  2:921970/149‑1 (MQ=255)                                                                                                                                       ttgtttttgttgGTTCGGGGACTACGGTGATCTACACCTTATTCTCGCTTCTACAGCAACAAACTCCACTTTCTATTTCTATCTTCGATCAGGCTGACGAGGCCGGTGTCGGGATGCCATACAGTGATGAGGAAAACTCAAAAatgatg                <  2:1304264/143‑1 (MQ=255)                                                                                                                                           tATTGTGGGTGCCGGGCCTACGGGGATCTACACCTTATTCTCGCTTCTACAGCAACAAACTCCACTTTCTATTTCTATCTTCGAGCAGGCTGACGAGGCCGGTGTCGGGATGCCATACAGTGATGAGGAAAACTCAAAAATGATGCTgg            <  2:1481433/149‑1 (MQ=255)                                                                                                                                            aTTGTGGGTGCCGGGCCTACGGGGATCTACACCTTATTCTCGCTTCTACAGCAACAAACTCCTCTTTCTATTTCTATCTTCGAGCAGGCTGACGAGGCCGGTGTCGGGATGCCATACAGTGATGAGGAAAACTCAAAAATGATGCTGGc           >  1:598297/1‑149 (MQ=255)                                                                                                                                                     gCCGGGCCTACGGGGATCTACACCTTATTCTCGCTTCTACAGCAACAAACTCCACTTTCTATTTCTATCTTCGAGCAGGCTGACGAGGCCGGTGTCGGGATGCCATACAGTGATGAGGAAAACTCAAAAATGATGCTGGCAAATATTGc  <  1:768016/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                     AGAATTTAATCAAGGGCGGTTAGCGCCCTTTTCATCCCTGTCTGAAATTTCTCAAATTCTAAAAATCTCAACCAAACTTATCTGATAACACTAAATTCGAAAGAATGCGTACAGGTAAGTAACAATGAAAAAAATTGCTATTGTGGGTGCCGGGCCTACGGGGATCTACACCTTATTCTCGCTTCTACAGCAACAAACTCCACTTTCTATTTCTATCTTCGAGCAGGCTGACGAGGCCGGTGTCGGGATGCCATACAGTGATGAGGAAAACTCAAAAATGATGCTGGCAAATATTGC  >  NC_000913/1747007‑1747303

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑