Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,606,215	(G)8→7	100%	noncoding (72/87 nt)	leuP ←	tRNA‑Leu

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,606,208	0	G	.	55.1%	18.5 / 185.8	89	noncoding (79/87 nt)	leuP	tRNA‑Leu
Reads supporting (aligned to +/- strand):  ref base G (23/16);  new base . (27/22);  total (50/39)
Fisher's exact test for biased strand distribution p-value = 8.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.76e-01

GTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTT‑GATATCGCCTCGTTTGGTG‑CGA‑‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATAT  >  NC_000913/4606109‑4606302                                                                                                        |                                                                                               gTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTT‑GATATCGCCTCGTTTGGTG‑CGA‑‑g                                                                                                <  1:1944596/100‑1 (MQ=255)    cGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTAAACTTACGTTTTGT‑GATATAGTGTCGTTTGGTG‑CGA‑‑aggg                                                                                             >  1:239230/1‑100 (MQ=17)     gTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTT‑GATATCGCCTCGTTTGGTG‑CGA‑‑ggggg                                                                                            <  2:991010/100‑1 (MQ=255)     gTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTT‑GATATCGCCTCGTTTGGTG‑CGA‑‑ggggg                                                                                            <  1:2508010/100‑1 (MQ=255)      tAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTT‑GATATCGCCTCGTTTGGTG‑CGA‑‑gggggg                                                                                           >  2:1014339/1‑100 (MQ=255)      tAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTT‑GATATCGCCTCGTTTGGTG‑CGA‑‑gggggg                                                                                           <  2:2564766/100‑1 (MQ=255)         ggACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTT‑GATATCACATAGTTTGGTG‑CGA‑‑ggggggggc                                                                                        >  2:1509100/1‑99 (MQ=255)          gACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTT‑GATATCGCCTCGTTTGGTG‑CGA‑‑gggggggact                                                                                       >  2:1527645/1‑97 (MQ=255)           acacTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTT‑GATATCGCCTCGTTTGGTG‑CGA‑‑gggggggact                                                                                       <  2:1144990/99‑4 (MQ=255)              cTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTT‑GATATCGCCTCGTTTGGTG‑CGA‑‑gggcgtgaact                                                                                      >  1:2383308/1‑89 (MQ=255)               tAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTT‑GATATCGCCTCGTTTGGTG‑CGA‑‑GGGGGGGGACTTGaa                                                                                  >  2:2338602/1‑100 (MQ=255)               tAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTAATCTTACTTTTTTT‑GATATCGCCTCGTTTGGTG‑CGA‑‑GGGGGGGGGCttca                                                                                   >  2:2573336/1‑97 (MQ=255)                 acacCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTT‑GATATCGCCTCGTTTGGTG‑CGA‑‑GGGGGGGGACTTGAAcc                                                                                >  2:2031060/1‑100 (MQ=255)                 acacCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTT‑GATATCGCCTCGTTTGGTG‑CGA‑‑‑GGGGGGGACTTGAAcc                                                                                <  2:992067/99‑1 (MQ=255)                  cacCTGAAGCTAGCGCGTCTAACAATTCCGCCACCTTCGCACAGTCGTCTTACTTTTTTT‑GATATCGCCTCGTTTGGGG‑CGAGGGGGGGGGGACTTGAAc                                                                                 <  1:599250/100‑1 (MQ=25)                   acCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTT‑GATATCGCCTCGTTTGGTG‑CGA‑‑‑GGGGGGGACTTGAAccccc                                                                             <  2:1581131/100‑1 (MQ=255)                    ccggcagatAGCGCGTCTACCAATCCCGCCTCCTGCGCACAGTCCTCTTACTTTTTTTGGAAATCGCCTCGTTTGGTG‑CGAG‑GGGGGGGGACTTGAAcc                                                                                <  2:1447999/91‑1 (MQ=255)                    ccTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTT‑GATATCGCCTCGTTTGGTG‑CGA‑‑‑GGGGGGGACTTGAAccccc                                                                             <  2:238674/99‑1 (MQ=255)                    ccTGAAGCTAGCGCGGCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTT‑GATATCGCCTCGTTTGGTG‑CGA‑‑‑GGGGGGGACTTGAACCCCCa                                                                            <  2:852143/100‑1 (MQ=255)                     cTGACGCTAGCCCGTCTAACAATTCCGCAACCTTCGCCCAGTCACCTTACTTTTTTT‑GAACTAGCATAGTTTGGTG‑CGA‑‑‑GGGGGGGACTTGAACCCCCa                                                                            <  1:1814421/99‑1 (MQ=19)                      tGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATGTTACTTTTTTT‑GATATCGCCTCGTTTGGTG‑CGA‑‑‑GGGGGGGACTTGAACCCCCACg                                                                          <  1:260912/100‑1 (MQ=255)                      tGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTT‑GATATCGCCTCGTTTGGTG‑CGA‑‑GGGGGGGGACTTGAACCCCCAc                                                                           >  2:1055771/1‑100 (MQ=255)                      tGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTT‑GATATCGCCTCGTTTGGTG‑CGA‑‑‑GGGGGGGACTTGAACCCCCACg                                                                          <  1:218388/100‑1 (MQ=255)                      tGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTT‑GATATCGCCTCGTTTGGTG‑CGA‑‑‑GGGGGGGACTTGAACCCCCACg                                                                          <  2:214754/100‑1 (MQ=255)                       ggaGCTAGCGCGGCTAGCAATGCTGCGACCTGGGCGCAGGGATGTTACGGTTGTG‑GAGATCGCCTCGTTGGGTG‑CGA‑‑GGGGGGGGACTTGAACCCCCACg                                                                          <  1:418040/98‑1 (MQ=12)                       gAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTT‑GATATCGCCTCGTTTGGTG‑CGA‑‑‑GGGGGGGACTTGAACCCCCACGt                                                                         <  2:1410001/100‑1 (MQ=255)                            tAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTAACTTTTTTT‑GATATCGCATCGTTTGGTG‑CGA‑‑GGGGGGGGACTTGAACCCCCACGTCCg                                                                      <  2:2573872/99‑1 (MQ=255)                             aGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTT‑GATATCGCCTCGTTTGGTG‑CGA‑‑‑GGGGGGGACTTGAACCCCCACGTCCGTaa                                                                   >  1:2469601/1‑100 (MQ=255)                               cgcgTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTT‑GATATCGCCTCGTTTGGTG‑CGA‑‑GGGGGGGGACTTGAACCCCCACGTCCGTAAg                                                                  <  2:769474/100‑1 (MQ=255)                                 cgTCTACCAATTCCGGCACCTTCGCACAGTCATCTTACTTTTATT‑GATATCGCCTCGTTTGGTG‑CGA‑‑‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAac                                                               <  1:2585277/100‑1 (MQ=38)                                 cgTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTT‑GATATCGCCTCGTTTGGTG‑CGA‑‑GGGGGGGGACTTGAACCCCCACGTCCGTa                                                                    <  1:2456939/96‑1 (MQ=255)                                 cgTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTT‑GATATCGCCTCGTTTGGTG‑CGA‑‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGaa                                                                >  1:555047/1‑100 (MQ=255)                                 cgTCTACCAATTCCGACACCTTCGCACAGTCATCTTATTTTTTTT‑AATATCGCCTCGGTTGGTG‑CGA‑‑‑GGGGGGGACTTGAAACAACAAGACAGTAAGaa                                                                <  1:1411279/99‑1 (MQ=255)                                  gTCTACCAATTCCGCCACCTTCGCACAGGGATCTTACTTTTTTT‑GATATCGCCTCGTTTGGTG‑CGA‑‑‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAaca                                                              <  1:172508/100‑1 (MQ=38)                                    cTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTT‑GATATCTCCTCGTTTGGTG‑CCA‑‑‑GGGGGGGACTTTTAACCCCACGTCCGTAAGAACACt                                                            >  1:462579/1‑100 (MQ=38)                                    cTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTT‑GATATCGCCTCGTTTGGTG‑CGA‑‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAacac                                                             >  1:1268/1‑100 (MQ=255)                                    cTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTT‑GATATCGCCTCGTTTGGTG‑CGA‑‑‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACt                                                            >  1:1403709/1‑100 (MQ=255)                                      actaaTTCCGGCACCATCGAACAGTCACAATACTTGTATT‑GAAATACCCTCATTTGTTG‑CGA‑‑‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTa                                                           <  1:1598373/96‑1 (MQ=12)                                      aCCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTT‑GATATCGCCTCGTTTGGTG‑CGA‑‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTa                                                           <  2:753737/100‑1 (MQ=255)                                      aCCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTT‑GATATCGCCTCGTTTGCTG‑CGA‑‑‑GGGGGGGACTTGAACCCCCACGTGCGTa                                                                    >  1:1278715/1‑90 (MQ=38)                                       ccAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTT‑GATATCGTCTCGTTTGGGG‑CGA‑‑‑GGGGGTGACTTGAACACCCACGTACGTAAGTacaatca                                                          >  1:659070/1‑95 (MQ=255)                                       ccAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTT‑GATATCGCCTCGTTTGGTG‑CGA‑‑‑GGGGGGGACTTGAACCCCCACGTCCGTCAGGACACTAac                                                         >  1:2501095/1‑100 (MQ=37)                                       ccAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTT‑GATATCGCCTCGTTTGGTG‑CGA‑‑‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAac                                                         >  1:2022981/1‑100 (MQ=255)                                       ccAATTCCGCCACCTTCGCACAGTCATCTTACTATTTTT‑GATATCGCCTCGTTTGGTG‑CGA‑‑GGGGGGGGAGTTGAACCCCCACGTCCGTAAGAACACTaa                                                          >  1:1557491/1‑100 (MQ=255)                                         aaTTCCGCCACCTTCGCACAGTCATCTTACTTTTTTT‑GATATCGCCTCGTTTGGTG‑CGA‑‑‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAacac                                                       >  2:19043/1‑100 (MQ=255)                                           ttCCGCCACCTTCGCACAGTCATCTTACTTTTTTT‑GATATCGCCTCGTTTGGTG‑CGA‑‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACAcc                                                      <  1:2031048/100‑1 (MQ=255)                                             ccGCCACCTTCGCACAGTCATCTTACTTTTTTT‑GATATCGCCTCGTTTGGTG‑CGA‑‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCt                                                     >  2:910911/1‑99 (MQ=255)                                             ccGCCACCTTCGCACAGTCATCTTAC‑TTTTTT‑GATATCGCCTCGTTTGGTG‑CGA‑‑‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGaa                                                  >  1:125503/1‑100 (MQ=255)                                             cagcgacgttCGCACAGGCATCATACATTTTAT‑GATATCGCCACGTTTGGTG‑CGA‑‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACAcc                                                      <  2:170210/90‑1 (MQ=18)                                              cGCCACCTTCGCACAGTCATCTTACTTTTTTT‑GATATCGCCTCGTTTGGTG‑CGA‑‑‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGaa                                                  >  1:256370/1‑100 (MQ=255)                                                 cACCTTCGCACAGTCATCTTACTTTTTTT‑GATATCGCCTCGTTTGGTG‑CGA‑‑‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCt                                               >  1:2443081/1‑100 (MQ=255)                                                   ccTTCGCACAGTCATCTTACTTTTTTT‑GATATCGCCTCGTTTGGTG‑CGA‑‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTa                                              >  1:25525/1‑100 (MQ=255)                                                    cTTCGCACAGTCATCTTACTTTTTTT‑TTTTTCGCTTCGTTTGGTG‑CGA‑‑GGGGGGGGACTTGCACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAg                                             >  1:448877/1‑100 (MQ=18)                                                    cTTCGCACAGTCATCTTACTTTTTTT‑GATATCGCCTCGTTTGGTG‑CGA‑‑‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAgc                                            >  2:2296983/1‑100 (MQ=255)                                                        gCACAGTCATCTTACTTTTTTT‑GATATCGCCTCGTTTGGTG‑CGAG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAgcgc                                          >  2:1460431/1‑100 (MQ=255)                                                        gCACAGTCATCTTACTTTTTTT‑GATATCGCCTCGTTTGGTG‑CGA‑‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGcgcg                                         >  2:12931/1‑100 (MQ=255)                                                         cacaGTCATCTTACTTTTTTT‑GATATCGCCTCGTTTGGTG‑CGA‑‑‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAAAACTAACACCTGAAGCTAGCGCGTc                                       >  2:1258162/1‑100 (MQ=255)                                                          acaGTCATCTTACTTTTTTT‑GATATCGCCTCGTTTGGTG‑CGA‑‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGcgcg                                         >  1:1781614/1‑98 (MQ=255)                                                           caGTCATCTTACTTTTTTT‑GATATCGCCTCGTTTGGTG‑CGA‑‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCt                                      <  2:2227439/100‑1 (MQ=255)                                                           caGTCATCTTACTTTTTTT‑GATATCGCCTCGTTTGGTG‑CGA‑‑‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTa                                     >  1:1152355/1‑100 (MQ=255)                                                             gTCATCTTACTTTTTTT‑TATATCTCCTCGTTTGGTG‑CGA‑‑‑GGGGGGGACTTGAACCCCCACGTTCGTAATAACACTAACACCTGAAGCTAGCGCGTCTAc                                    >  1:2260271/1‑99 (MQ=18)                                                             gTCATCTTACTTTTTTT‑GATATCGCCTCGTTTGGTG‑CGA‑‑‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTAcc                                   <  2:1286905/100‑1 (MQ=255)                                                                aTCTTACTTTTTTT‑GATATCGCCTCGTTTGGTG‑CGAG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCa                                  <  2:152369/100‑1 (MQ=255)                                                                aTCTTACTTTTTTT‑GATATCGCCTCGTTTGGTG‑CGAG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCa                                  >  2:2470462/1‑100 (MQ=255)                                                                aTCTTACTTTTTTT‑GATATCGCCTCGTTTGGTG‑CGA‑‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCaa                                 >  2:2276874/1‑100 (MQ=255)                                                                 tCTTACTTTTTTT‑GATATCGCCTCGTTTGGTG‑CGA‑‑‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAAtt                               <  2:548150/100‑1 (MQ=255)                                                                 tCTTACTTTTTTT‑GATATCGCCTCGTTTGGTG‑CGA‑‑‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAAtt                               >  1:321038/1‑100 (MQ=255)                                                                  cTTACTTTTTTT‑GATATCGCCTCGTTTGGTG‑CGA‑‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAAtt                               >  1:1354041/1‑100 (MQ=255)                                                                  cTTACTTTTTTT‑GATATCGCCTCGTTTGGTG‑CGA‑‑‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTc                              >  2:237803/1‑100 (MQ=255)                                                                  aggacggttttg‑GATATCGTCTAGTTTGGCG‑‑GT‑‑CGGGGGGTACTTGAACCCCCACGTCCGTAAGACCACGAACACCTGAAGCTAGCGCGTCTACCa                                  <  1:2369026/89‑1 (MQ=2)                                                                    tACTTTTTTT‑GATATCGCCTCGTTTGGTG‑CGA‑‑‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCg                            >  1:424561/1‑100 (MQ=255)                                                                    tACTTTTTTT‑GATATCGCCTCGTTTGCTG‑CGA‑‑‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCg                            <  1:337753/100‑1 (MQ=255)                                                                      cTTTTTTT‑GATATCGCCTCGTTTGGTG‑CGAG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCCATTCCg                            >  2:235415/1‑100 (MQ=255)                                                                      cTTTTTTT‑GATATCGCCTCGTTTGGTG‑CGA‑‑‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGcc                          >  2:290399/1‑100 (MQ=255)                                                                      cTTTTTTT‑GATATCGCCTCGTTTGGTG‑CGA‑‑‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGcc                          >  2:1918982/1‑100 (MQ=255)                                                                       ttttttt‑GATATCGCCTCGTTTGGTG‑CGA‑‑‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCa                         >  1:2563972/1‑100 (MQ=255)                                                                       ttttttt‑GATATCGCCTCGTTTGGTG‑CGA‑‑‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCa                         >  2:198856/1‑100 (MQ=255)                                                                       ttctttt‑GATATCGCCTCGTTTGGTG‑CGA‑‑‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCa                         <  2:1602286/97‑1 (MQ=255)                                                                          tttt‑GATATCGCCTCGTTTGGTG‑CGA‑‑‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCt                      >  1:1248978/1‑100 (MQ=255)                                                                            tt‑GATAACCCCGCGTTTGGTGCCGG‑‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCa                         <  1:2470444/97‑1 (MQ=14)                                                                             t‑GATATCGCCTCGTTTGGTG‑CGA‑‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACCCTAACACCTGAAGCTAGCGCGTCCACCAATTCCGCCACCTTc                    >  1:999382/1‑100 (MQ=255)                                                                                 tatCGGCTCGTTTGGTG‑CG‑‑‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGccc                <  2:2528869/100‑1 (MQ=255)                                                                                 tatCGCCTCGTTTGGTG‑CGA‑‑‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGccc                <  2:55060/100‑1 (MQ=255)                                                                                 taacacatCGTTTGGTG‑CGA‑‑‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGccc                <  1:771447/93‑1 (MQ=14)                                                                                  atCGCCTCGTTTGGTG‑CGA‑‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGccc                >  1:184261/1‑100 (MQ=255)                                                                                  atCGCCTCGTTTGGTG‑CGA‑‑‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCa               >  2:225925/1‑100 (MQ=255)                                                                                     gCCTCGTTTGGTG‑CGAGGGGGGGGGGACTTGAACCCCCACGTCCGTAAGAAGACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCa               >  1:280390/1‑100 (MQ=14)                                                                                     gCCTCGTTTGGTG‑CGA‑‑‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTGGCACCTGAAGCTAGCGCGTCTACCAATTGCGCCACCTTCGCCCAGt             >  1:2366689/1‑99 (MQ=14)                                                                                      ccTCGTTTGGTG‑CGA‑‑‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGc           <  2:1012178/100‑1 (MQ=255)                                                                                       cTCGTTTGGTG‑CGA‑‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGc                  >  1:777851/1‑93 (MQ=38)                                                                                       cTCGTTTGGTG‑CGA‑‑‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCg          >  2:2077830/1‑100 (MQ=255)                                                                                          gTTTGGTG‑CGA‑‑‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGc       <  2:701802/100‑1 (MQ=255)                                                                                            ttGGTG‑CGA‑‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAg        <  1:354424/98‑1 (MQ=255)                                                                                            ttGGGG‑CGA‑‑GGGTTGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGc       <  2:476835/99‑1 (MQ=255)                                                                                             tGGTG‑CGA‑‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGCACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCCTCGCCCacccaccctac     >  1:1782338/1‑90 (MQ=11)                                                                                             tGGTG‑CGA‑‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCaa     <  2:184262/100‑1 (MQ=255)                                                                                              ggTG‑CGA‑‑‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCGACCAATTCCGCCACCTTCGCGCAGTGCGAGCAat    >  2:514815/1‑99 (MQ=14)                                                                                               gTG‑CGA‑‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAata   >  1:1147969/1‑100 (MQ=255)                                                                                               gTG‑CGA‑‑‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAatat  <  2:1248982/100‑1 (MQ=255)                                                                                                        |                                                                                               GTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTT‑GATATCGCCTCGTTTGGTG‑CGA‑‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATAT  >  NC_000913/4606109‑4606302

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