Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	780,929	C→T	intergenic (+11/‑139)	lysW → / → valZ	tRNA‑Lys/tRNA‑Val

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	780,929	0	C	T	98.3%	203.0 / ‑6.0	60	intergenic (+11/‑139)	lysW/valZ	tRNA‑Lys/tRNA‑Val
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (28/31);  minor base G (1/0);  total (29/31)
Fisher's exact test for biased strand distribution p-value = 4.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

GGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCGGGTGA  >  NC_000913/780782‑781073                                                                                                                                                    |                                                                                                                                                 gggAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATTGaa                                                                                                                                               <  2:720514/151‑1 (MQ=18) gggAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATTGa                                                                                                                                                <  2:671479/150‑1 (MQ=18)    agagCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATTGAAGAc                                                                                                                                            >  2:731738/1‑151 (MQ=25)      aactccacccTTACAAGGAGGGGGTCGGCGATTCGAACCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATTGAAGACAg                                                                                                                                          <  1:653902/144‑1 (MQ=255)         aCCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCTCCCCCCACCGGGTCGTCAGCTCAGTTTTCCGAGCTGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATTGAAGACAGATg                                                                                                                                       <  1:231406/151‑1 (MQ=14)              ccTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATTGAAGACAGATGTTAAg                                                                                                                                  <  2:137387/151‑1 (MQ=255)                  ccgaggcggGGGTCGGCGGTTCGAACCCGTCATCACCCAACACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATTGAAGACAGAt                                                                                                                                        <  1:683713/138‑1 (MQ=255)                            gggaggcggTACGATCCCGTCATCACCCAACACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATTGAAGACAGATGTTAAGCGTGTAGGATAAc                                                                                                                     <  1:599430/146‑1 (MQ=255)                                ggcggTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATTGAAGACAGATGTTAAGCGTGTAGGATAACGTTGc                                                                                                                >  2:786550/1‑151 (MQ=255)                                ggcggTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGGTCGAATCCTGCACGACCCACCAGTTTTAAAATTGAAAACAGATGTTAAGCGTGTAGGATAACGtt                                                                                                                  >  2:799792/1‑149 (MQ=255)                                       cGATCCCGTCATCACCCAAAACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATTGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCa                                                                                                            <  1:718883/148‑1 (MQ=255)                                              gTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATTGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACg                                                                                                      <  2:400981/147‑1 (MQ=255)                                                catcaCCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATTGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCg                                                                                                 >  1:317551/1‑150 (MQ=255)                                                 atcaCCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATTTTAACATTGAACATAGATGTGAAGCTTTACGAATAACGTTGCGCAAGCAACGGCACGt                                                                                                >  1:109606/1‑150 (MQ=11)                                                 atcaCCCAACACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATTGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACg                                                                                                      <  1:831603/144‑1 (MQ=255)                                                 acaaccaaccaGCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATTGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCg                                                                                                 <  1:94625/146‑1 (MQ=255)                                                   cacccaccACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATTGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAgg                                                                                             <  2:847716/151‑1 (MQ=255)                                                    acccaccACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATTGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAggg                                                                                            <  1:235287/151‑1 (MQ=255)                                                    acccaccACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATTGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAggg                                                                                            >  2:199767/1‑151 (MQ=255)                                                     cccaccACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATGGAACACAGAGGTTAATTGTTTCGGTTAACGTTGCGTCAGCAACGGCCAGGAggg                                                                                            >  1:312304/1‑150 (MQ=14)                                                      ccaccaCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATTGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGgcg                                                                                          >  2:233453/1‑151 (MQ=255)                                                        accaccGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATTGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGgcgag                                                                                        <  2:57111/151‑1 (MQ=255)                                                          caccGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATTGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGgcgagcg                                                                                      >  1:144401/1‑151 (MQ=255)                                                            ccGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATTGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGagcgaa                                                                                    <  2:787397/151‑1 (MQ=255)                                                              gggTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATTGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGagcgaagc                                                                                  >  1:524526/1‑151 (MQ=255)                                                               ggTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATTGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGagcgaagcg                                                                                 <  2:363574/151‑1 (MQ=255)                                                                gTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATTGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGagcgaagcg                                                                                 >  1:389830/1‑150 (MQ=255)                                                                gTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATTGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGagcgaagcg                                                                                 >  1:108287/1‑150 (MQ=255)                                                                    ttAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATTGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTc                                                                             >  2:459340/1‑150 (MQ=255)                                                                      aGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATTGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATc                                                                          >  2:811104/1‑151 (MQ=255)                                                                          cAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATTGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCt                                                                        <  2:108310/149‑1 (MQ=255)                                                                            gTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATTGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCAc                                                                    <  1:199744/151‑1 (MQ=255)                                                                            gTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATTTTAACATTGAAGACGGATGATAAGCGGGTAGGATAACTTTGCGTCAGCAACGGCCGGTAGGGCAGGCGAAGCGAGTTATcc                                                                         >  1:3211/1‑146 (MQ=14)                                                                                     gCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATTGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACccacca                                                           <  1:311044/151‑1 (MQ=255)                                                                                       aGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATTGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACt                                                         >  2:348018/1‑151 (MQ=255)                                                                                         ttGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATTGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTa                                                        >  2:450467/1‑150 (MQ=255)                                                                                            aCTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATTGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGa                                                    <  2:245821/151‑1 (MQ=255)                                                                                               tttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATTGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACgg                                                 >  2:528544/1‑151 (MQ=255)                                                                                                  aaTCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATTGAAGACAGATGTTAAGCGTGTCGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTgg                                              >  2:129089/1‑151 (MQ=255)                                                                                                         tGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATTGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCgg                                        >  2:357854/1‑150 (MQ=255)                                                                                                          ggTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATTGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGt                                       <  2:771642/150‑1 (MQ=255)                                                                                                           gTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATTGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTCTTGAAGGGGGGCTCGGTg                                      >  1:593073/1‑150 (MQ=14)                                                                                                             cGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATTGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTATTAACGGGGGGTTCGGGGGa                                    >  1:41700/1‑150 (MQ=255)                                                                                                                aGGTTCGAATCCTGCACGACCCACCAGTTTTAACATTGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAagta                                <  2:517911/151‑1 (MQ=255)                                                                                                                  gTTCGAATCCTGCACGACCCACCAGTTTTAACATTGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAagtagt                              <  2:144424/151‑1 (MQ=255)                                                                                                                   ttCGAATCCTGCACGACCCACCAGTTTTAACATTGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGtt                             >  2:775405/1‑151 (MQ=255)                                                                                                                     cGAATCCTGCACGACCCACCAGTTTTAACATTGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGttt                            <  2:437248/150‑1 (MQ=255)                                                                                                                      gAATCCTGCACGACCCACCAGTTTTAACATTGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTgt                          >  2:666539/1‑151 (MQ=255)                                                                                                                       aaTCCTGCACGACCCACCAGTTTTAACATTGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTgta                         <  2:273830/151‑1 (MQ=255)                                                                                                                        aTCCTGCACGACCCACCAGTTTTAACATTGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGGGCGGAGGAAGTAGTTTgt                          >  1:235106/1‑149 (MQ=255)                                                                                                                           cTGCACGACCCACCAGTTTTAACATTGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGAccc                                                               >  2:403717/1‑109 (MQ=255)                                                                                                                           cTGCACGACCCACCAGTTTTAACATTGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGAccc                                                               <  1:403685/109‑1 (MQ=255)                                                                                                                           cTGCACGACCCACCAGTTTTAACATTGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTCATGACGGTTGGTTCGGAGGAATAAGCTTgtag                        >  1:35268/1‑148 (MQ=17)                                                                                                                               aCGACCCACCAGTTTTAACATTGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAg                 <  1:512846/151‑1 (MQ=255)                                                                                                                                 gACCCACCAGTTTTAACATTGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAgcg               >  2:469327/1‑151 (MQ=255)                                                                                                                                 gACCCACCAGTTTTAACATTGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAgc                <  1:249918/150‑1 (MQ=255)                                                                                                                                  aCCCACCAGTTTTAACATTGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAgcgc              <  2:584640/151‑1 (MQ=255)                                                                                                                                       ccaGTTTTAACATTGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTAt         <  1:551208/151‑1 (MQ=255)                                                                                                                                              tAACATTGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGAGGAAAAAGTTTGTAGTAACCAGCCCAGTAGCGGGt    >  1:335801/1‑149 (MQ=21)                                                                                                                                               aaCATTGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCGGGTGa  <  1:375357/150‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                 GGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCACTAATGACGGTGGGTTCGGTGGAAGTAGTTTGTAGTATCCAGCGCAGTATCGGGTGA  >  NC_000913/780782‑781073

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