Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	1,929,873	0	C	A	53.8%	6.1 / 111.2	65	intergenic (‑956/‑212)	USA300HOU_1794/USA300HOU_1796	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (14/16);  new base A (16/19);  total (30/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.72e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

ATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑T‑AA‑GGT‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAG  >  CP000730/1929735‑1930013                                                                                                                                               |                                                                                                                                             aTTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCGAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑T‑AA‑GGT‑GATGCtt                                                                                                                                            <  2:127767/141‑1 (MQ=255)  ttCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑T‑AA‑GGT‑GATGCTTa                                                                                                                                           <  1:260384/141‑1 (MQ=255)  ttCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑T‑AA‑GGT‑GATGCTTa                                                                                                                                           <  1:35926/141‑1 (MQ=255)  ttCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGGA‑T‑AA‑GGT‑GATGCTTa                                                                                                                                           >  1:349740/1‑141 (MQ=255)           tCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑T‑AA‑GGT‑GATGCTTATTTTCGTTc                                                                                                                                  >  1:365493/1‑141 (MQ=255)                tAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGGGAGGAGGTGTA‑T‑AA‑GGT‑GATGCTTATTTTCGTTCAcatc                                                                                                                             <  1:150342/141‑1 (MQ=255)                tAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑T‑AA‑GGT‑GATGCTTATTTTCGTTCAcatc                                                                                                                             <  2:202223/141‑1 (MQ=255)                tAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑T‑AA‑GGT‑GATGCTTATTTTCGTTCAcatc                                                                                                                             <  1:150344/141‑1 (MQ=255)                tAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑T‑AA‑GGT‑GATGCTTATTTTCGTTCAcatc                                                                                                                             <  1:28454/141‑1 (MQ=255)                      cATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑T‑AA‑GGT‑GATGCTTATTTTCGTTCACATCATAGCa                                                                                                                       >  2:62109/1‑141 (MQ=255)                            aCTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑T‑AA‑GGT‑GATGCTTATTTTCGTTCACATCATAGCACCAGtc                                                                                                                 >  2:249505/1‑141 (MQ=255)                                aTCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑T‑AA‑GGT‑GATGCTTATTTTCGTTCACATCATAGCACCAGtcatca                                                                                                             >  1:290419/1‑141 (MQ=255)                                    atcggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑TTT‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGc                                                                                                        <  2:87914/127‑1 (MQ=255)                                    atcggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑T‑TACTTT‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGc                                                                                                        <  2:279915/127‑1 (MQ=255)                                       gCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑T‑AA‑GGT‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtg                                                                                                      >  2:293589/1‑141 (MQ=255)                                         aagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑TTT‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGc                                                                                                   <  2:265279/132‑1 (MQ=255)                                           gtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TGAC‑TGT‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctg                                                                                                 <  1:164248/134‑3 (MQ=255)                                           gTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑T‑AA‑GGT‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGcc                                                                                                  >  2:186100/1‑141 (MQ=255)                                                 gtgCTTGTATATTTTCCGATTATATAAGAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑T‑AA‑GGT‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCg                                                                                            <  2:104132/141‑1 (MQ=255)                                                 gtgCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑T‑AA‑GGT‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCg                                                                                            <  2:159203/141‑1 (MQ=255)                                                   gCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑T‑AA‑GGT‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGtt                                                                                          <  1:155283/141‑1 (MQ=255)                                                   gCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑TTT‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtat                                                                                         <  1:198307/141‑3 (MQ=255)                                                    cTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑TTT‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtatt                                                                                        <  2:271620/141‑4 (MQ=255)                                                    cTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑TTT‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtatt                                                                                        <  2:94651/141‑4 (MQ=255)                                                    cTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑TTT‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtatt                                                                                        <  2:48027/141‑4 (MQ=255)                                                      tGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑T‑AA‑GGT‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGttttt                                                                                       >  2:279763/1‑141 (MQ=255)                                                      tGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑T‑AA‑GGT‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGttttt                                                                                       >  2:169786/1‑141 (MQ=255)                                                      tGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑T‑AA‑GGT‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGttttt                                                                                       <  1:408472/141‑1 (MQ=255)                                                        tataTTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑T‑AA‑GGT‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGttttttt                                                                                     >  1:79971/1‑141 (MQ=255)                                                            ttcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑TTT‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTAt                                                                                <  1:355912/135‑1 (MQ=255)                                                             tcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑TTT‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTAtt                                                                               >  1:253151/6‑141 (MQ=255)                                                             tcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAA‑‑TGGACGAAGAGAGGAGGTGCATT‑AC‑TTT‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTAtt                                                                               <  2:227860/136‑1 (MQ=255)                                                             tcacagatggtATAATAATAGAGTCGAATAGCTCTCTGGCGTCAATT‑‑GTCGAAGAGAGGACGTGCATT‑AC‑TTT‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCGTCAGTGGCTGTGCTGTTGCGTATTATACTTAtt                                                                               <  2:294776/136‑1 (MQ=255)                                                              cacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑TTT‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTg                                                                              >  1:392296/5‑141 (MQ=255)                                                              cacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑TTT‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTg                                                                              >  1:337127/5‑141 (MQ=255)                                                               acagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑TTT‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTgg                                                                             >  1:205007/4‑141 (MQ=255)                                                                cagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑TTT‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGc                                                                            >  1:150018/3‑141 (MQ=255)                                                                    tggtataATAATAGAGTCGCCTATCTCTCAGGCGTAAATT‑‑GACGCACAGCCGAGCTGCC‑T‑TT‑CTTCGCTAATCCTTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCttag                                                                        <  1:356325/138‑4 (MQ=255)                                                                        ataataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑TTT‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaaa                                                                    <  1:60162/141‑3 (MQ=255)                                                                        ataataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑TTT‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaaa                                                                    >  1:243561/1‑139 (MQ=255)                                                                        ataataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑TTT‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaaa                                                                    >  1:247734/1‑139 (MQ=255)                                                                        ataataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑TTT‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaaa                                                                    <  1:296840/141‑3 (MQ=255)                                                                        ataataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑T‑TACTTT‑GCTAATCATTTGCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaaa                                                                    <  1:104210/141‑3 (MQ=255)                                                                         taataaTAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑TTT‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaaac                                                                   >  2:346788/1‑138 (MQ=255)                                                                         taataaTAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑TTT‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaaac                                                                   >  2:344814/1‑138 (MQ=255)                                                                           ataatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑TTT‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGc                                                                 <  2:133806/137‑1 (MQ=255)                                                                            tagtagAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑T‑AA‑GGT‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGc                                                                 <  2:365493/141‑1 (MQ=255)                                                                                agagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑TTT‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaa                                                            >  2:354594/1‑138 (MQ=255)                                                                                agagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑TTT‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaa                                                            >  1:28047/1‑138 (MQ=255)                                                                                 gagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑TTT‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaa                                                           >  2:33102/1‑137 (MQ=255)                                                                                  agTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑TTT‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaat                                                          <  1:16954/141‑6 (MQ=255)                                                                                  agTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑TTT‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaat                                                          >  2:185088/1‑136 (MQ=255)                                                                                   gTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑T‑AA‑GGT‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACaa                                                          <  1:249505/141‑1 (MQ=255)                                                                                     cGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑T‑AA‑GGT‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACaaaa                                                        >  1:1374/1‑141 (MQ=255)                                                                                     cGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑T‑AA‑GGT‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACaaaa                                                        >  1:416504/1‑141 (MQ=255)                                                                                     cGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑TTT‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataat                                                       >  2:246334/1‑133 (MQ=255)                                                                                       ccTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑TTT‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatta                                                     >  1:177414/1‑131 (MQ=255)                                                                                       ccTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑TTT‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatta                                                     <  1:143421/141‑11 (MQ=255)                                                                                        cTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑TTT‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaa                                                    >  1:151526/1‑130 (MQ=255)                                                                                        cTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATT‑AC‑TTT‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaa                                                    <  2:249963/141‑12 (MQ=255)                                                                                              ctcAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑T‑AA‑GGT‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACa                                               <  2:151615/141‑1 (MQ=255)                                                                                                 aGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑T‑AA‑GGT‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACatat                                            <  1:293589/141‑1 (MQ=255)                                                                                                                      gaggagGTGTA‑T‑AA‑GGT‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATccc                       <  2:259647/141‑1 (MQ=255)                                                                                                                      gaggagGTGTA‑T‑AA‑GGT‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATccc                       >  2:415785/1‑141 (MQ=255)                                                                                                                                  t‑aA‑GGT‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGc            >  1:391145/1‑141 (MQ=255)                                                                                                                                               cTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGAAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAg  <  1:316249/141‑1 (MQ=255)                                                                                                                                               |                                                                                                                                             ATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑T‑AA‑GGT‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAG  >  CP000730/1929735‑1930013

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