Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,222,596	0	G	A	41.7%	85.6 / 43.8	72	F13F (TTC→TTT)	ymgD	uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (38/4);  new base A (0/30);  total (38/34)
Fisher's exact test for biased strand distribution p-value = 2.82e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGCTGGCATTACCCATTACCGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGT‑AATATAGGTTTCTGACGGCAAACCGTTTTT‑AATATTGTAGTCCTG‑CGCCCAGG‑TCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGAACGCCAGTGACAGCACC  >  NC_000913/1222450‑1222745                                                                                                                                                       |                                                                                                                                                      aCGCTGGCATTACCCATTACCGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGT‑AATATAGGTTTCTGACGGCAAACCGTTTTT‑AATATTGTAGTCCTG‑CGCCCAGG‑TCATTGGTGCGaaa                                                                                                                                                    >  2:1270523/1‑150 (MQ=255)     tGGCATTACCCATTACCGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGT‑AATATAGGTTTCTGACGGCAAACCGTTTTT‑AATATTGTAGTCCTG‑CGCCCAGG‑TCATTGGTGCGAaaacaa                                                                                                                                                >  1:218210/1‑150 (MQ=255)      ggCATTACCCATTACCGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGT‑AATATAGGTTTCTGACGGCAAACCGTTTTT‑AATATTGTAGTCCTG‑CGCCCAGG‑TCATTGGTGCGAaaacaaa                                                                                                                                               <  1:683096/150‑1 (MQ=255)       gCATTACCCATTACCGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGT‑AATATAGGTTTCTGACGGCAAACCGTTTTT‑AATATTGTAGTCCTG‑CGCCCAGG‑TCATTGGTGCGAaaacaaac                                                                                                                                              >  1:686263/1‑150 (MQ=255)          ttacccattaccGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGT‑AATATAGGTTTCTGACGGCAAACCGTTTTT‑AATATTGTAGTCCTG‑CGCCCAGG‑TCATTGGTGCGAAAACAAACAGGc                                                                                                                                          >  2:851218/1‑151 (MQ=255)                attaccGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGT‑AATATAGGTTTCTGACGGCAAACCGTTTTT‑AATATTGTAGTCCTG‑CGCCCAGG‑TCATTGGTGCGAAAACAAACAGGCCCGCCa                                                                                                                                    <  2:838629/151‑1 (MQ=255)                 ttaccGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGT‑AATATAGGTTTCTGACGGCAAACCGTTTTT‑AATATTGTAGTCCTG‑CGCCCAGG‑TCATTGGTGCGAAAACAAACAGGCCCGCCa                                                                                                                                    <  1:935400/150‑1 (MQ=255)                  taccGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGT‑AATATAGGTTTCTGACGGCAAACCGTTTTT‑AATATTGTAGTCCTG‑CGCCCAGG‑TCATTGGTGCGAAAACAAACAGGCCCGCCAGt                                                                                                                                  >  1:1107703/1‑151 (MQ=255)                   cccGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGT‑AATATAGTTTTCTGACGGAAACCCTTTTTT‑AATATTGTATTCCTG‑CCCCCAGG‑TCATTGGTGCAAAAACAAACAGGCCCGCCAGTa                                                                                                                                 <  1:827960/150‑1 (MQ=255)                    ccGCAAAAATTTCTGCACCCTGTGCGCTGTCAGTTTTTCCCTTTTCTTTGGTTTCTGCGCAAGT‑AATATGGTTTTCTGACGGCAAACCTTTTTT‑AATATTGTATTCCTG‑CGCCCAGG‑TCTTGGGTGCAAAAACAAAAAGCCCCGCCAGTa                                                                                                                                 <  1:1163700/150‑1 (MQ=255)                     cGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGT‑AATATAGGTTTCTGACGGCAAACCGTTTTT‑AATATTGTAGTCCTG‑CGCCCAGG‑TCATTGGTGCGAAAACAAACAGGCCCGCCAGTaa                                                                                                                                >  2:1608580/1‑150 (MQ=255)                          caaTTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGT‑AATATAGGTTTCTGACGGCAAACCGTTTTT‑AATATTGTAGTCCTG‑CGCCCAGG‑TCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCaa                                                                                                                           >  2:588691/1‑150 (MQ=255)                           attttCTGCACCCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGTTTCTGCGAAAGT‑AATATAGGTTTCGGACGGCAACCCTTTTTT‑AATATTGTAGCCCTG‑CCCCCAGG‑TCATGGGTGCAAAAACAAACAGGCCCGCCAGTAAAGCAAAt                                                                                                                         <  2:261536/149‑1 (MQ=255)                             tttCTGCACCCTGTGCGCTGTCAGTTTTTGCTTTTTCTTTGGTTTCTGCGAAAAA‑AATATGGTTTCCGGACGGAAAACCTTTTTT‑AATTTTGTAGCCCGG‑CGCCCGGG‑TCATGGGGGCAAAAAAAAACAGGCCCGCCAGTAAAg                                                                                                                              <  1:1242934/144‑1 (MQ=255)                                    acctgttccgcTGTTATTTTTTGCCATTTCGTTGGTTTCTGCGAAAGA‑AATATGGGTTTCTGACGGCAATCTTTTTTA‑AATTTGTTATTCCGG‑CCCCCAGT‑TTATGGGTGCAAAAACAAACGCGCCCGCCAGTAAAGCAAATTTTTtcatc                                                                                                                <  1:1699965/143‑1 (MQ=255)                                    aaCCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGT‑AATATAGGTTTCTGACGGCAAACCGTTTTT‑AATATTGTAGTCCTG‑CGCCCAGG‑TCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTtcatc                                                                                                                >  1:1093242/1‑151 (MQ=255)                                     aCCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGT‑AATATAGGTTTCTGACGGCAAACCGTTTTT‑AATATTGTAGTCCTG‑CGCCCAGG‑TCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTtcatc                                                                                                                >  1:721170/1‑150 (MQ=255)                                      ccggtgCGCTGTAATTTTTTGCCATTTCGTGGGCTTCT‑CGAAAATAAATAAAGTTTTCTGACGCCAATCCGTTTTT‑AATTTTGTGTTCCGA‑CCCCCGGT‑TCTTGGGGGAAAAACCAAACGGGCCCGCCAGTAAAGCAAATTTTTtcatca                                                                                                               <  1:539636/147‑1 (MQ=255)                                         gtgCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGT‑AATATAGGTTTCTGACGGCAAACCGTTTTT‑AATATTGTAGTCCTG‑CGCCCAGG‑TCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTcatcatcat                                                                                                           >  2:1275351/1‑151 (MQ=255)                                          tgCGCTGTCATTTTTGGCCTTTTCGTTGGTTTCTGCGAAATA‑AATAAAGTTTTCTGCCGGAAATCCTTTTTT‑AATTTTGTACCCCTC‑CGCCAAGT‑TCTTTGGTGGAAAAAAAAACAGGCCCGCCAGTAAAGCAAAttttt                                                                                                                     <  2:905457/140‑1 (MQ=255)                                           gcgcTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGT‑AATATAGGTTTCTGACGGCAAACCGTTTTT‑AATATTGTAGTCCTG‑CGCCCAGG‑TCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTc                                                                                                         >  1:1120842/1‑151 (MQ=255)                                             gcTGTCAGTTTTTGCCTTTTCGTTGGTTTCTGCGAAAGT‑AATATAGTTTTCTGACGGAAACCCTTTTTT‑AATATTGTAGCCCTG‑CCCCCAGG‑TCATTGGTGCAAAAACAAACGGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTc                                                                                                         <  2:622931/149‑1 (MQ=255)                                             gcTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGT‑AATATAGGTTTCTGACGGCAAACCGTTTTT‑AATATTGTAGTCCTG‑CGCCCAGG‑TCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCt                                                                                                       >  1:1373818/1‑151 (MQ=255)                                                 tttttttttCCCTTTTCGTTGTCTTCGCCACAAGT‑AATATAGTTTTCGGAC‑GAAAACCTTTTTTAAATTTTGTAGCCCGG‑CCCCCATT‑TTTTTGGTGAAAAAAAAAACAGCCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTAt                                                                                                    <  1:1082337/146‑1 (MQ=255)                                                 tCAGTTTTTGCCATTTCGTTGGCTTCTGCGCACGT‑AATATAGGTTTCTGACGGCAAACCGTTTTT‑AATATTGTAGTCCTG‑CGCCCAGG‑TCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTAtt                                                                                                   >  1:1698283/1‑151 (MQ=255)                                                    ttttttGCCTTTTCGTTGGTTTCTGCGAAAGT‑AATATAGTTTTCTGCCGGCAACCCTTTTTT‑AATATTGTATCCCGG‑CCCCCGGG‑TCATTGGTGCAAAAAAAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTAttt                                                                                                  <  2:1216632/148‑1 (MQ=255)                                                     tttttGCTTTTTCGTGGGTTTCGGCGAAAGA‑AATATGGTTTTCGGACGGCAAACTTTTTTT‑AATATTGTAGTCCGG‑CCCCCGGT‑TTATTGGTGCAAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTAttt                                                                                                  <  2:179950/148‑1 (MQ=255)                                                     tttttGCCTTTTGGTTGGTTTCTGCGAAAGT‑AATATAGGTTTCGGACGGCAACCCTTTTTT‑AATATTGTAGCCCTG‑CGCCCAGG‑TCATTGGTGAAAAAAAAAACAGCCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCa                                                                                                <  2:1341484/150‑1 (MQ=255)                                                     tttttGCCTTTTCGTTGGTTTCTGCGCAATT‑AATATAGTTTTCTGCCGGCAACCCTTTTTT‑AATATTGTAGTCCTG‑CGCCCAGG‑TCTTTGGTGCAAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCtt                                                                                                      <  2:858325/144‑1 (MQ=255)                                                      ttttGCCATTTCGTTGGCTTCTGCGCAAGT‑AATATAGGTTTCTGACGGCAAACCGTTTTT‑AATATTGTAGTCCTG‑CGCCCAGG‑TCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCAtt                                                                                              >  2:1283075/1‑151 (MQ=255)                                                       tttCCCTTTTCGTTGGTTTCTGCGAAAAT‑AATATAGTTTTCTGAC‑GCAAACCTTTTTTAAATATTGTACTCCTC‑CCCCCAGG‑TCTTTGGTGCAAAAAAAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCAtt                                                                                              <  1:389500/150‑1 (MQ=255)                                                        ttGCCTTTTCGTTGGTTTCTGCGAAATT‑AATATGGTTTTCTGACGGCAAACCTTTTTT‑AATATGGTAGTCCTG‑CGCCAAGG‑TCATTGGTGCAAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCAttt                                                                                             <  2:1293401/150‑1 (MQ=255)                                                           ccttttCTTTGTTTTCTGCGAAAGA‑AATATAGTTTTCTGAGGCCAACCCTTTTTT‑AATATTGTAGTCCGG‑CGCCCAGG‑TCATGGGTGCAAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTAcc                                                                                         <  2:9956/148‑1 (MQ=255)                                                           ccATTTCGTTGGCTTCTGCGCAAGT‑AATATAGGTTTCTGACGGCAAACCGTTTTT‑AATATTGTAGTCCTG‑CGCCCAGG‑TCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTAcc                                                                                         >  2:1099234/1‑151 (MQ=255)                                                              tttCGTGGGTTTCTGCGAAAGA‑AATATAGTTTTCTGACGGCAACCCGTTTTT‑AATTTTGTAGCCCTG‑CGCCCAGT‑TCTTTGGGGCAAAAAAAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTa                                                                                           <  2:1238828/146‑1 (MQ=255)                                                                   ttGGTTTCTGCGAAATA‑AATATAGGTTTCTGACGGCAAACCTTTTTT‑AATTTTGTAGTCCGG‑C‑CCCAGGTTTTTTGGTGCAAAAAAAAACAGCCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGa                                                                                     <  1:1029527/147‑1 (MQ=255)                                                                    tgttttCTGCGAAAGT‑AATATAGGTTTCTGACGGCAAACCTTTTTT‑AATTTTGTAGTCCTG‑CGCCCAGT‑TCTTTGGTGCAAAAAAAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTg                                                                                 <  1:819321/146‑1 (MQ=255)                                                                    tGGTTTCTGCGCAAGT‑AATATAGGTTTCTGACGACAACCCGTTTTT‑AATATTGTACTCCTG‑CGCCCAGG‑TCATTGGTGCAAAAACAAACATCCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTg                                                                                 <  1:1290668/150‑1 (MQ=255)                                                                    tGGCTTCTGCGCAAGT‑AATATAGGTTTCTGACGGCAAACCGTTTTT‑AATATTGTAGTCCTG‑CGCCCAGG‑TCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGc                                                                                >  2:679223/1‑151 (MQ=255)                                                                      gCTTCTGCGCAAGT‑AATATAGGTTTCTGACGGCAAACCGTTTTT‑AATATTGTAGTCCTG‑CGCCCAGG‑TCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGcac                                                                              >  1:58923/1‑151 (MQ=255)                                                                                aaGT‑AATATAGGTTTCTGACGGCAAACCGTTTTT‑AATATTGTAGTCCTG‑CGCCCAGG‑TCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTAcc                                                                    >  1:591719/1‑151 (MQ=255)                                                                                aaGA‑AATATAGGTTTCTAACGGCAACCCTTTTTA‑AATTTTGTAGCCCTG‑CCCCCAGG‑TCTTTGGTGCAAAAAAAAACAGCCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTa                                                                      <  1:1194676/149‑1 (MQ=255)                                                                                     aaTATAGGTTTCTGACGGCAAACCGTTTTT‑AATATTGTAGTCCTG‑CGCCCAGG‑TCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGcc                                                                 >  1:175752/1‑150 (MQ=255)                                                                                       tataGGTTTCTGACGGCAAACCGTTTTT‑AATATTGTAGTCCTG‑CGCCCAGG‑TCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAAt                                                              >  1:276118/1‑151 (MQ=255)                                                                                        ataGGTTTCTGACGGCAAACCGTTTTT‑AATATTGTAGTCCTG‑CGCCCAGG‑TCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATc                                                             >  2:958510/1‑151 (MQ=255)                                                                                             tttCTGCC‑GCAAACCGTTTTT‑AAATTTGGTGCTCTGCCGCCCGGG‑TATTTGGTGAAAAAACAAACGGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCa                                                            <  1:982894/147‑1 (MQ=255)                                                                                                    cGGCAACCCTTTTTT‑AATATTGTAGTCCTG‑CGCCCAGG‑TCATTGGTGCAAAAACAAACAGCCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCttt                                                  <  2:223113/150‑1 (MQ=255)                                                                                                     ggaaaaCCTTTTTA‑AATATTGTAGCCCGG‑CCCCCAGG‑TCTTTGGGGCAAAAAAAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCttt                                                  <  2:1492366/146‑1 (MQ=255)                                                                                                        aaaCCGTTTTT‑AATATTGTAGTCCTG‑CGCCCAGG‑TCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCg                                             >  1:1622089/1‑151 (MQ=255)                                                                                                            cGTTTTT‑AATATTGTAGTCCTG‑CGCCCAGG‑TCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGcccc                                         >  2:328279/1‑151 (MQ=255)                                                                                                               tttt‑AATATTGTAGTCCTG‑CGCCCAGG‑TCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTc                                      >  1:1298799/1‑151 (MQ=255)                                                                                                                ttt‑AATATTGTAGTCCTG‑CGCCCAGG‑TCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCg                                     >  1:438549/1‑151 (MQ=255)                                                                                                                ttt‑AATATTGTAGTCCGG‑CCCCCGGG‑TCTTTGGTGAAAAAAAAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTc                                      <  2:341717/150‑1 (MQ=255)                                                                                                                    aaTTTTGTAGTCCTG‑CGCCCAGT‑TCTTGGGTGCAAAAAAAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAgg                                  <  2:1026199/151‑1 (MQ=255)                                                                                                                      tatTGTAGTCCTG‑CGCCCAGG‑TCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGcc                                >  1:1057187/1‑151 (MQ=255)                                                                                                                       atTGTAGTCCTG‑CGCCCAGG‑TCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGcc                                >  1:1464571/1‑150 (MQ=255)                                                                                                                        ttttgtcgcccgcccccAGGTTCTTTGGTGAAAAAAAAAACGGCCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAgg                                  <  1:92778/135‑1 (MQ=255)                                                                                                                          tgtacccccgcccccAGG‑TTTTTGGTGAAAAAAAAAACGGCCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGcc                                <  1:628364/136‑1 (MQ=255)                                                                                                                                ccggcccccAGG‑TCATTGGTGCAAAAAAAACCAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCgg                        <  2:261936/144‑1 (MQ=255)                                                                                                                                 ccggcccccggattaTTGGTGAGAAAAAAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGt                       <  2:1350807/137‑1 (MQ=255)                                                                                                                                    cgcCCAGG‑TCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTc                      >  1:1057395/1‑148 (MQ=255)                                                                                                                                     gcCCAGG‑TCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGa                   >  1:1262710/1‑150 (MQ=255)                                                                                                                                      cccAGG‑TCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGAAc                 >  2:833768/1‑151 (MQ=255)                                                                                                                                       ccAGG‑TCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGAACg                >  2:1631055/1‑151 (MQ=255)                                                                                                                                       ccAGG‑TCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGAACg                >  2:946105/1‑151 (MQ=255)                                                                                                                                        cAGG‑TCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGAACGc               >  1:1186384/1‑151 (MQ=255)                                                                                                                                         aGG‑TCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGAACGcc              >  2:1672769/1‑151 (MQ=255)                                                                                                                                         aGG‑TCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGAACGcc              >  2:194732/1‑151 (MQ=255)                                                                                                                                            ttattgggTGCAAAAAAAACCGGCCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGAACGCCa             <  1:336615/144‑1 (MQ=255)                                                                                                                                             tCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGAACGCCAGt           >  1:84675/1‑151 (MQ=255)                                                                                                                                              cATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGAACGCCAGTg          >  1:970117/1‑151 (MQ=255)                                                                                                                                                     gCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGAACGCCAGTGACAGCAc   >  2:623862/1‑151 (MQ=255)                                                                                                                                                       gAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGAACGCCAGTGACAGCAcc  >  2:1512816/1‑150 (MQ=255)                                                                                                                                                       |                                                                                                                                                      ACGCTGGCATTACCCATTACCGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGT‑AATATAGGTTTCTGACGGCAAACCGTTTTT‑AATATTGTAGTCCTG‑CGCCCAGG‑TCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGAACGCCAGTGACAGCACC  >  NC_000913/1222450‑1222745

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