Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	1,929,917	C→T	65.9%	intergenic (‑1000/‑168)	USA300HOU_1794 ← / → USA300HOU_1796	hypothetical protein/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	1,929,917	0	C	T	65.9%	57.0 / 47.1	61	intergenic (‑1000/‑168)	USA300HOU_1794/USA300HOU_1796	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (8/13);  new base T (21/19);  total (29/32)
Fisher's exact test for biased strand distribution p-value = 4.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

GTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGA‑TGCTTAT‑‑TTTCGT‑TCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACA‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCT  >  CP000730/1929783‑1930051                                                                                                                                               |                                                                                                                                              ttgCTTGTATTCTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGA‑TGCTTAT‑‑TTTCGT‑TCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCg                                                                                                                                         <  2:245553/140‑1 (MQ=255) atgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GCCGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCAT‑‑TTTCGT‑TCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGt                                                                                                                                        >  1:291399/2‑141 (MQ=255) atgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCAT‑‑TTTCGT‑TCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGt                                                                                                                                        >  1:189774/2‑141 (MQ=255) atgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGC‑TAATCAT‑‑TTTCGT‑TCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGt                                                                                                                                        >  1:163342/2‑141 (MQ=255) atgCTTGTATATTCACAGAGGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTTAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGC‑TAATCAT‑‑TTTCGT‑TCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGt                                                                                                                                        >  1:232664/2‑141 (MQ=255)    cTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCAT‑‑TTTCGT‑TCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtatt                                                                                                                                     <  2:275754/141‑4 (MQ=255)    cTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCAT‑‑TTTCGT‑TCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtatt                                                                                                                                     <  2:216324/141‑4 (MQ=255)     ttGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCAT‑‑TTTCGT‑TCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtatta                                                                                                                                    <  1:251848/141‑5 (MQ=255)     ttGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCAT‑‑TTTCGT‑TCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtatta                                                                                                                                    <  2:173931/141‑5 (MQ=255)        tatattcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑C‑TTTGC‑TAATCAT‑‑TTTCGT‑TCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattatac                                                                                                                                 >  1:261721/11‑137 (MQ=255)        tataTTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑CTT‑TGC‑TAATCAT‑‑TTTCGT‑TCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattatac                                                                                                                                 <  1:96084/141‑8 (MQ=255)          tattcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑C‑T‑T‑‑‑TGCTAATGGTTACTTATCACATCATAGCACCAGTCATCAGGGGCTGGGCTGTTGCGTATTATACtt                                                                                                                               >  1:99374/9‑141 (MQ=255)              cacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCAT‑‑TTTCGT‑TCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTg                                                                                                                           >  2:197643/5‑141 (MQ=255)               acagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCAT‑‑TTTCGT‑TCACATCATAGCACCAGTCATCAGTGGCTGTGCTTTTGCGTATTATACTTATTgg                                                                                                                          <  1:217272/138‑1 (MQ=255)                ccGATGATATAATAGTAGAGGCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGA‑TGCTTAT‑‑TTTCGT‑TCACATCATAGCACCCGTTATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTgg                                                                                                                          <  2:58816/141‑1 (MQ=255)                cagatggtATAATAATAGAGGCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCAT‑‑TTTCGT‑TCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGc                                                                                                                         <  1:175042/139‑1 (MQ=255)                cagagggtataATAATAGAGTCGCCTATCTCTCCGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGC‑TAATCAT‑‑TTTCGT‑TCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGc                                                                                                                         <  2:357515/134‑1 (MQ=255)                 cGAGGATATAATAGGAGAGTCGCCGATGTCTCAGGCGGCAATGTAGGCGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGA‑TGCTTAT‑‑TTTCGT‑TCACATCAGAGCACCAGTCATCAGTGGCTGTGCCACTGCGTTTTTTTCTTATTGGc                                                                                                                         <  1:270609/141‑1 (MQ=255)                 agatggtATAATAAAAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCAT‑‑TTTCGT‑TCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCt                                                                                                                        >  2:243199/2‑141 (MQ=255)                  gatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGC‑TAATCAT‑‑TTTCGT‑TCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCtt                                                                                                                       >  1:87064/1‑140 (MQ=255)                   atggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCAT‑‑TTTCGT‑TCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCtta                                                                                                                      >  1:38217/5‑139 (MQ=255)                    tggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGC‑TAATCAT‑‑TTTCGT‑TCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCttag                                                                                                                     >  1:285321/4‑138 (MQ=255)                    tggtagaataagagagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCAT‑‑TTGCGT‑TCACAGCATAGCACCAGTCATCAGGGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAg                                                                                                                     <  2:305221/129‑1 (MQ=255)                    tggcataataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGC‑TAATCAT‑‑TTTCGT‑TCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCttag                                                                                                                     <  2:241046/137‑4 (MQ=255)                     ggtataATAATAAAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑CTT‑TGC‑TAATCAT‑‑TTTCGT‑TCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTGTACTTATTGGCTTAgg                                                                                                                    >  2:209550/3‑140 (MQ=255)                       tataATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGA‑TGCTTAT‑‑TTTCGT‑TCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTa                                                                                                                   >  1:336584/1‑141 (MQ=255)                       tataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCAT‑‑TTTCGT‑TCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaa                                                                                                                  <  2:173789/141‑2 (MQ=255)                       taaaataGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGA‑TGCTTAT‑‑TTTCGT‑TCGCATCATAGCACCAGTCATCAGTGGCTGTGCCGTTGCGTTTTTTTCTTCATGGCTAAGTa                                                                                                                   >  2:24364/4‑141 (MQ=255)                        ataataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGC‑TAATCAT‑‑TTTCGT‑TCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaaa                                                                                                                 <  2:56056/141‑3 (MQ=255)                           ataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTGT‑C‑C‑TGT‑CAATCAT‑‑TTTCGT‑TCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGc                                                                                                              <  1:82463/141‑1 (MQ=255)                           ataGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGA‑TGCTTAT‑‑TTTCGT‑TCACATCATAGCACCAGTCATCAGTGGCTGTGCCATCGCGTTTTTTTGTTATTGGCTCCGTAGACg                                                                                                               >  1:294474/1‑141 (MQ=255)                            tagtagAGTCGCCTATCTCTCCGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGA‑TGCTTAT‑‑TTTCGT‑TCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGc                                                                                                              <  2:282721/141‑1 (MQ=255)                            tagtagAGTCGCCTATCTCTCAGGCGTCAATTTAGTCGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGA‑TGCTTAT‑‑TTTCGT‑TCACATCATAGCACCAGTCATCAGTGGCTGTTCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGc                                                                                                              <  2:268317/141‑1 (MQ=255)                            taatagAGTCGCCTATCTCTCAGGCGTCGATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCAT‑‑TTTCGT‑TCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa                                                                                                             <  1:249006/138‑1 (MQ=255)                            aaagagagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCAT‑‑TTTCGT‑TCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa                                                                                                             <  1:37940/137‑1 (MQ=255)                                agagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCAT‑‑TTTCGT‑TCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGGGTATTATACTTATTGGCTTAGTAAACGCaacaa                                                                                                         >  1:159474/1‑138 (MQ=255)                                agagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑C‑TTTGC‑TAATCAT‑‑TTTCGT‑TCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTGTTATACTTATTGGCTTAGTGAACGCaacaa                                                                                                         >  1:33425/1‑138 (MQ=255)                                agagTCGCCTATCTCTCAGGCGTCAATT‑‑CACGAAGAGAGGAGGTGCA‑TTA‑C‑T‑TTACTAATCAT‑‑TTTCGT‑TCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAAGGCaacaa                                                                                                         >  1:231223/1‑138 (MQ=255)                                agagTCGCCTAGCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGC‑TAATCAT‑‑TTTCGT‑TCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaa                                                                                                         <  1:63764/141‑4 (MQ=255)                                 gagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCAT‑‑TTTCGT‑TCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaa                                                                                                        >  1:270385/1‑137 (MQ=255)                                       ccTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCAT‑‑TTTCGT‑TCACATCATAGCACCAGTCATCAGTGGATGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatta                                                                                                  >  1:172911/1‑131 (MQ=255)                                       ccTATCTCTCAGGCGTCAA‑‑TAGACGAAGGGAGGAGGTGCA‑TTA‑CTT‑TGC‑TAATCTT‑‑TTTCGT‑TCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatta                                                                                                  >  1:232189/1‑131 (MQ=255)                                         tATCTCTCAGGCGGCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCAT‑‑TTTCGT‑TCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaat                                                                                                <  1:222125/141‑13 (MQ=255)                                         tATCTCTAAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCAT‑‑TTTCGT‑TCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaat                                                                                                <  2:94746/141‑13 (MQ=255)                                          ttctctcAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGA‑TGCTTAT‑‑TTTCGT‑TCACATCATTGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGt                                                                                            <  2:313064/140‑1 (MQ=255)                                          aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCAT‑‑TTTCGT‑TCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaata                                                                                               >  1:32157/1‑128 (MQ=255)                                          aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCAT‑‑TTTCGT‑TCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaata                                                                                               >  1:221566/1‑128 (MQ=255)                                          aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCAT‑‑TTTCGT‑TCACATCATAGCACCAGTCAGCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaata                                                                                               <  1:346969/141‑14 (MQ=255)                                             tctcAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGA‑TGCTTAT‑‑TTTCGT‑TCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTCTTTTCTTATTCGCTAAGTGGACGCACT‑‑CCAC‑‑AATAGGCGAc                                                                                         >  2:230209/1‑141 (MQ=255)                                                    cGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGA‑TGCTTAT‑‑TTTCGT‑TCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGTCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACA‑TAT‑AGc                                                                                >  2:347322/1‑141 (MQ=255)                                                        aaTTTAGACGCAGAAAGGAGGTGTA‑TAA‑G‑G‑TGA‑TGCTTAT‑‑TTTCCT‑TCACATCAAAAAAAAAAAAAAAAGTGGCTGGGACAGTGCGTTTTTTTCTTACTGGCTACGTAGAAAAAAT‑‑ACAA‑‑AGCCGGTTTCA‑TAT‑Agaaaaa                                                                            <  2:155225/141‑6 (MQ=255)                                                                                  taA‑G‑G‑TGA‑TGCTTAT‑‑TTTCGT‑TCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACA‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGc                                                  <  2:212254/141‑1 (MQ=255)                                                                                                  tAT‑‑TTTCGT‑TCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACA‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAggg                                      <  1:249357/141‑1 (MQ=255)                                                                                                    t‑‑tttCGT‑TCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACA‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGa                                    <  1:171248/141‑1 (MQ=255)                                                                                                    t‑‑tttCGT‑TCACACCCTCGCACCAGTCATCAGTGGCTGTGCCATTGCGTTCTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACA‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGa                                    <  1:24364/141‑1 (MQ=255)                                                                                                                    catAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACA‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGt                       >  1:194722/1‑141 (MQ=255)                                                                                                                    catAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACA‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGGGAGGGGGTTGGTGTATAAGt                       >  1:91724/1‑141 (MQ=255)                                                                                                                     atAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCAACAAATAATTAATAGGCGACTTTATAAGCCAGCAACGA‑AAAAATCCCCACACTACTGC‑AATAGTGAGGGGATGGGTGcat                           >  2:125072/1‑141 (MQ=255)                                                                                                                         cACCAGTCATCAGGGGCTGTGCCATTGCGTTTTTTTCTTAGTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACA‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATa                  <  2:241423/141‑1 (MQ=255)                                                                                                                                     gTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACA‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGtt      >  2:309151/1‑141 (MQ=255)                                                                                                                                         cTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACA‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCt  <  1:148616/141‑1 (MQ=255)                                                                                                                                               |                                                                                                                                              GTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGA‑TGCTTAT‑‑TTTCGT‑TCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACA‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCT  >  CP000730/1929783‑1930051

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