Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	578,811	0	T	G	52.9%	‑5.6 / 76.1	52	Q28P (CAG→CCG)	borD	prophage lipoprotein BorD
Reads supporting (aligned to +/- strand):  ref base T (14/10);  new base G (14/13);  total (28/24)
Fisher's exact test for biased strand distribution p-value = 7.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.94e-02
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

CCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCTGTTT‑GTTTTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGCCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAAC  >  NC_000913/578575‑579041                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                        ccATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCTCCACAAATTTTAGCTGCAGCGACAGTGTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTT‑GTTTTGAACa                                                                                                                                                                                                                          <  1:277375/251‑1 (MQ=255)    tatCGGTGGACAACTCATGCACTTATTTTGAGCAATACACACGCGCTTCCAGCGGAGCATAAATGACTAACGTAAAAAAACCGAGCAATCCATTTACGACTGTTTGATGGGTTACTGTTTTAACAACATTTTCTGCGCCGCCCCACATTTTAGCTGCATCGCCAGATTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTAATGGTTTCCTTTGGTGCTACTGCTGTCTGTTT‑GTTTTGAACAGt                                                                                                                                                                                                                        >  2:379803/1‑250 (MQ=255)        gATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTT‑GTTTTGAACAGTaaacg                                                                                                                                                                                                                   >  1:50016/1‑249 (MQ=255)                      gCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCTGTTT‑GTTTTGAACAGTAAATGTCTGTTGAGCACAt                                                                                                                                                                                                     <  2:304610/251‑1 (MQ=255)                                         cacacGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCTGTTT‑GTTTTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCa                                                                                                                                                                                  <  1:146296/251‑1 (MQ=255)                                             cgcgCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGCGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCTGTTT‑GTTTTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAgcgc                                                                                                                                                                              >  2:1995/1‑251 (MQ=255)                                              gcgcTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTT‑GTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAgcgc                                                                                                                                                                              >  1:381538/1‑250 (MQ=255)                                                gcTTCCCGCGGAGTATACATGCCTCAAGTCATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTACCAACATTTTCTGCGCAGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTT‑GTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGt                                                                                                                                                                           >  2:188126/1‑251 (MQ=255)                                                gcTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTTCGAATTTTTGCTGGGTTTCTGTTTTAACAACATTTTCGGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTTTCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTT‑GTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGt                                                                                                                                                                           <  1:241752/251‑1 (MQ=255)                                                      aGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTT‑GTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAg                                                                                                                                                                     <  1:227203/251‑1 (MQ=255)                                                            gTATAAATGCCTAAAGAAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGAGCCGCCACAACTTTTAGCTGCATCGACAGTTTTCTTCTGCCAAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTT‑GTTTTGATCAGTAAACGTCTGTTGAGCACATCCTGCAATAAGCAGGGCCAGCGCCGTAGCGAGTAGCAt                                                                                                                                                               >  2:336058/1‑251 (MQ=255)                                                              ataAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCTGTTT‑GTTTTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCAttt                                                                                                                                                             >  2:160687/1‑251 (MQ=255)                                                                aaaTGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCTGTTT‑GTTTTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCAttttt                                                                                                                                                           <  1:357560/251‑1 (MQ=255)                                                                               taaAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTT‑GTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGGTATTc                                                                                                                                            >  2:334608/1‑251 (MQ=255)                                                                                aaaaCCGAGCAAGTCATTTACGAATGGTTGCTGGGTTGCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTT‑GTTTTGAACAGTAAACGTCTGTGGAGCACATCCTGTAATAAGCAGGTCCAGCGCAGTAGCGAGTAGCAGTTTTTTCATGTTGTTATTcc                                                                                                                                           <  2:200834/251‑1 (MQ=255)                                                                                 aaaCCGAGCAATCCATTTACGACTGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCTGTTT‑GTTTTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTccc                                                                                                                                          >  2:269652/1‑251 (MQ=255)                                                                                   aCCGAGCAAGCCATTTACGAATGTTTGCTGGGTGTCTTTTTTAACAACATTTTCTTCGCCGCCACAAATTGTATCTGTAGCGACAGTTTTCTTCTTCCCAATGCCAGAAACGAAGAAATGATGGGTGATGGTTTCGTTTGGTTCTAGTGCT‑TCCGTTTTGTTTTGAACAGTAAAGTTGTGTTGAGCACATCCTGTAATGAGCTGGGTCAGCGCAGTAGCGAGTAGCATTTTTTTCATGTTGTTATTCCCGa                                                                                                                                        <  1:219356/251‑1 (MQ=255)                                                                                        gCAATCCATTTACGAATGTTTGCTGGTTTTCTGTTTTAACAACATGTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCTGTTT‑GTTTTGAACAGGAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAgc                                                                                                                                                                                <  2:390109/206‑1 (MQ=255)                                                                                        gCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCTGTTT‑GTTTTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAgc                                                                                                                                                                                >  1:390109/1‑206 (MQ=255)                                                                                              cATTTACGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTT‑GTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGaa                                                                                                                             <  1:188126/251‑1 (MQ=255)                                                                                               aTTTACGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTT‑GTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGa                                                                                                                              >  1:273277/1‑249 (MQ=255)                                                                                                  tACGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTT‑GTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGt                                                                                                                                                                           <  2:125394/201‑1 (MQ=255)                                                                                                  tACGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTT‑GTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGt                                                                                                                                                                           >  1:125394/1‑201 (MQ=255)                                                                                                     gAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCAGAAATTTTAGCTGCATCGACAGTTTTCTTATGCCCAATTCCAGAAACGGAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTT‑GTTTTGAACAGAAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTTAAGTTCGc                                                                                                                       >  1:51189/1‑250 (MQ=255)                                                                                                          tttGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCTGTTT‑GTTTTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATc                                                                                                                 >  1:282319/1‑251 (MQ=255)                                                                                                             tttGGGTTTCTGGTTTAAAAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCTGTTT‑GTTTTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGt                                                                                                               <  2:4872/248‑1 (MQ=255)                                                                                                                 ggTTTCTGTTTTACCAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTT‑GTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGAAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTAtgtg                                                                                                          >  1:94489/1‑251 (MQ=255)                                                                                                                               caactTTTTCTGAGCCGGCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCTGTTT‑GTTTTGAACAGTAATTGTCTGTTGAGTACATCCTGTAATAAGCAGGGCCAGCTCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATGaaacc                                                                                            <  2:246818/251‑2 (MQ=255)                                                                                                                                          gcgcCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTT‑GTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAAc                                                                                  <  2:351351/250‑1 (MQ=255)                                                                                                                                          gcgcCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTT‑GTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACa                                                                                 <  1:355132/251‑1 (MQ=255)                                                                                                                                              cgccACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCTGTTT‑GTTTTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATg                                                                              >  1:294220/1‑250 (MQ=255)                                                                                                                                                    aaaTTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTAGTGCTAATGCTGTCTGTTT‑GTTTTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGa                                                                        >  1:406051/1‑250 (MQ=255)                                                                                                                                                                cATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTT‑GTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGt                                                           <  2:381538/251‑1 (MQ=255)                                                                                                                                                                    gACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCTGTTT‑GTTTTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAat                                                       <  1:269652/251‑1 (MQ=255)                                                                                                                                                                        gTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTT‑GTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGa                                                                             >  2:181116/1‑225 (MQ=255)                                                                                                                                                                        gTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTT‑GTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTGCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGa                                                                             <  1:181116/225‑1 (MQ=255)                                                                                                                                                                              ttctGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTT‑GTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTACACAAACCCTACACAATGAGTTGAACTTTCATATTGTTAATATTTATTAAt                                             >  2:39121/1‑251 (MQ=255)                                                                                                                                                                                ctGCCCAATTAAAGAAACGAAGAACTGATGGGTGATGGTTTCCTTTGGTACTACTGCTGCCGGTTT‑GTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCCTTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGt                                           >  2:252530/1‑251 (MQ=255)                                                                                                                                                                                     cAATTCAAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGAGCTACTGCTGTCTGTTT‑GCTTTGAACAGTAAATGAATGTTGAGCAAATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATGCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACACTAAACAATGAGATGAAATTTCATATTGTTAATATTTATTAATGTATGc                                       >  2:172584/1‑250 (MQ=255)                                                                                                                                                                                         tCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCTGTTT‑GTTTTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGCCAgg                                   <  1:160687/250‑1 (MQ=255)                                                                                                                                                                                            aGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTT‑GTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGTCAGGTGCg                               >  2:1651/1‑251 (MQ=255)                                                                                                                                                                                            aGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTT‑GTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGTCAGGTGCg                               <  1:334608/251‑1 (MQ=255)                                                                                                                                                                                               aaCGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTT‑GTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGTCAGGTGCGATg                            <  2:94489/251‑1 (MQ=255)                                                                                                                                                                                                  gaagaaATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTT‑GTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTCATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTCTTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGTCAGGTGCGATGAAt                         >  1:188009/1‑251 (MQ=255)                                                                                                                                                                                                    agaaATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCTGTTT‑GTTTTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAAtattt                                                   >  2:338129/1‑223 (MQ=255)                                                                                                                                                                                                    agaaATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCTGTTT‑GTTTTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAAtattt                                                   <  1:338129/223‑1 (MQ=255)                                                                                                                                                                                                                tGATGGTTTCCTTTGGTGCTACTGCTGTCTGTTT‑GTTTTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGCCAGGTGCGATGAATCGTCATTGTATTcc           >  1:396994/1‑251 (MQ=255)                                                                                                                                                                                                                    ggTTTCATTTGGTGCTACTGCTGTCTGTTT‑GTTTTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCAAGAGCAGCAGCGAGGAGCATTTTATTCAAGGGGTTAATCAAGATGCTATTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCATAAACAATGAGATGAAATTTAATATTGATAATATTTATTAATGTATGCCAGGTGAGATGAATAGTCATAGTaatcacgg        >  1:177268/1‑243 (MQ=255)                                                                                                                                                                                                                        tCCTTTGGTGCTACTGCTGCCGGTTT‑GTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGTCAGGTGCGATGAATCGTCATTGTATTCCCGGATTa    <  2:50016/250‑1 (MQ=255)                                                                                                                                                                                                                         ccTTTGGTGCTACTGCTGTCTGTTT‑GTTTTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGCCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAAc  >  1:186705/1‑251 (MQ=255)                                                                                                                                                                                                                                  cTACTGCTGTCTGTTT‑GTTTTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTaa                                                                                               <  1:264781/149‑1 (MQ=255)                                                                                                                                                                                                                                  cTACTGCTGTCTGTTT‑GTTTTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTaa                                                                                               >  2:264781/1‑149 (MQ=255)                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                        CCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCTGTTT‑GTTTTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGCCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAAC  >  NC_000913/578575‑579041

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