Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,113,401	0	G	T	14.3%	87.8 / 7.0	40	A26D (GCC→GAC)	yghG	secretin (GspDbeta) OM localization lipoprotein pilotin
Reads supporting (aligned to +/- strand):  ref base G (21/12);  new base T (1/4);  total (23/17)
Fisher's exact test for biased strand distribution p-value = 1.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.00e-01
Rejected: Not supported by required number of reads on both strands.

GATACTGGTTACCTGTACGTGTATCATTAATCGTTATGCTGTAATTAATTCCCTCAGTGATCATTAATTTCACCGTCGGGTCAGCGCACATTTGTCGTTGATAGCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGTTTTATCGACATGGTAAATAATCTCTAAGGTTATTAATAAGAGTTAAAATGTCACTTTGATAATGACGTGGTTATCATTAAAACAATGCCTGTAGATAAAGTGTTGCTATACCGCCCAGACTTAAACACGGTCCAAAAGGCAGTGTGGTTGATCCTCTTTTTGTAATAACGGCATATATCAGGCCGCAGCATGAGGCGATTAAAGCAACATTGGGCAGCGACAATGCCCCCA  >  NC_000913/3113111‑3113697                                                                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                                                                         ttcagacgtgtgctcttccgatctCATTAATCGTTATGCTGTAATTAATTCCCTCAGTGATCATTAATTTCACCGTCGGGTCAGCGCACATTTGTCGTTGATAGCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGTgccagccata                                                                                                                                                                                                                                                                                                <  2:414341/278‑3 (MQ=255)   tACTGGTTACCTGTACGTGTATCATTAATCGTTATGCTGTAATTAATTCCCTCAGTGATCATTAATTTCACCGTCGGGTCAGCGCACATTTGTCGTTGATAGCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGTgccagccataat                                                                                                                                                                                                                                                                                              <  2:622208/301‑5 (MQ=255)   tACTGGTTACCTGTACGTGTATCATTAATCGTTATGCTGTAATTAATTCCCTCAGTGATCATTAATTTCACCGTCGGGTCAGCGCACATTTGTCGTTGATAGCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGATATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGTgccagccataat                                                                                                                                                                                                                                                                                              <  2:17875/301‑5 (MQ=255)     ttGGTTACCTGTACGTGTTTTATTAATCGTTATGCTGTAAGTAATTCCCTCAGTGATCATTAATTTCACCGTCGGGTCAGCGCACATTTGTCGTTGATAGCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTtaat                                                                                                                                                                                                                                                                                            <  1:42595/300‑1 (MQ=255)       ggTTACCTGTACGTGTATCATTAATCGTTATGCTGTAATTAATTCCCTCAGTGATCATTAATTTCACCGTCGGGTCAGCGCACATTTGTCGTTGATAGCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCCCAGCCACTtaataa                                                                                                                                                                                                                                                                                          >  1:306424/1‑301 (MQ=255)          tACCTGTACGTTTATCATTAATCGTTATGCTGTAATTAATTCCCTCAGTGATCATTAATTTCACCGTCGGGTCAGCGCACATTTGTCGTTGATAGCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATcc                                                                                                                                                                                                                                                                                       <  2:472889/301‑1 (MQ=255)          tACCTGTACGTGTATCATTAATCGTTATGCTGTAATTAATTCCCTCAGTGATCATTAATTTCACCGTCGGGTCAGCGCACATTTGTCGTTGATAGCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGTTCATTTTTACCGTTATGCCTCTACTTTGCGCCAGCACTAAGGTCTATCCCGCAGATTTAATCTGCAGGTTTTGGCTGATCTTTTTCGCCTGTTTTTTCGCCCGTAACCTGGCATTTTCATTATGGCTTTCACAGCCACTTAATACTca                                                                                                                                                                                                                                                                                       <  2:306828/301‑2 (MQ=255)              ctcctacacgtcgctcttccgatctcTGTAATTAATTCCCTCATTGATCATTAATTTCACCGTCGGGTCAGCGCACATTTGTCGTTGATAGCTGGTTAAAAAGGCGTCAGTCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTa                                                                                                                                                                                                                                                                                   <  1:326011/276‑1 (MQ=255)                   ccctacacgacgctcttccgatctaattaatTCCCTCAGTGATCATTAATTTCACCGTCGGGTCAGCGCACATTTGTCGTTGATAGCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCt                                                                                                                                                                                                                                                                              <  1:300888/278‑1 (MQ=255)                         cATTAATCGTTATGCTGTAATTAATTCCCTCAGTGATCATTAATTTCACCGTCGGGTCAGCGCACATTTGTCGTTGATAGCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGTgccagccataagatcggaagagcgtcgtgaaggg                                                                                                                                                                                                                                                                        >  1:413778/1‑275 (MQ=255)                                       cTGTAATTAATTCCCTCAGTGATCATTAATTTCACCGTCGGGTCAGCGCACATTTGTCGTTGATAGCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCAATTAATAATCCTGTaagatcggaagagcacacgtctgaaa                                                                                                                                                                                                                                                          >  2:326415/1‑277 (MQ=255)                                       cTGTAATTAATTCCCTCAGTGATCATTAATTTCACCGTCGGGTCAGCGCACATTTGTCGTTGATAGCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAACCACTTAATAATCCTATAACGTCcaaccacacacatattaacc                                                                                                                                                                                                                                                          >  1:500620/1‑285 (MQ=255)                                           aattaattCCCTCAGTGATCATTAATTTCACCGTCGGGTCAGCGCACATTTGTCGTTGATAGCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCtagatcggaagagcacacgtctgaa                                                                                                                                                                                                                                                      >  2:301292/1‑277 (MQ=255)                                                                                 gtgggtgtaaTTTGTTGTTGTGAGGTGGTTAAAATGGCGTCTGGCGTTTGCGTGGTTGTGGTACCCGCCTCGTTGATAATGGTTATTTTTACCGTTGTGCCGCGACTTTTGGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATGTTGCGCCTGTTTGTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTATTACCCTCCCTGGCATTTGTTTTATCTACATGGTAAATAATctct                                                                                                                                                                                                                <  1:93841/292‑1 (MQ=255)                                                                                      tactcggctctcttccgttttGGTTAAAAAGTCGTCAGTCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCAGTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTGTCGCCAGTAAACTGGCATTTTCATTATTGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATAGTAATACCCTCCCTGGCATTTGTTTTATCGACATGGTAAATAATCTCTAAGGt                                                                                                                                                                                                           <  1:436458/281‑1 (MQ=255)                                                                                          aTTTGTCGTTGTTAGCTGGTTAAAAAGGCGTCAGGCTTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTGCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGTTTTATCGACATGGTAAATAATCTCTAAGGttatt                                                                                                                                                                                                       <  1:315181/301‑1 (MQ=255)                                                                                            tggagttcagacgtgtgctcttccgatctcaggcgtcTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGTTTTATCGACATGGTAAATAATCTCTAAGGTTATtaa                                                                                                                                                                                                     <  2:125488/273‑1 (MQ=255)                                                                                              gTCGTTGATAGCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATCATCCTGTAACGCTCAACAATAGCGATATTAATACCATCACTGGCATTTGTTTTATCGACATGATAAATAATCTCCAAGGTTATtaaaa                                                                                                                                                                                                   >  1:165929/1‑299 (MQ=255)                                                                                                           ggTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGTTTTATCGACATGGTAAATAATCTCTAAGGtagatcagaagagcaaaagtct                                                                                                                                                                                      >  2:437022/1‑280 (MQ=255)                                                                                                           ggTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGTTTTATCGACATGGTAAATAATCCCTAAGGATATTAATAAGAGTTAAAATGt                                                                                                                                                                                      >  1:157271/1‑301 (MQ=255)                                                                                                                         caggcgtcTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGTTTTATCGACATGGTAAATAATCTATAAGGTTATtaaagatcggaagagcgtcgtgaagggaaaga                                                                                                                                                                        >  1:125364/1‑272 (MQ=255)                                                                                                                                        tttGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGAGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGTTTTATCGACATGGTAAATAATCTCTAAGGTTATTAATAAGAGTTAAAATGTCACTTTGATAATGACGTAGTTATCATTaa                                                                                                                                                         >  2:447280/1‑301 (MQ=255)                                                                                                                                                                 ggTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCACCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGTTTTATCGACATGGTAAATAATCTCTAAGGTTATTAATAAGAGTTAAAATGTCACTTTGATAATGACGTGGTTATCATTAAAACAATGCCTGTAGATAAAGTGTTg                                                                                                                                >  2:280804/1‑301 (MQ=255)                                                                                                                                                                                 gggCCGCTACTTTGGGTCAGTACTAAGGTATATCCCGCTGATTTAATCGGCAGGGTTTGGCTGATATGTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATGTTCATGATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGTTTTATCGACATGGTAAATAATCTCTAAGGTTATTAATAAGAGTTAAAATGTCACTTTGATAATGACGTGGTTATCATTAAAACAATGCCTGTAGATAAAGTGTTGCTATACCGCCCAGACt                                                                                                                <  1:594747/299‑1 (MQ=255)                                                                                                                                                                                  tgCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGTTTTATCGACATGGTAAATAATCTCTAAGGTTATTAATAAGAGTTAAAATGTCACTTTGATAATGACGTGGTTATCATTAAAACAATGCCTGTAGATAAAGTGTTGCTATAACGCCCAGACtt                                                                                                               >  2:355376/1‑301 (MQ=255)                                                                                                                                                                                                    gCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCCTTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCCACAATAGCGATATTAATACCCTCCCTGGCATTTGTTTTATCGACATGGTAAATAATCTCTAAGGTTACTCATAAGAGTTACAATGTACCTTTGATAATGACGTGGTTATACTTAAAACAATGACTGTAAATAAAGTATTTCTATACCGACCAGACTCAAACAAggaaaaaaaagc                                                                                             >  1:464984/1‑291 (MQ=255)                                                                                                                                                                                                       agacgtgtgctcttccgatcttAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGTTTTATCGACATGGTAAATAATCTCTAAGGTTATTAATAAGAGTTAAAATGTCACTTTGATAATGACGTGGTTATCATTAAAACAATGCCTGTAGATAAAGTGTTGCTATACCGCCCAGACTTAAACACGGTCCAAAAGGCAgt                                                                                          <  2:617025/281‑1 (MQ=255)                                                                                                                                                                                                             tataTCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGTTTTATCGACATGGTAAATAATCTCTAAGGTTATTAATAAGAGTTAAAATGTCACTTTGATAATGACGTGGTTATCATTAAAACAATGCCTGTAGATAAAGTGTTGCTATACCGCCCAGACTTAAACACGGTCCAAAAGGCAGTGTGGtt                                                                                    >  2:296949/1‑301 (MQ=255)                                                                                                                                                                                                              atatCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGTTTTATCGACATGGTAAATAATCTCTAAGGTTATTAATAAGAGTTAAAATGTCACTTTGATAATGACGTGGTTATCATTAAAACAATGCCTGTAGATAAAGTGTTGCTATACCGCCCAGACTTAAACACGGTCCAAAAGGCAGTGTGGTTg                                                                                   >  2:610878/1‑301 (MQ=255)                                                                                                                                                                                                                 cgtgtgctcttccgatctcAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGCCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGTTTTATCGACATGGTAAATAATCTCTAAGGTTATTAATAAGAGTTAAAATGTCACTTTGATAATGACGTGGTTATCATTAAAACAATGCCTGTAGATAAAGTGTTGCTATACCGCCCAGACTTAAACACGGTCCAAAAGGCAGTGTGGTTGATc                                                                                <  2:534868/283‑1 (MQ=255)                                                                                                                                                                                                                  cccGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGTTTTATCGACATGGTAAATAATCTCTAAGGTTATTAATAAGAGTTAAAATGTCACTTTGATAATGACGTGGTTATCATTAAAACAATGCCTGTAGATAAAGTGTTGCTATACCGCCCAGACTTAAACACGGTCCAAAAGGCAGTGTGGTTGATcc                                                                               >  2:356198/1‑301 (MQ=255)                                                                                                                                                                                                                            tAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGTTTTATCGACATGGTAAATAATCTCTAAGGTTATTAATAAGAGTTAAAATGTCACTTTGATAATGACGTGGTTATCATTAAAACAATGCCTGTAGATAAAGTGTTGCTATACCGCCCAGACTTAAACACGGTCCAAAAGGCAGTGTGGTTGATCCTCTTTTTGta                                                                     >  2:326474/1‑301 (MQ=255)                                                                                                                                                                                                                            tAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGTTTTATCGACATGGTAAATAATCTCTAAGGTTATTAATAAGAGTTAAAATGTCACTTTGATAATGACGTGGTTATCATTAAAACAATGCCTGTAGATAAAGTGTTGCTATACCGCCCAAACTTAAACACGGTCCAAAAGGAAgtagatcggaaaagcgtcgtggt                                                                     >  1:616167/1‑280 (MQ=255)                                                                                                                                                                                                                                   cAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGCCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGTTTTATCGACATGGTAAATAATCTCTAAGGTTATTAATAAGAGTTAAAATGTCACTTTGATAATGACGTGGTTATCATTAAAACAATGCCTGTAGATAAAGTGTTGCTATACCGCCCAGACTTAAACACGGTCCAAAAGGCAGTGTGGTTGATcagatcggaagagcgtcgg                                                              >  1:534195/1‑283 (MQ=255)                                                                                                                                                                                                                                                       gccataatgaaaatgCCAGTTTACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGTTTTATCGACATGGTAAATAATCTCTAAGGTTATTAATAAGAGTTAAAATGTCACTTTGATAATGACGTGGTTATCATTAAAACAATGCCTGTAGATAAAGTGTTGCTATACCGCCCAGACTTAAACACGGTCCAAAAGGCAGTGTGGTTGATCCTCTTTTTGTAATAACGGCATATATCAGGCCgcagcaa                                          >  2:482263/15‑300 (MQ=255)                                                                                                                                                                                                                                                            tGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGTTTTATCGACATGGTAAATAATCTCTAAGGTTATTAATAAGAGTTAAAATGTCACTTTGATAATGACGTGGTTATCATTAAAACAATGCCTGTAGATAAAGCGTTGCTATACCGCCCAGACTTAAACACGGTCCAAAAGGCAGTGTGGTTGATCCTCTTTTTGTAATAACGGCATATATAAGGCCGCAGCATGAggg                                     >  2:330864/1‑300 (MQ=255)                                                                                                                                                                                                                                                               gaaaatgCCAGTTTACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGTTTTATCGACATGGTAAATAATCTCTAAGGTTATTAATAAGAGTTAAAATGTCACTTTGATAATGACGTGGTTATCATTAAAACAATGCCTGTAGATAAAGTGTTGCTATACCGCCCAGACTTAAACACGGTCCAAAAGGCAGTGTGGTTGATCCTCTTTTTGTAATAACGGCATATATCAGGCCGCAGCATGAAGCGaa                                  >  2:521646/7‑300 (MQ=255)                                                                                                                                                                                                                                                                             aCTGGCTTGTGCATGTTGGCTTGCACATCCACGTAATAATTCTGTATCGCTCATGTATAGCTTTATTAATATCCTCCCTGTCAGTTGTTTTATCGAAATGGTAAATAATCTCTAAGGTTATTAATAAGAGTTAAAATGTGACTTTGATAAGGACGTGGTTATCATTAAAACAATGCCTGTAGATAAAGTGTTGCTATACCGCCCAGACTTAAACACGGTGCAAAAGGCAGTGTGGTTGATCCTCTTTTTGTAATAACGGCATAGATCAGGCCGCAGCATGAGGCGATTAAAGCAACATTgg                    <  1:354892/301‑1 (MQ=255)                                                                                                                                                                                                                                                                                  cATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGTTTTATCGACATGGTAAATAATCTCTAAGGTTATTAATAAGAGTTAAAATGTCACTTTGATAATGACGTGGTTATCATTAAAACAATGCCTGTAGATAAAGTGTTGCTATACCGCCCAGACTTAAACACGGTCCAAAAGGCAGTGTGGTTGATCCTCTTTTTGTAATAACGGCATATAACAGGCCGCAGCATGAGGCGATTAAAGCAACATTGGGCaac               >  1:25757/1‑299 (MQ=255)                                                                                                                                                                                                                                                                                    ttttCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGTTTTATCGACATGGTAAATAATCTCTAAGGTTATTAATAAGAGTTAAAATGTCACTTTGATAATGACGTGGTTATCATTAAAACAATGCCTGTAGATAAAGTGTTGCTATACCGCCCAGACTTAAACACGGTCCAAAAGGCAGTGTGGTTGATCCTCTTTTTGTAATAACGGCATATATCAGGCCGCAGCATGAGGCGATTAAAGCAACATTGGGCAGCGa             <  1:481662/301‑1 (MQ=255)                                                                                                                                                                                                                                                                                            ttggctggcACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGTTTTATCGACATGGTAAATAATCTCTAAGGTTATTAATAAGAGTTAAAATGTCACTTTGATAATGACGTGGTTATCATTAAAACAATGCCTGTAGATAAAGTGTTGCTATACCGCCCAGACTTAAACACGGTCCAAAAGGCAGTGTGGTTGATCCTCTTTTTGTAATAACGGCATATATCAGGCCGCAGCATGAGGCGATTAAAGCAACATTGGGCAGCGACAATGccc     <  2:25808/300‑1 (MQ=255)                                                                                                                                                                                                                                                                                              ggctggcACAGCCACTTAATAATCCTGTAACGCTCAACAATTGCGATATTAATACCCGCCCTGGCATTTGTTTTATCGACATGGTAAATAATCTCTAAGGTTATTAATAAGAGTTAAAATGTCACTTTGATAATGACGTGGTTATCATTAAAACAATGCCTGTAGATAAAGTGTTGCTATACCGCCCAGACTTAAACACGGTCCAAAAGGCAGTGTGGTTGATCCTCTTTTTGTAATAACGGCATATATCAGGCCGCAGCATGAGGCGATTAAAGCAACATTGGGCAGCGACAATGccccc   <  1:520974/301‑1 (MQ=255)                                                                                                                                                                                                                                                                                               gctggcACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGTTTTATCGACATGGTAAATAATCTCTAAGGTTATTAATAAGAGTTAAAATGTCACTTTGATAATGACGTGGTTATCATTAAAACAATGCCTGTAGATAAAGTGTTGCTATACCGCCCAGACTTAAACACGGTCCAAAAGGCAGTGTGGTTGATCCTCTTTTTGTAATAACGGCATATATCAGGCCGCAGCATGAGGCGATTAAAGCAACATTGGGCAGCGACAATGCCccca  >  2:423009/1‑301 (MQ=255)                                                                                                                                                                                                                                                                                               actggcACAGCAATTTAATAATCCTGTAACGTTTAACTATAGCGTTATTAATACCCTCCCTGGCATTTGTTTTATCGACATGGTAAATAATCTCTAAGGTTATTAATAAGAGTTAAAATGTCACTTTGATAATGACGTGGTTATCATTAAAACAATGCCTGTAGATAAAGTGTTGCTATACCGCCCAGACTTAAACACGGTCCAAAAGGCAGTGTGGTTGATCCTCTTTTTGTAATAACGGCATATATCAGGCCGCAGCATGAGGCGATTAAAGCAACATTGGGCAGCGACAATGCCccca  <  1:326070/300‑1 (MQ=255)                                                                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                                                                         GATACTGGTTACCTGTACGTGTATCATTAATCGTTATGCTGTAATTAATTCCCTCAGTGATCATTAATTTCACCGTCGGGTCAGCGCACATTTGTCGTTGATAGCTGGTTAAAAAGGCGTCAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGTTTTATCGACATGGTAAATAATCTCTAAGGTTATTAATAAGAGTTAAAATGTCACTTTGATAATGACGTGGTTATCATTAAAACAATGCCTGTAGATAAAGTGTTGCTATACCGCCCAGACTTAAACACGGTCCAAAAGGCAGTGTGGTTGATCCTCTTTTTGTAATAACGGCATATATCAGGCCGCAGCATGAGGCGATTAAAGCAACATTGGGCAGCGACAATGCCCCCA  >  NC_000913/3113111‑3113697

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑