Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	76,180	C→T	19.7%	intergenic (+238/‑896)	USA300HOU_0072 → / → USA300HOU_0075	universal stress protein/lysophospholipase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	76,180	0	C	T	19.7%	194.1 / 31.0	72	intergenic (+238/‑896)	USA300HOU_0072/USA300HOU_0075	universal stress protein/lysophospholipase
Reads supporting (aligned to +/- strand):  ref base C (22/35);  new base T (7/7);  total (30/42)
Fisher's exact test for biased strand distribution p-value = 5.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.38e-01

TAATATAACTAATATTGCTGCACCTGTATCACTTAGAACAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATAT  >  CP000730/76046‑76308                                                                                                                                       |                                                                                                                                 tAATATAACTAATATTGCTGCACCTGTATCACTTAGAACAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGt                                                                                                                            <  1:107862/141‑1 (MQ=255) tAATATAACTAATATTGCTGCACCTGTATCACTTAGAACAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGt                                                                                                                            <  2:376884/141‑1 (MQ=255) tAATATAACTAATATTGCTGCACCTGTATCACTTAGAACAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGt                                                                                                                            <  1:29630/141‑1 (MQ=255) tAATATAACTAATATTGCTGCACCTGTATCACTTAGAACAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGt                                                                                                                            <  1:363240/141‑1 (MQ=255)         cTAATATTGCTGCACCTGTATCACTTAGAACAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAAtttt                                                                                                                    <  1:373727/141‑1 (MQ=255)         cTAATATTGCTGCACCTGTATCACTTAGAACAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAAtttt                                                                                                                    <  2:193542/141‑1 (MQ=255)         cTAATATTGCTGCACCTGTATCACTTAGAACAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAAtttt                                                                                                                    <  1:173299/141‑1 (MQ=255)         cTAATATTGCTGCACCTGTATCACTTAGAACAGCTAACCAAAGTGTTAGTAATCCAGGGAAAACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTaaaaaaaa                                                                                                                    <  2:88548/141‑5 (MQ=255)         atAATATTGCTGCACCTGTATCACTTAGAACAGCTAACCAAAGTATTAAAAAAAAAAAAAAAAATAACACAAAGGAAAATAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGTATCTATAgaaaaaaaaa                                                                                                                    <  2:331210/140‑10 (MQ=255)           aaTATTGCTGCACCTGTATCACTTAGAACAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTaa                                                                                                                  <  1:406125/141‑1 (MQ=255)                 gctgcACCTGTATCACTTAGAACAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACaaaaa                                                                                                            <  1:298340/141‑1 (MQ=255)                  ctgcACCTGTATCACTTAGAACAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATCATTATAGCCAAATATCGGTTATGATGCAGATGGGAGACAATAGTATAATTTTAACTAAAAg                                                                                                           >  1:89072/1‑141 (MQ=255)                     cACCTGTATCACTTAGAACAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAg                                                                                                        <  1:241368/141‑1 (MQ=255)                     cACCTGTATCACTTAGAACAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAg                                                                                                        <  1:255794/141‑1 (MQ=255)                     cACCTGTATCACTTAGAACAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTAAAAAAAAAAAAAAAAGGAg                                                                                                        <  2:420745/141‑1 (MQ=255)                          gTATCACTGAGAACAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAAAAAAAAAAAAAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGtcttc                                                                                                   <  2:5385/141‑1 (MQ=37)                                 ttAGAACAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGa                                                                                            <  2:31792/141‑1 (MQ=255)                                  tAGAACAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGaa                                                                                           <  2:171808/141‑1 (MQ=255)                                    gAACAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACt                                                                                         >  1:44193/1‑141 (MQ=255)                                     aaCAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTa                                                                                        <  2:73215/141‑1 (MQ=255)                                     aaCAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTAAAAAAAAAAAAAAAAGGAGTCTTCTGTATGAACTa                                                                                        <  2:319242/141‑1 (MQ=38)                                      aCAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTAt                                                                                       >  2:158814/1‑141 (MQ=255)                                      aCAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTAt                                                                                       >  2:397797/1‑141 (MQ=255)                                      aCAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTAt                                                                                       >  2:397815/1‑141 (MQ=255)                                              ccAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGata                                                                               <  1:347936/141‑1 (MQ=255)                                                 aaGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATaa                                                                            <  2:290475/141‑1 (MQ=255)                                                 aaGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTATAAAAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATaa                                                                            <  2:203014/141‑1 (MQ=255)                                                  tctgtTAGTAATCATTGGAATGCTAACACTAAGGCGATTAATTTAGTAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGTATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCTGCTATaaa                                                                           <  2:111027/139‑1 (MQ=21)                                                  aGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATaaa                                                                           <  1:283211/141‑1 (MQ=255)                                                  aGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATaaa                                                                           <  2:310557/141‑1 (MQ=255)                                                     gtTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACaa                                                                        >  1:198226/1‑141 (MQ=255)                                                           aaTCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAAc                                                                  >  1:138594/1‑141 (MQ=255)                                                                 gggAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGAt                                                            >  2:280606/1‑141 (MQ=255)                                                                 gggAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGAt                                                            >  2:443262/1‑141 (MQ=255)                                                                 gggAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGAt                                                            >  1:315520/1‑141 (MQ=255)                                                                        cTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCa                                                     <  1:229399/141‑1 (MQ=255)                                                                               aaaGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGc                                              >  1:268738/1‑141 (MQ=255)                                                                                    cAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGTCGTTg                                         <  2:44193/141‑1 (MQ=255)                                                                                         aatttaatAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGctac                                    <  2:365330/141‑1 (MQ=255)                                                                                         aatttaatAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGctac                                    <  2:151214/141‑1 (MQ=255)                                                                                         aatttaatAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGctac                                    >  2:192583/1‑141 (MQ=255)                                                                                         aatttaatAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTATAAAAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGctac                                    <  2:222222/141‑1 (MQ=255)                                                                                            ttaatAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTAt                                 >  2:404894/1‑141 (MQ=255)                                                                                               ataatTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTa                              <  2:173725/141‑1 (MQ=255)                                                                                                taatTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTaa                             <  2:236618/141‑1 (MQ=255)                                                                                                         agtaaatgaGGTTCTGTTGCAAAGTTGAATTTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACCCGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCAt                    <  1:221860/133‑1 (MQ=11)                                                                                                         agtaaatgaGGTTCTGTTGCAAAGTTGAATTTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCAt                    <  1:34237/133‑1 (MQ=1)                                                                                                         agtaaatgaGGTTCTGTTGCAAAGTTGAATTTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCAt                    >  1:48855/9‑141 (MQ=11)                                                                                                         agtaaatgaGGTTCTGTTGCAAAGTTGAATTTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCAt                    <  1:164694/133‑1 (MQ=255)                                                                                                         agtaaatgaGGTTCTGTTGCAAAGTTGAATGTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCAt                    >  2:433622/9‑141 (MQ=255)                                                                                                           taaatgaGGTTCTGTTGCAAAGTTGAATTTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTg                  <  1:23642/135‑1 (MQ=255)                                                                                                            aaatgaGGTTCTGTTGCAAAGTTGAATTTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGa                 >  1:316306/6‑141 (MQ=2)                                                                                                            aaatgaGGTTCTGTTGCAAAGTTGAATTTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGa                 <  1:36750/136‑1 (MQ=2)                                                                                                            aaatgaGGTTCTGTTGCAAAGTTGAATTTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGa                 >  2:256594/6‑141 (MQ=0)                                                                                                             atatAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAg                <  1:158814/141‑1 (MQ=39)                                                                                                             atatAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAg                <  1:397797/141‑1 (MQ=40)                                                                                                             atatAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAg                <  1:397815/141‑1 (MQ=39)                                                                                                             aatgaGGTTCTGTTGCAAAGTTGAATTTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAg                <  1:154003/137‑1 (MQ=2)                                                                                                              atgaGGTTCTGTTGCAAAGTTGAATTTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGt               >  1:254230/4‑141 (MQ=2)                                                                                                              atgaGGTTCTGTTGCAAAGTTGAATTTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGt               <  2:374503/138‑1 (MQ=2)                                                                                                              atgaGGTTCTGTTGCAAAGTTGAATTTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGt               >  1:340089/4‑141 (MQ=2)                                                                                                                taGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTa             >  1:402867/1‑141 (MQ=38)                                                                                                                taGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTa             >  1:60858/1‑141 (MQ=37)                                                                                                                taGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTa             >  2:270018/1‑141 (MQ=38)                                                                                                                taGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTa             >  1:402674/1‑141 (MQ=38)                                                                                                                gaGGTTCTGTTGCAAAGTTGAATTTATAATATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAgcta             >  2:108591/2‑138 (MQ=0)                                                                                                                gaGGTTCTGTTGCAAAGTTGAATTTATAATATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAgcta             >  1:426075/2‑138 (MQ=2)                                                                                                                   gTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTATCg          >  1:151437/1‑141 (MQ=35)                                                                                                                   gTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTATCg          >  1:151405/1‑141 (MQ=35)                                                                                                                     tCTGTTGCAAAGTTGAATCTATAATATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGCTATCGTg        >  2:228652/1‑141 (MQ=21)                                                                                                                           gCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGatatat  >  2:450226/1‑141 (MQ=35)                                                                                                                           gCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGatatat  >  2:159879/1‑141 (MQ=35)                                                                                                                                       |                                                                                                                                 TAATATAACTAATATTGCTGCACCTGTATCACTTAGAACAGCTAACCAAAGTGTTAGTAATCCAGGGAATACTAACACAAAGGCAATTAATTTAATAATTATAGCAAAATATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATAT  >  CP000730/76046‑76308

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑