Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,222,596	0	G	A	37.9%	83.2 / 32.1	66	F13F (TTC→TTT)	ymgD	uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (29/12);  new base A (0/25);  total (29/37)
Fisher's exact test for biased strand distribution p-value = 2.07e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCAACGCTGGCATTACCCATTACCGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAG‑TAATATAGGTTTCTGA‑CGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGG‑TCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGAACGCCAGTGACAGCACC  >  NC_000913/1222447‑1222745                                                                                                                                                         |                                                                                                                                                      gCAACGCTGGCATTACCCATTACCGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAG‑TAATATAGGTTTCTGA‑CGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGG‑TCATTGGTGCGa                                                                                                                                                      <  2:579841/151‑1 (MQ=255) gCAACGCTGGCATTACCCATTACCGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAG‑TAATATAGGTTTCTGA‑CGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGG‑TCATTGGTGCGa                                                                                                                                                      <  2:582481/151‑1 (MQ=255)   aaCGCTGGCATTACCCATTACCGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAG‑TAATATAGGTTTCTGA‑CGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGG‑TCATTGGTGCGaa                                                                                                                                                     >  1:65915/1‑150 (MQ=255)   aaCGCTGGCATTACCCATTACCGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAG‑TAATATAGGTTTCTGA‑CGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGG‑TCATTGGTGCGaa                                                                                                                                                     >  1:943038/1‑150 (MQ=255)      gCTGGCATTACCCATTACCGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAG‑TAATATAGGTTTCTGA‑CGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGG‑TCATTGGTGCGAaaac                                                                                                                                                  >  2:732832/1‑150 (MQ=255)         ggCATTACCCATTACCGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAG‑TAATATAGGTTTCTGA‑CGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGG‑TCATTGGTGCGAaaacaaac                                                                                                                                              >  2:403272/1‑151 (MQ=255)          gCATTACCCATTACCGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAG‑TAATATAGGTTTCTGA‑CGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGG‑TCATTGGTGCGAAaacaaaca                                                                                                                                             <  2:6806/151‑1 (MQ=255)          gCATTACCCATTACCGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAG‑TAATATAGGTTTCTGA‑CGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGG‑TCATTGGTGCGAAaacaaaca                                                                                                                                             <  1:922643/151‑1 (MQ=255)               ccccattaccGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAG‑TAATATAGGTTTCTGA‑CGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGG‑TCATTGGTGCGAAAACAAACAGGccc                                                                                                                                        >  2:659344/2‑151 (MQ=255)               acccattaccGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAG‑TAATATAGGTTTCTGA‑CGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGG‑TCATTGGTGCGAAAACAAACAGGccc                                                                                                                                        >  1:223818/1‑151 (MQ=255)               acccattaccGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAG‑TAATATAGGTTTCTGA‑CGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGG‑TCATTGGTGCGAAAACAAACAGGccc                                                                                                                                        <  1:919631/151‑1 (MQ=255)                     taccGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAG‑TAATATAGGTTTCTGA‑CGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGG‑TCATTGGTGCGAaaac                                                                                                                                                  >  2:676689/1‑135 (MQ=255)                     taccGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAG‑TAATATAGGTTTCTGA‑CGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGG‑TCATTGGTGCGAaaac                                                                                                                                                  <  1:677074/135‑1 (MQ=255)                     taccGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAG‑TAATATAGGTTTCTGA‑CGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGG‑TCATTGGTGCGAAAACAAACAGGCCCGCCAg                                                                                                                                   <  1:812250/150‑1 (MQ=255)                           aacatTTTCGGCACCTTGTGCGCTGTCATTTTTTGCCATTTCGTGGGTTTCTCCGCAAT‑TAATATAGGTTTCTGA‑CGGCAAACCTTTTTTAATATTGTAGTCCTCCGCCGAGT‑TCATTGGTGCAAAAAAAAACGGGCCCGCCAGTAAAGc                                                                                                                             <  2:945975/150‑1 (MQ=255)                              aaTTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAG‑TAATATAGGTTTCTGA‑CGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGG‑TCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAAt                                                                                                                         >  1:682394/1‑151 (MQ=255)                               aTTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAG‑TAATATAGGTTTCTGA‑CGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGG‑TCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAAtt                                                                                                                        >  1:356800/1‑151 (MQ=255)                                   cTGCAACCTGTGCGCTGTCATTTTTTGCCTTTTCGTTGGCTTCTGCGAAAG‑TAATATAGGTTTCTGA‑CGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGG‑TCATTGGTGCAAAAAAAAACAGGCCCGCCAGTAAAGCAAAtttt                                                                                                                      <  1:121434/149‑1 (MQ=255)                                   cTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAG‑TAATATAGGTTTCTGA‑CGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGG‑TCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAAtttttt                                                                                                                    >  2:554765/1‑151 (MQ=255)                                        aCCTGTGCGCTGTCATTTTTTGCCTTTTCGTTGGTTTCTGCGAAAG‑TAATATGGTTTTCTGA‑CGGCAAACCTTTTTTAATTTTGTAGTCCTGCGCCCAGT‑TTTTTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTtcatc                                                                                                                <  2:294311/150‑1 (MQ=255)                                        aCCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAG‑TAATATAGGTTTCTGA‑CGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGG‑TCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTtcatc                                                                                                                >  1:151869/1‑150 (MQ=255)                                              gcgcTGTCATTTTTTGCTATTTCGTTGGCTTCTGCGCAAG‑TAATATAGGTTTCTGA‑CGGCAACCCTTTTTTAATATTGTAGTCTTGCGCCAATT‑TTATTGGTGCAAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTtcatcatca                                                                                                            <  2:506247/148‑1 (MQ=255)                                                gcTGTCATTTTTTGCCTTTTCGTTGGTTTCTGCGAAAG‑TAATATAGTTTTCTGA‑CGGAAACCCTTTTTAAATATTGTAGTCCTGCCCCCAGG‑TATTGGGTGCAAAAACAAACGGGCCCGCCAGTAAAGCAAATTTTTTcatcatcat                                                                                                           <  1:573327/147‑1 (MQ=255)                                                  tGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAG‑TAATATAGGTTTCTGA‑CGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGG‑TCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCtt                                                                                                      >  1:442575/1‑150 (MQ=255)                                                       ttttttCCTATTTCTTTGTTTTCTGCCAAAG‑TAATATAGTTTTCTGA‑CGGCAACCCTTTTTTAATATTGTAGCCCTGCCCCCATT‑TCTTTGGTAGAAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTc                                                                                                 <  2:494924/149‑1 (MQ=255)                                                        tttttGCCTTTTCGTTGGTTTCTGCGCAAG‑TAATATAGGTTTCTGA‑CGGCAACCCTTTTTAAATATTGTAGTCCTGCGCCCAGG‑TCATTGGTGCAAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTc                                                                                                 <  1:925032/149‑1 (MQ=255)                                                        tttttGCCTTTTCGTTGGTTTCTGCGAAAG‑TAATATAGGTTTCTGA‑CGGCAAACCTTTTTTAATATTGTAGTCCTGCCCCCAGG‑TCATTGGTGAGAAAAAAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTAttt                                                                                                  <  1:733255/148‑1 (MQ=255)                                                          tttGCCATTTCGTTGGCTTCTGCGAAAG‑TAATATAGGTTTCGGA‑CGGAAACCCTTTTTTAATATTGTAGCCCTGCGCCCAGG‑TCATTGGTGCAAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCAtt                                                                                              <  1:946854/150‑1 (MQ=255)                                                            tGCCTTTTCGTTGGCTTCTGCGAAAA‑TAATAGAGGTTTCTGA‑CGAAAAACCTTTTTTAATATTGTACTCCTGCGCCCAGT‑TCTTTGGTGCGAAAAAAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCAttt                                                                                             <  2:359950/149‑1 (MQ=255)                                                              ccATTTGTTGGGTTTTTGCGCAAG‑AAATATGGTTTCCTGA‑CGGCAACCCTTTTTTAATTTTGTAGCCCTCCCCCAAGT‑TCTTTGGTGCAAAAAAAACCGGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTAc                                                                                          <  1:40120/150‑1 (MQ=255)                                                                 tttCTTTGGTTTCTGCGAAAG‑TAATATGGTTTTCTGA‑CGGAAAACCTTTTTAAATATTGTAGTCCGGCCCCCAGG‑TATTTGGTGAAAAAACAAACGGCCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCa                                                                                       <  1:5616/150‑1 (MQ=255)                                                                 tttCGTTGGTTTGTGCGAAAT‑TAATATAGTTTTCGGC‑GGGCAACCTTTTTTTAATTTTATAGCCCCGCCCCCGGG‑TCTTTGGTGCGAAAAAAACCGGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCAtt                                                                                              <  1:468359/143‑1 (MQ=255)                                                                 tttCGTTGGTTTCTGCGCAAG‑TAATATAGGTTTCTGG‑CGGCAACCCTTTTTTAATATTGTAGTCCTGCGCCCAGG‑TCATGGGTGCAAAAAAAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTAccc                                                                                        <  1:286782/149‑1 (MQ=255)                                                                    cGTTGTTTTCTGC‑CAAGATAATA‑AGGTTTCAGGCCGGCAAACCTTTTTAAATTTTTTAGTCC‑CCCCCCAGGTTCTTTGGTGAAAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGaa                                                                                    <  2:424217/150‑1 (MQ=255)                                                                              ctgcgaatgtggtataGGTTTCTGA‑CGGCAAACCGTTTTTAATATTGTAGTCCTGCGACCAGG‑TCATTGGCGCGAACACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCCTTCCTTATTTTATTTTCCCCAAACTTGcaccaa                                                                           >  2:252649/13‑149 (MQ=255)                                                                                cgcAAG‑TAATATAGGTTTCTGA‑CGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGG‑TCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGt                                                                       >  1:241981/1‑151 (MQ=255)                                                                                 caaaaatAATATAGGTTTCTGA‑CGGCAACCCTTTTTAAATATTGTAGTCCGCCCCCCAGG‑TCATTGGTGCAAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCaccacc                                                                          <  1:771146/142‑1 (MQ=255)                                                                                  cAAG‑TAATATAGGTTTCTGA‑CGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGG‑TCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTAc                                                                     >  1:594486/1‑151 (MQ=255)                                                                                    aG‑TAATATAGGTTTCTGA‑CGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGG‑TCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCg                                                                   >  1:290604/1‑151 (MQ=255)                                                                                        aaTATAGTTTTCGGA‑CGGCAACCCTTTTTTAATATTGTAGTCCTGCCCCCAGT‑TCTTGGGTGCAAAAAAAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCa                                                                <  1:140504/151‑1 (MQ=255)                                                                                               gTTTCTGA‑CGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGG‑TCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACgg                                                         >  1:693146/1‑151 (MQ=255)                                                                                                  ttccggcgaaaaaCCTTTTTTATGATTGCTGCCCCGGCCCCAGT‑TCATTGGGGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGc                                                        <  1:866314/140‑1 (MQ=255)                                                                                                     tgaccccAACCCTTTTTTAATATTGTAGCCCGGCCCCCAGT‑TCTTTGGTGCAAAAACAAACGGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCa                                                       <  2:293976/142‑1 (MQ=255)                                                                                                         ggCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGG‑TCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTa                                                 >  1:633897/1‑150 (MQ=255)                                                                                                         gacAACCCGTTTTAAATTTTGGAGTCCGGCCCCGGGG‑ACATTGGTGAAAAAAAAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTaa                                                <  1:485481/149‑1 (MQ=255)                                                                                                            aacccttttttAATTTTGTGGCCCGGCGCCCAGG‑TATTTGGTGAAAAAAAAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTaa                                                <  1:401137/142‑1 (MQ=255)                                                                                                                   ttttAATATTGTAGTCCTGCGCCCAGG‑TCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTc                                      >  1:335954/1‑151 (MQ=255)                                                                                                                   ttttAATATTGTAGCCCGGCCCCCAGG‑TCATTGGTGCAAAAAAAACCAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCtt                                       <  1:555906/150‑1 (MQ=255)                                                                                                                       atttttGTACTCCTCCCCCCGGT‑TCATTGGTCCAAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAg                                   <  2:878053/146‑1 (MQ=255)                                                                                                                       aaTATTGTAGCCCGGCGCCCGGG‑TCATTGGTGCAAAAACAAACGGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCg                                     <  1:440344/148‑1 (MQ=255)                                                                                                                        atatTGTAGTCCTGCGCCCAGG‑TCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAACTTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACAGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGc                                 >  2:271877/1‑151 (MQ=255)                                                                                                                        atatTGTAGTCCTGCGCCCAGG‑TCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGc                                 >  1:871895/1‑151 (MQ=255)                                                                                                                         tatTGTAGCCCTGCGCCCAGT‑TCATTGGTGCAAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGc                                 <  1:119642/150‑1 (MQ=255)                                                                                                                            tGTAGTCCTGCGCCCAGG‑TCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCAtt                             >  1:230336/1‑151 (MQ=255)                                                                                                                              tAGCCCGGCCCCCAGG‑TCATGGGTGAAAAAAAAAACGGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCa                               <  1:495362/147‑1 (MQ=255)                                                                                                                                  cctccccccAGG‑TCTTTGGTGCAAAAAAAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCgg                        <  1:61446/144‑1 (MQ=255)                                                                                                                                  ccTGCGCCCAGG‑TCATTGGTGTGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTAATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCg                         >  2:360369/1‑149 (MQ=255)                                                                                                                                    tGCGCCCAGG‑TCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTc                      >  1:401109/1‑150 (MQ=255)                                                                                                                                      cgcCGAGG‑TCATTGGTCCAAAAACAACCGGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGt                       <  1:268007/147‑1 (MQ=255)                                                                                                                                      cgcCCAGG‑TCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTc                      >  2:749592/1‑148 (MQ=255)                                                                                                                                      cccccAGG‑TCATTGGTGCAAAAACAAACAGCCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTc                      <  1:607650/146‑1 (MQ=255)                                                                                                                                        cccAGG‑TCTTGGGTGCAAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGa                   <  1:445345/149‑1 (MQ=255)                                                                                                                                             g‑TCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGAACGcc              >  1:96644/1‑149 (MQ=255)                                                                                                                                               tCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGAACGCCAg            >  2:262445/1‑150 (MQ=255)                                                                                                                                                 aTTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGAACGCCAGTg          >  2:648992/1‑150 (MQ=255)                                                                                                                                                      tGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGAACGCCAGTGACAGc     >  1:179437/1‑150 (MQ=255)                                                                                                                                                         gAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGAACGCCAGTGACAGCAcc  >  2:61726/1‑150 (MQ=255)                                                                                                                                                         |                                                                                                                                                      GCAACGCTGGCATTACCCATTACCGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAG‑TAATATAGGTTTCTGA‑CGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGG‑TCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGAACGCCAGTGACAGCACC  >  NC_000913/1222447‑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: ‑