Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	925,912	0	G	C	38.9%	144.6 / 34.4	74	noncoding (60/88 nt)	serW	tRNA‑Ser
Reads supporting (aligned to +/- strand):  ref base G (7/37);  new base C (28/0);  total (37/37)
Fisher's exact test for biased strand distribution p-value = 8.90e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTCGAACCCCCGATACG‑TTGCCGTAT‑ACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACCAAACCTCAT‑GGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGAAGAATTTTAATGATAAT  >  NC_000913/925762‑926062                                                                                                                                                         |                                                                                                                                                              cGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTCGAACCCCCg                                                                                                                                                               <  1:252153/151‑1 (MQ=255) cGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTCGAACCCCCg                                                                                                                                                               <  2:458743/151‑1 (MQ=255) cGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTCGAACCCCCg                                                                                                                                                               <  2:167642/151‑1 (MQ=255)    aCAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTCGAACCCCCGATa                                                                                                                                                            <  1:535478/151‑1 (MQ=255)                gAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGGGTCGAACCCCCCCTACG‑TTGCCGttt                                                                                                                                                >  2:781764/1‑148 (MQ=255)                  aaTTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTCGAACCCCCGATACG‑TTGCCGTAT‑aca                                                                                                                                            <  2:1173465/151‑1 (MQ=255)                  aaTTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTCGAACCCCCGATACG‑TTGCCGTAT‑aca                                                                                                                                            <  2:282449/151‑1 (MQ=255)                      cacaGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGGGTTGGACCCCCCACACG‑TTGCCGTAT‑acac‑a                                                                                                                                         >  2:50954/1‑149 (MQ=25)                        caGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTCGAACCCCCGATACG‑TTGCCGTAT‑ACAC‑A‑Cttt                                                                                                                                    <  2:615061/151‑1 (MQ=255)                          ggTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTCGAACCCCCGATACG‑TTGCCGTAT‑ACAC‑A‑CTTTcc                                                                                                                                  <  2:269127/151‑1 (MQ=255)                             aaCATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTCGAACCCCCGATACG‑TTGCCGTAT‑ACAC‑A‑CTTTCCAgg                                                                                                                               <  2:4665/151‑1 (MQ=255)                             aaCATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTCGAACCCCCGATACG‑TTGCCGTAT‑ACAC‑A‑CTTTCCAgg                                                                                                                               <  2:811856/151‑1 (MQ=255)                                  aCTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GG‑GGGGGGGGGTTGAACCCCCCAAACC‑CTTCCCTAT‑ACAC‑ACCTTTCCAg                                                                                                                                >  2:379838/1‑146 (MQ=21)                                   cTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTCGAACCCCCGATACG‑TTGCCGTAT‑ACAC‑A‑CTTTCCAGGCGTGc                                                                                                                          <  2:530494/150‑1 (MQ=255)                                     ccACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGGTTCGAAACCCCCCTAAG‑TTGTCGTGT‑ACAA‑A‑CCTTTCACGCGGGCTcc                                                                                                                       >  1:494476/1‑151 (MQ=21)                                            cccACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTCGAACCCCCGATACG‑TTGCCGTAT‑ACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGc                                                                                                                 <  1:344557/150‑1 (MQ=255)                                            cccACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGA‑GA‑GGGGGGGGGTCGAACCCCCCGTACC‑TTGCCGGAT‑ACAC‑A‑CTTTCCATGCGGGGTCCCTCCGcc                                                                                                                >  2:433086/1‑151 (MQ=18)                                             ccaccaTGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGGTTGGAAACCCCCCGATG‑TTGCCGTAT‑TAAC‑A‑ATTTCCACGCGGGCttc                                                                                                                       >  2:849157/1‑141 (MQ=18)                                              caccaTGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTCGAACCCCCGATACG‑TTGCCGTAT‑ACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCAc                                                                                                              <  2:1204694/151‑1 (MQ=255)                                               accaTGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGGGTTGAACCCCCCGTACG‑ATGCCGCAT‑ACAC‑A‑CCTTTCTTGCGGGCTCCTTCcaccac                                                                                                              >  1:111762/1‑150 (MQ=18)                                                 caTGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTCGAACCCCCGATACG‑TTGCCGTAT‑ACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCg                                                                                                           <  2:1090080/151‑1 (MQ=255)                                                   tGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGGTTCGAACCCCCCATACG‑TTGCCGTAT‑AAAC‑A‑CCTTCCAGGCGGGCTCCTTCAGCCAccccg                                                                                                          >  1:663014/1‑146 (MQ=18)                                                   tGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGA‑TA‑GGGGGGGGGTTTAACCCCCCACAAG‑TTTCCGGAT‑ACAC‑A‑GTTTTCAAACGTGCCCCTTcccccgcccagg                                                                                                         >  1:801741/1‑140 (MQ=16)                                                    gTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTCGAACCCCCGATACG‑TTGCCGTAT‑ACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGac                                                                                                        <  1:620516/151‑1 (MQ=255)                                                      tCACTGCCGTACAAACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTCGAACCCCCGATACG‑TTGCCGTAT‑ACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGaca                                                                                                       <  1:1152603/150‑1 (MQ=255)                                                         cTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTCGAACCCCCGATACG‑TTGCCGTAT‑ACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCt                                                                                                   <  2:235128/150‑1 (MQ=255)                                                          tGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGGGTTGAAACCCCCATACG‑CTGCCGTAT‑ACAC‑A‑CTTTTCAGGCGGGGGCCCTCAGCCCCACGGAGA‑cccc                                                                                                  >  2:184592/1‑148 (MQ=17)                                                                acagacagATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GAGGGGGGGGGTTTAAACCCCCCA‑ACG‑TTGCCCGAT‑ACAC‑A‑CTTTCCAAGGGTGCTCCCTCCACCCCCCGGAAA‑CCCCA‑CCaaa                                                                                           >  1:459452/1‑150 (MQ=14)                                                                acagacagATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTCGAACCCCCGATACG‑TTGCCGTAT‑ACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCaaa                                                                                           <  1:1422937/150‑1 (MQ=255)                                                                  agacagaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTCGAACCCCCGATACG‑TTGCCGTAT‑ACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTg                                                                                        <  1:685775/151‑1 (MQ=255)                                                                  agacagaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGA‑GG‑GGGGGGGGGTCTAAACCCCCCTACG‑TTGCCCTAT‑AAAC‑A‑ACTTCCAGGGGGGCTGCCTCTGCCACCCGGGCA‑CCCCC‑CCCAAtac                                                                                        >  2:273699/1‑149 (MQ=9)                                                                     cagaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTCGAACCCCCGATACG‑TTGCCGTAT‑ACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGttt                                                                                     <  2:572619/151‑1 (MQ=255)                                                                        aTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGGTTCGGACCCCCCCTAAG‑TTTCCGGCT‑AAAC‑A‑CTTTCCCAGCGGGCTCCTTCCGCCACCCGGGCA‑CCCCA‑CCCAAAtttttt                                                                                    >  1:275121/1‑149 (MQ=14)                                                                           aaaTGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTCGAACCCCCGATACG‑TTGCCGTAT‑ACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTAc                                                                               <  1:31949/151‑1 (MQ=255)                                                                            aaTGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGGGTTGAACCCCCCCCACG‑ATTCCGTAT‑AAAC‑A‑ATTTTCAGGCGGGGTCCTTCATCCCCCCCGACA‑CCCCC‑CCAAAATGTTTTTCt                                                                                 >  1:950888/1‑148 (MQ=9)                                                                            aaTGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGA‑GA‑AGGGGGGGGTTGAAACCCCGATACG‑TTTCCCTAT‑ACAC‑A‑CCTTCCCGGGGGGGTCCTTTAGCCCCCCCGACA‑CCCCA‑CCCAAATGTTTTTCTAcc                                                                              >  2:360765/1‑151 (MQ=9)                                                                             aTGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGG‑GA‑GGGGGGGGTTTGGAAACCCCCTAAG‑CTGGCGTAT‑ACAA‑A‑CTTTTCAGGGGGGCTCCTTTCAACACCCGGACA‑CCCCC‑CCAAAATTTTTTTCTatc                                                                              >  1:1261593/1‑148 (MQ=2)                                                                               gCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTCGAACCCCCGATACG‑TGGTCGTAT‑ACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGTCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACCa                                                                             <  1:902409/149‑1 (MQ=255)                                                                               gCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGA‑GA‑GGGGGGGGTTTGAACCCCCCCCACG‑TTGCCGTAT‑ATAA‑A‑CTTTTCACGCGGGCCCCCTCAGCCACCCGGGCA‑CCCCA‑CCCAAATAttttttttcc                                                                              >  2:97725/1‑142 (MQ=2)                                                                                cGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GG‑GGGGGGGGTTCGAACCCCCGCTACG‑TTGCCGGAT‑AAAC‑A‑CTTTCCAGGCGCGCCCCCTCACCCCCCCGGGCA‑CCTCACCCAAATttttttttacccaa                                                                            >  1:1141226/1‑142 (MQ=11)                                                                                 gAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTCGAACCCCCGATACG‑TTGCCGTAT‑ACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACCAAAcc                                                                         <  1:369280/151‑1 (MQ=255)                                                                                  aaaaaaaaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑TA‑GGGGGGGAATCGGACCCCCCCTACG‑TTGCCGGAT‑TAAC‑A‑CCTTTCCGGCGGGCCCCCTCCGCCCCTCCGACA‑CCCCC‑CCAAATTTTTTTTCTTCCAAAccc                                                                        >  1:1328054/1‑150 (MQ=9)                                                                                   aaaaaaaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GG‑GGGGGGGGGTTGGACCCCCCCTACG‑TTGCCGGTT‑AAACAA‑CTTTCC‑GGCGGGGTCCCTTTGCCACCCCGACA‑CCCCACCCAAATATTTTTTTTACCAAAc                                                                          >  2:1383954/1‑149 (MQ=2)                                                                                    aaaaaaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGGGTTGCACCCCCCATACG‑CTGCCCTAT‑ACAC‑C‑CCTTTCCGGCGTGCTCCTTCAGCCACTCGGCCA‑CCCCA‑CCCAATTTTTTTTGTACCcaaac                                                                         >  1:473322/1‑146 (MQ=11)                                                                                     aaaaaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGA‑GA‑AGGGGGGGGTCGAAACCCCCCAATG‑TTGCCCTTT‑AACC‑A‑ATTTCTCGGGGGGCTTCTCCTGCCCCACGCACACCCCCA‑CCAAATttttttctcccacacccc                                                                        >  2:955226/1‑137 (MQ=2)                                                                                       aaaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GG‑GGGGGGGGGTTGAACCCCCCCCACG‑TTGCCGTAT‑ACAC‑A‑CTTTCCAGGCGGGCTCCCTCAGCCACCCCGACA‑CCCCA‑CCCAAATGTTTTTTTACCAAAACTCCT‑gg                                                                  >  2:571660/1‑151 (MQ=14)                                                                                        aaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGGATTGAAACCCCCAAACC‑TTGCCCGATAAAAC‑A‑TTTTCCA‑GCGCGGTCCTTCCGCCCCCCGGGCACCCCCA‑CCAAATT‑TTTTTCTAACAAACCCCAT‑ggg                                                                 >  1:1276919/1‑151 (MQ=9)                                                                                         aGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTCGAACCCCCGATACG‑TTGCCGTAT‑ACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACCAAACCTCAT‑GGGt                                                                <  1:1220307/151‑1 (MQ=255)                                                                                              gCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGGG‑GGGGGGGGTTGGAACCCCC‑CTAAG‑TTTCCGTAT‑ACAC‑A‑CTTTTCCAGCGGGGTCCCTCCGCCCCCCCGGCA‑CCCCC‑CCCAAATGTTTTTGTTCCAAACCCCAT‑ggggggg                                                             >  1:461317/1‑148 (MQ=2)                                                                                              gCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGGGTCGGACCCCCCATACG‑TTGCCGTAT‑ACAC‑A‑CATTCCAGGGGGGCTCCTTCATCCACCCGGACA‑CCCCC‑CCAAAATGTTTTTCTTCCAAAACCCAT‑gggggg                                                              >  1:91099/1‑148 (MQ=11)                                                                                                 cTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GG‑GGGGGGGGTTTGAAACCCCCATACG‑TTTGCGTAT‑ACAA‑A‑CCTTTCCGGGGGGGTCCTTCCGCCACCCCGACA‑ACCCC‑CCCAAATGTTTTTTTACCAAAACCCAT‑TGGGGGGAAccg                                                        >  1:1172657/1‑149 (MQ=2)                                                                                                 cTTTCGTACGAGCTCTTCTTTAAATATGGCGGG‑GAGGGGGGGGGTTGGAACCCCC‑CCACG‑TTGCCGTAT‑AAAC‑A‑CTTTTCCGGCGGGCTCCCTCCGCCCCTCGGACA‑CCCCC‑CCAAAATTTTTTTCTACCCAAACCCAT‑GGGTgggga                                                           >  1:869644/1‑145 (MQ=2)                                                                                                     cGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGGTTTGCACCCCCCCTACG‑GTTCCGGAT‑ACAC‑A‑CTTTCCCGGCGGGCTCCTTTAGCCACTCGGACA‑CCCCA‑CCAAAATGTTTTTCTACCAAAACTCATGGGGGGGCAACgggcg                                                     >  1:1401267/1‑149 (MQ=14)                                                                                                       tACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGGTTTGGACCCCCCCGACGTTTGCCG‑AT‑ACAA‑C‑CTTTTCCGGCGGGCTCCCTCAGCCACCCCGGCA‑CCCCC‑CCAAAATGTTTTTTTAACAAAACCCAC‑ggggggggaaaggg                                                      >  1:147624/1‑139 (MQ=2)                                                                                                       tACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTCGAACCCCCGATACG‑TTGCCGTAT‑ACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACCAAACCTCAT‑GGGTGGCAACGGGgcgc                                                   <  1:1074055/150‑1 (MQ=255)                                                                                                         cGAGCTCTTCTTTAAATATGGCGGG‑GC‑GGGGGGGGCTCTAACCCCCCCAACG‑TTTACGCAT‑ACACAA‑CTTTCCA‑GCGTGCTCCTTTTCCCACACGGGGA‑CCTCA‑CCAAAAAGCTTTTCTACCCAACCTCAT‑TGGGGGaaaaaggggggt                                                  >  1:17781/1‑140 (MQ=2)                                                                                                          gAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTCGAACCCCCGATACG‑TTGCCGTAT‑ACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACCAAACCTCAT‑GGGTGGCAACGGGGCGctact                                               <  1:466007/151‑1 (MQ=255)                                                                                                            gCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTCGAACCCCCGATACG‑TTGCCGTAT‑ACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACCAAACCTCAT‑GGGTGGCAACGGGGCGCTACtat                                             <  1:50954/151‑1 (MQ=255)                                                                                                               cttcttTAAATATGGCGGT‑GA‑GGGGGGGATTCGAACCCCCGATACG‑TTGCCGTAT‑ACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACCAAACCTCAT‑GGGTGGCAACGGGGCGCTACtata                                            <  2:1113192/149‑1 (MQ=255)                                                                                                                ttcttTAAATATGGCGGT‑GA‑GGGGGGGAGTCGAAACCCCCATACG‑TTGTCGTAT‑ACAC‑A‑ATTTCCAGGCGTGCTCCTTCCGCCACCCGGGCC‑CCCCA‑CCCAAATGTTTTTCTATCAAAACTCAT‑GGGGGGGAACCGGGGGCTACTATAggg                                         >  1:539065/1‑151 (MQ=17)                                                                                                                 tcttTAAATATGGCGGT‑GA‑GGGGGGGATTCGAACCCCCGATACG‑TTGCCGTAT‑ACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACCAAACCTCAT‑GGGTGGCAACGGGGCGCTACTATAgg                                          <  1:55397/149‑1 (MQ=255)                                                                                                                          atGGCGGT‑GA‑GGGGGGGATTCGAACCCCCGATACG‑TTGCCGTAT‑ACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACCAAACCTCAT‑GGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAg                                 >  1:113316/1‑149 (MQ=255)                                                                                                                          atGGCGGT‑GA‑GGGGGGGATTCGAACCCCCGATACG‑TTGCCGTAT‑ACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACCAAACCTCAT‑GGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGt                                <  1:247080/150‑1 (MQ=255)                                                                                                                             gCGGT‑GA‑GGGGGGGATTCGAACCCCCGATACG‑TTGCCGTAT‑ACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACCAAACCTCAT‑GGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTaaaa                            <  1:952058/151‑1 (MQ=255)                                                                                                                              cGGT‑GA‑GGGGGGGATTCGAACCCCCGATACG‑TTGCCGTAT‑ACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACCAAACCTCAT‑GGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAAc                           <  2:120239/151‑1 (MQ=255)                                                                                                                                gT‑GA‑GGGGGGGATTCGAACCCCCGATACG‑TTGCCGTAT‑ACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACCAAACCTCAT‑GGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACgg                         <  2:958431/151‑1 (MQ=255)                                                                                                                                    a‑GGGGGGGATTCGAACCCCCGATACG‑TTGCCGTAT‑ACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACCAAACCTCAT‑GGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTc                       >  1:46042/1‑150 (MQ=255)                                                                                                                                    a‑GGGGGGGATTAGAACCCCCGATACG‑TTGCCGTAT‑ACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACCAAACCTCAT‑GGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGGAAAACGGTc                       >  2:209283/1‑150 (MQ=255)                                                                                                                                       ggggggATTCGAACCCCCGATACG‑TTGCCGTAT‑ACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACCAAACCTCAT‑GGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCa                      >  2:1008577/1‑149 (MQ=255)                                                                                                                                         ggggATTCGAACCCCCGATACG‑TTGCCGTAG‑ACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACCAAACCTCAT‑GGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCaagaa                  <  1:819274/151‑1 (MQ=255)                                                                                                                                           ggATTCGAACCCCCGATACG‑TTGCCGTAT‑ACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACCAAACCTCAT‑GGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAagaaga                <  1:1150196/151‑1 (MQ=255)                                                                                                                                               tCGAACCCCCGATCCG‑TTGCCGTAT‑ACAC‑C‑CTTTCCAGGCGTGCCCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACCAAACCTCAT‑GCGTGGCAACGGGGCGCTACTATAGGGCGTTGGAGTAAAACGGTCAAGAAGAAttt            >  2:393836/1‑151 (MQ=255)                                                                                                                                                  aaCCCCCGATACG‑TTGCCGTAT‑ACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACCAAACCTCAT‑GGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGAAGAATTTTaa         <  2:459446/151‑1 (MQ=255)                                                                                                                                                  aaCCCCCGATACG‑TTGCCGTAT‑ACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACCAAACCTCAT‑GGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGAAGAATTTTaa         <  1:1468623/151‑1 (MQ=255)                                                                                                                                                     ccccGATACG‑TTGCCGTAT‑ACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACCAAACCTCAT‑GGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGAAGAATTTTAATGa      <  2:664096/151‑1 (MQ=255)                                                                                                                                                     ccccGATACG‑TTGCCGTAT‑ACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACCAAACCTCAT‑GGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGAAGAATTTTAATGa      <  1:293643/151‑1 (MQ=255)                                                                                                                                                        cGATACG‑TTGCCGTAT‑ACAC‑A‑CTTTCCAGGCGTGCTCCGTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACCAAACCTCAT‑GGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGAAGAATTTTAATGATaa   <  2:1427770/151‑1 (MQ=255)                                                                                                                                                         gATACG‑TTGCCGTAT‑ACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACCAAACCTCAT‑GGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGAAGAATTTTAATGATAat  <  1:881867/151‑1 (MQ=255)                                                                                                                                                         |                                                                                                                                                              CGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑GGGGGGGATTCGAACCCCCGATACG‑TTGCCGTAT‑ACAC‑A‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCA‑CCAAATTGTTTTGCTACCAAACCTCAT‑GGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGAAGAATTTTAATGATAAT  >  NC_000913/925762‑926062

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