Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	925,900	0	A	G	51.7%	63.9 / 87.6	58	noncoding (72/88 nt)	serW	tRNA‑Ser
Reads supporting (aligned to +/- strand):  ref base A (4/24);  new base G (30/0);  total (34/24)
Fisher's exact test for biased strand distribution p-value = 1.60e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.37e-11

GAGTTCACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTC‑GAACCCCCGATACGTT‑GCCGTATACA‑CA‑C‑TTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCA‑CCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGAAG  >  NC_000913/925751‑926047                                                                                                                                                        |                                                                                                                                                          gAGTTCACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGa                                                                                                                                                           <  2:1103827/150‑1 (MQ=255)       aCTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGGTTT‑Gaa                                                                                                                                                    >  2:501545/1‑150 (MQ=255)       aCTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTC‑GAAc                                                                                                                                                   <  1:1121288/151‑1 (MQ=255)       aCTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTC‑GAAc                                                                                                                                                   <  2:1138195/151‑1 (MQ=255)       aCTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGG‑GG‑GGGGGGGGTTC‑Gaa                                                                                                                                                    >  1:775754/1‑150 (MQ=255)        cTGCCGTACAGGCAGCTTAGAAATTCACAGGTACCATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGGGTT‑GAAcc                                                                                                                                                  >  2:733576/1‑151 (MQ=255)             gaaCAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACCGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGGGTT‑GAAccccccc                                                                                                                                             >  1:1165472/3‑149 (MQ=255)               aCAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTC‑GAACCCCCGAt                                                                                                                                            <  2:1218957/150‑1 (MQ=255)                cAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTC‑GAACCCCCGATAc                                                                                                                                          <  2:1098174/151‑1 (MQ=255)                   gCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GT‑GGGGGGGGTTT‑GAAccccccccaagtt                                                                                                                                       >  1:945746/1‑143 (MQ=37)                     aGCTTAGAAATTCACAGGTAACATACTACACCCGCCCCCCATGTTCACTGCCGTACCGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GG‑GGGGGGGGTTT‑GCAACCCCGacaagct                                                                                                                                       >  1:1003404/1‑143 (MQ=25)                       cTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGGTTT‑GAACCCCCCGTACGTT‑gtc                                                                                                                                   >  1:331497/1‑148 (MQ=255)                                tCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GG‑GGGGGGGGGTC‑GAACCCCCCCTAAGTT‑TGCGGATaca‑ca                                                                                                                         >  2:1333904/1‑150 (MQ=25)                                tCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTC‑GAACCCCCGATACGTT‑GCCGTATaca‑ca                                                                                                                         <  2:723197/150‑1 (MQ=255)                                     ggTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTC‑GAACCCCCGATACGTT‑GCCGTATACA‑CA‑C‑TTTcc                                                                                                                 <  1:388889/151‑1 (MQ=255)                                        aaCATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTC‑GAACCCCCGATACGTT‑GCCGTATACA‑CA‑C‑TTTCCAgg                                                                                                              <  2:592429/151‑1 (MQ=255)                                          cATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTC‑GAACCCCCGATACGTT‑GCCGTATACA‑CA‑C‑TTTCCAGGCg                                                                                                            <  1:159293/151‑1 (MQ=255)                                            tACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGGTTC‑GAACCCCCCACACGCT‑GCCCTATACA‑CA‑C‑ATTCCCGGCg                                                                                                            >  1:899086/1‑149 (MQ=21)                                             aCTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTC‑GAACCCCCGATACGTT‑GCCGTATACA‑CA‑C‑TTTCCAGGCGTg                                                                                                          <  1:977697/150‑1 (MQ=255)                                              cTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GC‑TGGGGGGGGTT‑GAAACCCCCATCCCTT‑GCCGTTTTCA‑CACC‑TTTCCggggg                                                                                                            >  1:908708/1‑146 (MQ=21)                                              cTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGGTTT‑GAACCCCCCATACGTT‑TCCGTATTCA‑CA‑C‑CTTTCAGGCggg                                                                                                          >  1:193732/1‑147 (MQ=21)                                                 cACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGG‑GA‑GGGGGGGGGTC‑GAAACCCCGACACGTT‑TCCCTGTACA‑CC‑C‑CTTCCCCGCGGGCTcc                                                                                                      >  1:1286228/1‑150 (MQ=21)                                                  aCCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGA‑GA‑GGGGGGGGGTC‑GACACCCCCCTACGTT‑TCCGCATACA‑AA‑C‑CTTTCCCGCGTGCTcc                                                                                                      >  1:735393/1‑149 (MQ=21)                                                      gCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGGTTC‑GAAACCCCCCTATGTT‑GTCGTATACA‑CA‑C‑TTTCCCGGCGGGCTccct                                                                                                    >  1:1009956/1‑145 (MQ=18)                                                         caccaTGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGGTTT‑GAACCCCCGATACGTT‑GCCGTATACA‑CA‑C‑CTTTCCGGCGTGCCCCTTCAGcccc                                                                                             >  2:509395/1‑149 (MQ=18)                                                               gTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTC‑GAACCCCCGATACGTT‑GCCGTATACA‑CA‑C‑TTTCCAGGCGTGCTCCTTCAGCCACTCGGa                                                                                        <  1:386857/150‑1 (MQ=255)                                                                ttCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGGTTC‑GAACCCCCCACACGTT‑TCCGTATACA‑AA‑C‑TTTTCCGGCGGGCTCCCTCAGCCCCTCCGa                                                                                        >  1:880958/1‑149 (MQ=17)                                                                 tCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTC‑GAACCCCCGATACGTT‑GCCGTATACA‑CA‑C‑TTTCCAGGCGTGCTCCTTCAGCCACTCGGaca                                                                                      <  2:522571/150‑1 (MQ=255)                                                                  cACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTC‑GAACCCCCGATACGTT‑GCCGTATACA‑CA‑C‑TTTCCAGGCGTGCTCCTTCAGCCACTCGGacac                                                                                     <  2:1402633/150‑1 (MQ=255)                                                                     tGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GC‑GGGGGGGGGTCGGAACCCCCCA‑ACCCT‑GTCGTTTACACCC‑C‑TTTCCGGGCGCGC‑CCTTCAGcccccccggcacccc                                                                                  >  1:1116883/1‑137 (MQ=12)                                                                       ccGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGGTTC‑GAAACCCCGATAAGTT‑GTCGTATAAA‑CA‑C‑TTTTCAGGCGTGCTCCTTCAGCCACCCGGACAACCCA‑cc                                                                              >  2:811867/1‑151 (MQ=17)                                                                        cGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGGGTT‑GGACCCCCCCTACGTT‑GTCGTATACA‑AA‑C‑TTTTCAGGCGGGCTCCCTCAGCCACCCCGACAccccccc                                                                               >  1:995228/1‑145 (MQ=17)                                                                         gTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGA‑GA‑AGGGGGGGTTT‑GAACCCCCCATACGTT‑TCCGTATACA‑CA‑C‑ATTCCAGGGGTGGTCCCTCAGCCACCCCGACAccccccc                                                                               >  1:720190/1‑144 (MQ=17)                                                                           acagacagATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGGTTC‑GAACCCCCCATACGTT‑GCCGTATAAA‑CA‑C‑TTTTCCGGCGGGCTCCTTCCGCCACTCGGACAccccccccaa                                                                            >  2:1389860/1‑142 (MQ=14)                                                                           acagacagATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTC‑GAACCCCCGATACGTT‑GCCGTATACA‑CA‑C‑TTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCA‑CCAAAt                                                                          <  1:988049/151‑1 (MQ=255)                                                                              gacagaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTC‑GAACCCCCGATACGTT‑GCCGTATACA‑CA‑C‑TTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCA‑CCAAATTGt                                                                       <  1:1370492/151‑1 (MQ=255)                                                                                 agaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTT‑GAACCCCCGCTACGTT‑GTCGTATACA‑CA‑C‑ATTCCAGGCGTGGTCCTTCAGCCACTCCGACACCCCA‑CCCAAATGtttt                                                                    >  1:764807/1‑151 (MQ=14)                                                                                  gaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTC‑GAACCCCCGATACGTT‑GCCGTATACA‑CA‑C‑TTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCA‑CCAAATTGTTTTg                                                                   <  2:61928/151‑1 (MQ=255)                                                                                           cGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTC‑GAACCCCCGATACGTT‑GCCGTATACA‑CA‑C‑TTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCA‑CCAAATTGTTTTGCTACCAAAc                                                          <  2:183894/151‑1 (MQ=255)                                                                                                 aaaaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGGTTT‑TAACCCCCCCTACGTT‑TCCGTATACA‑CA‑C‑TTTTCAGGCGTGCTCCTTCAGCCACCCGGGCACCCCC‑CCAAATTGTTTTGTTAACCAACCTc                                                       >  1:480438/1‑148 (MQ=14)                                                                                                    aGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGGTTC‑GAACCCCCCCTACTTT‑TTCCTATAAA‑CA‑C‑TTTCCAGGCGGGCTCCCTCCGCCACCCCGACACCCCC‑CCAAAATGTTTTTCTTCCAAACATccttg                                                   >  1:1017517/1‑145 (MQ=9)                                                                                                    aGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGGGTC‑GAACCCCCCCAACGTT‑GCCGTATACA‑CA‑C‑TTTCCAGGCGTGCTCCTTCAGCCCCTCCGGCACCCCA‑CCCAATTGTTTTTCTACCAAACCTCCTggg                                                  >  2:1277805/1‑150 (MQ=14)                                                                                                           aCTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGGTTT‑GGAACCCCCATAAGGT‑TGCGTATACA‑CA‑A‑TTTTCAGGCGTGGTCCCTCACCCACTCGGACACCCCC‑CCCAAATTTTTTTCTTCCAAAACTCCTgtgggggaa                                            >  1:1261685/1‑141 (MQ=2)                                                                                                           aCTTTCGTACGAGCTCTTCTTTAAATATGGCGGC‑GA‑GGGGGGGATTT‑GAAACCCCCATACGTT‑GGCGCTTACA‑AA‑C‑TTTCCAGGCGGGCGCCTTCCGCCACCCCGGCACCCCA‑CCAAATTTTTTTTTTTCCCAACCACATggggg                                                >  1:639582/1‑143 (MQ=2)                                                                                                                cGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGGTTC‑GAACCCCCCATACGTT‑GCCGTATACA‑CA‑C‑TTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCCCA‑CCAAAATGTTTTTCTACCAAACCTCCTGGGGGGCAACgggggg                                     >  2:1081375/1‑149 (MQ=14)                                                                                                                cGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTC‑GAACCCCCGATACGTT‑GCCGTATACA‑CA‑C‑TTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCA‑CCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGgcg                                     <  2:981152/151‑1 (MQ=255)                                                                                                                    cGAGCTCTTCTTTAAATATGGCGGT‑GAGGGGGGGGGTTG‑GAACCCCCG‑TTCGTT‑GTCGTATAAA‑AA‑C‑CTTTCCGGCGTGCCCCTTCTCCCCCCCCGCCACCCCA‑CCCAAATCTTTTGCTCCCCAACCTCCTGGGGGGGCACgggg                                       >  1:419577/1‑145 (MQ=2)                                                                                                                     gAGCTCTTCTTTAAATATGGCGGTGGG‑GGGGGGGGTTG‑GAACCCCC‑CTACCTTGGCCGT‑TACA‑CA‑C‑CTTTTAGGGGGGTTCCTTTACCCCCTCCGGCAACTCC‑CCCAATTTTTTTTTTTCCCAAACTCCTGGGTGGCAACGGGGGGctac                                 >  2:286895/1‑150 (MQ=255)                                                                                                                      aGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGGTTT‑GAACCCCCCATACGTT‑GCCGTATAAA‑CA‑C‑TTTTCAGGCGGGCTCCCTCCGCCACCCGGACACCCCA‑CCCAAATGTTTTTCTACCAAACCTCCTGGGGGGGAACCGGGGGctact                                >  2:208492/1‑150 (MQ=14)                                                                                                                           ttcttTAAATATGGCGGT‑GA‑GGGGGGGGGTT‑TAACACCCCATACTTT‑TTCGCATATA‑CA‑CTTTTCCGGGGGGGCCCCTCCA‑CCACCCCGGCACCTCACCCAAAATTTTTTTTTCCAAAACCCCTTGGGGGGCAACgggcg                                      >  1:511504/1‑138 (MQ=255)                                                                                                                                tAAATATGGCGGT‑GA‑GGGGGGGATTC‑GAACCCCCGATACGTT‑GCCGTATACA‑CA‑C‑TTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCA‑CCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGtt                      <  2:899238/150‑1 (MQ=255)                                                                                                                                              gA‑GGGGGGGATTC‑GAACCCCCGATACGTT‑GCCGTATACA‑CA‑C‑TTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCA‑CCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTc        >  2:764325/1‑151 (MQ=255)                                                                                                                                               a‑GGGGGGGATTC‑GAACCCCCGATACGTT‑GCCGTATACA‑CA‑C‑TTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCA‑CCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTc        >  2:229227/1‑150 (MQ=255)                                                                                                                                                   gggggATTC‑GAACCCCCGATACGTT‑GCCGTATACA‑CA‑C‑TTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCA‑CCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCaaga    <  2:574730/151‑1 (MQ=255)                                                                                                                                                   gggggATTC‑GAACCCCCGATACGTT‑GCCGTATACA‑CA‑C‑TTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCA‑CCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCaaga    <  2:1361049/151‑1 (MQ=255)                                                                                                                                                   gggggATTC‑GAACCCCCGATACGTT‑GCCGTATACA‑CA‑C‑TTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCA‑CCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCaaga    <  2:857548/151‑1 (MQ=255)                                                                                                                                                     gggATTC‑GAACCCCCGATACGTT‑GCCGTATACA‑CA‑C‑TTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCA‑CCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCaagaag  <  1:393901/151‑1 (MQ=255)                                                                                                                                                     gggATTC‑GAACCCCCGATACGTT‑GCCGTATACA‑CA‑C‑TTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCA‑CCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCaagaag  <  1:885304/151‑1 (MQ=255)                                                                                                                                                     gggATTC‑GAACCCCCGATACGTT‑GCCGTATACA‑CA‑C‑TTGCCAGGCGTGCTCCTTCAGCCGCGCGGACACCTCA‑CCAAACTGTTTTGCCATCAAAGCTCATGGGTGGCACCGGGGCGCTACTATAGGGAGTTGGATTAAAACGGTCaagaag  <  1:989757/151‑1 (MQ=21)                                                                                                                                                        |                                                                                                                                                          GAGTTCACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTC‑GAACCCCCGATACGTT‑GCCGTATACA‑CA‑C‑TTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCA‑CCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGAAG  >  NC_000913/925751‑926047

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