Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	887,271	C→T	L23L (TTG→TTA)	ybjJ ←	inner membrane protein YbjJ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	887,271	0	C	T	96.9%	188.6 / ‑2.4	64	L23L (TTG→TTA)	ybjJ	inner membrane protein YbjJ
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base T (25/37);  total (25/39)
Fisher's exact test for biased strand distribution p-value = 5.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.71e-01

GAGCAGGGGCGATGTCAGCCAGAGTGCCAGACTTAATATCATCATCCCGATCAATGCGCAGGACATCGTGACCAGGATGACATTACGTGTCCCAAAGCGTTTCACTAACCACGCCGAGCAGAGAATACCGCTCATCGAACCGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGAGAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATACATGCTGTTACGCAC  >  NC_000913/887025‑887515                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                     gagCAGGGGCGATGTCAGCCAGAGTGCCAGACTTAATATCATCATCCCGATCAATGCGCAGGACATCGTGACCAGGATGACATTACGTGTCCCAAAGCGTTTCACTAACCACGCCGAGCAGAGAATACCGCTCATCGAACCGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGAGAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGTaaaa                                                                                                                                                                                                                                                  <  1:410344/251‑1 (MQ=255)       gggCGATGTCAGCCAGAGTGCCAGACTTAATATCATCATCCCGATCAATGCGCAGGACATCGTGACCAGGATGACATTACGTGTCCCAAAGCGTTTCACTAACCACGCCGAGCAGAGAATACCGCTCATCGAACCGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGAGAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGTAAAAagaaga                                                                                                                                                                                                                                            <  2:456333/251‑1 (MQ=255)           gATGTCAGCCAGAGTGCCAGACTTAATATCATCATCCCGATCAATGCGCAGGACATCGTGACCAGGATGACATTACGTGTCCCAAAGCGTTTCACTAACCACGCCGAGCAGAGAATACCGCTCATCGAACCGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGAGAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGTAAAAAGAAGAACAt                                                                                                                                                                                                                                        >  2:536735/1‑251 (MQ=255)                   ccAGAGTGCCAGACTTAATATCATCATCCCGATCAATGCGCAGGACATCGTGACCAGGATGACATTACGTGTCCCAAAGCGTTTCACTAACCACGCCGAGCAGAGAATACCGCTCATCGAACCGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGAGAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGTAAAAAGAAGAACATAAACAgcg                                                                                                                                                                                                                                <  1:244346/251‑1 (MQ=255)                            cataCTTAATATCATCATCCCGATCAATGCGCAGGACAGCGGGACCAGGATGACATTACGTGTCCCAAAGCGTTCCACTAACCACGCCGAGCAGAGAATACCGCTCATCGAACCGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGAGAATATCGCGGATAGCAGACGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGTAAAAAGAAGAACATAAACAGCGCCCAGGTTc                                                                                                                                                                                                                       <  2:51660/248‑1 (MQ=255)                             aGACTTAATATCATCATCCCGATCAATGCGCAGGACATCGTGACCAGGATGACATTACGTGTCCCAAAGCGTTTCACTAACCACGCCGAGCAGAGAATACCGCTCATCGAACCGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGAGAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGTAAAAAGAAGAACATAAACAGCGCCCAGGTTCg                                                                                                                                                                                                                      >  2:142741/1‑251 (MQ=255)                             aGACTTAATATCATCATCCCGATCAATGCGCAGGACATCGCGACCAGGATGACATTACGTGTCCCAAAGCGTTTCACTAACCACGCCGAGCAGAGAATACCGCTCATCGAACCGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGAGAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGTAAAAAGAAGAACATAAACAGCGCCCAGGTTCg                                                                                                                                                                                                                      >  2:119075/1‑251 (MQ=255)                                    atatCATCATCCCGATCAATGCGCAGGACATCGTGACCAGGATGACATTACGTGTCCCAAAGCGTTTCACTAACCACGCCGAGCAGAGAATACCGCTCATCGAACCGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGAGAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGTAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTc                                                                                                                                                                                                               >  1:500904/1‑251 (MQ=255)                                        catcatcCCGATCAATGCGCAGGACATCGTGACCAGGATGACATTACGTGTCCCAAAGCGTTTCACTAACCACGCCGAGCAGAGAATACCGCTCATCGAACCGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGAGAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATg                                                                                                                                                                                                           <  1:320835/251‑1 (MQ=255)                                           catcCCGATCAATGCGCAGGACATCGTGACCAGGATTACATTACGTGTCCCAAAGCGTTTCACTAACCACGCCGAGCAGAGAATACCGCTCATCGAACCGATCGACAGACCAAAGAGATCACCGCCCATTTCAGCGATCGAGACAGAGAGAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGTAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCAt                                                                                                                                                                                                        <  1:457881/251‑1 (MQ=255)                                               ccGATCAATGCGCAGGACATCGTGACCAGGATGACATTACGTGTCCCAAAGCGTTTCACTAACCACGCCGAGCAGAGAATACCGCTCATCGAACCGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGAGAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGTAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACg                                                                                                                                                                                                    <  1:434708/251‑1 (MQ=255)                                                    cAATGAGCAGGACATCGGGACCAGGATGACATAACGTGTCCCAAAGCGTTTCACTAACCACGCCGAGCAGAGAATACCGCTCATCGAACCGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGAGAATATCGCGGATAGCAGGCGTACGGGTGGCCCAGGACGCCATTAACAGGCCTGGTAAAAAGAAGAACATAAACAGCGACCAGGTTCGGCGGTTCAATGCATTACGTgaa                                                                                                                                                                                                <  2:462726/250‑1 (MQ=255)                                                     aaTGCGCAGGACATCGTGACCAGGATGACATTACGTGTCCCAAAGCGTTTCACTAACCACGCCGAGCAGAGAATACCGCTCAACGAACCGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGAGAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGTAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTgaaga                                                                                                                                                                                              <  1:265813/251‑1 (MQ=255)                                                         cgcAGGACATCGTGACCAGGATGACATTACGTGTCCCAAAGCGTTTCACTAACCACGCCGAGCAGAGAATACCGCTCATCGAACCGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGAGAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGTAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAAttt                                                                                                                                                                                          >  2:260037/1‑251 (MQ=255)                                                                    gTGACCAGGATGACATTACGTGTCCCAAAGCGTTTCACTAACCACGCCGAGCAGAGAATACCGCTCATCGAACCGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGAGAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGTAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGa                                                                                                                                                                               <  1:270636/251‑1 (MQ=255)                                                                          aGGATGACATTACGTGTCCCAAAGCGTTTCACTAACCACGCCGAGCAGAGAATACCGCTCATCGAACCGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGAGAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGTAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGt                                                                                                                                                                         <  1:174100/251‑1 (MQ=255)                                                                                  aTTACGTGTCCCAAAGCGTTTCACTAACCACGCCGAGCAGAGAATACCGCTCATCGAACCGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGAGAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGTAAAAAGAAGAACAt                                                                                                                                                                                                                                        <  2:450643/180‑1 (MQ=255)                                                                                  aTTACGTGTCCCAAAGCGTTTCACTAACCACGCCGAGCAGAGAATACCGCTCATCGAACCGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGAGAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGTAAAAAGAAGAACAt                                                                                                                                                                                                                                        >  1:450643/1‑180 (MQ=255)                                                                                  aTTACGTGTCCCAAAGCGTTTCACTAACCACGCCGAGCAGAGAATACCGCTCATCGAACCGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGAGAATATCGCGGATAGCAGGCGTACGGGATGCCCAGGACGCCATTAACAGGCCTGGTAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATa                                                                                                                                                                  >  1:374816/1‑250 (MQ=255)                                                                                      cGTGTCCCAAAGCGTTTCACTAACCACGCCGAGCAGAGAATACCGCTCATCGAACCGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGAGAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGTAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATaagaag                                                                                                                                                             <  1:188447/251‑1 (MQ=255)                                                                                                    tttCACTAACCACGCCGAGCAGAGAATACCGCTCATCGAACCGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGAGAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGTAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGacaaca                                                                                                                                                >  2:10028/1‑250 (MQ=255)                                                                                                      tCACTAACCACGCCGAGCAGAGAATACCGCTCATCGAACCGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGAGAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACAtt                                                                                                                                              <  1:449476/250‑1 (MQ=255)                                                                                                           aaCCACGCCGAGCAGAGAATACCGCTCATCGAACCGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGAGAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGTAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAg                                                                                                                                        <  1:145366/251‑1 (MQ=255)                                                                                                            aCCACGCCGAGCAGAGAATACCGCTCATCGAACCGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGAGAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGTAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAgg                                                                                                                                       <  1:104987/251‑1 (MQ=255)                                                                                                                     aGCAGAGAATACCGCTCATCGAACCGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGAGAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGTAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTAc                                                                                                                               >  1:384057/1‑250 (MQ=255)                                                                                                                               aCCGCTCATCGAACCGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGAGAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGTAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTacac                                                                                                                     >  1:520664/1‑250 (MQ=255)                                                                                                                                    tCATCGAACCGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGAGAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGTAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGc                                                                                                               >  2:471008/1‑251 (MQ=255)                                                                                                                                        cGAACTGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGAGAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGTAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCaga                                                                                                            <  1:260037/250‑1 (MQ=255)                                                                                                                                        cGAACCGATCGACATACCAAATAGAACACCGCCCATTGCAGCGAGCGAGACAGAGAGAATATCGCGGAGAGCATGCGTACGGGTTGCCCAGGACGCCATTATCAGGCCTGGTAAAAAGAAGAACATAAACAGCACCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAAGTTACGGTCATAGATCACGTCGAAATAATAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCaga                                                                                                            <  1:52325/250‑1 (MQ=255)                                                                                                                                         gAACCGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGAGAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGTAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCaa                                                                                                                                                                                                             <  1:351649/152‑1 (MQ=255)                                                                                                                                         gAACCGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGAGAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGTAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCaa                                                                                                                                                                                                             >  2:351649/1‑152 (MQ=255)                                                                                                                                             cGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGAGAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGTAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAgg                                                                                                      >  2:53517/1‑251 (MQ=255)                                                                                                                                             cGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGAAAGAGAGAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGTAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAgg                                                                                                      >  1:562254/1‑251 (MQ=255)                                                                                                                                              gATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGAGAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGTAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAAttt                                                                                                                                                                                          <  1:118629/166‑1 (MQ=255)                                                                                                                                              gATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGAGAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGTAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAAttt                                                                                                                                                                                          >  2:118629/1‑166 (MQ=255)                                                                                                                                                   aCAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGAGAATATCACGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGTAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATg                                                                                                 >  1:292693/1‑250 (MQ=255)                                                                                                                                                         cAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGAGAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGTAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCg                                                                                          <  2:384057/251‑1 (MQ=255)                                                                                                                                                           aaGAGAACACCGCCCATTTCAGCGATCGAGACAGAGAGAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGTAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTc                                                                                        <  1:53517/251‑1 (MQ=255)                                                                                                                                                                  cacCGCCCATTTCAGAGATCGAGACAGAGAGAATATAGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGTAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTa                                                                                  >  2:360454/1‑250 (MQ=255)                                                                                                                                                                        ccATTTCAGCGATCGAGACAGAGAGAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGTAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATcc                                                                           <  2:67546/251‑1 (MQ=255)                                                                                                                                                                        ccATTTCAGCGATCGAGACAGAGAGAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGACTGGTAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATcc                                                                           <  2:602568/251‑1 (MQ=255)                                                                                                                                                                         cctttCAGCGATCGAGACAGAGAGAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGTAAAAAGAAGAACCTAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCg                                                                          >  2:527680/3‑251 (MQ=255)                                                                                                                                                                         cATTTCAGCGATCGAGACAGAGAGAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGTAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCg                                                                          <  1:578113/251‑1 (MQ=255)                                                                                                                                                                          aTTTCAGCGATCGAGACAGAGAGAATATCGCGGCTAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGTAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGa                                                                                                     >  2:509353/1‑223 (MQ=255)                                                                                                                                                                          aTTTCAGCGATCGAGACAGAGAGAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGTAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGa                                                                                                     <  1:509353/223‑1 (MQ=255)                                                                                                                                                                              cAGCGATCGAGACAGAGAGAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGTAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACg                                                                     <  1:142741/251‑1 (MQ=255)                                                                                                                                                                                     cGAGACAGAGAGACAATCGCGGGTAGCAGGCGAACGGGTTGCCCAGGACGCCATTAACAGGCATGGTACAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGACGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATATGTACACAATTGCAGAGAGAGGAAATGACATGCGTCCCGCTAACGATCCGCCACGACgagc                                                               >  2:120885/1‑249 (MQ=255)                                                                                                                                                                                        gaCAGAGAGAATAGCGCGGATAGCAGTCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGTAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAAGTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACTGTCCGCAACGGCGAGaaaaa                                                           <  2:317111/251‑1 (MQ=255)                                                                                                                                                                                            gagagaATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGTAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAAtt                                                        >  1:346111/1‑250 (MQ=255)                                                                                                                                                                                                gaATATCTCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGTAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCa                                                   <  1:339803/251‑1 (MQ=255)                                                                                                                                                                                                    atCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGTAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGcc                                               <  1:135385/251‑1 (MQ=255)                                                                                                                                                                                                           aTAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGTAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTgg                                        <  2:500904/251‑1 (MQ=255)                                                                                                                                                                                                             aGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGTAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGa                                       <  1:87613/250‑1 (MQ=255)                                                                                                                                                                                                               cAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGTAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGc                                    <  2:629049/251‑1 (MQ=255)                                                                                                                                                                                                               cAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCAGGTAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGAGCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGc                                    <  2:73759/251‑1 (MQ=255)                                                                                                                                                                                                                  gCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGTAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGATGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGt                                 <  1:72462/251‑1 (MQ=255)                                                                                                                                                                                                                   cGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGTAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTGCGTGAAGAATTTACGGTCATAGATCGTGTCAAAATAAGAAGAGGGAAGACAACATAAGCAAGGTTGTGTACATTTGTACACAATTGCAGTGATAGGAAATGACATGCGTCGCGCGAACTATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCGGGAGGCGGTg                                <  1:398230/251‑1 (MQ=255)                                                                                                                                                                                                                             gCCCAGGACGCCATTAACAGGCCTGGTAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACg                      <  1:5680/251‑1 (MQ=255)                                                                                                                                                                                                                                     cGCCATTAACAGGCCTGGTAAAACGAAGCACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATACAt              >  2:210771/1‑251 (MQ=255)                                                                                                                                                                                                                                         aTTAACAGGCCTGGTAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTAAATGCATTACGTGAAGAATTTACGGTCCTAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATAAATGCTg          >  2:426614/1‑251 (MQ=255)                                                                                                                                                                                                                                         aTTAACAGGCCTGGTAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTCTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAAATTTACGGAATACctgctg          >  1:562170/1‑251 (MQ=255)                                                                                                                                                                                                                                           tAACAGGCCTGGTAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACAAAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTCCGGAAAACATGCGGtt        >  2:566476/1‑251 (MQ=255)                                                                                                                                                                                                                                           tAACAGGCCTGGTAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAAGTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCGGGAGGCGGTGAAACTGTACGGAATACATGCTGt         <  2:375024/250‑1 (MQ=255)                                                                                                                                                                                                                                                 gCCTGGTAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATACATGCTGTTACGCAc  <  1:426614/251‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                     GAGCAGGGGCGATGTCAGCCAGAGTGCCAGACTTAATATCATCATCCCGATCAATGCGCAGGACATCGTGACCAGGATGACATTACGTGTCCCAAAGCGTTTCACTAACCACGCCGAGCAGAGAATACCGCTCATCGAACCGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGAGAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATACATGCTGTTACGCAC  >  NC_000913/887025‑887515

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