Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	888,014	0	T	G	32.6%	95.4 / 11.3	49	noncoding (44/79 nt)	rybB	small regulatory RNA RybB
Reads supporting (aligned to +/- strand):  ref base T (9/24);  new base G (15/1);  total (24/25)
Fisher's exact test for biased strand distribution p-value = 1.06e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AAACCGCTATTAAAAACGGTAATGCAAAACTGGATGCAGCGCAGTCAGCAAACGCTCGAACAATGGTTTGAACCCGGAACCGCCCGCGCGCTTGATGCGTTTATTGAGGGGATGACGCTGCATTTTGTCACCGACCGTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGG‑AAAAGTTCTGCGGTTGTGCAAAAAAATTTCATTTTCAGGGCAACTTCAGTTTTATCCTAATCCTGGCCATACCATGACGATGATTGTCCCTGCCAGCGTCAGCAGGACGTTGGCGATTGCATAGGTGCCCGCATAGCCCAGCGCCGGGATGTTACTGCGAGCTGTATCACTGATGATCTCC  >  NC_000913/887768‑888263                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                           aaaCCGCTATTAAAAACGGTAATGCAAAACTGGATGCAGCGCAGTCAGCAAACGCTCGAACAATGGTTTGAACCCGGAACCGCCCGCGCGCTTGATGCGTTTATTGAGGGGATGACGCTGCATTTTGTCACCGACCGTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGAt                                                                                                                                                                                                                                                        <  1:252626/251‑1 (MQ=255) aaaCCGCTATTAAAAACGGTAATGCAAAACTGGATGCAGCGCAGTCAGCAAACGCTCGAACAATGGTTTGAACCCGGAACCGCCCGCGCGCTTGATGCGTTTATTGAGGGGATGACGCTGCATTTTGTCACCGACCGTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGAt                                                                                                                                                                                                                                                        <  1:260306/251‑1 (MQ=255)   aCCGCTATTAAAAACGGTAATGCAAAACTGGATGCAGCGCAGTCAGCAAACGCTCGAACAATGGTTTGAACCCGGAACCGCCCGCGCGCTTGATGCGTTTATTGAGGGGATGACGCTGCATTTTGTCACCGACCGTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATgg                                                                                                                                                                                                                                                      <  1:161054/251‑1 (MQ=255)         aTTAAAAACGGTAATGCAAAACTGGATGCAGCGCAGTCAGCAAACGCTCGAACAATGGTTTGAACCCGGAACCGCCCGCGCGCTTGATGCGTTTATTGAGGGGATGACGCTGCATTTTGTCACCGACCGTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTcc                                                                                                                                                                                                                                                 <  2:33263/250‑1 (MQ=255)                 cGGTAATGCAAAACTGGATGCAGCGCAGTCAGCAAACGCTCGAACAATGGTTTGAACCCGGAACCGCCCGCGCGCTTGATGCGTTTATTGAGGGGATGACGCTGCATTTTGTCACCGACCGTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGGTGATAGGGGCCACaaaaa                                                                                                                                                                                                                                          >  2:30026/1‑248 (MQ=255)                  ggTAATGCAAAACTGGATGCAGCGCAGTCAGCAAACGCTCGAACAATGGTTTGAACCCGGAACCGCCCGCGCGCTTGATGCGTTTATTGAGGGGATGACGCTGCATTTTGTCACCGACCGTAAGCCGCTATCGCGCGAGGAGATTTTGAGGCTGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGGTTATGGGCTCCACAAAATggg                                                                                                                                                                                                                                       >  1:394916/1‑251 (MQ=255)                      aTGCAAAACTGGATGCAGCGCAGTCAGCAAACGCTCGAACAATGGTTTGAACCCGGAACCGCCCGCGCGCTTGATGCGTTTATTGAGGGGATGACGCTGCATTTTGTCACCGACCGTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGGTGATGGGCTCCACAAAAAGGGGAc                                                                                                                                                                                                                                    >  1:327590/1‑250 (MQ=255)                        gCAAAACTCGATGCAGCGCAGTCAGCAAACGCTCGAACAATGGTTTGAACCCGGAACCGCCCGCGCGCTTGATGCGTTTATTGAGGGGATGACGCTGCATTTTGTCACCGACCGTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGGTGAGGGGGTCCCCAAAAAggggggac                                                                                                                                                                                                                                  >  2:350402/1‑246 (MQ=255)                         cAAAACTGGATGCAGCGCAGTCAGCAAACGCTCGAACAATGGTTTGAACCCGGAACCGCCCGCGCGCTTGATGCGTTTATTGAGGGGATGACGCTGCATTTTGTCACCGACCGTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCa                                                                                                                                                                                                                                <  2:432578/251‑1 (MQ=255)                          aaaaCTGGATGCAGCGCAGTCAGCAAACGCTCGAACAATGGTTTGAACCCGGAACCGCCCGCGCGCTTGATGCGTTTATTGAGGGGATGACGCTGCATTTTGTCACCGACCGTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGGTGATGGGCTCCCCAACAAGGGGACATCa                                                                                                                                                                                                                                >  2:385169/1‑250 (MQ=255)                                              cAGCAAACGCTCGAACAATGGTTTGAACCCGGAACCGCCCGCGCGCTTGATGCGTTTATTGAGGGGATGACGCTGCATTTTGTCACCGACCGTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAAt                                                                                                                                                                                                           <  1:571564/251‑1 (MQ=255)                                                                          ccGGAACCGCCCGCGCGCTTGATGCGTTTATTGAGGGGATGACGCTGCATTGTGTCACCGACCGTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGAGTGATGCCCTGCGG‑AAAAGTTCt                                                                                                                                                                              <  2:467672/251‑1 (MQ=255)                                                                              aaCCGCCCGCGCGCTTGATGCGTTTATTGAGGGGATGACGCTGCATTTTGTCACCGACCGTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGGTGATGGGCTCCACAAAAAGGGGAAATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGG‑AAAAGATCTGCgg                                                                                                                                                                          >  1:263491/1‑251 (MQ=255)                                                                               aCCGCCCGCGCGCTTGATGCGTTTATTGAGGGGATGACGCTGCATTTTGTCACCGACCGTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGGCTGATGCCCTGCGG‑AAAAGTTCTGCGGt                                                                                                                                                                         <  1:78506/251‑1 (MQ=255)                                                                                          gcTTGATGCGTTTATTGAGGGGATGACGCTGCATTTTGTCACCGACCGTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGG‑AAAAGTTCTGCGGTTGTGCaaaaaa                                                                                                                                                              <  1:324874/251‑1 (MQ=255)                                                                                             tGATGCGTTTATTGAGGGGATGACGCTGCATTTTGTCACCGACCGTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGGGGGACTAATTAAGACTGATGCCCTGCGGGAAAAGT‑CTGCGGGTGTGCAAAAAAAtt                                                                                                                                                           >  1:257123/1‑251 (MQ=255)                                                                                               aTGCGTTTATTGAGGGGATGACGCTGCATTTTGTCACCGACCGTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGG‑AAAAGTTCTGCGGTTGTGCAAAAAAATTTc                                                                                                                                                         <  1:158639/251‑1 (MQ=255)                                                                                                     ttATTGAGGGGATGACGCTGCATTTTGTCACCGACCGTAAGCCCCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGATGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGAGGGGATCCACAAAAAGGGGACATACAAGAAAAGCAATGGCATTAATTAAGACTGATGCCCTGCGG‑AAAAGTTCTGCGGGTGTGCAAAAAAAATTCAttttt                                                                                                                                                   >  2:174376/1‑250 (MQ=255)                                                                                                      tATTGAGGGGATGACGCTGCATTTTGTCACCGACCGTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGG‑AAAAGTTCTGCGGTTGTGCAAAAAAATTTCATTTTc                                                                                                                                                   <  2:558302/250‑1 (MQ=255)                                                                                                            ggggATGACGCTGCATTTTGTCACCGACCGTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGG‑AAAAGTTCTGCGGTTGTGCAAAAAAATTTCATTTTCAGGGCaa                                                                                                                                            <  2:69394/251‑1 (MQ=255)                                                                                                                 tGACGCTGCATTTTGTCACCGACCGTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGGTGATGGGCTCCACAAAATGGGGGCACCAAAGAAAAGCAGGGGCACTAATTAAGACTGATGGCCTGCGG‑AAAAGTTTTGCGGGTGTTCAAAAAAATTTCCTTTTTAGGGGAAcctc                                                                                                                                        >  2:544634/1‑247 (MQ=255)                                                                                                                                   ccgaccgTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGG‑AAAAGTTCTGCGGTTGTGCAAAAAAATTTCATTTTCAGGGCAACTTCAGTTTTATCCTAATCCTgg                                                                                                                     <  1:174552/251‑1 (MQ=255)                                                                                                                                       ccgTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGG‑AAAAGTTCTGCGGTTGTGCAAAAAAATTTCATTTTCAGGGCAACTTCAGTTTTATCCTAATCCTGGCCAt                                                                                                                 <  1:350402/251‑1 (MQ=255)                                                                                                                                           aaGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGG‑AAAAGTTCTGCGGTTGTGCAAAAAAATTTCATTTTCAGGGCAACTTCAGTTTTATCCTAACCCTGGCCATAcc                                                                                                              <  1:431643/250‑1 (MQ=255)                                                                                                                                            aGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAAAGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGG‑AAAAGTTCTGCGGGTGTGCAAAAAAATTTCATTTTCAGGGCAACTTCAGTTTTATCCTAATCCTGGCCATACCAt                                                                                                            >  1:10840/1‑251 (MQ=255)                                                                                                                                                              agaTTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCCCAAAATGGGGGCATAAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGG‑AAAAGATCTGCGGTTGTGCAAAAAAATTTCATTTTCAGGGCCACTTCAGTTTTATCCTAATCCTGGCCATACCATGACGATGATTGTCCCTGc                                                                                          >  1:342040/1‑251 (MQ=255)                                                                                                                                                              agaTTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGG‑AAAAGTTCTGCGGTTGTGCAAAAAAATTTCATTTTCAGGGCAACTTCAGTTTTATCCTAATCCTGGCCATACCATGACGATGATTGTCCCTGc                                                                                          <  1:236923/251‑1 (MQ=255)                                                                                                                                                                               gagaGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGG‑AAAAGTTCTGCGGTTGTGCAAAAAAATTTCAtttt                                                                                                                                                    <  1:421361/176‑1 (MQ=255)                                                                                                                                                                               gagaGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGGTGATGGGCTCCACAAAATTGGGACATCAAAGAAAAGCAGTGGCACTAATTAAAACTGATGCCCTGCGG‑AAAAGTTCTGCGGGTGTGCAAAAAAATTTCATTTTCAGGGCAACTTCAGTTTTATCCTAATCCTGGCCATACCATGACGATGATTTTCCCCGCCACCGTCAGCCGGACGt                                                                         >  1:235465/1‑251 (MQ=255)                                                                                                                                                                               gagaGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGGTGATGGGCCCCACAAAAAGGGGACATCAAAGAAAAGAAGTGGCACTAATTAAGACTGATGCCCTGCGG‑AAAAGTTCTGCGGTTGTGCAAAAAAATTTCAtttt                                                                                                                                                    >  2:421361/1‑176 (MQ=255)                                                                                                                                                                                      ttGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGGTGATGGGCTCCACAAAAAAGGGACAACAAAAAAAAGCAGTGGCACTAATTAAGACTGATGGCCTGCGG‑AAAAGTTCTGCGGGTGTGCAAAAAAAATTTTTTTTCAGGGCAACTTTCGTTTTATTCTAATCCTTGCCATACCATTATGGTGATTGTCCCTGCCAGCGTCATCAGGAAGTTTGCg                                                                    >  2:371889/1‑249 (MQ=255)                                                                                                                                                                                         cAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGGTGATGGGCTCCCCAAAATGGGGACAACAAAGAAAAGCAGGGGCACTAATTAAGACTGATGCCCTGCGG‑AAAAGTTCTGAGGTTGTGCAAAAAAATTTTTTTTTCAGGGAAACTTCAGTTTTATTCTAATCCTGTCCCTACCATGACGATGATTGTCCCTGCCAGCGTCAGCCAGACGTTGGCGATTGc                                                               >  1:49512/1‑251 (MQ=255)                                                                                                                                                                                                               aTAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGAACCACAAAAATGGGACACCAAAGAAAAACAAGGGCCCTAATTAAAACTGATGCCCTGCGGGAAAAATT‑TGCGGTTGGGCAAAAAAAATTTATTTTTAGGGGAACTTCAGTTTTATTCTAATCCAGGCCATCACATTACGATGATTGTCCCCGCCAGCGTCAAAAGGACGTTGGCGATTGGATAGGGGCCCGCaaaaccccgc                                         >  2:437781/1‑242 (MQ=255)                                                                                                                                                                                                                     aaaaaCCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCCCCCCAAAAAAGGGGCATCCAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCCGCGG‑AAAAAATTTGCGGTTGGGCAAAAAAATTTTATTTTTAGGGCCACTTCAGTTTTATTCCAATCCTGGCCATACCATGACGCTGATTGTCCCTGCCCGCCTCAGCAAGACGGTGGCGGTTTCATAGGGGCCCTCATAGCCCAGCCCCggg                                   >  2:469964/1‑251 (MQ=255)                                                                                                                                                                                                                         aCCCATCAACCTTGAACCGAAATGGCGGGGGTGATGGGCTCCACAAAATGGGGGCATCAAAGAAAAACAGTGGCACTAATTAAGACTGATGCCCTGCGG‑AAAAGTTCTGCGGTTGTGCAAAAAAATTTCATTTTCAGGGCAACTTCAGTTTTATCCTAATCCTGGCCATACCATGACGATGATTGTCCCTGCCAGCGTCAGCAGGACGTTGGCGATTGCATAGGTGCCCGCATAGCCCAGCGCCGGGATGt                               >  1:167657/1‑251 (MQ=255)                                                                                                                                                                                                                          cccATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGG‑AAAAGTTCTGCGGTTGTGCAAAAAAATTTCATTTTCAGGGCAACTTCAGTTTTATCCTAATCCTGGCCATACCATGACGATGATTGTCCCTGCCAGCGTCAGCAGGACGTTGGCGATTGCATAGGTGCCCGCATAGCCCAGCGCCGGGATGtt                              <  1:385169/251‑1 (MQ=255)                                                                                                                                                                                                                          cccATCAACCTTGAACCGAAATGGCGGGGGTGAAGGGCTCCCCAAAATGGGGGCATCAAAAAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGG‑AAAAGGTCTGCGGTTGTGCAAAAAAAATTCATTTTTAGGGGAAATCCAGTTTTATCCTAATCCTGGCCCAACCATGACGATGAATGTCCCTGCCAGCGTCAGCAGGACGTTGGCGATTTCATAAGTTCCCCCATAGACCAACGCCGGGATGtt                              >  2:423875/1‑251 (MQ=255)                                                                                                                                                                                                                            cATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGG‑AAAAGTTCTGCGGTTGTGCAAAAAAATTTCATTTTCAGGGCAACTTCAGTTTTATCCTAATCCTGGCCATACCATGACGATGATTGTCCCTGCCAGCGTCAGCAGGACGTTGGCGATTGCATAGGTGCCCGCATAGCCCAGCGCCGGGATGTTAc                            <  2:342040/251‑1 (MQ=255)                                                                                                                                                                                                                               naaCCTTGAACCGAAATGGCGGGGGTGATGGGCTCCACAAAATGGGGACATCAAAAAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGG‑AAAAGGTCTGCGGGTGTGCAAAAAAATTTTATTTTCAGGGCAACTTTAGTTTTATCCTAATCCTGGCCATACCATTACGATGATTTTCCCTGCCAGCGTCAGCAGGACGTTGGCGATTGCATAGGTGCCCGCATAGCCCCGCGCCGGGGTGTTACttc                         >  2:276545/2‑249 (MQ=255)                                                                                                                                                                                                                                   cTTGAACCGAAATGGCGGGGGTGGTGGGCTCCACAAAATGGGGACATCAAAAAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGG‑AAAAGATCTGCGGGTGGGCAAAAAAATTTCATTTTCAGGGCAACTTCAGTTTTATCCTAATCCTGGCCATACCATGACGATGATTGTCCCTGCCCGCGGCAGCACGACGTTGGCGATTGCATAGGGGCCCGCATATCCCAGCGCCGGGATGTTACTGCGAgg                     >  2:178050/1‑250 (MQ=255)                                                                                                                                                                                                                                       aaCCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGG‑AAAAGTTCTGCGGTTGTGCAAAAAAATTTCATTTTCAGGGCAACTTCAGTTTTATCCTAATCCTGGCCATACCATGACGATGATTGTCCCTGCCAGCGTCAGCAGGACGTTGGCGATTGCATAGGTGCCCGCATAGCCCAGCGCCGGGATGTTACTGCGAGCTGTa                 <  1:30026/251‑1 (MQ=255)                                                                                                                                                                                                                                             aaTGGCGGGGGTGATGGGCTCCACAAAAAGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGG‑AAAGGTTCTGCGGTTGGGCAAAAAAATTTCATTTTCAGGGCAACTTCAGTTTTATCCTAGTCCTGGCCATACCATGACGATGATTGTCCCTGCCAGCGTCAGCAGGACGTTGGCGATTGCATAGGTGCCCGCATAGCCCAGCGCCGGGATGTTACTGCGAGCTGTATCAcgg           <  2:147911/251‑3 (MQ=255)                                                                                                                                                                                                                                                ggCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGG‑AAAAGTTCTGCGGTTGTGCAAAAAAATTTCATTTTCAGGGCAACTTCAGTTTTATCCTAATCCTGGCCATACCATGACGATGATTGTCCCTGCCAGCGTCAGCAGGACGTTGGCGATTGCATAGGTGCCCGCATAGCCCAGCGCCGGGATGTTACTGCGAGCTGTATCACtgatg        <  2:185221/251‑1 (MQ=255)                                                                                                                                                                                                                                                ggCGGGGTTGATGGGCTCCACAAAATAGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGG‑AAGAGTTCTGCGGTTGGGCAAAAAAATTTAATTTTCAGGGCAACTTCAGTTTTATCCTAATCCTGGCCATACCATGACGATGATTGTCCATGCCAGCGTCAGCAGGACGTTGGCGATTGCATAGGTGCCCGCATAGCCCAGCGCCGGGATGTTACTGCGAGCTGTATCCCtgatg        >  2:103590/1‑251 (MQ=255)                                                                                                                                                                                                                                                  cGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGG‑AAAAGTTCTGCGGTTGTGCAAAAAAATTTCATTTTCAGGGCAACTTCAGTTTTATCCTAATCCTGGCCATACCATGACGATGATTGTCCCTGCCAGCGTCAGCAGGACGTTGGCGATTGCATAGGTGCCCGCATAGCCCAGCGCCGGGATGTTACTGCGAGCTGTATCACTgatgat      >  2:188194/1‑251 (MQ=255)                                                                                                                                                                                                                                                   ggggTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGG‑AAAAGTTCTGCGGTTGTGCAAAAAAATTTCATTTTCAGGGCAACTTCAGTTTTATCCTAATCCTGGCCATACCATGACGATGATTGTCCCTGCCAGCGTCAGCAGGACGTTGGCGATTGCATAGGTGCCCGCATAGCCCAGCGCCGGGATGTTACTGCGAGCTGTATCACTGATGAtc     <  1:409102/251‑1 (MQ=255)                                                                                                                                                                                                                                                    gggTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGG‑AAAAGTTCTGCGGTTGTGCAAAAAAATTTCATTTTCAGGGCAACTTCAGTTTTATCCTAATCCTGGCCATACCATGACGATGATTGTCCCTGCCAGCGTCAGCAGGACGTTGGCGATTGCATAGGTGCCCGCATAGCCCAGCGCCGGGATGTTACTGCGAGCTGTATCACTGATGAtct    >  1:333057/1‑251 (MQ=255)                                                                                                                                                                                                                                                    gggTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGG‑AAAAGTTCTGCGGTTGTGCAAAAAAATTTCATTTTCAGGGCAACTTCAGTTTTATCCTAATCCTGGCCATACCATGACGATGATTGTCCCTGCCAGCGTCAGCAGGACGTTGGCGATTGCATAGGTGCCCGCATAGCCCAGCGCCGGGATGTTACTGCGAGCTGTATCACTGATGAtct    <  1:686/251‑1 (MQ=255)                                                                                                                                                                                                                                                      gTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGG‑AAAAGTTCTGCGGTTGTGCAAAAAAATTTCATTTTCAGGGCAACTTCAGTTTTATCCTAATCCTGGCCATACCATGACGATGATTGTCCCTGCCAGCGTCAGCAGGACGTTGGCGATTGCATAGGTGCCCGCATAGCCCAGCGCCGGGATGTTACTGCGAGCTGTATCACTGATGATCTcc  <  1:469964/251‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                           AAACCGCTATTAAAAACGGTAATGCAAAACTGGATGCAGCGCAGTCAGCAAACGCTCGAACAATGGTTTGAACCCGGAACCGCCCGCGCGCTTGATGCGTTTATTGAGGGGATGACGCTGCATTTTGTCACCGACCGTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGG‑AAAAGTTCTGCGGTTGTGCAAAAAAATTTCATTTTCAGGGCAACTTCAGTTTTATCCTAATCCTGGCCATACCATGACGATGATTGTCCCTGCCAGCGTCAGCAGGACGTTGGCGATTGCATAGGTGCCCGCATAGCCCAGCGCCGGGATGTTACTGCGAGCTGTATCACTGATGATCTCC  >  NC_000913/887768‑888263

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